《MATLAB程序设计与应用(第二版)实验参考答案(可编辑)》由会员分享,可在线阅读,更多相关《MATLAB程序设计与应用(第二版)实验参考答案(可编辑)(38页珍藏版)》请在金锄头文库上搜索。
1、MATLAB程序设计与应用(第二版)实验参考答案MATLAB程序设计与应用(第二版)实验参考答案%实验一 MATLAB运算基础%第一题%(1)z1 2*sin 85*pi/180 / 1+exp 2 %(2)x 2,1+2i;-0.45,5;z2 0.5*log x+sqrt 1+x.2 % 3 a -3.0:0.1:3.0;z3 exp 0.3*a -exp -0.3*a /2.*sin a+0.3 +log 0.3+a /2 % 4 ;z4 t.2.* t 0&t 1 + t.2-1 .* t 1&t 2 + t.2-2*t+1 .* t 2&t 3 %第二题A 12 34 -4;34 7
2、 87;3 65 7;B 1 3 -1;2 0 3;3 -2 7;A+6*BA-B+eye size A A*BA.*BA3A.3A/BBAA,BA 1,3,: ;B2%第三题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 25B 3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11C A*BF size C D C F 1 -2:F 1 ,F 2 -1:F 2 whos%第四题% 1 :A 100:999;B rem A,21 ;C length find B 0 % 2 :A lsdhKS
3、DLKklsdkl;k find A A&A Z ;A k %实验二第一题E eye 3 ;R rand 3,2 ;O zeros 2,3 ;S diag 2,3 ;A E,R;O,S;A2B E, R+R*S ;O,S2%第二题H hilb 5 P pascal 5 Hh det H Hp det P Th cond H Tp cond P %第三题:A fix 10*rand 5 H det A Trace trace A Rank rank A Norm norm A 第四题:A -29,6,18;20,5,12;-8,8,5V,D eig A %数学意义略第五题:% 1 :A 1/2,
4、1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;b 0.95,0.67,0.52;x inv A *b% 2 :B 0.95,0.67,0.53;x inv A *B% 3 :cond A %第五题:A hilb 4 A :,1 A 4,: B 0.95,0.67,0.52;X inv A *BB1 0.95,0.67,0.53;X1 inv A *B1N cond B N1 cond B1 Na cond A %矩阵A为病态矩阵第题A 1,4,9;16,25,36;49,64,81B sqrtm A C sqrt A %sqrtm函数是以矩阵为单位进行计算,sqrt函数是以矩
5、阵中的元素进行计算%实验三 选择程序结构设计%第一题x -5.0,-3.0,1.0,2.0,2.5,3.0,5.0;y ; %建立存放所有y值的矩阵for x0 x if x0 0&x0 -3 y y,x0*x0+x0-6; elseif x0 0&x0 5&x0 2&x0 3 y y,x0*x0-5*x0+6; else y y,x0*x0-x0-1; endendx %输出所有xy %输出所有y%第一题程序二x -5,-3,1,2,2.5,3,5;y ;for a 1:7if x a 0&x a -3y y, x a 2+x a -6;elseif x a 0&x a 5&x a 2&x
6、a 3y y, x a 2-5*x a +6;elsey y,x a *x a -x a -1;endend%第二题:x input 请输入一个百分制成绩: ;if x 100|x 0disp 您输入的成绩不是百分制成绩,请重新输入。 ;elseif x 100&x 90disp A ;elseif x 89&x 80disp B ;elseif x 79&x 70disp C ;elseif x 69&x 60disp D ;elsedisp E ;endend%第二题s input 请输入一个成绩(0分到100分之间): ; %s用于存放成绩while 1 %判断输入成绩的合理性 if s
7、 0|s 100 disp 输入的成绩需在0到100之间,请重新输入: s input 请输入一个成绩(0分到100分之间): ; else break; endendswitch fix s/10 %对成绩做出等级判断 case 9,10 disp A case 8 disp B case 7 disp C case 6 disp D otherwise disp E end%第三题n input 请输入员工工号: ;h input 该员工工作时数是: ;if h 120x h-120 *84* 1+0.15 +120*84;elseif h 60x h*84-700;elsex h*84;
8、enddisp num2str n ,号员工的应发工资为x ;%第四题(还可以用switch语句实现)a fix 10+ 99-10 *rand 1,2 %产生两个随机整数x a 1 ;y a 2 ;t input 请输入运算符号:,s ;if t + z x+y;elseif t - z x-y;elseif t * z x*y;elseif t / z x/y;enddisp num2str x ,t,num2str y , ,num2str z %输出运算结果%第五题a rand 5,6 %产生5x6的随机矩阵n input 请输入您要输出矩阵的第几行: ;if n 5 disp 超出了
9、矩阵的行数,矩阵的最后一行为: a 5,: else disp 矩阵的第,num2str n ,行为: a n,: end实验四第一题s 0;n input n ? ;for i 1:n s s+1/i/i;endPI sqrt 6*s pi%第一题n input n ? ;a 1:n;b 1./a.2;PI sqrt 6*sum b pi%第二题y 0;n 1;while y 3 y y+1/ 2*n-1 ;n n+1;endy y-1/ 2* n-1 -1 n n-2%第题 input a ? ;b input b ? ;Xn 1;Xn1 a/ b+Xn ;n 0;while abs Xn
10、1-Xn 1e-5 Xn Xn1; Xn1 a/ b+Xn ; n n+1; if n 500 break; endendnXn1r1 -b+sqrt b*b+4*a /2r2 -b-sqrt b*b+4*a /2%第四题for i 1:100if i 1f i 1;elseif i 2f i 0;elseif i 3f i 1;elsef i f i-1 -2*f i-2 +f i-3 ;endend f min f sum f length find f 0 length find f 0 length find f 0 %第五题:s 0;n 0;for i :49b i* i+1 -1;
11、m fix sqrt b ;for j 2:mif rem b,j 0breakendendif j mn n+1;s s+b;endendns%实验五 函数文件%第一题y exp x ,log x ,sin x ,cos x ;%在命令窗口调用上述函数文件:y mat1 1+i %第二题function a,b,N,M shiyanwu2 m,n,t A m*cos t*pi/180 ,-m,-sin t*pi/180 ,0;m*sin t*pi/180 ,0,cos t*pi/180 ,0;0,n,-sin t*pi/180 ,0;0,0,-cos t*pi/180 ,1;B 0,9.8*
12、m,0,9.8*n;C inv A *B;a C 1 ;b C 2 ;N C 3 ;M C 4 ;%在命令窗口调用该函数文件m1 input m1 ;m2 input m2 ;theta input theta ;a1,a2,N1,N2 shiyanwu2 m1,m2,theta %第二题function X mat2 m1,m2,t g 9.8;A m1*cos t*pi/180 ,-m1,-sin t*pi/180 ,0;m1*sin t*pi/180 ,0,cos t*pi/180 ,0;0,m2,-sin t*pi/180 ,0;0,0,-cos t*pi/180 ,1;B 0;m1*g;0;m2*g;X inv A *B;%在命令窗口调用该函数文件:X mat2 1,1,60 %第三题function flag mat3 x flag 1;for i 2:sqrt x if rem x,i 0 flag 0; break;endend%在命令窗口调用该函数文件:for i 10:99j 10*rem