MATLAB仿真作业..docx

Matlab语言及其在控制与仿真中的应用仿真实验与作业任课教师 杨 惠 珍 学生姓名 项 志 杰 班 级 11205-6 学 号 2011200872a=rand(5,6)b=a(1,3)c=a(4,2)d=a(:,1)e=a(4,:)[i,j]=find(a0.5)运行结果：a = 0.9501 0.7621 0.6154 0.4057 0.0579 0.2028 0.2311 0.4565 0.7919 0.9355 0.3529 0.1987 0.6068 0.0185 0.9218 0.9169 0.8132 0.6038 0.4860 0.8214 0.7382 0.4103 0.0099 0.2722 0.8913 0.4447 0.1763 0.8936 0.1389 0.1988b = 0.6154c = 0.8214d = 0.9501 0.2311 0.6068 0.4860 0.8913e = 0.4860 0.8214 0.7382 0.4103 0.0099 0.2722i = 1 3 5 1 4 1 2 3 4 2 3 5 3 3j = 1 1 1 2 2 3 3 3 3 4 4 4 5 6A=[1 -1;2 5];I=eye(2,2);B=2*A^4-12*A^3+19*A^2-29*A+37*I;F=inv(B)运行结果：F = 0.3043 0.0435 -0.0870 0.1304syms x1 x2 x3 tdiff_equ1=Dx1=-x1+x2;diff_equ2=Dx2=-4*x1+3*x2;diff_equ3=Dx3=-8*x1+8*x2-x3;[x1,x2,x3]=dsolve(diff_equ1,diff_equ2,diff_equ3,t)运行结果：x1 = exp(t)*(C2+C3*t) x2 = exp(t)*(2*C2+2*C3*t+C3)x3 = 4*C2*exp(t)+4*C3*exp(t)*t+2*C3*exp(t)+exp(-t)*C1AA=[1 -1;0 2];BB=[-3 4;1 0];CC=[1 3;-2 2];XX=lyap(AA,BB,-CC)运行结果：XX = 0 2.0000 1.0000 -1.0000t1=-10:0.01:10;t2=10:0.001:20;y1=-t1+5;y2=(1+t2).*sin(t2);plot(t1,y1);hold on;plot(t2,y2);hold off运行结果：自定义函数建立一个lorenz.m文件：function dx=lorenz(t,x)dx=zeros(3,1);dx(1)=-8*x(1)/3+x(2)*x(3);dx(2)=-10*x(2)+10*x(3);dx(3)=-x(1)*x(2)+28*x(2)-x(3);end求解如下：[t,x]=ode45(@lorenz,[0,100],[0;0;0.001]);figure(1)plot(t,x(:,1))figure(2)plot(t,x(:,2))figure(3)plot(t,x(:,3))figure(4)plot3(x(:,1),x(:,2),x(:,3))运行结果：A=[-1 1 0;0 -1 -3;-1 -5 -3];B=[0 0;1 0;0 1];C=[0 -1 0];D=[1 -5];sys1=ss(A,B,C,D);sys2=tf(sys1)运行结果：Transfer function from input 1 to output:s^3 + 4 s^2 - 12 s - 18-----------------------s^3 + 5 s^2 - 8 s - 15 Transfer function from input 2 to output:-5 s^3 - 25 s^2 + 43 s + 78--------------------------- s^3 + 5 s^2 - 8 s – 15s=tf(s)Gc=0.48*(1+1/(1.814*s)+0.4353*s/(1+0.04353*s));G=10/(s+1)^3;s1=Gc*Gsys=feedback(s1,1)sys1=ss(sys)sys2=zpk(sys)运行结果：Transfer function:s Transfer function: 4.169 s^2 + 8.916 s + 4.8-------------------