《Matlab與數(shù)學(xué)實(shí)驗(yàn)》期末作業(yè)(共5頁)

1、深圳大學(xué)考試答題紙(以論文、報(bào)告等形式考核專用)二一 四 二一 五 學(xué)年度第 一 學(xué)期課程編號(hào)1601600001課程名稱Matlab與數(shù)學(xué)實(shí)驗(yàn)主講教師仇高新評(píng)分學(xué) 號(hào)2013800394姓名周卓鵬專業(yè)年級(jí)2013級(jí) 電子科學(xué)與技術(shù)教師評(píng)語：題目：Matlab與數(shù)學(xué)實(shí)驗(yàn)期末作業(yè)1. 用數(shù)值積分公式計(jì)算 (結(jié)果保留小數(shù)點(diǎn)后8位)：(1) 取積分步長, 用梯形公式計(jì)算S= 6.26659757 。(2) 要求相對(duì)誤差為10-6, 用Simpson公式S= 6.26873546 ，相應(yīng)的Matlab命令是： s=quad(1-0.152*sin(h.2),0,2*pi,10e-6) 2.設(shè) 用數(shù)值解

2、法算出y(1)= 0.2711 （精確到4位小數(shù)）, 你用的方法是 先用ode45函數(shù)求出數(shù)值解，再用interp1函數(shù)插值求y（1） 調(diào)用的 Matlab命令是:function dy=hwfun2(x,y)dy =y(2);y(2)*sin(x)-y(1)*exp(x);format shortx,y=ode45(hwfun2,0,2,1,0);y1=interp1(x,y(:,1),1)3. 用電壓V=14伏的電池給電容器充電，電容器上t時(shí)刻的電壓滿足：，其中是電容器的初始電壓，是充電常數(shù)。試用下列數(shù)據(jù)確定和。t（秒）0.30.51.02.04.07.0v(t)5.68736.14347

3、.16338.862611.032812.6962 你用的方法是 lsqcurvefit ，結(jié)果是= 4.9711 ，= 3.5869 。4. 小型火箭初始質(zhì)量為900千克，其中包括600千克燃料?；鸺Q直向上發(fā)射時(shí)燃料以15千克/秒的速率燃燒掉，由此產(chǎn)生30000牛頓的恒定推力。當(dāng)燃料用盡時(shí)引擎關(guān)閉。設(shè)火箭上升的整個(gè)過程中，空氣阻力與速度平方成正比，比例系數(shù)為0.4（千克/米）。重力加速度取9.8米/秒2.A. 建立火箭升空過程的數(shù)學(xué)模型（微分方程）；解：引擎關(guān)閉前： h=(30000-0.4*h2)/(900-15*t)-9.8 所以：function dh=hwfun41(t,h)dh=

4、h(2);(30000-0.4*h(2)2)/(900-15*t)-9.8;關(guān)閉引擎后：h=-0.4*h2/300-9.8 所以：function dh=hwfun42(t,h)dh=h(2);(-0.4*h(2)2)/300-9.8;B. 求引擎關(guān)閉瞬間火箭的高度、速度、加速度，及火箭到達(dá)最高點(diǎn)的時(shí)間和高度。解：t1=600/15 h0=0; v0=0;t,h=ode45(hwfun41,0,t1,h0,v0);h1=h(end,1)v1=h(end,2)a1=(30000-0.4*v12)/300-9.8a11=-0.4*v12/300-0.98t,h=ode45(hwfun42,0,20

5、,h1,v1);h2=interp1(h(:,2),h(:,1),0)t2=interp1(h(:,2),t,0)結(jié)果為：t1 = 40 h1 =8.3230e+003 v1 = 258.9822a1 =0.7709 a11 = -90.4091h2 =9.1916e+003 t2 = 10.9511即：引擎關(guān)閉瞬間火箭的高度為8.32km,速度為258.98m/s,前一瞬間加速度為0.77m/s2,后一瞬間加速度為-90.41m/s2,火箭到達(dá)最高點(diǎn)時(shí)時(shí)間為t1+t2=50.95s,高度為9.19km。5. 冰淇淋的下部為椎體，上部為半球。設(shè)它由錐面和球面圍成，用蒙特卡羅方法計(jì)算它的體積。解

6、：易知椎體在（x,y,z）|0=x=1,0=y=1,0=z=1范圍內(nèi)，半球在（x,y,z）|0=x=1,0=y=1,1=zsqrt(x2+y2) n=n+1; end;end;s1=n/tmax;n=0;for t=1:tmax x=rand(); y=rand(); z=rand()+1; if(z=0 i=i+1;end;ip=0;for n=1:i p=p+13-n-3*n(2/3);end;p=p-20結(jié)果為：i = 4 p =0.4380即第4年時(shí)獲得最大利潤43.8萬元。8. 編寫能實(shí)現(xiàn)如下功能的主程序與子函數(shù)文件。屏幕提示從鍵盤輸入一個(gè)四位整數(shù)，輸入的整數(shù)按如下規(guī)則加密后輸出。加

7、密規(guī)則：每位數(shù)字都加上7，然后除以10所得的余數(shù)取代該數(shù)字。解：function hwfun8x=input(請(qǐng)輸入一個(gè)四位整數(shù)： );while x9999 fprintf(格式錯(cuò)誤!n) x=input(請(qǐng)輸入一個(gè)四位整數(shù)： );end;x=enc(x);fprintf( 加密后： %4dn,x) function a=enc(b)a=0;for i=1:4 x(i)=fix(rem(b/10(i-1),10); x(i)=rem(x(i)+7,10); a=a+x(i)*10(i-1);end運(yùn)行結(jié)果：hwfun8請(qǐng)輸入一個(gè)四位整數(shù): 1234 加密后: 89019. 結(jié)合你所學(xué)的某一門

8、課程，如高等數(shù)學(xué)、線性代數(shù)、大學(xué)物理等必修課程，或電路分析、信號(hào)與系統(tǒng)等專業(yè)選修課程，以課程中某一單元的知識(shí)點(diǎn)為例說明Matlab語言在其中的應(yīng)用的例子。要求首先對(duì)該知識(shí)點(diǎn)做背景介紹，然后通過具體實(shí)例的方式加以說明。解：電路分析中的節(jié)點(diǎn)電壓法及網(wǎng)孔電流法是求解復(fù)雜線性電路的一種簡單可靠的方法，該方法法原理簡單，方程對(duì)稱，不易出錯(cuò)。但用該方法所列的方程都是多元方程組，人工運(yùn)算解答繁瑣耗時(shí)，特別當(dāng)未知量較多時(shí)，人工運(yùn)算很不現(xiàn)實(shí)。而如果用matlab進(jìn)行方程組求解，將會(huì)變得非常簡單，極大提高效率，還能減小出錯(cuò)的概率。例如下題：以u(píng)0為參考電壓，用節(jié)點(diǎn)電壓法求u1,u2,u3:容易列出方程組：（1/1+1/0.2+1/0.5）*u1-2/1-i-3*u2-1/0.5*u3=0(1/1+1/0.5)*u2+i-1/1*u3=0(1/0.5+1/1+1/1)*u3-3*u2-1/0.5*u1-1/1*u2-3=0u1+1=u2整理得：8*u1+3*u2-2*u3-i=2 3*u2- u3+i=0-2*u1-4*u2+4*u3 =3 u1- u2 = -1matlab命令：u1,u2,u3,i0=solve(8*u1+3*u2-2*u3-i0=2,3