刘卫国matlab课后习题答案 .pdfVIP

刘卫国matlab课后习题答案

【篇一：matlab程序设计与应用第二版刘卫国课后实验

答案】

%第一小题

z1=2*sin(85*pi/180)/(1+exp(2))%%第二小题x=[2,1+2i;-0.45,5];

z2=1/2*log(x+sqrt(1+x.^2));z2

%%第三小题a=-3.0:0.1:3.0;

z3=1/2*(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)+log((0.3+a)/2)

%%第四题t=0:0.5:2.5

z4=(t=0t1).*(t.^2)+(t=1t2).*(t.^2-1)+(t=2t3).*(t.^2-2*t+1)

t2：

a=[12,34,-4;34,7,87;3,65,7]b=[1,3,-1;2,0,3;3,-2,7]disp(a+6*b=);

disp(a+6*b);disp(a-b+i=);disp(a-b+eye(3));disp(a*b=);

disp(a*b);disp(a.*b=);disp(a.*b);disp(a^3=);disp(a^3);

disp(a.^3=);disp(a.^3);disp(a/b=);disp(a/b);disp(b\a=);

disp(b\a);disp([a,b]=);disp([a,b]);

disp([a([1,3],:);b^2]=);disp([a([1,3],:);b^2]);t3:z=1:25;

a=reshape(z,5,5);

b=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];

1

c=a*b

d=c(3:5,2:3)t4-1:

a=100:999;

b=find(rem(a,21)==0);c=length(b)t4-2:

a=input(请输入一个字符串：,s);b=find(a=aa=z);a(b)=[];disp(a);

实验二：t1：

e=eye(3),r=rand(3,2),o=zeros(2,3),s=diag([1,2]);

a=[e,r;o,s]disp(a^2=);disp(a^2);

disp([e,r+rs;o,s^2]);b=[e,r+r*s;o,s^2]t2：h=hilb(5)p=pascal(5)

hh=det(h)hp=det(p)th=cond(h)tp=cond(p)a=abs(th-1);

b=abs(tp-1);ifab

disp(帕萨卡矩阵p性能更好);elseifab

disp(希尔伯特矩阵h性能更好);else

disp(两个矩阵性能相同);endt3:a=1:25;

a=reshape(a,5,5)

disp(行列式的值：);disp(det(a));

disp(矩阵的秩：);disp(rank(a));

disp(矩阵的迹：);

disp(trace(a));

disp(矩阵的范数：);disp(norm(a));t4:

a=[-29,6,18;20,5,12;-8,8,5][v,d]=eig(a)t5:

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]x1=a\b

b(3)=0.53x2=a\b

disp(系数矩阵a的条件数：);disp(cond(a));t6:a=1:25;

a=reshape(a,5,5)disp(b1=sqrtm(a));b1=sqrtm(a)

disp(b2=sqrt(a));b2=sqrt(a)disp(b1*b1);b1*b1

disp(b2.*b2);b2.*b2

实验三：t1:

x=-5.0:2:5.0fori=1:length(x);

if(x(i)0)(x(i)~=-3)y(i)=x(i)^2+x(i)-6;elseif

(x(i)=0)(x(i)5)(x(i)~=2)(x(i)~=3);y(i)=x(i)^2-5*x(i)+6;else

y(i)=x(i)^2-x(i)-1;endendyt2-if:

s=input(pleaseenterthescore:);while(s0||s100)

di