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)函數