MATLAB 實(shí)習(xí)報(bào)告

1、1.計(jì)算矩陣 與 之和程序：cleara=5 3 5;3 7 4;7 9 8;b=2 4 2;6 7 9;8 3 6;c=b+a運(yùn)行結(jié)果：c =7 7 -79 14 1315 12 142. 求的共軛轉(zhuǎn)置。程序：x=4+8i 3+5i 2-7i 1+4i 7-5i ;3+2i 7-6i 9+4i 3-9i 4+4i;x運(yùn)行結(jié)果：ans = 4.0000- 8.0000i 3.0000 - 2.0000i3.0000-5.0000i 7.0000+6.0000i2.0000+7.0000i 9.0000-4.0000i1.0000-4.0000i 3.0000+9.0000i3. 對(duì)于，如果，求

2、解X。程序：clearA=4 9 2;7 6 4;3 5 7;B=37;26;28;X=AB運(yùn)行結(jié)果：X =-0.51184.04271.33184.已知：，分別計(jì)算a的數(shù)組平方和矩陣平方，并觀察其結(jié)果。程序：cleara=1 2 3;4 5 6;7 8 9;a.2a2運(yùn)行結(jié)果：ans =1 4 916 25 3649 64 81ans =30 36 4266 81 96102 126 1505. 計(jì)算多項(xiàng)式乘法(x2+2x+2)(x2+5x+4)程序：a=1 2 2;b=1 5 4;h=conv(a,b)運(yùn)行結(jié)果：h =1 7 16 18 86. 對(duì)下式進(jìn)行部分分式展開(kāi)：程序：a=1 3 4

3、 2 7 2;b=3 2 5 4 6;r,s,k=residue(b,a)運(yùn)行結(jié)果：r =1.1274 + 1.1513i1.1274 - 1.1513i-0.0232 - 0.0722i-0.0232 + 0.0722i0.7916 s =-1.7680 + 1.2673i-1.7680 - 1.2673i0.4176 + 1.1130i0.4176 - 1.1130i-0.2991 k = 7.計(jì)算多項(xiàng)式的微分和積分。程序：clearf=sym(4*x4-12*x3-14*x2+5*x+9)diff(f)clearf=sym(4*x4-12*x3-14*x2+5*x+9);int(f)運(yùn)行

4、結(jié)果：ans =16*x3-36*x2-28*x+5ans =4/5*x5-3*x4-14/3*x3+5/2*x2+9*x8. 因式分解： 程序：syms x;f=x4-5*x3+5*x2+5*x-6;factor(f)運(yùn)行結(jié)果： ans = (x-1)*(x-2)*(x-3)*(x+1)9. ，用符號(hào)微分求df/dx。程序：syms a x;f=a,x2,1/x;exp(a*x),log(x),sin(x);df=diff(f) 運(yùn)行結(jié)果：df = 0, 2*x, -1/x2 a*exp(a*x), 1/x, cos(x)10. 求代數(shù)方程組關(guān)于x,y的解。程序：S=solve(a*x2+b

5、*y+c=0,b*x+c=0,x,y);disp(S.x=),disp(S.x)disp(S.y=),disp(S.y)運(yùn)行結(jié)果：S.x=-c/bS.y=-c*(a*c+b2)/b311. 有一組測(cè)量數(shù)據(jù)滿足，t的變化范圍為0-10，用不同的線型和標(biāo)記點(diǎn)畫出a=0.1、a=0.2和a=0.5三種情況下的曲線。程序:t=0:0.5:10;y1=exp(-0.1*t);y2=exp(-0.2*t);y3=exp(-0.5*t);plot(t,y1,-.b,t,y2,:*r,t,y3,-g)運(yùn)行結(jié)果：12. 在5.1題結(jié)果圖中添加標(biāo)題和圖例框。程序：t=0:0.5:10;y1=exp(-0.1*t)

6、;y2=exp(-0.2*t);y3=exp(-0.5*t);plot(t,y1,-.b,t,y2,:*r,t,y3,-g)legend(a=0.1,a=0.2,a=0.5)title(y=exp(-a*x)運(yùn)行結(jié)果：13. x= 66 49 71 56 38，繪制餅圖，并將第五個(gè)切塊分離出來(lái)。程序：x=66 49 71 56 38;y=0 0 0 0 1;pie(x,y)運(yùn)行結(jié)果：14. ，當(dāng)x和y的取值范圍均為-2到2時(shí)，用建立子窗口的方法在同一個(gè)圖形窗口中繪制出三維線圖、網(wǎng)線圖、表面圖。程序：x,y=meshgrid(-2:0.2:2);z=x.*exp(-x.2-y.2);mesh(x

7、,y,z)subplot(2,2,1),plot3(x,y,z)title(plot3(x,y,z)subplot(2,2,2),mesh(x,y,z)title(mesh(x,y,z)subplot(2,2,3),surf(x,y,z)title(surf(x,y,z)運(yùn)行結(jié)果：15. 有傳遞函數(shù)如下的控制系統(tǒng)，用Simulink建立系統(tǒng)模型，并對(duì)系統(tǒng)的階躍響應(yīng)進(jìn)行仿真。程序： 運(yùn)行結(jié)果：16. 建立一個(gè)簡(jiǎn)單模型，用信號(hào)發(fā)生器產(chǎn)生一個(gè)幅度為2V、頻率為0.5Hz的正弦波，并疊加一個(gè)0.1V的噪聲信號(hào)，將疊加后的信號(hào)顯示在示波器上并傳送到工作空間。程序：運(yùn)行結(jié)果：17. 已知五個(gè)數(shù)據(jù)點(diǎn): (1

8、,5.5), (2,43.1), (3,128), (4,290.7), (5,498.4)(1)用三次曲線擬合上述數(shù)據(jù)點(diǎn)(2)在同一圖形中繪出數(shù)據(jù)點(diǎn)和擬合曲線程序：clearx=1 2 3 4 5;y=5.5 43.1 128 290.7 498.4;n=2;p=polyfit(x,y,n);x1=linspace(0,6,20);y1=polyval(p,x1);plot(x,y,o-b,x1,y1,r-.,linewidth,2);運(yùn)行結(jié)果：18. 解線性方程組a1x+b1y+c1z=d1a2x+b2y+c2z=d2a3x+b3y+d3z=d3，其中ai, bi, ci, di為常數(shù),

9、x, y, z為變量程序：g1=sym(a1*x+b1*y+c1*z=d1); g2=sym(a2*x+b2*y+c2*z=d2);g3=sym(a3*x+b3*y+c3*z=d3); f=solve(g1,g2,g3)x,y,z=solve(a1*x+b1*y+c1*z=d1,a2*x+b2*y+c2*z=d2,a3*x+b3*y+c3*z=d3)運(yùn)行結(jié)果：x =-(-b1*c2*d3+b1*d2*c3-c1*b3*d2+c1*b2*d3-d1*b2*c3+d1*b3*c2)/(a1*b2*c3-a1*b3*c2-a2*b1*c3+b3*a2*c1-b2*a3*c1+a3*b1*c2)y =

10、(-a1*c2*d3+a1*d2*c3-d2*a3*c1+c2*a3*d1+a2*c1*d3-a2*d1*c3)/(a1*b2*c3-a1*b3*c2-a2*b1*c3+b3*a2*c1-b2*a3*c1+a3*b1*c2)z =(b3*a2*d1-a1*b3*d2+a1*b2*d3-b2*a3*d1-a2*b1*d3+a3*b1*d2)/(a1*b2*c3-a1*b3*c2-a2*b1*c3+b3*a2*c1-b2*a3*c1+a3*b1*c2)19. 繪制下列各種函數(shù)圖形(1)繪制下列極坐標(biāo)圖形 r=3(1-cosq) r=2(1+cos q) r=2(1+sin q) r=cos3 q(

11、2)求函數(shù)z的三維圖形.定義區(qū)間與z的函數(shù)表達(dá)式如下:-5 x 5 , -5 y 5程序：(1)x=linspace(0,2*pi);r=3*(1-cos(x);polar(x,r);x=linspace(0,2*pi);r=2*(1+cos(x);polar(x,r);x=linspace(0,2*pi);r=2*(1+sin(x)polar(x,r);x=linspace(0,2*pi);r=cos(3*x);polar(x,r);運(yùn)行結(jié)果：(2)x=linspace(-5,5);y=linspace(-5,5); x,y=meshgrid(x,y); z=1(x+1).2+(y+1).2

12、+1-1(x-1).2+(y-1).2+1; mesh(x,y,z);運(yùn)行結(jié)果：20. 用MATLAB命令繪制下列數(shù)學(xué)函數(shù)的圖形(1) 繪制如下函數(shù)圖形y(t)=1-2exp(-t)sin(t) (0 t 8)且在x軸寫上“Time”,y軸寫上“Amplitude”標(biāo)號(hào),圖形標(biāo)題為“Decaying-oscillating Exponential”(2) 繪制如下圖形y(t)=5exp(-0.2t)cos(0.9t-30O)+0.8exp(-2t) (0 t 30)(3)在0 t10區(qū)間內(nèi)繪制下面圖形y(t)=1.23cos(2.83t+240O)+0.625x(t)=0.625t(4)在0

13、t20區(qū)間內(nèi),且在同一圖中繪制下述函數(shù)圖形y1(t)=2.62exp(-0.25t)cos(2.22t+174O)+0.6 y2(t)=2.62exp(-0.25t)+0.6 y3(t)=0.6程序：(1) t=0:0.2:8;y=1-2.*exp(-t).*sin(t);plot(t,y,linewidth,3);title(Decaying-oscillating Exponential);xlabel(Time)ylabel(Amplitude)(2) t=0:0.2:30;y=5*exp(-0.2*t).*cos(0.9*t-300)+0.8.*exp(-2*t);plot(t,y,linewidth,3);(3) t=0:0.2:10;y=1.23*cos(2.83*t+2400)+0.625;x=0.625*t;plot(t,y,-or,t,x,b*);(4) t=