2次遅れ系にPI制御系を入れた時の閉ループ伝達関数Tc
s=%s G=1/(s^2+2*s+1) C=1+1/s Tc=G*C/(1+G*C) nyquist(syslin("c",Tc) #単位円プロット→show_margins(SISO,'nyquist')関数を使うときは不要 a = linspace(0, 2*%pi, 100); x=cos(a) y=sin(a) plot2d(x,y)
show_margins関数
show_margins(syslin("c",Tc)) show_margins(syslin("c",Tc),"nyquist")
p_margin, g_margin関数
-->p_margin(syslin("c",Tc)) ans = 25.841933 -->g_margin(syslin("c",Tc)) ans = Inf