Matlab第五章答案_原创精品文档.pdf

第一题

（1）a=[198;725;3-27]%产生矩阵

det(a)%检验是否可逆ans=-442,非0，可逆

div(a)%求逆矩阵

（2）b=[10-75;0-2672;7435;8-3215]

det(b)

div(b)

第二题

(1)A=[123;225;351];

B=[111231];

X=B/A

(2)A=[3105;0673;0430;2-126];

B=[2478];

X=B/A

第三题

（1）t=[12345678910];

y=[4.8424.3623.7543.3683.1693.0833.0343.016

3.0123.005];

A=[ones(size(t))exp(-t)];

T=[0:.1:10];

Y=[ones(size(T))exp(-T)]*C;

plot(T,Y,-,t,y,o)

title(采用y(t)≈c1+c2e–t的拟合)

xlabel(),ylabel()

（2）t=[0.3.81.11.62.3];

y=[.82.72.63.60.55.50];

A=[ones(size(t))t.*exp(-t)];

T=[0:.1:2.5];

Y=[ones(size(T))T.*exp(-T)]*C;

plot(T,Y,-,t,y,o)

title(采用y(t)≈d1+d2te–t拟合)

xlabel(),ylabel()

第四题

A=[11.5912.8115.66;15.24.1813.61;10.59

7.599.22];

[L,U]=lu(A)

[QR]=qr(A)

B=[16.004.41-10.37-21.61;0.88-20.0412.86

8.56;-1.4310.7118.81-5.99;-12.4824.35

-23.910.34];

[C,D]=lu(B)

[EF]=qr(B)

第五题

（1）A=[5-5-6;3-25;2-1-4];

x0=[1;-4;5];

X=[];

fort=0:.01:1

X=[Xexpm(t*A)*x0];

end

plot3(X(1,:),X(2,:),X(3,:),-o)

gridon

(2)A=[12-31;301-2;1-205;2301];

x0=[1;-1;2;1];

X=[];

fort=0:.01:1

X=[Xexpm(t*A)*x0];

end

plot3(X(1,:),X(2,:),X(3,:),-o)

gridon

第六题

（1）A=[11.5912.8115.66;15.24.1813.61;

10.597.599.22];

lambda=eig(A)

[V,D]=eig(A)

（2）B=[16.004.41-10.37-21.61;0.88-20.04

12.868.56;-1.4310.7118.81-5.99;-12.48

24.35-23.910.34];

lambda=eig(B)

[V,D]=eig(B)

第七题

（1）x=[12345678910];

y=[15.039.566.085.589.067.512.0-86.4-236.9

-448.4];

p=polyfit(x,y,2);

x2=1:.1: