matlab習(xí)題與答案

1、 Matlab及工程應(yīng)用-作業(yè)胡章芳 電氣工程與自動化（1）班0513080109燕如龍重慶郵電大學(xué)電工程學(xué)院習(xí)題一 MATLAB運(yùn)算基礎(chǔ)一、目的1、熟悉啟動和退出MATLAB的方法；2、熟悉MATLAB命令窗口的組成；3、掌握建立矩陣的方法；4、掌握MATLAB各種表達(dá)式的書寫規(guī)則以及常用函數(shù)的使用。二、內(nèi)容1 先求下列表達(dá)式的值，然后顯示MATLAB工作空間的使用情況并保存全部變量。提示：利用冒號表達(dá)式生成a向量，求各點(diǎn)的函數(shù)值時(shí)用點(diǎn)乘運(yùn)算。>> z1=2*sin(85*pi/180)/(1+exp(2)z1 = 0.2375>> x=2 1+2i;-0.45 5;

2、>> z2=1/2*log(x+sqrt(1+x2)z2 = 0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i>> a=-3.0:0.1:3.0;z3=1/2*(exp(0.3*a)-exp(0.3*a).*sin(a+0.3)z3 = Columns 1 through 6 0 0 0 0 0 0 Columns 7 through 12 0 0 0 0 0 0 Columns 13 through 18 0 0 0 0 0 0 Columns 19 through 24 0 0 0

3、 0 0 0 Columns 25 through 30 0 0 0 0 0 0 Columns 31 through 36 0 0 0 0 0 0 Columns 37 through 42 0 0 0 0 0 0 Columns 43 through 48 0 0 0 0 0 0 Columns 49 through 54 0 0 0 0 0 0 Columns 55 through 60 0 0 0 0 0 0 Column 61 0>> A=12 34 -4;34 7 87;3 65 7;>> B=1 3 -1;2 0 3;3 -2 7;>> A+6

4、*Bans = 18 52 -10 46 7 105 21 53 49>> A-B+eye(3)ans = 12 31 -3 32 8 84 0 67 1>> A*Bans = 68 44 62 309 -72 596 154 -5 241>> A.*Bans = 12 102 4 68 0 261 9 -130 49>> A3ans = 37226 233824 48604 247370 149188 600766 78688 454142 118820>> A.3ans = 1728 39304 -64 39304 343 658

5、503 27 274625 343>> A/Bans = 16.4000 -13.6000 7.6000 35.8000 -76.2000 50.2000 67.0000 -134.0000 68.0000>> BAans = 109.4000 -131.2000 322.8000 -53.0000 85.0000 -171.0000 -61.6000 89.8000 -186.2000>> A,Bans = 12 34 -4 1 3 -1 34 7 87 2 0 3 3 65 7 3 -2 7>> A(1,3,:);B2ans = 12 34

6、-4 3 65 7 4 5 1 11 0 19 20 -5 40（1）求它們的乘積C；（2）將矩陣C的右下角3x 2子矩陣賦給D；（3）查看MATLAB工作空間的使用情況。>> A=1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25 ;>> B=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11;>> C=A*BC = 93 150 77 258 335 237 423 520 397 588 705 557 753 890 717>> D=A(

7、3:5,4:end)D = 14 15 19 20 24 25>> whoYour variables are:A B C D >> whos Name Size Bytes Class A 5x5 200 double array B 5x3 120 double array C 5x3 120 double array D 3x2 48 double arrayGrand total is 61 elements using 488 bytes4完成下列操作：（1）求100，999之間能被21整除的數(shù)的個數(shù)。提示；先利用冒號表達(dá)式，再利用find和Length函數(shù)。

8、>> a=100:999; c=rem(a/21,1); d=find(c=0); n=length(d)n = 43習(xí)題二程序設(shè)計(jì)及函數(shù)和數(shù)據(jù)的可視化一、目的1、掌握建立和執(zhí)行M文件的方法；2、掌握利用if、switch、for、While語句實(shí)現(xiàn)的方法；3、掌握定義和調(diào)用MATLAB函數(shù)的方法；4、掌握繪制二維、三維圖形的常用函數(shù)以及繪制圖形的輔助操作。二、內(nèi)容1求下列分段函數(shù)的值 要求：（1）用if語句實(shí)現(xiàn)，分別輸出x-5.0，-3.0，1.0，2.0，2.5，3.0，5.0時(shí)的y值。>> x=input('please input x:'); i

9、f x<0&x=-3 y=x2+x-6 elseif 0<=x<10,x=2&x=3 y=x2-5*x+6 else y=x2-x-1 endplease input x:-5.0y =14please input x:-3.0y = 30please input x:1.0y = 2please input x:2.0y = 0please input x:2.5y = -0.2500please input x:3.0y = 0please input x:5.0y = 62輸入一個百分制成績，要求輸出成績等級A,B,C,D,E.其中9099分為A，808

10、9分為B，7079分為C，6069分為D，60分以下為E。要求： 分別用if語句和switch語句實(shí)現(xiàn)。function switch_case(number); g=input('please input number:','s'); switch(g) case('9099') disp('A'); case('8089') disp('B'); case('7079') disp('C'); case('6069') disp('D

11、9;); case('<60') disp('E'); otherwise disp('輸入錯誤!') end 3根據(jù)，求p的近似值，當(dāng)n分別取100、1000、10000時(shí)，結(jié)果是多少？n=input('please input n:');s=0;for m=1:ns=s+1.0/(m*m); t=sqrt(6*s); endtplease input n:100t = 3.1321please input n:1000t = 3.1406please input n:10000t = 3.14154一個三位整數(shù)各位數(shù)字的

12、立方和等于該數(shù)本身則稱該數(shù)為水仙花數(shù)。試輸出全部水仙花數(shù)。 for a=1:9; for b=0:9; for c=0:9; i=100*a+10*b+c; if a3+b3+c3=i; disp(i); end endendend 153 370 371 4075定義三個函數(shù)文件，分別求給定雙曲正弦、余弦和正切函數(shù)（如下），并在命令文件中調(diào)用該函數(shù)文件。 并用你定義的函數(shù)分別繪制上述三條曲線。>> function imean stdev=stat(x) mean=1/2*(exp(x)-exp(-x); stdev=sinh(x); mean=stdev;function im

13、ean stdev=stat(x) mean=1/2*(exp(x)+exp(-x); stdev=sinh(x); mean=stdev;>> function imean stdev=stat(x) mean=(exp(x)-exp(-x)/(exp(x)+exp(-x); stdev=sinh(x); mean=stdev; 6已知y1x2，y2cos(2x)，y3=y1*y2，完成下列操作：(1)在同坐標(biāo)系下用不同的顏色和線型繪制3條曲線； x=0:2*pi; y1=x.2; plot(x,y1,'r-'); hold on; y2=cos(2*x); pl

14、ot(x,y2,'g-.'); hold on; y3=y1.*y2; plot(x,y3); hold off; (2)以子圖形式繪制3條曲線。x=0:2*pi; y1=x.2; plot(x,y1) ; x=0:2*pi; y2=cos(2*x); plot(x,y2) x=0:2*pi; y1=x.2; y2=cos(2*x); y3=y1.*y2; plot(x,y3); (3)分別用條形圖、階梯圖、桿圖和填充圖繪制3條曲線。 x=0:2*pi; y1=x.2; stairs(x,y1,'r-'); hold on; y2=cos(2*x); stair

15、s(x,y2,'g-.'); hold on; y3=y1.*y2; stairs(x,y3); hold off; x=0:2*pi; y1=x.2; stem(x,y1,'r-'); hold on; y2=cos(2*x); stem(x,y2,'g-.'); hold on; y3=y1.*y2; stem(x,y3); hold off; x=0:2*pi; y1=x.2; hist(x,y1,'r-'); hold on; y2=cos(2*x); hist(x,y2,'g-.'); hold on;

16、y3=y1.*y2; hist(x,y3); hold off; 7繪制極坐標(biāo)曲線Pasin(b+n)，并分析參數(shù)a，b，n對曲線形狀的影響. t = 0 : .01 : 2 * pi; polar(t,2*sin(4+5*t),'-r') t = 0 : .01 : 2 * pi; polar(t,22*sin(44+55*t),'-r') 8分別用plot和fplot函數(shù)繪制函數(shù)ysin(1/x)的曲線，分析兩曲線的差別。alpha=1./x; y=sin(alpha); plot(x,y) x=0:1:50;y=sin(1./x);figure(2)fpl

17、ot(x) sin(1./x),0,50)習(xí)題三Matlab在工程中的應(yīng)用一、習(xí)題目的1、掌握生成特殊矩陣和矩陣分析的方法；2、掌握線性方程組的求解方法；3、了解統(tǒng)計(jì)和分析的方法以及多項(xiàng)式的常用運(yùn)算；4、掌握求數(shù)值導(dǎo)數(shù)、數(shù)值積分、常微分方程數(shù)值求解、非線性代數(shù)方程數(shù)組求解的；5、掌握符號表達(dá)式的運(yùn)算法則、符號矩陣運(yùn)算、符號函數(shù)極限及導(dǎo)數(shù)、符號函數(shù)定積分和不定積分的方法。二、習(xí)題內(nèi)容1分別用3種不同的數(shù)值方法解線性方程組。>> A=6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 0 2;>> B=-4 13 1 11'>> x=ABx

18、 = 0.6667 -1.0000 1.5000 -0.0000>> A=6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 0 2; B=-4 13 1 11' C=inv(A); x=C*Bx = 0.6667 -1.0000 1.5000 -0.0000>> x,y,z,u=solve('6*x+5*y-2*z+5*u+4','9*x-y+4*z-u-13','3*x+4*y+2*z-2*u-1','3*x-9*y+2*u-11') x = 0 y = 2/3 z = -1 u

19、= 3/22將20個學(xué)生5門功課的成績存入矩陣P中，進(jìn)行如下處理：(1) 分別求每門課的最高分、最低分及相應(yīng)學(xué)生序號。(2) 分別求每門課的平均分和標(biāo)推方差。(3) 5門課總分的最高分、最低分及相應(yīng)學(xué)生序號。(4) 將5門課總分按從大到小順序存入zcj中，相應(yīng)學(xué)生序號存入xsxh.p=1,67,78,65,87,89;2,98,78,67,89,78;.3,99,78,98,67,78;4,78,98,45,67,89;.5,64,65,78,99,86;6,70,71,73,87,88;.7,82,73,67,56,89;8,61,47,87,56,0;.9,46,97,67,87,56;1

20、0,34,68,56,87,88;.11,56,78,78,74,75;12,56,67,86,76,98;.13,78,44,56,43,56;14,98,32,56,65,16;.15,67,76,57,87,78;16,45,56,56,67,87;.17,67,67,56,83,84;18,87,89,89,90,69;.19,9,78,78,66,87;20,56,78,78,78,89;a=p(:,2);b=p(:,3);c=p(:,4);d=p(:,5);e=p(:,6);a1=max(a) a2=min(a)b1=max(b)b2=min(b)c1=max(c)c2=min(c)d1=max(d)d2=min(d)e1=max(e)e2=min(e)mean(p)std(p) a1 = 99a2 = 9b1 = 98b2 = 32c1 = 98c2 = 45d1 = 99d2 = 43e1 = 98e2 = 0ans = 10.5000 65.9000 70.9000 69.6500 75.5500 74.0000ans = 5.9161 22.4591 16.5112 14.1022 14.3361 25.1438 3有3個多項(xiàng)式Pl(x)x4十2x3十4x2十5，P2(x)=x十2，P3(x)x2十2x十