matlab练习3.doc

上机练习(三) 用Matlab语言完成下列问题 (1) 已知矩阵 运用相关运算指令，回答下列问题 (a) 分别是多少？（其中I为3×3矩阵） (b) 的结果是多少？它们是否相同？为什么？ 解： (a) A=[4 12 20;12 45 78;20 78 90]; B=[1 2 3;4 5 6;7 8 9]; I=eye(3); J=inv(A); A+5*B-I A+J*B ans = 8 22 35 32 69 108 55 118 134 ans = 3.9167 12.8333 21.7500 12.1111 44.8889 77.6667 20.0000 78.0000 90.0000 (b)B=[1 2 3;4 5 6;7 8 9]; B^3 B.^3 ans = 468 576 684 1062 1305 1548 1656 2034 2412 ans = 1 8 27 64 125 216 343 512 729 前者：3个矩阵相乘。后者：矩阵里每元素三次方。 (2) 分别用for和while循环结构编写程序，求出 的结果，并写出源程序。 解： A=[1 2;2 1]; sum1=0; for i=1:5 sum1=sum1+A^i; end sum1 i=1; sum2=0; while i=5 sum2=sum2+A^i; i=i+1; End sum2 sum1 = 181 182 182 181 sum2 = 181 182 182 181 (3) 求出以下矩阵A的行列式、迹、秩、特征多项式和特征值。 解： A=[7 3 -2;3 4 -1;-2 -1 3]; det(A) %行列式； trace(A) %迹； rank(A) %秩； D=eig(A) %特征值； AA=sym(A); %特征多项式； poly(AA) ans = 46 ans = 14 ans = 3 D = 2.0000 2.3944 9.6056 ans = x^3-14*x^2+47*x-46 2．考虑如下Lorenz模型，其微分方程为： 其中参数为： ，试写出该方程的S函数表示，并分别通过Matlab 编程和Simulink对其进行仿真分析并绘制图形。 解： （1）S函数的表示 function [sys,x0]=vdp(t,x,u,flag) if flag==0 sys=[3;0;3;0;0;0];x0=[1 1 1]; elseif abs(flag)==1 sys(1)=-8/3*x(1)+x(2)*x(3); sys(2)=-10*x(2)+10*x(3); sys(3)=-x(2)*x(1)+28*x(2)-x(3); elseif flag==3 sys=x else sys=[] End 保存文件名vdp； (2)matlab 编程 平面图形 function dx=rs(t,x) dx=[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(2)*x(1)+28*x(2)-x(3); 保存文件名rs； t_end=10; x0=[1;1;1]; [t,x]=ode45(rs,[0,t_end],x0); plot(t,x(:,1),-,t,x(:,2),r,t,x(:,3),k) xlabel(时间（t）); ylabel(); legend(X1,X2,X3); title(Lorenz模型) 三维图形 x0=[1;1;1]; %初始状态 [t,y]=ode45(rs,[0,100],x0); plot3(y(:,1),y(:,2),y(:,3)) Simulink