6
Feb
轻微的扰动——摄动法简介(2)
By 苏剑林 | 2013-02-06 | 34126位读者 | 引用为了让大家更加熟悉摄动法的基本步骤，本文再讲一个用摄动法解代数方程的例子。这是从实际研究中出来的：
$$\begin{eqnarray*} x=\frac{k(1+k^2+k^4+l^2)}{2(1+k^2)^2} \\ k=\frac{dy}{dx}\end{eqnarray*} $$
这是一道微分方程。要求解这道方程，最好的方法当然是先从第一式解出$k=k(x)$的形式然后再积分。但是由于五次方程没有一般的显式解，所以迫使我们要考虑近似解。当然，一般来说熟悉mathematica的人都会直接数值计算了。我这里只考虑摄动法。
我们将原方程变为下面的形式：
$$x=\frac{k}{2}[1+\frac{l^2}{(1+k^2)^2}]$$
21
Feb
[问题解答]有多少位数字？
By 苏剑林 | 2013-02-21 | 13972位读者 | 引用解决完上一题《有多少个5？》后，子瑞表示看到一道类似的题目，当然，这道题比上一道难一些：
一个数，各个数字加起来等于900，乘以2后各个数字加起来还是等于900，已知这个数字只有3、4、5、6组成，请问满足条件的最大数与最小数的积有多少位数？
要解答这个问题，我们只需要知道最大数和最小数分别有多少位即可。因为最大数必然是6...3的形式，而最小数只能是3...6的形式，它们的位数之和就是所求的位数。
怎样比较两个数的大小呢？显然，在不同位数的数时，位数多的数要大，同样位数才从高到低逐位比较。因此，我们应当考虑位数的最大与最小。
7
Mar
高斯型积分的微扰展开（二）
By 苏剑林 | 2015-03-07 | 20760位读者 | 引用为什么第二篇姗姗来迟？
其实要写这系列之前，我已经构思好了接下来几篇的内容，本来想要自信地介绍自己想到的一些积分展开的技巧；而且摄动法我本身就比较熟悉，所以正常来说不会这么迟才有第二篇。然而，在我写完第一篇，准备写第二篇的期间，我看到了知乎上的这篇回复：
http://www.zhihu.com/question/24735673
这篇文章大大地拓展了我对级数的认识。里边谈及到了积分的展开是一个渐近级数。这让我犹豫了，怀疑这系列有没有价值，因为渐近级数意味着不管怎样的展开技巧，得到的级数收敛半径都是0。
后来再想想，就算是渐近级数，也有改进的空间，有加速收敛的方法，所以我想我这几篇文章，应该还有一点点意义吧，还可以顺便介绍一下渐近级数和奇点的相关理论。嗯，就这么办吧。
曾经我会一字不差地看完你的日志，
一点蛋疼的破事都会当成宝贝一样。
和你分享，
跟你在一起,
笑点低的莫名其妙。
你知道我所有的事，
我也收藏着你太多的秘密。
我们可以一直聊到凌晨，
好像从来不缺话题。
可是...
可是...
后来，我们慢慢失去了联系。
等我们发现
时间是贼了，
它早己把我们
说不完的话
偷光了。
偶尔遇见，
也只能尴尬一笑，
寒暄几句，
便再无联络。
你一定以为无情的我把过去都忘记了，
你以为我把你看得不再重要。
那么，你肯定不知道，
我常梦见我们一起仰望过的那片天空呢。
亲爱的老朋友，
和亲爱的曾经心心相印的人。
不联系不是因为你不重要，
而是我好怕，
我不再重要。
13
Mar
单摆运动级数解：初试同伦分析
By 苏剑林 | 2013-03-13 | 18469位读者 | 引用开始之初，我偶然在图书馆看到了一本名为《超越摄动：同伦分析方法导论》，里边介绍了一种求微分方程近似解的新方法，关键是里边的内容看起来并不是十分难懂，因此我饶有兴致地借来研究了。果然，这是一种非常有趣的方法，在某种意义上来说，还是非常简洁的方法。这解决了我一直以来想要研究的问题：用傅里叶级数来近似描述单摆运动的近似解。当然，它带给我的冲击不仅仅是这些。为了得出周期解，我又同时研究了各种摄动方法的技巧，如消除长期项的PL(Poincaré–Lindstedt)方法。这同时增加了我对各种近似解析方法的了解。从开学到现在快三周的时间，我一直都在研究这些问题。
这已经是去年写的稿件了，刊登在今年二月份的《天文爱好者》上，本文的标题还登载了该期天爱的封面上，当时甚是高兴呢！在此与大家分享、共勉。
相信许多天文爱好者都知道第一、第二、第三宇宙速度的概念，也会有不少的天爱自己动手计算过它们。我们道，只要发射速度达到7.9km/s，宇宙飞船就可以绕地球运行了；超过11.2km/s，就可以抛开地球，成为太阳系的一颗“人造行星”；再大一点，超过16.7km/s，那么就连太阳也甩掉了，直奔深空。
16.7km/s，咋看上去并不大，因为地球绕太阳运行的速度已经是30km/s了，这个速度在宇宙中实在是太普通了。但是对于我们目前的技术来说，它大得有点可怕。维基百科上的资料显示，史上最强劲的火箭土星五号在运送阿波罗11号到月球时，飞船最终也只能加速到接近逃逸速度，即11.2km/s，而事实上第三宇宙速度已经是是目前人造飞行器的速度极限了。可是没有速度，我们就不能发射探测器去探索深空，那些科幻小说中的“星际移民”，就永远只能停留在小说上了。
很早以前我就对这个问题感兴趣了，但是一直搁置着，没有怎么研究。最近在阅读《引力与时空》的“潮汐力”那一节时重新回到了这个问题上，决定写点什么东西。在这里不深究流体静力平衡的定义，顾名思义地理解，它就是流体在某个特定的力场下所达到的平衡状态。流体静力学告诉我们：
达到流体静力平衡时，流体的面必定是一个等势面。
这是为什么呢？我们从数学的角度来简单分析一下：只考虑二维情况，假如等势面方程是$U(x,y)=C$，那么两边微分就有
$$0=dU=\frac{\partial U}{\partial x}dx+\frac{\partial U}{\partial y}dy=(\frac{\partial U}{\partial x},\frac{\partial U}{\partial y})\cdot (dx,dy)$$
这意味着向量$(\frac{\partial U}{\partial x},\frac{\partial U}{\partial y})$和向量$(dx,dy)$是垂直的，前者便是力的函数，后者就是一个切向量（三维就是一个切平面）。也就是说合外力必然和流体面垂直，这样才能提供一个相等的方向相反的内力让整个结构体系处于平衡状态！
14
Apr
费曼积分法(8)：求高斯积分
By 苏剑林 | 2013-04-14 | 48250位读者 | 引用自从了解了费曼积分法之后，我就一直想着用费曼积分法来求高斯积分$\int_0^{\infty} e^{-x^2}dx=\frac{\sqrt{\pi}}{2}$这个神奇的积分，但一直无果。在《数学桥》里边，作者是通过将其转变为二重积分来解决的，简洁而巧妙。但是为了显示费曼积分法的威力，我一直想找到高斯积分的其他求法。上星期在《数学物理方法》中看到作者用拉普拉斯变换求出了该积分，眼睛不禁为之一亮，不过这属于积分变换内容，属于“积分符号内取积分”的技巧，在此不作讨论。今天在网上查找资料时，在“赵洁”的一篇论文《含参变量积分》中，看到了一种属于费曼积分法范畴内的方法，特与大家分享。
从“事后分析”来看，高斯积分的结果涉及到了$\sqrt{\pi}$这个量，一般来说我们常见的公式出现$\pi$的不少，可是几乎没有出现$\sqrt{\pi}$的，所以一般来说我们都将它平方。我们引入
$$f(x)=(\int_0^x e^{-t^2}dt)^2$$
最近评论