東北石油大學(xué)測(cè)控專業(yè)Matlab報(bào)告.doc

東 北 石 油 大 學(xué)MATLAB環(huán)境與應(yīng)用院 系 電氣信息工程學(xué)院 專業(yè)班級(jí) 測(cè)控技術(shù)與儀器班 學(xué)生姓名 學(xué)生學(xué)號(hào) 指導(dǎo)教師 2011年 月 日第一部分 作業(yè)1. 在MATLAB中如何建立矩陣，并將其賦予變量a？程序：a=5 7 3;4 9 1運(yùn)行結(jié)果：a = 5 7 3 4 9 12計(jì)算矩陣與之和。程序：a=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 -7 9 14 13 15 12 143 求的共軛轉(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.0000i 3.0000-5.0000i 7.0000+6.0000i 2.0000+7.0000i 9.0000-4.0000i 1.0000-4.0000i 3.0000+9.0000i 7.0000+5.0000i 4.0000-4.0000i4對(duì)于，如果，求解X。程序：clearA=4 9 2;7 6 4;3 5 7;B=37;26;28;X=AB運(yùn)行結(jié)果：X = -0.5118 4.04271.33185 已知：，分別計(jì)算a的數(shù)組平方和矩陣平方，并觀察其結(jié)果。程序：a=1 2 3;4 5 6;7 8 9;a.2a2運(yùn)行結(jié)果：ans = 1 4 9 16 25 36 49 64 81ans = 30 36 42 66 81 96 102 126 1506. ，觀察a與b之間的六種關(guān)系運(yùn)算的結(jié)果。程序：a=1 2 3;4 5 6;b=8 -7 4;3 6 2;abans = 0 1 0 1 0 1a=bans = 0 1 0 1 0 1abans = 1 0 1 0 1 0a=bans = 1 0 1 0 1 0a=bans = 1 1 1 1 1 17. 將矩陣、和組合成兩個(gè)新矩陣：（1）組合成一個(gè)43的矩陣，第一列為按列順序排列的a矩陣元素，第二列為按列順序排列的b矩陣元素，第三列為按列順序排列的c矩陣元素，即 程序：a=4 2;5 7;b=7 1;8 3;c=5 9;6 2;% (1) d=a(:)b(:)c(:)運(yùn)行結(jié)果：d=4 7 55 8 62 1 97 3 2（2）按照a、b、c的列順序組合成一個(gè)行矢量，即 程序：e=a(:)b(:)c(:)運(yùn)行結(jié)果：e= 4 5 2 7 7 8 1 3 6 9 28. 將(x-6)(x-3)(x-8)展開為系數(shù)多項(xiàng)式的形式。程序：a=1 -6;b=1 -8;c=1 -30;y=conv(conv(a,b),c)運(yùn)行結(jié)果：y =1 -44 468 -14409. 求解多項(xiàng)式x3-7x2+2x+40的根。程序：a=1,-7,2,40;roots(a)運(yùn)行結(jié)果：ans = 5.0000 4.0000 -2.000010. 求解在x=8時(shí)多項(xiàng)式(x-1)(x-2) (x-3)(x-4)的值。程序：a=1 -1;b=1 -2;c=1 -3;d=1 -4;f=conv(conv(conv(a,b),c),d);x=8;y=polyval(f,x)運(yùn)行結(jié)果：y =84011. 計(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 812. 計(jì)算多項(xiàng)式除法(3x3+13x2+6x+8)/(x+4)。程序：h=3 13 6 8;a=1 4;f,r=deconv(h,a)運(yùn)行結(jié)果：f = 3 1 2r = 0 0 0 013. 對(duì)下式進(jìn)行部分分式展開： 程序：a=1 3 4 2 7 2;b=3 2 5 4 6;r,s,k=residue(b,a)運(yùn)行結(jié)果：r = 1.1274 + 1.1513i 1.1274 - 1.1513i -0.0232 - 0.0722i -0.0232 + 0.0722i 0.7916 s = -1.7680 + 1.2673i -1.7680 - 1.2673i 0.4176 + 1.1130i 0.4176 - 1.1130i -0.2991 k = 14. 計(jì)算多項(xiàng)式的微分和積分。程序：f=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)行結(jié)果：ans =16*x3-36*x2-28*x+5ans =4/5*x5-3*x4-14/3*x3+5/2*x2+9*x15. 解方程組。程序：a=2 9 0;3 4 11;2 2 6;b=13;6;6;x=ab運(yùn)行結(jié)果：x = 7.4000 -0.2000 -1.400016. 用符號(hào)函數(shù)法求解方程at2+b*t+c=0。程序：r=solve(a*t2+b*t+c=0,t)運(yùn)行結(jié)果：r = -1/2*(b-(b2-4*a*c)(1/2)/a -1/2*(b+(b2-4*a*c)(1/2)/a17. 因式分解：程序：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)18. ，用符號(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)19. 求代數(shù)方程組關(guān)于x,y的解。程序：S=solve(a*x2+b*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)/b320. 繪制曲線，x的取值范圍為-5,5。程序；x=-5:0.2:5;y=x.3+x+1;plot(x,y)運(yùn)行結(jié)果：21 有一組測(cè)量數(shù)據(jù)滿足，t的變化范圍為010，用不同的線型和標(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,-ob,t,y2,:*r,t,y3,-.g)運(yùn)行結(jié)果：22. 在5.1題結(jié)果圖中添加標(biāo)，并用箭頭線標(biāo)識(shí)出各曲線a的取值。程序：title(ityrm=e-itat)title(ityrm=e-itat,FontSize,12)text(t(6),y1(6),leftarrowitarm=0.1,FontSize,11)text(t(6),y2(6),leftarrowitarm=0.2,FontSize,11)text(t(6),y3(6),leftarrowitarm=0.5,FontSize,11)運(yùn)行結(jié)果：23. 在5.1題結(jié)果圖中添加標(biāo)題和圖例框。程序：title(ityrm=e-itat,FontSize,12)legend(a=0.1,a=0.2,a=0.5)運(yùn)行結(jié)果：24. x= 66 49 71 56 38，繪制餅圖，并將第五個(gè)切塊分離出來。程序：x=66 49 71 56 38;y=0 0 0 0 1;pie(x,y)運(yùn)行結(jié)果：25. ，當(dāng)x和y的取值范圍均為-2到2時(shí)，用建立子窗口的方法在同一個(gè)圖形窗口中繪制出三維線圖、網(wǎng)線圖、表面圖。程序：x,y=meshgrid(-2:.2:2);z=x.*exp(-x.2-y.2);mesh(x,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é)果：26. 有傳遞函數(shù)如下的控制系統(tǒng)，用Simulink建立系統(tǒng)模型，并對(duì)系統(tǒng)的階躍響應(yīng)進(jìn)行仿真。 程序：simulink運(yùn)行結(jié)果：27. 建立一個(gè)簡(jiǎn)單模型，用信號(hào)發(fā)生器產(chǎn)生一個(gè)幅度為2V、頻率為0.5Hz的正弦波，并疊加一個(gè)0.1V的噪聲信號(hào)，將疊加后的信號(hào)顯示在示波器上并傳送到工作空間。程序：simulink運(yùn)行結(jié)果：28編寫M腳本文件，分別使用for和while循環(huán)語(yǔ)句計(jì)算的程序。程序：(1)for語(yǔ)句clearsum=0;for i=1:1:10; sum=sum+ii;endsum(2)while語(yǔ)句cleari=1;sum=0;while i11 sum=sum+ii;endsum運(yùn)行結(jié)果：sum = 1.0405e+01029編寫M函數(shù)文件，將某班學(xué)生某門課的成績(jī)?yōu)椋?0,75,85,96,52,36,86,56,94,84,77，用switch結(jié)構(gòu)統(tǒng)計(jì)各分段的人數(shù)，并將各人的成績(jī)變?yōu)閮?yōu)、良、中、及格和不及格表示，統(tǒng)計(jì)人數(shù)和成績(jī)變換都用子函數(shù)實(shí)現(xiàn)。程序：clearx=60 75 85 96 52 36 86 56 94 84 77; a=0;b=0;c=0;d=0;e=0; x=fix(x/10);for i=1:1:11switch x(i)case 9a=a+1;disp()case 8b=b+1;disp()case 7c=c+1;disp()case 6d=d+1;disp()otherwisee=e+1;disp()endenddisp()adisp()bdisp()cdisp()ddisp()e運(yùn)行結(jié)果：及格中良優(yōu)不及格不及格良不及格優(yōu)良中得優(yōu)人數(shù)a = 2得良人數(shù)b = 3得中人數(shù)c = 2及格人數(shù)d = 1不及格人數(shù)e = 3第二部分 自選習(xí)題1角度x=30 45 60，求x的正弦、余弦、正切和余切。程序：x=30 45 60;x1=x/180*pi;運(yùn)行結(jié)果：sin(x1)ans = 0.5000 0.7071 0.8660cos(x1)ans = 0.8660 0.7071 0.5000tan(x1)ans = 0.5774 1.0000 1.7321cot(x1)ans = 1.7321 1.0000 0.57742用四舍五入的方法將數(shù)組2.4568 6.3982 3.9375 8.5042取整。程序：a=2.4568 6.3982 3.9375 8.5042;round(a)運(yùn)行結(jié)果：ans = 2 6 4 93,x從0到2， =0.02，求y的最大值、最小值、均值和標(biāo)準(zhǔn)差。程序：x=0:0.02*pi:2*pi;y=sin(x);ymax=max(y)ymin=min(y)ymean=mean(y)ystd=std(y)運(yùn)行結(jié)果：ymax = 1ymin = -1ymean = 2.2995e-017ystd =0.70714x=1 3 5 7 9,y=2 4 6 8 10,計(jì)算x的協(xié)方差、y的協(xié)方差、x與y的互協(xié)方差。程序：x=1 3 5 7 9;y=2 4 6 8 10;cx=cov(x)cy=cov(y)cxy=cov(x,y)運(yùn)行結(jié)果：cx = 10cy = 10cxy = 10 10 10 105將(x-6)(x-3)(x-8)展開為系數(shù)多項(xiàng)式的形式。程序：a=6 3 8;pa=poly(a);ppa=poly2sym(pa)運(yùn)行結(jié)果：ppa =x3-17*x2+90*x-1446用符號(hào)計(jì)算驗(yàn)證三角等式：程序：syms phi1 phi2;y=simple(sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2)運(yùn)行結(jié)果： y =sin(phi1-phi2)7.符號(hào)函數(shù)繪圖法繪制函數(shù)x=sin(3t)cos(t)，y=sin(3t)sin(t)的圖形，t的變化范圍是0,2程序：syms t;ezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),0,pi)運(yùn)行結(jié)果：8矩陣a=，求其行列式和逆矩陣。程序：a=2 9 0;3 4 11;2 2 6;ad=det(a)ai=inv(a)運(yùn)行結(jié)果：ad = 40ai = 0.0500 -1.3500 2.4750 0.1000 0.3000 -0.5500 -0.0500 0.3500 -0.47509建立一個(gè)模擬系統(tǒng)，將攝氏溫度轉(zhuǎn)換為華氏溫度()結(jié)果：10利用函數(shù)編程，判斷兩個(gè)正數(shù)的大小。程序：function minmax(a,b)if (a=0|bbfprintf(the larger number is %f,a)else if abfprintf(the larger number is %f,b)elsefpintf(%f equals %f,a,b)endendend運(yùn)行窗口的程序：cleara=input(input a number:);input a number:10b=input(input a number:);input a number:15minmax(a,b)運(yùn)行結(jié)果：the larger number is15.0000001011繪制peak函數(shù)的表面圖，用colormap函數(shù)改變預(yù)置的色圖，觀察色彩的分布情況。程序：surf(peaks(30);colormap