东北大学自动化控制系统计算机辅助设计实验

东北大学11控制系统计算机辅助设计——MATLAB 语言与应用控制系统计算机辅助设计第一部分1 function>> x0=[0;0;0];>> [t,y]=ode45('lorenzeq',[0,100],x0);>> plot(t,y)>> figure;plot3(y(:,1),y(:,2),y(:,3)),grid,2function>> Aeq=[];Beq=[];A=[];B=[];>> xm=[0;0]; xM=[];x0=[0;0];>> f1=inline('x(1)^2-2*x(1)+x(2)');>> [x,f]=fmincon(f1,x0,A,B,Aeq,Beq,xm,xM,'f2a');x',fans = 1.0000 0f =-13（a）> >s=tf('s');G=(s^3+4*s+2)/s^3/(s^2+2)/[(s^2+1)^3+2*s+5]Transfer function: s^3 + 4 s + 2------------------------------------------------------s^11 + 5 s^9 + 9 s^7 + 2 s^6 + 12 s^5 + 4 s^4 + 12 s^3（b）>> z=tf('z',0.1);>> H=(z^2+0.568)/(z-1)/(z^2-0.2*z+0.99)Transfer function:-----------------------------4>> A=[0 1 0;0 0 1;-5 -4 -13];>> B=[0;0;2];>> C=[1 0 0;0 0 0;0 0 0];>> D=[0];精品文档---下载后可任意编辑>> G=ss(A,B,C,D);>> Ga = x1 x2 x3 x1 0 1 0 x2 0 0 1 x3 -5 -4 -13b = u1 x1 0x2 0 x3 2c = x1 x2 x3 y1 1 0 0 y2 0 0 0 y3 0 0 0d = u1 y1 0 y2 0 y3 0Continuous-time model.>> G=tf(G)Transfer function from input to output... 2 #1: ---------------------- s^3 + 13 s^2 + 4 s + 5 #2: 0 #3: 0>> GG=zpk(G)Zero/pole/gain from input to output... 2 #1: ---------------------------------- (s+12.72) (s^2 + 0.2836s + 0.3932) #2: 0 #3: 0根据微分方程也可以直接写出传递函数模型：>> num=[2];>> den=[1,13,4,5];>> G=tf(num,den);>> GTransfer function: 2----------------------s^3 + 13 s^2 + 4 s + 5>> GG=zpk(G)Zero/pole/gain: 2----------------------------------(s+12.72) (s^2 + 0.2836s + 0.3932)5>> num=[1,2];>> den=[1,1,0.16];>> H=tf(num,den,'Ts',1);>> HTransfer function:z + 2--------------Sampling time: 16functionifend>> syms J Kp Ki s;>> gc=(Kp*s+Ki)/s;>> g=(s+1)/(J*s^2+2*s+5);>> gg=feedback(g*gc,1)>> gg=feedback(g*gc,1)gg =((Ki + Kp*s)*(s + 1))/(J*s^3 + (Kp + 2)*s^2 + (Ki + Kp + 5)*s + Ki)7（a）>> s=tf('s');>> G=(211.87*...