利用計算機編程模擬質(zhì)點在平方反比引力場中的運動_第1頁
利用計算機編程模擬質(zhì)點在平方反比引力場中的運動_第2頁
利用計算機編程模擬質(zhì)點在平方反比引力場中的運動_第3頁
利用計算機編程模擬質(zhì)點在平方反比引力場中的運動_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

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

文檔簡介

1、利用計算機編程模擬質(zhì)點在平方反比引力場中的運動題目：研究質(zhì)點在平方反比引力場中的運動，例如行星繞太陽的運動。設(shè)質(zhì)量為m0的質(zhì)點位于力心且固定不動，質(zhì)量為m 的質(zhì)點在m0產(chǎn)生的引力場中運動，當(dāng)m 與m0相距r 時，質(zhì)點所受萬有引力為F= Gm m0/r2 ，G 為引力常量。要求：(1)當(dāng)質(zhì)點總能量大于、等于和小于零時，畫出質(zhì)點在平方反比引力場中的運動軌跡。(2)當(dāng)質(zhì)點總能量小于零且保持不變時，改變角動量的大小，畫出質(zhì)點相應(yīng)的運動軌跡。(3)學(xué)習(xí)將極坐標(biāo)變換成直角坐標(biāo)的方法以及利用對稱性畫曲線的方法。function guidaoclear; clc;global G MG=1; M=5000;

2、r0=100;E=-18,0,22;v0=sqrt(2*E+2*G*M/r0);figureaxis(-400,200,-1000,1000);xlabel(x); ylabel(y);title(不同E 時的不同曲線);s1=E0;hold onsun=line(0,0,color,r,marker,.,markersize,20);for k=1:3t u=ode45(guidaofun,0:0.05:120,r0,0,0,v0(k)/r0);x,y=pol2cart(u(:,3),u(:,1);X=flipud(x);x; Y=-flipud(y);y;text(-100-(k-1)*1

3、00,-(k-1)*300,sk);comet(X,Y);plot(X,Y)endE=-10;a=-G*M/2/E;r0=a,50,15;v0=sqrt(2*E+2*G*M./r0);figureaxis(-600,300,-300,300);xlabel(x); ylabel(y);title(E 相同L 不同時的不同曲線);hold onaxis equalsun=line(0,0,color,r,marker,.,markersize,20);for k=1:3t u=ode45(guidaofun,0:0.05:180,r0(k),0,0,v0(k)/r0(k);x,y=pol2cart(u(:,3),u(:,1);X=flipud(x);x; Y=-flipud(y);y;text(-100-(k-1)*100,-(k-1)*300,sk);comet(X,Y);plot(X,Y)endfunction ydot=guidaofun(t,y)global G Myd

溫馨提示

  • 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)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論