地球衛(wèi)星三維運(yùn)行軌道MATLAB仿真_第1頁(yè)
地球衛(wèi)星三維運(yùn)行軌道MATLAB仿真_第2頁(yè)
地球衛(wèi)星三維運(yùn)行軌道MATLAB仿真_第3頁(yè)
地球衛(wèi)星三維運(yùn)行軌道MATLAB仿真_第4頁(yè)
地球衛(wèi)星三維運(yùn)行軌道MATLAB仿真_第5頁(yè)
全文預(yù)覽已結(jié)束

下載本文檔

版權(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ì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論