北理工自动控制实验报告

Pole-Zero Map1.50.9450.9760.890.810.680.50.31Imaginary Axis (seconds-1)0.50.994504321-0.50.994-10.9760.9450.89-3.5-3-2.50.81-2-10.68-1.5-10.50.3-0.50-1.5-5-4.5-4Real Axis (seconds)

2.已知控制系统的状态空间方程如下

试用MATLAB 建立系统的传递函数模型、零极点增益模型及系统的状态空间方程模型，并绘制系统零极点图。

代码：

A=[0 1 0 0;0 0 1 0;0 0 0 1;-1 -2 -3 -4]; B=[0;0;0;1]; C=[10 2 0 0]; D=0;

Gss=ss(A,B,C,D); Gtf=tf(Gss); Gzpk=zpk(Gss); pzmap(Gzpk); grid on; 实验结果：

零极点图：

Pole-Zero Map0.80.9840.9640.9350.870.740.450.60.9930.4Imaginary Axis (seconds-1)0.9990.2504321-0.20.999-0.4-0.60.9930.984-0.8-5-4.5-4-3.50.964-3-2.50.935-20.87-1.50.74-10.45-0.50Real Axis (seconds-1)

3已知三个系统的传递函数分别为

试用MATLAB 求上述三个系统串联后的总传递函数。 代码： num1=[2 6 5]; den1=[1 4 5 2]; G1=tf(num1,den1); num2=[1 4 1]; den2=[1 9 8 0]; G2=tf(num2,den2);

num3=conv(conv([1 3],[1 7]),[0 5]); den3=conv(conv([1 1],[1 4]),[1 6]); G3=tf(num3,den3);

G=G1*G2*G3; 实验结果：

4.已知如图E2-1 所示的系统框图

试用MATLAB 求该系统的闭环传递函数。 代码： num1=1; den1=[1 1]; G1=tf(num1,den1); num2=1; den2=[ 1]; G2=tf(num2,den2); num3=3; den3=[1 0]; G3=tf(num3,den3); G=(G1+G2)*G3; T=feedback(G,G2); U=G/(1+G2*G); 实验结果：