ex525.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 10 行

dnum=[3.4 -1.5];dden=[1 -1.6 0.8]; Ts=0.05;sysd=tf(dnum,dden,Ts);
num1=0.25;den1=[1 1];num2=1;den2=[1 2];
sys1=tf(num1,den1);sys2=tf(num2,den2);sysc2=sys1*sys2;
syscd=c2d(sysc2,Ts,'zoh');dsysc=sysd*syscd; [dnumc,ddenc]=tfdata(dsysc,'v');
%以上求出开环Z域脉冲传递函数模型：dsysc；dnumc,ddenc；Ts=0.05
[Gm,Pm,Wcg,Wcp]=margin(dsysc),margin(dsysc)   %求出系统频率特性参数，绘出Bode图
w=0.01:0.01:100;                              %对频率w向量赋值
figure(2);dnyquist(dnum,dden,Ts,w);             %绘制Nyquist图，或用nyquist(sysd,w)函数
figure(3);nichols(dsysc);                %绘制Nyquist图，或用nichols(dnum,dden,Ts,w)函数