6 Nov

警察捉贼,追牛问题,导弹跟踪

王二小的牛跑了,当他发现时,牛在他正南方300米。且一直向正西方向匀速的跑,王二小立即追牛,他不是朝着一个固定的方向,而是每时每刻都朝着牛的方向跑,且速度是牛速度的4/3倍。当他追上牛时王二小共跑了多远?

问题分析

米拉斯反潜导弹

米拉斯反潜导弹

咋看起来,追牛和导弹是风牛马不相及的两件事:一个是生活小事,一个是物理问题,怎么能够扯到一块呢?

回想一下平时警察抓小偷的过程。警察不是物理学家,不会也可不能先去研究小偷的逃走路线函数,然后设计最小追赶时间的路程吧?那么,在不能预知小偷逃跑路线的前提下,警察要怎样捉小偷呢?很简单,盯死他!是的,只要你以更快的速度,一直朝着他跑,总能够追到的。继续联想下:要想用导弹跟踪摧毁一首敌舰,不也是只能够采用这种方式吗?回看文章开始的“追牛问题”,本质上不是一样的吗?以下是上海交大提出的导弹跟踪问题:

点击阅读全文...

14 May

双固定引力中心问题

我们在研究地球附近的小天体运动时,如果把天体和地球看作一个二体系统,那最多只能算上一个零级近似,如果使用“地球+月球+小天体”组成的圆形限制性三体问题模型,那可以算上一个二级近似了。那么,一级近似又是什么了。BoJone认为,它就是本文将要讲的“双固定引力中心问题”了,也叫“双不动中心问题”,英文名是two fixed-center problem。这是一种特殊的限制性三体问题。在这个三体系统中,两个主天体(或称有限质量天体)固定不动,第三个小天体在两个固定的主天体吸引下运动。欧拉、拉格朗日、勒让德、雅可比等人很早就研究过这个问题。其中,欧拉最先成功地求出了这个系统的积分。[引用]

另外,双固定引力中心问题还有另外一个应用,在研究人造卫星的运动时,可以只考虑地球引力,但是由于地球不是完美的球体,把其看成一个质点其实不十分精确,要是把它拆分为两个引力源,就可以很大程度上提高精确度。毕竟双固定引力中心问题是完全可以积分的,可以作为一个比较好的中间轨道(介乎圆锥曲线和精确轨道之间的)。

点击阅读全文...

8 Jul

一道比较函数大小的题目

前几天刚结束的云浮高二期末考数学试卷中,有一道题目让我比较深刻。因为在当时我无法去证明它,只是用了举例子的方法得出了答案。刚才思考了一下,在此给出证明过程。题目如下:

定义在(0,+∞)的函数f(x)满足$x f'(x) \leq f(x)$,对于任意的0 < a < b,比较$a f(b)$和$b f(a)$的大小。

点击阅读全文...

26 Sep

数学基本技艺之23、24(上)

23、求解拟齐次方程$\frac{dy}{dx}=x+\frac{x^3}{y}$
24、求解拟齐次方程$\ddot{x}=x^5+x^2\dot{x}$

把这两道题目放在一起说是因为我觉得这两道题目本质上是一样的,当然,不管怎样,24题更复杂一些。在24题中,设$\dot{x}=y$,则$\ddot{x}=y\frac{dy}{dx}$,于是原方程就变成:
$$\frac{dy}{dx}=x^2+\frac{x^5}{y}$$
这样就跟23题的形式差不多了。

点击阅读全文...

6 Oct

从马尔科夫过程到主方程(推导过程)

主方程(master equation)是对随机过程进行建模的重要方法,它代表着马尔科夫过程的微分形式,我们的专业主要工具之一就是主方程,说宏大一点,量子力学和统计力学等也不外乎是主方程的一个特例。

然而,笔者阅读了几个著作,比如《统计物理现代教程》,还有我导师的《生物系统的随机动力学》,我发现这些著作对于主方程的推导都很模糊,他们在着力解释结果的意义,但并不说明结果的思想来源,因此其过程难以让人信服。而知乎上有人提问《如何理解马尔科夫过程的主方程的推导过程?》但没有得到很好的答案,也表明了这个事实。

马尔可夫过程

主方程是用来描述马尔科夫过程的,而马尔科夫过程可以理解为运动的无记忆性,说通俗点,就是下一刻的概率分布,只跟当前时刻有关,跟历史状态无关。用概率公式写出来就是(这里只考虑连续型概率,因此这里的$p$是概率密度):
$$\begin{equation}\label{eq:maerkefu}p(x,\tau)=\int p(x,\tau|y,t) p(y,t) dy\end{equation}$$
这里的积分区域是全空间。这里的$p(x,\tau|y,t)$称为跃迁概率,即已经确定了$t$时刻来到了$y$位置后、在$\tau$时刻达到$x$的概率密度,这个式子的物理意义是很明显的,就不多做解释了。

点击阅读全文...

22 Oct

分享:孟岩的《理解矩阵》一文

数学演算

数学演算

之前已经提到我要自学相对论和量子力学。作为现代物理的两大支柱,所用的数学也是很“现代”的,不能总是用高中那套简单的模式来计算,所以线性代数是我要熟悉的一门课程之一。现在大一还没开设线性代数课程,但是我所持的观点是:“任何东西只要你需要它,你就应该去学,而且能够学会。”其实我初三暑假的时候就开始接触了线性代数,我看的那本教材,跟国内其他线性代数教材一样,采用了一种只要求记忆和计算的方式来教授,先讲从线性方程组引出行列式,再到矩阵。我那时也在背诵,知道了了行列式怎么算的,行列式可以用来解方程组,矩阵是怎么相乘的等等。但我完全不知道为什么,我甚至不懂为什么这门课程叫“线性代数”。(当然,也有可能是那时的数学水平不够)国外很多教程都讲的很好,很规范地教,但是对于国内像我这样平庸的学生又显得过于专业。我一直期待有这样的一个平衡点,可惜一直没有找到,所以只能从各种渠道摸索。

点击阅读全文...

28 Feb

行列式的导数

在讨论曲线坐标系的积分时,通常都会出现行列式这个东西,作为“体积元”的因子。在广义相对论中,爱因斯坦场方程的作用量就带有度规的行列式,而在对其进行变分时,自然也就涉及到了行列式的求导问题。我参考了朗道的《场论》以及《数理物理基础--物理需用线性高等数学导引》,了解到相关结果,遂记录如下。

推导


\begin{equation}\boldsymbol{A}(t)=\left(a_{ij}(t)\right)_{n\times n}\end{equation}
是一个n阶矩阵,其中每个矩阵元素都是t的函数。其行列式为$|\boldsymbol{A}|$,自然地,考虑
\begin{equation}\frac{d}{dt}|\boldsymbol{A}|\end{equation}

点击阅读全文...

15 Jul

《新理解矩阵6》:为什么只有方阵有行列式?

学过线性代数的朋友都知道,方阵和非方阵的一个明显不同是,对于方阵我们可以计算它的行列式,如果不是方阵的话,就没有行列式这个概念了。在追求统一和谐的数学系统中,为什么非方阵却没有行列式?也许对于这个问题最恰当的回答是——因为不够美。对于非方阵,其实也可以类似地定义它的行列式,定义出来的东西,跟方阵的行列式具有同样的性质,比如某行乘上一个常数,行列式值也就乘以一个常数,等等;而且还可以把其几何意义保留下来。但是,非方阵的行列式是不够美的,因为对于一个一般的整数元素的方阵,我们的行列式是一个整数;而对于一个一般的整数元素的非方阵,却导致了一个无理数的行列式值。另外,一个也比较重要的原因是,单单是方阵的行列式也够用了。综合以上两个理由,非方阵的行列式就被舍弃不用了。

非方阵的行列式不够漂亮

$n$阶方阵的行列式是每个向量的线性函数,它代表着向量之间的线性相关性;从几何上来讲,它就是向量组成的平行n维体的(有向)体积。我们当然期望非方阵的行列式也保留这些性质,因为只有这样,方阵行列式的那些运算性质才得以保留,比如上面说的,行列式的一行乘上一个常数,行列式值也乘上一个常数。我们考虑$m\times n$的矩阵,其中$ m < n $,我们将它看成是$m$个$n$维向量的组合。最简单的,我们先考虑$1\times 2$矩阵的行列式,也就是二维向量$(a,b)$的行列式。

点击阅读全文...