實驗五行星的軌道和位置

實驗五行星的軌道和位置第一頁，共三十二頁，2022年，8月28日

數(shù)學(xué)實驗行星的軌道 和位置第二頁，共三十二頁，2022年，8月28日樂經(jīng)良

他以幾乎神一般的思維力,最先說明了行星的運動和圖象,慧星的軌道和大海的潮汐. ─Newton墓志銘第三頁，共三十二頁，2022年，8月28日樂經(jīng)良

哥白尼(波蘭，1473-1543)日心說

地球－我們的家園

46億歲,赤道半徑6378.14公里，比極半徑長21公里

金星－看起來最亮的行星半徑約為6073公里,表面溫度高達(dá)465至485度,自轉(zhuǎn)方向與其它行星相反

背景介紹16世紀(jì)前，人們認(rèn)為太陽只有6大行星托勒密（古希臘）地心說第四頁，共三十二頁，2022年，8月28日樂經(jīng)良

水星－距太陽最近的行星半徑為2440公里,較小，難以觀察

火星－離地球最近、人們最關(guān)注的行星：火星上有無生命？

土星－最美麗的行星衛(wèi)星數(shù)目最多,23顆.光環(huán)由無數(shù)塊冰狀物組成的

木星－行星中的巨無霸赤道半徑約為71400公里，是地球的11.2倍第五頁，共三十二頁，2022年，8月28日樂經(jīng)良

行星運行三大規(guī)律

開普勒(1571－1630)(觀察分析數(shù)據(jù))

在第谷·布拉赫(1546-1601)的基礎(chǔ)上提出1.行星運行的軌道是以太陽為一個焦點的橢圓;2.從太陽指向某一行星的線段在單位時間內(nèi) 掃過的面積相同;3.行星運行周期的平方與其運行軌道橢圓長軸 的立方之比值是不隨行星而改變的常數(shù).第六頁，共三十二頁，2022年，8月28日樂經(jīng)良(Adams1845,Leverrier1846)

萬有引力定律

“自然哲學(xué)的數(shù)學(xué)原理”(1687牛頓)天王星－樂師(Herschel)發(fā)現(xiàn)的行星(1781)海王星－筆尖上的行星冥王星－離太陽最遠(yuǎn)、未知數(shù)最多的行星

2006年8月24日國際天文學(xué)聯(lián)合大會決定： 冥王星降級為“矮行星”(大行星的定義)

太陽只有八大行星！第七頁，共三十二頁，2022年，8月28日樂經(jīng)良

二實驗?zāi)康?/p>

本實驗主要涉及常微分方程。通過實驗復(fù)習(xí)：微分方程的建模和解法，數(shù)值積分的計算。另外，還將介紹：建立數(shù)學(xué)模型時復(fù)坐標(biāo)系的選取，基于壓縮映像的求根方法，微分方程的Runge-Kutta法第八頁，共三十二頁，2022年，8月28日樂經(jīng)良

三實際問題

地球距太陽最遠(yuǎn)處(遠(yuǎn)日點)距離為

1.521×1011m,此時地球繞太陽運動(公轉(zhuǎn))的速度為2.929×104m/s,試求: 1)地球距太陽的最近距離

2)地球繞太陽運轉(zhuǎn)的周期

3)在從遠(yuǎn)日點開始的第100天結(jié)束時,

地球的位置與速度第九頁，共三十二頁，2022年，8月28日樂經(jīng)良加速度

數(shù)學(xué)模型

在運動學(xué)中常采用復(fù)坐標(biāo)系 設(shè)太陽中心所在位置為復(fù)平面之原點,在時刻t,

行星位于以下復(fù)數(shù)代表的點

速度為

四數(shù)學(xué)模型第十頁，共三十二頁，2022年，8月28日根據(jù)Newton第二定律比較虛實部導(dǎo)出微分方程組四數(shù)學(xué)模型第十一頁，共三十二頁，2022年，8月28日問題（5.4）~（5.9）就是行星繞太陽運行的軌跡的數(shù)學(xué)模型。四數(shù)學(xué)模型第十二頁，共三十二頁，2022年，8月28日四數(shù)學(xué)模型第十三頁，共三十二頁，2022年，8月28日于是我們得到了行星運動的形式較為簡單的數(shù)學(xué)模型：四數(shù)學(xué)模型第十四頁，共三十二頁，2022年，8月28日rdC1tp樂經(jīng)良

I行星的軌跡請嘗試推導(dǎo)出行星的軌道方程？

p1ecosr2(p,e是常數(shù),根據(jù)相關(guān)已知數(shù)據(jù)導(dǎo)出)

改寫前面積分表達(dá)式成為210dT1C1 2

1(1-ecos)給出時間T1，要求位置即求出θ1與r，較難！五解析方法第十五頁，共三十二頁，2022年，8月28日五解析方法第十六頁，共三十二頁，2022年，8月28日II行星的周期五解析方法第十七頁，共三十二頁，2022年，8月28日III行星的位置五解析方法第十八頁，共三十二頁，2022年，8月28日hnthtrjsdv0.05331.65001.49401.99592.98190.011681.68001.49331.99792.98340.0053371.685001.49321.99822.98360.00116861.686001.49311.99832.9837

六數(shù)值方法第十九頁，共三十二頁，2022年，8月28日

六數(shù)值方法第二十頁，共三十二頁，2022年，8月28日

六數(shù)值方法第二十一頁，共三十二頁，2022年，8月28日

六數(shù)值方法第二十二頁，共三十二頁，2022年，8月28日

六數(shù)值方法第二十三頁，共三十二頁，2022年，8月28日

五數(shù)值方法

六數(shù)值方法第二十四頁，共三十二頁，2022年，8月28日樂經(jīng)良艾薩克?煉金術(shù)士造幣廠總監(jiān)較之科學(xué)他更多致力于《圣經(jīng)》的研究

牛頓SirIsaacNewton(英格蘭1643－1727年)

科學(xué)史上最有影響力的人 物理學(xué)家數(shù)學(xué)家天文學(xué)家 哲學(xué)家

專心于科學(xué)研究到癡情第二十五頁，共三十二頁，2022年，8月28日樂經(jīng)良性格內(nèi)向獨身一生牛頓的一句名言

WhatDescartesdidwasagoodstepyouhave addedmuchseveralways&especiallyintaking thecoloursofthethinplatesintophilosophcal consideration IfIhaveseenfurtheritisbystandingonyeshouldersofGiants.第二十六頁，共三十二頁，2022年，8月28日hh樂經(jīng)良III微分方程的Runge-kutte方法以一元為例dx dtf(t,x),xt0x0設(shè)步長為h,則tkkh,記xkx(tk)

xk1x(tkh)（Taylor展開）

33!

22!x(tk)x(tk)x(tk)hx(tk)其中x(tk)f(tk,xk)x(tk)ft(tk,xk)fx(tk,xk)x(tk)

ft(tk,xk)fx(tk,xk)f(tk,xk)可以求出各階導(dǎo)數(shù)第二十七頁，共三十二頁，2022年，8月28日xkh(K1樂經(jīng)良

計算高階導(dǎo)數(shù)代之以f在一些點的值的組合當(dāng)Taylor展開到四階項，可取16262616K4)K2K3xk1其中K1f(tk,xk)K2f(tkh/2,xkhK1/2)K3f(tkh/2,xkhK2/2)K4f(tkh,xkhK3)Runge-Kutte迭代格式

在Matlab可 以直接調(diào)用第二十八頁，共三十二頁，2022年，8月28日樂經(jīng)良

使用Matlab

先定義一階微分方程組函數(shù)組functiondy=m5_2_fun(t,y)C1=4.455e15;MG=1.989e30*6.672e-11;dy=zeros(3,1);dy(1)=C1^2/y(2)^3-MG/y(2)^2;dy(2)=y(1);dy(3)=C1/y(2)^2;第二十九頁，共三十二頁，2022年，8月28日樂經(jīng)良

再調(diào)用Runge‐kutte方法專用程序：functionT=m5_2(h)[t,y]=ode45(@m5_2_fun,[0:h:400*24*3600],[0,1.521e11,0]);n=max(find(y(:,3)<2*pi));%查找小于

2pi所對應(yīng)的最大n值T=t(n);r=y(round(n/2),2);polar(y(:,3),y(:,2))第三十頁，共三十二頁，2022年，8月28日樂經(jīng)良

軌道圖形將計算所得的數(shù)據(jù)，利用Matlab作圖可得第三十一頁，共