matlab习题与答案

1、 Matlab及工程应用-作业胡章芳 电气工程与自动化（1）班0513080109燕如龙重庆邮电大学电工程学院习题一 MATLAB运算基础一、目的1、熟悉启动和退出MATLAB的方法；2、熟悉MATLAB命令窗口的组成；3、掌握建立矩阵的方法；4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。二、内容1 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。提示：利用冒号表达式生成a向量，求各点的函数值时用点乘运算。>> 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整除的数的个数。提示；先利用冒号表达式，再利用find和Length函数。

8、>> a=100:999; c=rem(a/21,1); d=find(c=0); n=length(d)n = 43习题二程序设计及函数和数据的可视化一、目的1、掌握建立和执行M文件的方法；2、掌握利用if、switch、for、While语句实现的方法；3、掌握定义和调用MATLAB函数的方法；4、掌握绘制二维、三维图形的常用函数以及绘制图形的辅助操作。二、内容1求下列分段函数的值 要求：（1）用if语句实现，分别输出x-5.0，-3.0，1.0，2.0，2.5，3.0，5.0时的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语句实现。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根据，求p的近似值，当n分别取100、1000、10000时，结果是多少？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一个三位整数各位数字的

12、立方和等于该数本身则称该数为水仙花数。试输出全部水仙花数。 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定义三个函数文件，分别求给定双曲正弦、余弦和正切函数（如下），并在命令文件中调用该函数文件。 并用你定义的函数分别绘制上述三条曲线。>> 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)在同坐标系下用不同的颜色和线型绘制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绘制极坐标曲线Pasin(b+n)，并分析参数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函数绘制函数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)习题三Matlab在工程中的应用一、习题目的1、掌握生成特殊矩阵和矩阵分析的方法；2、掌握线性方程组的求解方法；3、了解统计和分析的方法以及多项式的常用运算；4、掌握求数值导数、数值积分、常微分方程数值求解、非线性代数方程数组求解的；5、掌握符号表达式的运算法则、符号矩阵运算、符号函数极限及导数、符号函数定积分和不定积分的方法。二、习题内容1分别用3种不同的数值方法解线性方程组。>> 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个学生5门功课的成绩存入矩阵P中，进行如下处理：(1) 分别求每门课的最高分、最低分及相应学生序号。(2) 分别求每门课的平均分和标推方差。(3) 5门课总分的最高分、最低分及相应学生序号。(4) 将5门课总分按从大到小顺序存入zcj中，相应学生序号存入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个多项式Pl(x)x4十2x3十4x2十5，P2(x)=x十2，P3(x)x2十2x十