num=0.692;
den=[1,-1.758,0.375];
[z,p,k]=tf2zp(num,den);
p
figure(1)
dnyquist(num,den,0.1)
title('離散Nyquist 曲線圖');
xlabel('實數軸');
ylabel('虛數軸');
figure(2)
[num1,den1]=cloop(num,den);
dimpulse(num1,den1);
title('離散沖激響應');
xlabel('時間');
ylabel('振幅');