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

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進行舉報或認領(lǐng)

文檔簡介

1、地球衛(wèi)星三維運行軌道MATLAB 仿真1、問題的描述軌道上運行的地球衛(wèi)星，根據(jù)牛頓第二定律F=ma 以及萬有引力定律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)星的三個方向的初始位置和速度如下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)星軌道求解過程實際上是求解一個二階常微分方程，可首先將該方程轉(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)星軌道與地球的相對位置關(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. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論