matlab程序设计与应用第二版刘卫国课后实验解答

《matlab程序设计与应用第二版刘卫国课后实验解答》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、Matlab程序设计与应用第二版刘卫国课后实验答案实验一:C=A*BT1:D=C(3:5,2:3)T4-1:%%第一小题z1=2*sin(85*pi/180)/(1+exp(2))a=100:999;%%第二小题b=find(rem(a,21)==0);x=[2,1+2i;-0.45,5];c=length(b)z2=1/2*log(x+sqrt(1+x.^2));T4-2:z2a=input('请输入一个字符串:','s');%%第三小题b=find(a>='A'&a<='Z');a=-3.0:0.1:3.0;a(b)=[];z3=1

2、/2*(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)+ldisp(a);og((0.3+a)/2)实验二:%%第四题T1:t=0:0.5:2.5E=eye(3),R=rand(3,2),O=zeros(2,3),S=diag([1,2z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t]);>=2&t<3).*(t.^2-2*t+1)A=[E,R;O,S]disp('A^2=');T2:disp(A^2);A=[12,34,-4;34,7,87;3,65,7]disp('[

3、E,R+RS;O,S^2]');B=[1,3,-1;2,0,3;3,-2,7]B=[E,R+R*S;O,S^2]disp('A+6*B=');T2:disp(A+6*B);H=hilb(5)disp('A-B+I=');P=pascal(5)disp(A-B+eye(3));Hh=det(H)disp('A*B=');Hp=det(P)disp(A*B);Th=cond(H)disp('A.*B=');Tp=cond(P)disp(A.*B);a=abs(Th-1);disp('A^3=');b=abs(Tp-1);disp(A^3);

4、ifa>bdisp('A.^3=');disp('帕萨卡矩阵P性能更好');disp(A.^3);elseifa

5、p(det(A));z=1:25;disp('矩阵的秩:');A=reshape(z,5,5)';disp(rank(A));B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];disp('矩阵的迹:');1disp(trace(A));enddisp('矩阵的范数:');if(s>=90&s<=100);disp('A');disp(norm(A));T4:elseif(s>=80&s<90);A=[-29,6,18;20,5,12;-8,8,5]disp('B');[V,D]=eig(A)elseif(

6、s>=70&s<80);T5:disp('C');A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]elseif(s>=60&s<70);B=[0.95,0.67,0.52]'disp('D');X1=ABelseB(3)=0.53disp('E');X2=ABenddisp('系数矩阵A的条件数:');T2-switch:disp(cond(A));s=input('pleaseenterthescore:');T6:while(s<0

7、

8、s>100)a=1:25;disp('thescoreisnot

9、reasonable');A=reshape(a,5,5)s=input('pleaseenterthescore:');disp('B1=sqrtm(A)');endB1=sqrtm(A)switchfix(s/10)disp('B2=sqrt(A)');case{9,10}B2=sqrt(A)disp('A');disp('B1*B1');case{8}B1*B1disp('B');disp('B2.*B2');case{7}B2.*B2disp('C');实验三:case{6}T1:disp('D');x=-5.0:2:5.0ca

10、se{0,1,2,3,4,5}fori=1:length(x);disp('E');if(x(i)<0)&(x(i)~=-3)endy(i)=x(i)^2+x(i)-6;T3:elseift=input('