




下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、地球衛(wèi)星三維運(yùn)行軌道MATLAB 仿真1、問(wèn)題的描述軌道上運(yùn)行的地球衛(wèi)星,根據(jù)牛頓第二定律F=ma 以及萬(wàn)有引力定律3E F=-GmM*r/r ,可得3E a=-GM*r/r ,即 3E 3E 3E x''= -GM*x/ry''= -GM*y/r; z''= -GM*z/r(1)式中,(x ,y ,z )表示衛(wèi)星的三維坐標(biāo),為-1122G=6.672*10(Nm /kg 引力常數(shù),24E M =5.97*10(kg是地球的質(zhì)量。假定衛(wèi)星的三個(gè)方向的初始位置和速度如下x(0,y(0,z(0,x'(0,y'(0,z'(0=2
2、043922.17 8186504.63 4343461.71 -5379.54 -407.10 3516.05。該衛(wèi)星軌道求解過(guò)程實(shí)際上是求解一個(gè)二階常微分方程,可首先將該方程轉(zhuǎn)換為一階常微分方程,令TX=x,y,z,x',y',z',故公式(1)可轉(zhuǎn)化為 3E (4(5(6X'( =-GM/r(1(2(3X X X t A A X A X A X * (2)初始條件即為X(0=2043922.17 8186504.63 4343461.71 -5379.54 -407.10 3516.05。2、MATLAB 仿真代碼分兩段程序:(1)子程序?qū)⒍A微分方程轉(zhuǎn)換
3、為一階微分方程,代碼如下function fy=vdp(t,xr=x(12+x(22+x(32;G=3.986005e14;A=-G/r(3/2;fy=x(4x(5x(6A*x(1A*x(2A*x(3;End(2)主程序如下,注意:為更好地查看衛(wèi)星軌道與地球的相對(duì)位置關(guān)系,此處將地球模型圖的繪制代碼一并給出clear allclose allclcy0=2043922.166765 8186504.631471 4343461.714791 -5379.544693 -407.095342 3516.052656;t,result=ode45(vdp,0:1:9000,y0;x=result(:,1;y=result(:,2;z=result(:,3;X,Y,Z=sphere(200;RE=0.64e7;X=RE*X;Y=RE*Y;Z=RE*Z;figure(1hold ongrid onmesh(X,Y,Z%繪制地球plot3(x,y,z%繪制衛(wèi)星軌道仿真結(jié)果如下(給出兩張圖): -10-8-6-
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 倉(cāng)儲(chǔ)合同里免責(zé)合同范本
- 產(chǎn)品保本合同范本
- 小學(xué)道法教學(xué)中翻轉(zhuǎn)課堂的應(yīng)用
- 包包寄售合同范本
- 二手車市場(chǎng)發(fā)票合同范本
- 醫(yī)療器械 廣告合同范本
- 《JJG196-2006-常用玻璃量器檢定規(guī)程》
- 加工合同范本包括些費(fèi)用
- 化糞池清淤合同范本
- 借貸股權(quán)質(zhì)押合同范本
- 水環(huán)境綜合治理服務(wù)方案(技術(shù)標(biāo))
- 【原創(chuàng)】頭腦特工隊(duì)開(kāi)的那些心理學(xué)腦洞
- 美甲藝術(shù)全套教學(xué)課件
- 高等數(shù)學(xué)上冊(cè)目錄同濟(jì)第七版
- 中國(guó)古代餐具
- 電動(dòng)執(zhí)行機(jī)構(gòu)安裝施工工藝標(biāo)準(zhǔn)
- 施工日志模板
- 粗原料氣的凈化-二氧化碳的脫除(合成氨生產(chǎn))
- Agilent7820A氣相色譜儀操作規(guī)程知識(shí)講解
- 中醫(yī)適宜技術(shù)模擬試題(附答案)
- 加涅的信息加工理論-課件
評(píng)論
0/150
提交評(píng)論