clear all
syms s w
G = 1/((s)*(s+1)*(s+2)); %transfer function
G_w = subs(G,s,j*w);
W= [-100:0.01:100]; %[min_range:step size:max_range]
nyq = eval(subs(G_w,w,W));
x = real(nyq)
y = imag(nyq
有沒有一種方法來確定nyquist圖的穩定性，以檢查它是穩定還是不穩定？我目前需要手動輸入到另一個m.file傳遞函數來查看PZ位置，以確定 syms s w
G = (2*s+1)/((s^2)+(5*s)+6); %transfer function
G_w = subs(G,s,j*w);
W= [-100:0.1:100]; %[min_range:step
未定義變量I寫檢查奈奎斯特准則的封閉系統的穩定性的算法（http://en.wikipedia.org/wiki/Nyquist_stability_criterion） function answear=stability(re,im)
%% Function check stability of system
%re is real part of transmitation
%im