MATLAB程序设计与应用(第二版)实验参考答案.doc

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) t=0:0.5:2.5; z4=t.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3) %第二题 A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; A+6*B A-B+eye(size(A)) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] %第三题 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*B F=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=lsdhKSDLKklsdkl; k=find(A=AA=Z); A(k)=[] %实验二 MATLAB矩阵分析与处理 %第一题 E=eye(3); R=rand(3,2); O=zeros(2,3); S=diag([2,3]); A=[E,R;O,S]; A^2 B=[E,(R+R*S);O,S^2] %第二题 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,5] [V,D]=eig(A) %数学意义略 %第五题方法一: %(1): A=[1/2,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)*B B1=[0.95,0.67,0.53]; X1=inv(A)*B1 N=cond(B) N1=cond(B1) Na=cond(A) %矩阵A为病态矩阵 %第六题 A=[1,4,9;16,25,36;49,64,81] B=sqrtm(A) C=sqrt(A) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算 %实验三 选择程序结构设计 %第一题程序一 x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0]; y=[]; %建立存放所有y值的矩阵 for x0=x if x00x0~=-3 y=[y,x0*x0+x0-6]; elseif x0=0x05x0~=2x0~=3 y=[y,x0*x0-5*x0+6]; else y=[y,x0*x0-x0-1]; end end x %输出所有x y %输出所有y %第一题程序二 x=[-5,-3,1,2,2.5,3,5]; y=[]; for a=1:7 if x(a)0x(a)~=-3 y=[y,(x(a))^2+x(a)-6]; elseif x(a)=0x(a)5x(a)~=2x(a)~=3 y=[y,(x(a))^2-5*x(a)+6]; else y=[y,x(a)*x(a)-x(a)-1]; end end %第二题程序一: x=input(请输入一个百分制成绩：); if x100|x0 disp(您输入的成绩不是百分制成绩，请重新输入。); else if x=100