matlab答案1 8章,MATLAB基础与应用教程 习题答案 作者 蔡旭晖 刘卫国 蔡立燕 第1-8章答案...

8.

x=linspace(-3,3,100);y=linspace(-3,3,100);

[x y]=meshgrid(x,y);

fxy=-5./(1+x.^2+y.^2);

i=find(abs(x)<=0.8 & abs(y)<=0.5);

fxy(i)=NaN;

surf(x,y,fxy)

9.

u=linspace(1,10,100);v=linspace(-pi,pi,100);

[u,v]=meshgrid(u,v);

x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.^2;

x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2;

surf(x,y,z);shading interp;

light('position',[1,0,1]);

10.

t=0:pi/100:2*pi;

y=sin(t);

comet(t,y)

第4章

一、思考题

5.

(1)A=eye(3);

(2)C=100+(200-100)*rand(5,6);

(3)D=1+sqrt(0.2)*randn(1,500);

(4)E=ones(size(A));

(5)A=A+30*eye(size(A));

(6)B=diag(diag(A))

二、实验题

1.

P=pascal(5);H=hilb(5);

Dp=det(P);Dh=det(H);

Kp=cond(P);Kh=cond(H);

P矩阵的性能更好，因为Kp较小

2.

A=[1,-1,2,3;0,9,3,3;7,-5,0,2;23,6,8,3]

B=[3,pi/2,45;32,-76,sqrt(37);5,72,4.5e-4;exp(2),0,97]

A1=diag(A);B1=diag(B);

A2=triu(A);B2=triu(B);

A3=tril(A);B3=tril(B);

rA=rank(A);rB=rank(B);

nA=norm(A);nb=norm(B);

cA=cond(A);cB=cond(B);

3.

A=[31,1,0;-4,-1,0;4,-8,-2]；

[V,D]=eig(A);

4.

A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([2,2,2,2,2]) b=[1,0,0,0,0]';

x1=inv(A)*b;

x2=A\b;

[L,U]=lu(A);

x3=U\(L\b);

[Q,R]=qr(a);

[Q,R]=qr(A);