用MATLAB画零极点图

1、用matlab将传递函数化为零极点增益模型并绘制零极点图将传递函数化为零极点增益模型并绘制零极点图 num=3 2 8; den=1 3 8 4 2; g=tf(num,den) transfer function: 3 s2 + 2 s + 8-s4 + 3 s3 + 8 s2 + 4 s + 2 g1=zpk(g); z=g1.z; p=g1.p; z=z:; p=p:; k=g1.k; pzmap(g); pzmap(g1); grid on将传递函数化为零极点增益模型并绘制零极点图 transfer function: 3 s2 + 2 s + 8 - s4 + 3 s3 + 8 s2

2、 + 4 s + 2 将上式化为零极点增益模型，并绘制零极点图。 num=3 2 8; den=1 3 8 4 2; g=tf(num,den)求传递函数。命令后无“;”则屏幕立即显示，否则不显示。transfer function: 3 s2 + 2 s + 8-s4 + 3 s3 + 8 s2 + 4 s + 2 g1=zpk(g);%化为零极点增益形式 z=g1.z;%将g1零点存入z p=g1.p; %将g1极点点存入p k=g1.k;%将g1增益存入k z=z:z = -0.3333 + 1.5986i -0.3333 - 1.5986i p=p:p = -1.2496 + 2.20

3、82i -1.2496 - 2.2082i -0.2504 + 0.4980i -0.2504 - 0.4980i显示零点显示极点%下面语句也可获得零极点。 z,p,k=tf2zp(num,den) pzmap(g1);%绘制零极点图 grid on %打开绘图网络鼠标指向某零点或极点会显示该方框。极点、阻尼比、超调量、固有频率wn.21nw-1.4-1.2-1-0.8-0.6-0.4-0.20-2.5-2-1.5-1-0.500.511.522.5pole-zero mapreal axisimaginary axiswnwn由斜边和实部求出。超调量由求出。21/ empwnwnwn化传递函

4、数为零极点形式 num=1 2 1; den=1 5 3 8 9; t=0.2; g2=tf(num,den,t) %有采样时间 t，则显示为脉冲传递函数，自变量为 ztransfer function: z2 + 2 z + 1- ，sampling time: 0.2 z4 + 5 z3 + 3 z2 + 8 z + 9 transfer function: z2 + 2 z + 1-z4 + 5 z3 + 3 z2 + 8 z + 9 sampling time: 0.2 g22=zpk(g2) %求离散系统g2零极点zero/pole/gain: (z+1)2-(z+4.635) (z

5、+1) (z2 - 0.6347z + 1.942) sampling time: 0.2 num=1 2 1; den=1 5 3 8 9; t=0.2; g1=tf(num,den) %没有采样时间 t，则显示为传递函数,自变量为 stransfer function: s2 + 2 s + 1-s4 + 5 s3 + 3 s2 + 8 s + 9 g11=zpk(g1) %求连续系统g1零极点zero/pole/gain: (s+1)2-(s+4.635) (s+1) (s2 - 0.6347s + 1.942) num=1 2 1; den=1 5 3 8 9; z,p,k=tf2zp

6、(num,den)%化传递函数为零极点增益形式z = -1 -1p = -4.6347 0.3174 + 1.3569i 0.3174 - 1.3569i -1.0000 k = 1化零极点增益形式为传递函数 num1,den1=zp2tf(z,p,k) %化零极点增益形式为传递函数num1 = 0 0 1 2 1den1 = 1.0000 5.0000 3.0000 8.0000 9.0000系统的串、并联的传递函数 2 s2 + 6 s + 5 g1= - s3 + 4 s2 + 5 s + 2 s2 + 7 s + 12 g2= - s2 + 3 s + 2 分别求g1、g2串、并联后的

7、传递函数。 num2=1 7 12; den2=1 3 2; g2=tf(num2,den2)； num1=2 6 5; den1=1 4 5 2; g1=tf(num1,den1)； g=series(g1,g2) 或 g=g1*g2 %串联transfer function: 2 s4 + 20 s3 + 71 s2 + 107 s + 60-s5 + 7 s4 + 19 s3 + 25 s2 + 16 s + 4 gp=parallel(g1,g2) 或 gp=g1+g2 %并联transfer function:s5 + 13 s4 + 57 s3 + 112 s2 + 101 s + 34- s5 + 7 s4 + 19 s3 + 25 s2 + 16 s + 4用零阶保持器将连续系统离散化)()()(?)(1)()1(1)(sgsgzzgzgsesgsssghsth求零阶保持器： num=1; den=conv(1 0,1 1); g=tf(num,den); g1=c2d(g,1,zoh