# 19.4 行星运动
> 原文: [http://math.mit.edu/~djk/calculus_beginners/chapter19/section04.html](http://math.mit.edu/~djk/calculus_beginners/chapter19/section04.html)
行星和太阳之间的引力相互作用由反平方中心力定律描述。
为方便起见,我们将太阳置于坐标的原点,并在![](https://img.kancloud.cn/ac/9f/ac9f945b8c4c3cae0b8fa394c83a71e5_55x18.gif)点开始我们的行星,其位置的初始一阶导数由![](https://img.kancloud.cn/66/f8/66f8d848224b7b696a14fc7346be64cb_54x18.gif)给出。
我们将假设行星比太阳轻得多(因为地球与太阳相比),太阳不会移动。 (实际上,行星运动中固定的是系统的质心。木星和土星足够大,当它们在我们天空的同一部分时,所有行星的质心都不在太阳内,这样太阳就会移动,但不是很多。)
利用这些坐标和这个假设![](https://img.kancloud.cn/d7/72/d7726585e7eb772c859b54b5445030df_12x14.gif)的运动方程,行星的位置矢量服从方程
![](https://img.kancloud.cn/1a/b1/1ab1f0a0328db58e75a2b77b2fa789ad_83x39.gif)
由于行星上的力指向太阳,我们在![](https://img.kancloud.cn/e3/52/e3525a37737641aedf64a87c4d730002_40x18.gif)平面上开始行星,我们的![](https://img.kancloud.cn/f5/cc/f5cc9e47883d3f86d5795d82e36323c4_9x8.gif)坐标将始终为![](https://img.kancloud.cn/fb/02/fb0294de83ba79f08ad4634c597b108c_9x12.gif),我们可以忽略它。
这是一个二阶微分方程,有两个因变量,![](https://img.kancloud.cn/c1/1d/c11d457f49fddd6913fa6a0ec01e3fa8_29x18.gif)和![](https://img.kancloud.cn/e4/30/e430e5ba9b0a679a55c1670ad4204a03_29x18.gif)。我们可以在一个电子表格上设置它,每个列都有![](https://img.kancloud.cn/77/b7/77b767c523b51de233c43edf9f7b9582_42x16.gif)和![](https://img.kancloud.cn/77/90/7790dd0efb4a03a4c876741804d9b559_10x8.gif)和![](https://img.kancloud.cn/6c/70/6c704047d3148fd7a8b563aaf79dd7f4_9x12.gif)的导数。在坐标方面,运动方程是
![](https://img.kancloud.cn/d2/cf/d2cf49d287770539322e16dae01d7f07_228x41.gif)
由于![](https://img.kancloud.cn/5e/b9/5eb9d72c1ec656dd6de588089727a8e6_8x8.gif)出现在这两个方程中且![](https://img.kancloud.cn/5e/b9/5eb9d72c1ec656dd6de588089727a8e6_8x8.gif)为![](https://img.kancloud.cn/ec/f9/ecf9fd5a0ed14ce4911c1475c8b27277_77x23.gif),因此也可以方便地将列专用于![](https://img.kancloud.cn/5e/b9/5eb9d72c1ec656dd6de588089727a8e6_8x8.gif)。设置![](https://img.kancloud.cn/52/65/52650bcac8b74b381d8d20c612d56c49_41x13.gif)定义![](https://img.kancloud.cn/5e/b9/5eb9d72c1ec656dd6de588089727a8e6_8x8.gif)的比例,但不定义![](https://img.kancloud.cn/64/9c/649cbb0dc497a73d3cb2a543eaaf7933_6x12.gif)。这意味着我们可以选择我们的时间单位,以便![](https://img.kancloud.cn/b5/46/b546a89f36ed0a48be4457149ac399aa_8x8.gif)为![](https://img.kancloud.cn/55/95/559537f1e11c68d8ba3d9f6d540de6b0_7x13.gif)。
有了这个选择,我们可以按如下方式设置电子表格:
我们将时间变量![](https://img.kancloud.cn/64/9c/649cbb0dc497a73d3cb2a543eaaf7933_6x12.gif)放在 A 列中,并从第 7 行开始,A7 设置为 0.我们必须为![](https://img.kancloud.cn/da/b8/dab8f201ba10fb5a14e991ab157a9c7c_10x13.gif)选择一个增量,您可以确定最喜欢的那个。它必须足够小,以便![](https://img.kancloud.cn/dd/9c/dd9c7cdbcbb952f01fcc4ea284d83799_21x39.gif)很小,但足够大,你可以绘制轨道。您可以从![](https://img.kancloud.cn/3c/88/3c88b61b3e8aee20cc2f1be8f22df58f_68x17.gif)开始,如果不能正常工作则更改它。我们可以将字母 d 放在 A2 中,将其值放在 B2 中。我们需要指定的其他参数是![](https://img.kancloud.cn/77/90/7790dd0efb4a03a4c876741804d9b559_10x8.gif)和![](https://img.kancloud.cn/6c/70/6c704047d3148fd7a8b563aaf79dd7f4_9x12.gif)的导数的初始值。因此,在 A3 中输入“初始 x 速度”,在 B3 中输入其值(比如 0),在 A4 中输入“初始 y 速度”,在 B4 中输入其值(比如 1)。
把 t 放在 A6 中,x 放在 B6 中,y 放在 C6 中,r 放在 D6 中,x'放在 E6 中,y'放在 F6 中。我们将![](https://img.kancloud.cn/77/90/7790dd0efb4a03a4c876741804d9b559_10x8.gif)和![](https://img.kancloud.cn/6c/70/6c704047d3148fd7a8b563aaf79dd7f4_9x12.gif)放在 B 和 C 列中,因此在 B7 中放置 1,在 C7 中放置 0。我们将![](https://img.kancloud.cn/5e/b9/5eb9d72c1ec656dd6de588089727a8e6_8x8.gif)放在 D 列中,将 D7 设置为=(B7 ^ 2 + C7 ^ 2)^ 0.5。我们将![](https://img.kancloud.cn/15/9c/159c6fe833282b916d1d000eff6881dd_20x38.gif)(称之为![](https://img.kancloud.cn/65/63/6563134763fd723bef4e4897d6012eb0_49x12.gif))放在 E 列中,将 E7 设置为= B3 并将![](https://img.kancloud.cn/d6/5d/d65defeb43a657e28ed81daa0bf869b9_19x38.gif)(称之为![](https://img.kancloud.cn/22/1d/221d532f58a41787e4fd7846bbc4fbad_49x12.gif))放入 F 列 F7 设置为= B4。
我们接下来将 A8 设置为= A7 + $ B $ 2
B8 到= B7 + $ B $ 2 * E7
C8 到= C7 + B $ 2 * F7(你可以将 B8 复制到 C8)
将 D7 复制到 D8
将 E8 设为= E7- $ B $ 2 * B7 / $ D7 ^ 3
将 E8 复制到 F8
现在从列中复制 A8 到 F8
这将为您的参数值提供最粗略的解决方案。
**完成后,列 B 和 C 的![](https://img.kancloud.cn/61/41/6141ff672360e1d599330a2b6f77b1a5_27x12.gif)图将给出空间轨道。根据需要调整参数。**
<button aria-controls="planetary-motion-spreadsheet" aria-expanded="false" class="btn bg-light border-secondary" data-target="#planetary-motion-spreadsheet" data-toggle="collapse" id="toggle-spreadsheet-table" type="button">显示表</button>[](../download/planetary-motion.xlsx)
Number of steps<button aria-expanded="false" aria-haspopup="true" class="btn btn-sm bg-light border-secondary dropdown-toggle" data-toggle="dropdown" id="nbr-steps-btn" type="button" value="25">25</button>[10](#) [25](#) [50](#) [100](#)Number of digits after decimal point<button aria-expanded="false" aria-haspopup="true" class="btn btn-sm bg-light border-secondary dropdown-toggle" data-toggle="dropdown" id="nbr-digits-btn" type="button" value="10">10</button>[5](#) [10](#) [15](#)
**练习 19.5:设置它。 ![](https://img.kancloud.cn/15/9c/159c6fe833282b916d1d000eff6881dd_20x38.gif)和![](https://img.kancloud.cn/d6/5d/d65defeb43a657e28ed81daa0bf869b9_19x38.gif)的哪些值在这些坐标中给出圆周运动?**
在过去,以数字方式处理这些方程式是非常可怕的。相反,来自牛顿的物理学家通过引入数量来解决方程式,即能量和角动量,这些量不会随着这种运动而改变,并且通过推理而不是数值计算推导出轨道。
几个世纪以来,天文学家一直在仔细观察行星的实际行为,并在开普勒的三个定律中进行了清晰的总结,如下:
**1.受相同力量影响的行星和其他物体的运动位于“圆锥截面”的轨道上:椭圆或双曲线或非常特殊的抛物线(全部以太阳为焦点)或直线。**
**2.在任何轨道上每单位时间扫出的面积是不变的。**
3.椭圆轨道的周期与其半径的度量之间存在某种特定的关系,我们不再进一步讨论这种关系。
**最后注释:**最后几章包含许多未包含在任何正常单变量微积分课程中的材料。这些材料的目的是为了您的享受而不是恐吓您。问题在于,这里的 applet 和方法可以让你比定期的微积分课程更快地学习微积分。但是你学到和保留的东西很大程度上取决于你花多少时间去做。如果最终的结果是你花了很少的时间学习微积分,那对你来说就不好了。因此,您可能花费相同的时间,并了解更多!
- 第 0 章:为何学习微积分?
- 0.1 你应该知道什么
- 0.2 什么是微积分?我们为什么要研究它?
- 第 1 章:数字
- 1.1 什么是数字?有理数
- 1.2 小数和实数
- 1.3 复数
- 复数运算
- 1.4 可数集(消遣)
- 第 2 章:使用电子表格
- 2.1 什么是电子表格?
- 2.2 斐波纳契数
- 2.3 帕斯卡的三角形
- 2.4 与电子表格集成
- 第 3 章:线性函数
- 3.1 什么是函数?
- 3.2 线性函数
- 3.3 线性
- 第四章:函数的二次型和导数
- 4.1 更复杂的函数
- 4.2 二次函数的斜率
- 第 5 章:有理函数和导数的计算
- 5.1 有理函数的导数
- 第 6 章:指数函数,替换和链规则
- 6.1 最有用函数的导数
- 第 7 章:三角函数及其导数
- 7.1 二维数学
- 7.2 三角学和导数以及加法定理
- 第 8 章:反函数及其导函数
- 8.1 反函数
- 8.2 微分反函数
- 8.3 更多规则
- 第 9 章:数值微分和不可微函数
- 9.1 数值微分
- 9.2 绘制导数图
- 9.3 不可微函数
- 第 10 章:微分的回顾
- 10.1 复习
- 第 11 章:微分在求解方程中的应用
- 11.1 求解方程
- 第 12 章:反导数
- 12.1 反导数
- 第 13 章:曲线下面积;定积分
- 13.1 区域:定义,名称和符号
- 13.2 微积分和确定区域的基本定理
- 13.3 积分的诀窍
- 第 14 章:数值积分
- 14.1 数值积分计划
- 14.2 积分的“规则”
- 14.3 为什么这些规则有效?
- 第 15 章:平行数字的面积和体积;行列式
- 15.1 有符号面积和体积
- 15.2 表示平行边的图形
- 15.3 行列式的属性
- 15.4 求解行列式
- 15.5 用于求解电子表格中的行列式的爱丽丝梦游仙境方法
- 第 16 章一些纯数学
- 16.1 极限和点集拓扑简介
- 16.2 紧集
- 16.3 杂注
- 16.4 Lebesgue 积分
- 第 17 章:物理的建模应用
- 17.1 垂直运动建模
- 17.2 弹簧建模(谐波振荡器)
- 17.3 受迫振荡
- 17.4 简单电路
- 第 18 章捕食者猎物模型
- 18.1 捕食者猎物模型
- 第 19 章:求解微分方程
- 19.1 计划
- 19.2 一阶微分方程
- 19.3 二阶微分方程
- 19.4 行星运动