MATLAB程序設(shè)計(jì)應(yīng)用試卷與答案(共3頁(yè))

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