matlab习题与答案

1、Matlab及工程应用作业胡章芳电气工程与自动化（1）班0513080109燕如龙重庆邮电大学电工程学院习题一 MATLAB运算基础一、目的1、熟悉启动和退出MATLAB的方法；2、熟悉MATLAB命令窗口的组成；3、掌握建立矩阵的方法；4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。二、内容1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存 全部变量。八、I 2sin 85m三F1 + e(2) z2 = -yin (x 十)，其中第=：I2 -0.455 J03。 0山(3) z3二 sin (a +0.3),a = -3.0, -2,9, -2,89282

2、9,3,0提示：利用冒号表达式生成J向量，求各点的函数值时用点乘运算。 z 1 =2*sin(85*pi/l 80)/( 1 +exp(2) zl =0.2375 x=2 l+2i;-0.45 5; z2=l/2*log(x+sqrt( l+xA2)z2 =0.7114 - O.O253i0.8968 + O.3658i0.2139 + O.9343i1.1541 - 0.0044i a=-3.O:O:3O;z3=l/2*(exp(O.3*a)-exp(O.3*a).*sin(a+O.3)z3 =Columns 1 through 60 0 0000Columns 7 through 120

3、0 0000Columns 13 through 180 0 0000Columns 19 through 240 0 0000Columns 25 through 300 0 0000Columns 31 through 360 0 0000Columns 37 through 420 0 0000Columns 43 through 480 0 0000Columns 49 through 540 0 0000Columns 55 through 600 0 0000Column 6102.已知T234-4,T3-rA =34787,B二203-3657 -3-27 求下列表达式的值：(1

4、) A + 6* B和A-B + H其中1为单位矩阵)。(2) A* B 和 A. * Bo(3) A3 和 A.3C(4) A/B 及 BAoA,B和A(l,3, A=l 2 34 -4;34 7 87;3 65 7; B=l 3-l;2 0 3;3 -2 7; A+6*B ans =1852-10467105215349 A-B+eye(3)ans =1231-3328840671A*B ans =684462309-72596154-5241A.*B1210246802619-13049 AA3ans =372262338244860424737014918860076678688454

5、142118820 A.A3ans =172839304273930434327462564658503343 A/Bans =16.4000 -13.60007.600035.8000 -76.200050.200067.0000-134.000068.0000 BAans =109.4000 -131.2000 322.8000-53.000085.0000-171.0000-61.600089.8000-186.2000 A,Bans =1234-413-13478720336573-27 A(L3,:);BA2ans =3 6574 511101920-540-3170-6169B 0

6、23-4970-413lb（1）求它们的乘积C；（2）将矩阵C的右下角3x 2子矩阵赋给D：（3）查看MATLAB工作空间的使用情况。 A= 1 2 3 45;67 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25 ; B二3 0 16;176 9;0 23 4;9 7 0;4 13 11; C=A*B9315077258335237423520397588705557753890717 D=A(3:5,4:end)D =141519202425 whoYour variables are: whosNameSizeBytesClassA5

7、x5200double arrayB5x3120double arrayC5x3120double arrayD3x248double arrayGrand total is 61 elements using 488 bytes4. 完成下列操作：(1) 求100, 999之间能被21整除的数的个数。提示；先利用冒号表达式，再利用find和Length函数。 a= 100:999;c=rem(a/21 J);d=find(c=0);n=length(d)n =43习题二 程序设计及函数和数据的可视化一、目的1、掌握建立和执行M文件的方法；2、掌握利用if、switch、for、While语句

8、实现的方法;3、掌握定义和调用MATLAB函数的方法；4、掌握绘制二维、三维图形的常用函数以及绘制图形的辅助操作。二、内容（久2十X 6 , it2 - 5% + 6, X2 - X - I ,1. 求下列分段函数的值算0且策萨一3Ow兀 10,x2且戈工3其他要求：（1）用 if 语句实现，分别输出 x=-5.0, -3.0, 1.0, 2.0, 2.5, 3.0, 5.0 时 的y值。 x=input（*please input x:r）;if xvO&x=3y=xA2+x-6elseif Ov二xvlO,x=2&x=3y=xA2-5*x+6elsey=xA2-x-lendplease i

9、nput x:-5.0y =14please input x:-3.0y =30please input x:1.0y =2please input x:2.0please input x:2.5-0.2500please input x:3.0y =oplease input x:5.0y =62. 输入一个百分制成绩，要求输出成绩等级A,B,C,D,E.其中90-99分为A, 80-89 分为B, 70-79分为C, 6069分为D, 60分以下为E。要求：分别用if语句和switch语句实现。function switch_case(number);g=input(please inpu

10、t number:,s);switch(g)case(90-99)disp(A);case(K089)disp(B);case(7079)disp(C);case(6069)dispCD);case(60)disp(E);otherwisedispC输入错误!)end3根据冬=丄+丄+丄+丄，求的近似值，当n分别取100、1000、10000 6 r 23 汇时，结果是多少？n=input(*please input n:9);s=0;for m=l:ns=s+1.0/(m*m);t=sqrt(6*s);endtplease input n:1003.1321please input n:10

11、003.1406please input n: 100003.14154. 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。试输出 全部水仙花数。for a=l:9;for b=0:9;for c=0:9;i=100*a+10*b+c;if aA3+bA3+cA3=i;disp(i);endendendend1533703714075. 定义三个函数文件，分别求给定双曲正弦、余弦和正切函数（如下），并在命 令文件中调用该函数文件。并用你定义的函数分别绘制上述三条曲线。 function iniean stdev=stat(x)mean= 1 /2*(exp(x)-exp(-x);

12、stdev=sinh(x);mean=stdev;function imean 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.已知 yl=x2, y2=cos(2x), y3=yl*y2,完成下列操作:(1)在同一一坐标系下用不同的颜色和线型绘制3条曲线;x=0:2*pi;yl=x、2;plot(xrylf *r);hold

13、on;y2=cos(2*x);plot(xzy2f *g-. 1);hold on;y3=yl.*y2;plot(xzy3);hold off;以子图形式绘制3条曲线。 x=0:2*pi;yl=x;plot(xzyl)；x=0:2*pi; y2=cos(2*x);plot(xry2)x=0:2*pi; yl=x; y2=cos(2*x); y3=yl.*y2; plot(xzy3);(3)分别用条形图、阶梯图、杆图和填充图绘制3条曲线。 x=0:2*pi;yl=x;stairs(x z ylrr-1);hold on;y2=cos(2*x);stairs(xzy2 r g-)；hold on;

14、y3=yl.*y2;stairs(x,y3);hold off;403020100-2030x=0:2*pi;yl=x入2;stem(x z ylr r1); hold on;y2=cos(2*x); stem(xry2 f 9g- . 9 ); hold on;y3=yl.*y2;stem(x r y3);hold off;403020100:-10-20x=O:2*pi;yl=x入2;hist(xzylrr);hold on;y2=cos(2*x);hist(xry2,g-);hold on;y3=yl.*y2;hist(xzy3)；hold off;7.绘制极坐标曲线P=asin(b+n

15、0),并分析参数a, b, n对曲线形状的影响.t = 0 :01 : 2 * pi;polar(tz 2*sin(4+5*t) ,1-r1)t = 0 :01:2* pi;polar(t,22*sin(44+55*t)t 1-r1)2708.分别用plot和fplot函数绘制函数y=sin(l/x)的曲线，分析两曲线的差别。alpha=l/x;y=sin(alpha);plot(x,y)x=0:1:50;y=sin(1/x)；figure(2)fplot(x) sin(l/x),0,50)习题三 Matlab在工程中的应用一、习题目的1、掌握生成特殊矩阵和矩阵分析的方法；2、掌握线性方程组的

16、求解方法；3、了解统计和分析的方法以及多项式的常用运算；4、掌握求数值导数、数值积分、常微分方程数值求解、非线性代数方程数组求 解的；5、掌握符号表达式的运算法则、符号矩阵运算、符号函数极限及导数、符号函 数定积分和不定积分的方法。二、习题内容1.分别用3种不同的数值方法解线性方程组。6彩+ 5尹-2旷+5杉=-49久 一 y + 4z u = 133% + 2z - 21/ = 13% - 9y + 2u = 11 A=6 5 -2 5;9-1 4-l;342 -2;3 -9 0 2; B=-4 13 1 11; x=AB0.6667-1.00001.5000-0.0000 A=6 5 -2

17、 5;9-1 4-l;342 -2;3 -9 0 2;B=-4 13 1 11;C=inv(A);x=C*B0.6667-1.00001.5000X =-0.0000 x,y,z,u=solveC6*x+5*y-2*z+5*u4479*x-y44*z-u-l 33*x+4*y+2*z-2*u-1 3*x-9* y+2*u-ir)x =0y =2/3z =-1 u =3/22.将20个学生5门功课的成绩存入矩阵P中，进行如下处理：(1) 分别求每门课的最高分、最低分及相应学生序号。(2) 分别求每门课的平均分和标推方差。(3) 5门课总分的最高分、最低分及相应学生序号。 将5门课总分按从大到小顺

18、序存入zcj中，相应学生序号存入xsxh.p=l,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;9z46z97,67,87,56;10,34,68,56z87,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,

19、87;17,67,67,56,83,84;18,87,89,89,90,69;19,9,78,78,66,87;20,56z78z78,78,89;a=p(:,2)；b=p(:,3)；a=p(：,4);d=p(:,5);e=p(:,6)；al=max(a)a2=min(a)bl=max (b) b2=min (b) cl=max (c) c2=min (c) dl=max (d) d2=min (d) el=max (e) e2=min (e) mean(p) std(p)al =99a2 =9bl =98b2 =32cl =98c2 =45dl =99d2 =43el =98e2 =0ans

20、 =10.500065 900070.900069.650075.550074.0000ans =5 916122.459116.511214 102214336125.14383有 3 个多项式 Pl(x) = x4 十 2x?十 4x2 十 5, P2(x)=x 十 2, P3(x)=x2 十 2x 十 3, 试进行下列操作：(1)求 P(x)=Pl(x)*P2(x)*P3(x)o syms x pl p2 p3 pp 1 =xA4+2*xA3+4*xA2+5;p2=x+2;p3=xA2+2*x+3;p=pl*p2*p3P =(xA4+2*xA3+4*xA2+5)*(x+2)*(xA2+2*x+3)(2)求P(x)的根。 x=solve(,(xA4+2*xA344*xA2+5)*(x+2)*(xA2+2*x+3)=0,)x =-20.315579229403120539036910574680410.9977863722378007037546503988143*10.9