8 Jul

均匀球状星团内恒星的运动

我们考虑一个球状的星团,并假设它是各向同性的,即距离球心r处的物质密度ρ只与r有关,ρ=ρ(r)。那么,在半径为r的球形区域内的总质量为:
$$M(r)=\int_0^r 4\pi x^2 \rho(x) dx$$

想象有一颗质量比较小的恒星(其实相对于星团总质量,每一颗恒星的质量都很小)在星团的引力作用下运动(就好像太阳系绕着银河系运动一样),且恒星并没有受到其他物质(如星际尘埃等)的阻力。我们之前已经证明过,各向同性的球壳内部的引力是为0的,那么这种情况下的运动就相当于恒星只受到它到球心处的一个球形区域内的质量的引力吸引。根据万有引力定律,选择星团球心为参考系,可以得出
$$\ddot{\vec{r}}=-GM(r)\frac{\vec{r}}{r^3}$$

点击阅读全文...

25 Jul

收谷问题(1)

在农村,7月是忙碌的月份,农民们要忙着收割稻谷,收割完后要晒谷,同时还得准备“下秧”,准备新一轮的耕,BoJone家自然也不例外。不过我家田比较少(1亩左右),收割机几分钟搞定,谷也三两天就晒完了。不过在晒谷的时候,BoJone在考虑一个“收谷”问题:

晒谷时得先把成堆的谷子摊开,薄薄地平铺在平地上,等到傍晚或即将下雨时(这是最惨的情况,搞不好会淋谷)就将其收起来。问题就源于这里,一般来说我们会把谷均匀地铺成矩形,要把所有的谷都推到矩形里或外的哪一点上,才使得我们做功做小?

这个问题还可以推广开来,例如对于一地任意形状的谷子(如三角形),把它集中堆到哪个点最“轻松”?一堆固定质量的谷子,要把它平铺成什么形状,才使得收谷时最“轻松”?当然,这个问题的解不仅仅用于“收谷”,在很多规划建设中也可以应用到,例如要在一个人口大致均匀的城市中建设一个服务中心,这个服务中心应该建在哪里?这有点类似于我们之前讨论过的费马点问题 ,都是费马点只考虑了三个点的距离,而这个问题得考虑所有点的距离。

点击阅读全文...

21 Aug

有理直角三角形的面积能否为整数?

这是一个古老而有趣的问题,但在引入这个问题之前,我们首先来看一个简单的问题:

整数边直角三角形的面积能否为一个完全平方数?

答案是不能。我们可以举一些例子来检验一下,例如边长为3,4,5的直角三角形面积为6,6不是一个平方数;再如边长为5,12,13的直角三角形面积为30,30也不是一个平方数...当然,数学的最近目的是要求严格证明,而不是简单举例,否则就只得称为不完全归纳,这样得出来的是一个猜想,而不是“定理”,就好象著名的“哥德巴赫猜想”...本文我们将试图证明这个命题。

我们稍后还会发现,这个问题和以下问题是等价的:
是否存在一个面积为1的三边长都是有理数的直角三角形?

更让人意外的是,这个问题也等价于方程$x^4+y^4=z^4$并没有整数解,换句话说,我们要证明n=4时的“费马大定理”

点击阅读全文...

14 Jan

诡异的Dirac函数

量子力学中有一个很诡异的函数——Dirac函数,它似乎在物理的不少领域都有很大作用,它也具有明显的物理意义,但认真地看它却又感觉它根本就不是函数!这个“似而非是”的东西究竟是什么呢?让我们从一个物理问题引入:

设想一条质量为1,长度为$2l$的均匀直线,很显然直线的密度为$\rho=\frac{1}{2l}$;将直线的中点放置于坐标轴的原点,我们就有
$$\rho(x)=\left\{ \begin{array}{c}\frac{1}{2l} (-l \leq x \leq l)\\0 (x < -l , x > l)\end{array}\right.$$

所以有
$$\int_{-\infty}^{+\infty} \rho(x)dx=1$$

点击阅读全文...

24 Jun

为方轮自行车铺路

方轮自行车

方轮自行车

你见过正方形轮子的自行车吗?一般认为,只有圆形的车轮才能使我们的车子平稳向前移动,但这只是针对平直道路而言的。谁规定路一定是平的?只要铺好一条适当的道路,正方形车轮的自行车照样可以平稳前行!本文就让我们为方轮自行车铺一条路。

其实,方轮自行车已经不是新鲜玩意了,它早已出现在不少科技馆中。从图片中可以看到,它的特殊轨道是有许多段弧组成的,每一段弧的长度等于正方形的边长。车轮前行时,正方形会保持与弧形相切(确保不会打滑)。这样的路的形状是什么曲线呢?很幸运,它并不十分复杂,而且让人意外的是,它就是我们之前已经研究过的“悬链线”!原来,要设计这样的一个曲线的轨道,不需要多么高深的设计师,只需要我们手拿一条铁链,让它自由垂下......

点击阅读全文...

28 Sep

开始学习数学软件Scilab

其实很早之前我就想学习一款数学软件的使用,以前很感兴趣的是mathematica,也玩弄过一阵子,但毕竟在高中没有多大需要,也就没有坚持下来。更重要的是,这些软件都是要收费的。上了大学后,听了师兄姐对数学建模的讲述,发现他们基本上也是用mathematica或者matlab的,但这两个软件都是要收费的,我不大想用破解版本。既然我都已经用上了ubuntu了,那么我就该好好利用它。据说命令跟matlab很相似的软件是scilab,还有octave,不同的是这些都是开源免费的。

出于熟悉代码操作和数学软件编程的目的,我选择了学习scilab。虽然网上说octave与matlab的相似程度更高,但是我感觉scilab比octave用的更广一些,所以就用它。所谓“一理通百理明”,先专心学好一个。

下面是我编写的第一个scialb程序,利用威尔逊方法来进行素性测试。这个代码的主要目的是练习条件语句和循环语句,以及一些输出输入的技巧而已。程序本身比较丑陋。

//我的第一个scilab程序
//完成于2012.09.27

label1=['p:';];  //定义标签
B=x_mdialog(['本程序使用威尔逊方法判断进行素数测试。';'请输入要判断的数'],label1,['127';]);  //输入框
p=evstr(B(1));  //提取输入框里边的数字进行赋值
i=1;
j=1;
q=p-1;
while i<q
    j=j*i;
    j=modulo(j,p);//这个是模函数。
    i=i+1;
end
if j==1
    messagebox(['这是一个素数';],['测试结果']);  //输出,其中后边的“测试结果”是输入框的标题
else
    messagebox(['这是一个合数';],['测试结果']);
end

点击阅读全文...

18 Oct

证明光速不变的一个理想实验??

在狭义相对论发表之前和之后,都有不少实验从不同角度论证了它的正确性。这些实验大多数是实际测量得出结果的,当然也存在着一些“理想实验”,这些实验只需要一定的逻辑推理,而实际上是无法完成的。下面就是我很久之前在某本书(很抱歉,我真的忘记书名了)看到的一个用来推翻光速可叠加的伽利略变换的理想实验。它只用寥寥几句,就好像已经证明了“c+c=c”(c是真空中的光速)的事实。可是“c+c=c”在狭义相对论上是作为原理出现的,是不可能通过逻辑推理来证明的。事实究竟如何?我们先来看这个实验。

光速不变的理想实验

光速不变的理想实验

任意选定一个坐标原点。设想原点的正北方$c\cdot t_0$处有一架以光速$c$朝南运行的飞机1;原点的正西方$c\cdot t_0$处有一架以光速$c$朝东运行的飞机2。假设就这样匀速运动着,显然,$t_0$时间后,将会发生惨剧(飞机相撞)。

点击阅读全文...

18 Dec

黑洞融合的简单模拟

在天文爱好者眼中,黑洞是一个球体,其半径为$\frac{2GM}{c^2}$;这是广义相对论的施瓦兹黑洞的结果,也从经典力学推导推导出来,虽然用经典力学是错误的,但是对于多数的天文爱好者(包括笔者)来说,这是目前唯一的一种可行的理解方法(广义相对论那些复杂推导会让我们很崩溃的)。当然,事实上,黑洞不是一个球体,它只是一个密度很大的点。至于密度有多大,目前公认的说法是无穷大,但是严格的物理是不接受这个说法的,或者说,物理是不会接受任何无穷大的说法,所以现在积极发展量子引力理论来统一相对论和量子力学,不过这是另话了。$\frac{2GM}{c^2}$只不过是黑洞的视界,视界之内,我们就什么也不知道了。本文主要就从经典力学的角度探讨一下两个黑洞的合并过程中其视界的变化。读者将会发现,这些视界的形状相当有趣。

经典力学中的黑洞是这样定义的:天体表面的逃逸速度超过了光速,于是连光都无法逃脱,所以这个“洞”就很黑。也就是说,光子的总能量(引力势能与动能之和,经典力学意义下的)要为负,负数表示受到束缚。用数学公式来讲,就是:

$$\frac{1}{2}mc^2 - \frac{GM_1 m}{r_1}-\frac{GM_2 m}{r_2}-...-\frac{GM_n m}{r_n} \leq 0$$

点击阅读全文...