matlab计算零极点分布图

最近看书过程中看到一个matlab内容和大家分享分享

通过编程来计算单位冲击响应，单位阶跃y(t)和频率响应H(jw)

程序如下

clc%清空命令行 clear%从工作区中删除项目、释放系统内存 clf%清空当前图窗窗口 num=[1]; den=[1 2 2 1]; sys=tf(num,den);%构成传递函数 poles=roots(den);%求极点 subplot(211);%图1 pzmap(sys);%零极点分布图显示 [r,p,k]=residue(num,den)%实现部分分式展开的 t=0:0.02:10; h1=impulse(num,den,t);%单位冲击响应 subplot(222) plot(t,h1,'linewidth',4) title('冲击响应','fontsize',16) xlabel('时间\t','fontsize',16) h2=step(num,den,t);%单位阶跃响应 subplot(223) plot(t,h2,'linewidth',4) title('阶跃响应','fontsize',16) xlabel('时间\t','fontsize',16) [H,w]=freqs(num,den); subplot(224) plot(w,abs(H),'linewidth',4) xlabel('频率\omega','fontsize',16) title('幅频特性','fontsize',16) ylabel('幅值H(jw)','fontsize',16)

2、零极点分布图