轻微的扰动——摄动法简介(1)
By 苏剑林 | 2013-01-16 | 48248位读者 | 引用为了计算实际问题,我们总会采用各种各样的理想模型。一般而言,一个模型越接近实际现象,它往往会越复杂。而忽略掉多数微小的干扰,只保留一些主要的项,这通常可以得到一个相当简单、能够精确解出的模型。以这样的一个可以精确解出的近似模型为基础,逐渐地把微小项的影响添加进去,使得我们的答案越来越准确,这就是摄动法的思想,也称作“微扰理论”。这种方法源于求解天体力学的N体问题,而现在已经发展成为一门相当系统的学科,并应用到了相对多的领域,如量子力学、电子理论等。
其实不难发现,实际问题中存在不少这样的例子,即当我们要计算某个现象时,先考虑最突出的,然后再考虑细节。比如说,要计算地球的轨道,先把它看成一个与太阳组成的纯粹的二体系统,然后把各种微小效应加进去,比如月球的影响、各大行星的影响甚至由于地球的不规则形状所产生的影响等。当然,不仅仅是这一类复杂的“大问题”,我们平常可能会遇到的一些“小问题”有可能也让摄动法派上用场。本文试图将摄动法介绍给各位读者。
摄动法的主要步骤是先忽略微小影响(令小参数为0),求出精确解;然后把所要求的解表达为关于小参数的幂级数。这个方法可以用于解答代数方程、微分方程等等各种领域。下面先以一个简单的代数方程来说明:
一、求解方程:$\varepsilon x^3+x^2=p^2$
网友:椭圆定长弦中点轨迹的一种解法
By 苏剑林 | 2013-02-02 | 34789位读者 | 引用大概在半年前,我曾用“化圆法”解决了椭圆内定长弦中点轨迹问题,求出了轨迹方程。前几天,我收到了网名为“理想”的网友的Email,他提出了自己对这个问题的解法,并得到了形式不同的轨迹方程,因此对两者的等价性表示疑惑。经过检验,我跟他的轨迹方程基本上是等价的,不过,他求出的轨迹方程总包括了原点,这是一点不足之处。但是看起来,他的轨迹方程却感觉好看一些。这的确很让人意外,因为从他的化简过程来看,有种“化简为繁”的味道,却得出了相当简洁的答案,着实有趣。
经过网友的同意,将他的过程贴在这里与大家分享!后面附有pdf文档,欢迎下载阅读。希望在科学空间可以看到更多的读者留下的痕迹。
椭圆定长弦中点轨迹的一种解法
作者:理想
本文介绍了一种计算椭圆定长弦中点轨迹的方法。设椭圆长、短轴分别为$2a$、$2b$,弦长为$2r$,随着弦的两端在椭圆上滑动,弦的中点形成的轨迹为:
$$(\frac{x^2}{a^2} + \frac{y^2}{b^2} - 1)(\frac{x^2}{a^4} + \frac{y^2}{b^4} + \frac{r^2}{a^2b^2}) + \frac{r^2}{a^2b^2} = 0$$
它不是一个椭圆,而是一个高次曲线。
[问题解答]双曲线上的最短距离
By 苏剑林 | 2013-02-04 | 26998位读者 | 引用本文我们来探讨下列积分的极值曲线:
$$S=\int f(x,y)\sqrt{dx^2+dy^2}=\int f(x,y)ds$$
这本质上也是一个短程线问题。但是它形式比较简答,物理含义也更加明显。比如,如果$f(x,y)$是势函数的话,那么这就是一个求势能最小的二维问题;如果$f(x,y)$是摩擦力函数,那么这就是寻找摩擦力最小的路径问题。不管是哪一种,该问题都有相当的实用价值。下面将其变分:
$$\begin{aligned} \delta S =&\int \delta[f(x,y)\sqrt{dx^2+dy^2}] \\ =&\int [ds\delta f(x,y)+f(x,y)\frac{\delta (dx^2+dy^2)}{2ds}]\\ =&\int ds(\frac{\partial f}{\partial x}\delta x+\frac{\partial}{\partial y}\delta y)+f \frac{dx d(\delta x)+dy d(\delta y)}{ds} \\=&\int ds(\frac{\partial f}{\partial x}\delta x+\frac{\partial}{\partial y}\delta y)+f \frac{dx}{ds} d(\delta x)+\frac{dy}{ds} d(\delta y) \end{aligned}$$
力学系统及其对偶性(二)
By 苏剑林 | 2013-11-14 | 18359位读者 | 引用如果仅仅从牛顿第二定律的角度来进行变换推导,那么关于力学定律的对偶性的结果无疑仅仅是初等的。对于理论分析来说,更方便的是从做小作用量原理的形式出发,事实上,这种形式计算量也是很少的,甚至比直接代入运动方程变换更加便捷。
上一篇文章中我们讲到,变换$z \mapsto z^2$将一个原点为几何中心的椭圆映射为一个原点为焦点的椭圆,并且相信这种变换可以将胡克定律跟牛顿万有引力定律联系起来。然后就立即给出了变换$w=z^2,d\tau=|z^2|dt$。但是这个变换本身并不显然的,假如我们仅仅发现了$z \mapsto z^2$的几何意义,如何相应地得出$d\tau=|z^2|dt$这个变换呢?本文初步地解决这个问题。
几何作用量
让我们回顾力学的最小作用量原理:
$$ S = \int_{{t_1}}^{{t_2}} L dt = \int_{{t_1}}^{{t_2}} {(T - U)} dt $$
最近的我的主要学习是在研究路径积分,在推导路径积分的一种新的变换方法(或者是一个新的视角吧),但是有道坎还是迈不过去,因此blog中也一直更新寥寥。说到积分与微分,这两个本是互逆的东西,但是在复数的统一之下,它们两个去可以相互转化。比如说,薛定谔方程是量子力学的微分形式,而路径积分实际上可以说是量子力学的积分形式,这让我有些想法,是不是任何微分形式的数学都存在一个积分形式的版本呢?如果是,是微分版本优还是积分版本优?
在数学分析中,我们会感觉到求导会比求积分容易很多,求导有现成的公式等等。但是微分有个最大的缺点,它是多分量的,比如,势函数是一个标量,但是微分(求梯度)之后就变成了三分量的矢量(即作用力),多分量事实上是不好处理了,为了处理这类问题,又引入了大量的算符。积分的特点在于它的标量性,也许计算很复杂,但是思想确实容易把握的,我更喜欢积分形式的理论(比如作用量原理、路径积分等。)
说到数学分析中常见而又著名的定积分,不得不提到以下三角函数积分了。
$$\int_0^{\pi/2} \sin^{2n} \theta d\theta$$
不难证明,它也等于
$$\int_0^{\pi/2} \cos^{2n} \theta d\theta$$
一维弹簧的运动(下)
By 苏剑林 | 2014-03-13 | 27111位读者 | 引用在上一篇文章中,我们得到了一维弹簧运动的方程
$$m\frac{\partial^2 X}{\partial t^2}=k\frac{\partial^2 X}{\partial \xi^2}$$
并且得到了通解
$$X=F(u)+H(v)=F(\xi+\beta t)+H(\xi-\beta t)$$
或者
$$X(\xi,t)=\frac{1}{2}\left[X_0(\xi+\beta t)+X_0(\xi-\beta t)\right]+\frac{1}{2\beta}\int_{\xi-\beta t}^{\xi+\beta t} X_1 (s)ds$$
在文章的末尾,提到过这个解是有些问题的。现在让我们来详细分析它。
Project Euler 454 :五天攻下“擂台”
By 苏剑林 | 2014-06-27 | 28959位读者 | 引用进入期末了,很多同学都开始复习了,这学期我选的几门课到现在还不是很熟悉,本想也在趁着这段时间好好看看。偏生五天前我在浏览数学研发论坛的编程擂台时看到了这样的一道题目:
设对于给定的$L$,方程
$$\frac{1}{x}+\frac{1}{y}=\frac{1}{n}$$
满足$0 < x < y \leq L$的正整数解共有$f(L)$种情况。比如$f(6)=1,f(12)=3,f(1000)=1069$,求$f(10^{12})$。
这道题目的来源是Project Euler的第454题:Diophantine reciprocals III(丢潘图倒数方程),题目简短易懂,但又不失深度,正符合我对理想题目的定义。而且最近在学习Python学习得不亦乐乎,看到这道题目就跃跃欲试。于是乎,我的五天时间就没有了,而且过程中几乎耗尽了我现在懂的所有编程技巧。由于不断地测试运行,我的电脑发热量比平时大了几倍,真是辛苦了我的电脑。最后的代码,自我感觉已经是我目前写的最精彩的代码了。在此与大家共享和共勉~
上述表达式是分式,不利于编程,由于$n=\frac{xy}{x+y}$,于是上述题目也等价于求$(x+y)|xy$(意思是$x+y$整除$xy$)的整数解。
最近评论