科学空间|Scientific Spaces

在上一些关于限制性三体问题的探讨中，我们得出了在平面上的方程：
$$\ddot{R}+2i\omega \dot{R}=\omega^2 R-GM\frac{R-l_1}{|R-l_1|^3}-Gm\frac{R-l_2}{|R-l_2|^3}\tag{32}$$
能量积分为：
$$\frac{1}{2}|\dot{R}|^2=\frac{1}{2} \omega^2 |R|^2+\frac{GM}{|R-l_1|}+\frac{Gm}{|R-l_2|}-C\tag{33}$$
下面就以这两个方程为基础，再说说限制性三体问题的那些事儿...

点击阅读全文...

非抛物面望远镜的校正镜方程求解
The Corrector Plate of Non-parabola Telescope

本文在牧夫天文论坛的讨论：
http://www.astronomy.ac/bbs/thread-160257-1-1.html

为了克服折射望远镜的色差问题，1670年，牛顿制造了第一台实用的反射式望远镜，将望远镜的主镜由玻璃透镜换成了抛物反射面，从而消除了色差。然而，相比球面镜，大口径的抛物面并不容易磨制。因为制作大球面镜只需要将曲率相等的小镜片相对自由组合在一起就行了，而抛物线每点的曲率并不相等，所以需要逐个磨制曲率不等的小镜片，并按照严格的顺序组合起来。这无疑大大增加了磨制难度。

Lamost是目前世界最大的施密特望远镜

为了解决这一难题，天文学家们想到了一个折衷的办法：以球面为主镜，并配以校正镜来校正球差。迎着这一思路，施密特望远镜随之而生。而当代的大望远镜基本上都是沿用这一思路。然而，校正镜是一个比抛物面更加复杂的四次曲面，磨制工艺要求更高，因此，校正镜也不宜过大。

点击阅读全文...

通常来说，选取惯性系为参考系，列出的三体问题方程为
$$\ddot{\vec{r}}_k=\sum_{i=1,i != k}^{n} Gm_i\frac{\vec{r}_i-\vec{r}_k}{|\vec{r}_i-\vec{r}_k|^3}$$

历史上出现过很多不同形式的变换，使得三体问题的运动方程有了各样的形式，如Lagrange形式、Jacobi形式、Hamilton形式等。这些变换形式都各有特点，都能够在一定程度上化简三体问题。BoJone在研究摆弄等质量型三体问题的运动方程时，也发现了一种很有趣的变换，在此贴出与大家分享。

设$\vec{R}_1=\vec{r}_1-\vec{r}_2,\vec{R}_2=\vec{r}_2-\vec{r}_3,\vec{R}_3=\vec{r}_3-\vec{r}_1$，则三体问题的运动方程变为

点击阅读全文...

拉格朗日

BoJone在之前的《自然极值》系列已经花了一定篇幅来讲述“极值”在自然界中是多么的普遍，它能够引导我们进行某些问题的思考，从而获得简单快捷的解答。接下来，我要说的一个更加令人惊讶的“事实”：“极值”不仅仅在某些数学或物理问题上给予我们创造性的思考，它甚至构建了整个经典力学乃至于整个物理学！这不是夸大其辞，这是物理学中被称为“最小作用量原理”的一个原理，很多物理学家（如费恩曼）被它深深吸引着，甚至认为它就是“上帝创造世界的终极公式”！（关于做小作用量原理，大家不妨看一下范翔所写的《最小作用量原理与物理之美》系列文章）

话说在18世纪，欧拉和拉格朗日开创了一条独特的道路，即用变分法来研究经典力学，从而使经典力学焕发出了新的活力，也由此衍生出了一个叫“理论力学”或“分析力学”的分支。用变分法研究力学有很多的好处，变分的对象一般都是标量函数，我们只需要写出动力系统的动能与势能表达式，就可以进行一系列的研究，比如列出质点的运动方程、判断平衡点的稳定性、求周期轨道等等（由于BoJone对理论力学研究还不够深入，无法举太多例子，但请相信，其作用远远不止这些），省去了不少繁琐的矢量性分析，这些都是在变分法发明前难以研究的。

点击阅读全文...

转载自：eaglefantasy.com

有感于Matrix67神牛的这篇文章（强烈建议大家去读一读），我也发表一下自己对于教材编写的一点看法。

1.对线性代数的吐槽

（没学过线性代数的同学请忽略下面3段往后接着看。）

我一直觉得线性代数用那种严格公理化的语言写成课本根本不适合初学者学习，一开始学习线性代数的时候，我本人对很多概念的直观意义根本就是完全不知道。我们的课本是丘维声的《简明线性代数》，我在此毫不掩饰的表示对这本教材的鄙视：这本教材居然是按照这样的顺序讲线性代数的：线性方程组->行列式->线性方程组的进一步讨论->矩阵的运算->一大堆东西->线性空间->线性映射->一大堆东西。这个狗屁顺序直接导致我前半个学期一直以为线性代数就是研究怎么解线性方程组的，我心想，这么简单的问题，具体问题谁都会解，值得这么大动干戈的定义出这么大堆东西么。。。一直到线性空间那一个章节以前，我完全就不知道线性代数整个是在干什么..后来学的多了我才知道，其实线性代数就是研究线性空间和线性映射的嘛，什么线性方程组，根本没那么重要。一个更加合理的顺序是：先讲线性空间、线性映射，其中明确说明矩阵就是线性映射，然后再讲行列式，然后线性方程组只作为一个例子出现就可以了。

点击阅读全文...

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

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

点击阅读全文...

看完了“双不动中心”问题，我们不妨再来看一个貌似简单一点的力学问题，在一个固定质点的引力吸引的基础上，增加一个恒力作用，研究这样的力场中小天体的运动。

咋看上去这个问题比“双不动中心”简单多了，至少运动方程也显得更简单：
$$\ddot{vec{r}}=-GM\frac{\vec{r}}{|\vec{r}|^3}+\vec{F}$$

其中$\vec{F}$是一个常向量。不过让人比较意外的是，这个问题本质上和“双不动中心”是一样的，它可以看作是双不动中心问题的一个极限情况。而且它们的解法也是惊人地相似，下面我们就来分析这一个过程。

首先很容易写出这个方程的能量守恒积分：
$$1/2 \dot{vec{r}}^2-GM\frac{1}{|\vec{r}|}-\vec{F}\cdot \vec{r}=h$$

点击阅读全文...

我们考虑一个球状的星团，并假设它是各向同性的，即距离球心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}$$

点击阅读全文...