下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、 專業(yè)班級(jí): 學(xué)號(hào): 姓名: 命題教師: 該項(xiàng)由出卷人填寫 裝 訂 線一二三四五總 分標(biāo)準(zhǔn)分2020202020100得 分一、給出迭代方程 先編寫求解方程的函數(shù)文件,然后調(diào)用該函數(shù)文件求30 000個(gè)點(diǎn)上的x,y,最后在所有的()坐標(biāo)處標(biāo)記一個(gè)點(diǎn)(不要連線)會(huì)出圖形。這種圖形又稱為埃農(nóng)(Henon)引力線圖,它將迭代出來(lái)的隨機(jī)點(diǎn)吸引到一起,最后得出貌似連貫的引力線圖。二、分別利用數(shù)值積分法、符號(hào)積分法和Simulink仿真求 。三、已知阿波羅(Apollo)衛(wèi)星的運(yùn)動(dòng)軌跡()滿足下列微分方程:試在以上初值下進(jìn)行數(shù)值求解,并繪制出阿波羅衛(wèi)星位置()的軌跡。 提示:先選擇一組狀態(tài)變量,寫出一階常
2、微分方程組,并定義相應(yīng)的函數(shù)文件,然后求方程的數(shù)值解。四、實(shí)驗(yàn)圖4所示是一個(gè)蹺蹺板,兩板夾角為,左邊板長(zhǎng)為1.5m,上面的小孩重500N,右邊板長(zhǎng)為2m,小孩重400N。求當(dāng)蹺蹺板平衡時(shí),左邊木板與水平方向夾角的大小。要求先求解析解,然后給出兩種解決方案。 提示:這是一個(gè)力矩平衡問題,可列方程求解析解。可以考慮的兩種方案是用迭代法解方程和分別繪制兩個(gè)小孩所產(chǎn)生力矩隨變化的曲線,兩曲線的交點(diǎn)即是蹺蹺板平衡時(shí)的。 蹺蹺板示意圖五、某公司投資2000萬(wàn)元建成一條生產(chǎn)線。投產(chǎn)后,在時(shí)刻t的追加成本和追加收益分別為(百萬(wàn)元/年),(百萬(wàn)元/年)。試確定該生產(chǎn)線在何時(shí)停產(chǎn)可獲最大利潤(rùn)?最大利潤(rùn)是多少?提示
3、:利用函數(shù)(百萬(wàn)元),由于H(t)-G(t)單調(diào)下降,所以H(t)=G(t)時(shí),R(t)取得最大利潤(rùn)。1解: 源程序如下: x(1)=0;y(1)=0;for i=1:30000;x(i+1)=1+y(i)-1.4*x(i)2;y(i+1)=0.3*x(i);hold onplot(x(i),y(i),*b)end埃農(nóng)(Henon)引力線圖如下:2解:數(shù)值積分法源程序代碼: X=0:0.001:1;Y=(1/sqrt(2*pi)*exp(-X.2/2);trapz(X,Y) 運(yùn)行結(jié)果如下:ans = 0.3413符號(hào)積分法源程序代碼: x=sym(x);f=(1/sqrt(2*pi)*exp(
4、-x.2/2);I=int(f,0,1);double(I)運(yùn)行結(jié)果如下:ans = 0.3413Simulink仿真3. function dx=appollo(t,x)mu=1/82.45;mustar=1-mu;r1=sqrt(x(1)+mu)2+x(3)2);r2=sqrt(x(1)-mustar)2+x(3)2);dx=x(2)2*x(4)+x(1)-mustar*(x(1)+mu)/r13-mu*(x(1)-mustar)/r23x(4)-2*x(2)+x(3)-mustar*x(3)/r13-mu*x(3)/r23;-x0=1.2;0;0;-1.04935751;%x0(i)對(duì)應(yīng)
5、與xi的初值 options=odeset(reltol,1e-8); tic t,y=ode45(appollo,0,20,x0,options); toc plot(y(:,1),y(:,3) title(Appollo衛(wèi)星運(yùn)動(dòng)軌跡) xlabel(X) ylabel(Y)運(yùn)行結(jié)果如下:Elapsed time is 0.137341 seconds.圖如下:4.解:方案一:迭代法源程序如下 for alpha=0:0.001:pi/3;if 750*cos(alpha)-800*cos(pi/3-alpha)10e-6continue 專業(yè)班級(jí): 學(xué)號(hào): 姓名: 該項(xiàng)由出卷人填寫裝 訂
6、線 專業(yè)班級(jí): 學(xué)號(hào): 姓名: 該項(xiàng)由出卷人填寫裝 訂 線endbreakendalphacos(alpha)運(yùn)行結(jié)果如下:alpha = 0.4680ans = 0.8925方案二:曲線相交法源程序如下alpha=0:pi/1000:pi/3;y1=750*cos(alpha);y2=800*cos(pi/3-alpha);plot(alpha,y1,g-,alpha,y2,r)grid on運(yùn)行得到的效果圖如下:該圖放大后可大致得到alpha和cos(alpha)的值,見下圖:5.解:構(gòu)造函數(shù)f(t)=H(t)-G(t)=13-t-3t2/3=0;令t1/3=x,則f(t)=-t3-3t2
7、+13可得矩陣P=-1,-3,0,13求最佳生產(chǎn)時(shí)間的源程序如下:p=-1,-3,0,13;x=roots(p);t=x.3運(yùn)行結(jié)果如下:t = 3.6768 +21.4316i 3.6768 -21.4316i 4.6465 再分別將t的三個(gè)值帶入函數(shù)f(t),比較大小后,得到最大利潤(rùn)與最佳生產(chǎn)時(shí)間。求最大利潤(rùn)的程序代碼如下: t=3.6768 +21.4316i;x=0:0.01:t;y=13-x-3*x.(2/3);trapz(x,y)運(yùn)行結(jié)果: ans = 25.2583 t=3.6768 -21.4316i;x=0:0.01:t;y=13-x-3*x.(2/3);trapz(x,y)運(yùn)行結(jié)果: ans = 25.2583 t=4.64
溫馨提示
- 1. 本站所有資源如無(wú)特殊說明,都需要本地電腦安裝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ù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 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ì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 藝術(shù)展覽館樣板房裝修合同
- 2024年姿態(tài)敏感器項(xiàng)目建議書
- 美發(fā)店裝修保修期規(guī)定
- 糧食收購(gòu)運(yùn)輸承攬合同
- 武漢商場(chǎng)裝修合同模板
- 旅游景區(qū)裝修保密合同
- 游戲開發(fā)項(xiàng)目居間投資協(xié)議
- 液化氣罐運(yùn)輸安全管理合同
- 陶瓷制品倉(cāng)儲(chǔ)配送合同
- 通信基站運(yùn)輸協(xié)議范本
- 小學(xué)北師大版三年級(jí)數(shù)學(xué)上冊(cè)期中復(fù)習(xí)(課件)
- 梁祝(梁山伯與祝英臺(tái))克萊德曼(原版)鋼琴雙手簡(jiǎn)譜 鋼琴譜
- 公共關(guān)系學(xué)-實(shí)訓(xùn)項(xiàng)目1:公關(guān)三要素分析
- 房地產(chǎn)廣告圍擋施工投標(biāo)文件范本
- 人教版2022年四年級(jí)上冊(cè)數(shù)學(xué)期中考試考點(diǎn)檢查試卷
- 腫瘤患者的營(yíng)養(yǎng)治療及全程管理
- 慢性腎小球腎炎英文教學(xué)課件
- 智慧礦山建設(shè)規(guī)劃技術(shù)方案
- 課件《“多元一體”視域下的中國(guó)古代民族關(guān)系》
- 三大秘訣解決鋰電池不一致的問題
- 更換限矩型液力偶合器作業(yè)流程
評(píng)論
0/150
提交評(píng)論