Matlab解自控.doc

1、描述系统的微分方程：，用MATLAB命令求： （1）以分子、分母多项式形式建立系统的传递函数； （2）以零点、极点表达式形式建立系统传递函数； （3）求系统的零点、极点、增益； （4）画出零点、极点分布图； （5）在时间轴为20s的范围画出系统的单位阶跃响应曲线和脉冲响应曲线，并从图中求出等参数。 （1）num=[4 36 32];den=[1 10 30 40 24]; G=tf(num,den) Transfer function: 4 s^2 + 36 s + 32 --------------------------------- s^4 + 10 s^3 + 30 s^2 + 40 s + 24 （2）G1=zpk(G) Zero/pole/gain: 4 (s+8) (s+1) --------------------------- (s+6) (s+2) (s^2 + 2s + 2) （3）[zz,pp,kk]=zpkdata(G,v) zz = -8 -1 pp = -6.0000 -2.0000 -1.0000 + 1.0000i -1.0000 - 1.0000i kk = 4 (4) pzmap(G) （5）figure(2) step(G,20) 见图二 Tr=1.5s Tp=2.41s Ts=7.36s 超调量=14.3% figure(3) impulse(G,20) 见图三 Tr=0.411s Tp=0.64s Ts=6.8s 超调量=16526% 2、分析系统 的稳定性，并说明哪个系统是稳定的，造成不稳定的原因是什么？说明零点、极点对系统稳定性的影响。 num0=[1 3 4];den0=[1 8 30 76 80]; num1=[1 3 4];den1=[1 8 30 76 -80]; num2=[1 3 -4];den2=[1 8 30 76 80]; G1=tf(num0,den0) Transfer function: s^2 + 3 s + 4 -------------------------------- s^4 + 8 s^3 + 30 s^2 + 76 s + 80 G2=tf(num1,den1) Transfer function: s^2 + 3 s + 4 -------------------------------- s^4 + 8 s^3 + 30 s^2 + 76 s - 80 G3=tf(num2,den2) Transfer function: s^2 + 3 s - 4 -------------------------------- s^4 + 8 s^3 + 30 s^2 + 76 s + 80 pzmap(G1) pzmap(G2) pzmap(G3) step(G1,20) step(G2,20) step(G3,20) 1 3 稳定 2 不稳定 1 3 收敛 2 发散 3、设系统传递函数分析、绘制系统的单位阶跃响应和脉冲响应曲线，求出系统的极点，指出该系统响应是收敛的、发散的还是等幅振荡。 num=[1 1 1 1];den=[1 6 17 54 72];G=tf(num,den) Transfer function: s^3 + s^2 + s + 1 -------------------------------- s^4 + 6 s^3 + 17 s^2 + 54 s + 72 zz,pp,kk]=zpkdata(G,v) zz = -1.0000 -0.0000 + 1.0000i -0.0000 - 1.0000i pp = -0.0000 + 3.0000i -0.0000 - 3.0000i -4.0000 -2.0000 kk = 1 阶越响应 step(G,20) 等幅振荡 脉冲响应 impulse(G,20) 等幅振荡 4、系统结构图图中： 求系统传递函数，系统零点、极点和增益，系统的单位阶跃响应。 num0=[3 2 1];den0=[4 5 8]; num1=[4 0];den1=[1 2 10]; num2=[2 8];den2=[1