[adsense:responsibe:9545213979]
clf; s1=%s/(2*%pi); den1=s1*(s1+1)*(s1+2); num1=s1^4+1; s2=%s; den2=s2*(s2+1)*(s2+2); num2=s2^4+1; g1=syslin('c',num1/den1); g2=syslin('c',num2/den2); nyquist(g1); xgrid; mtlb_axis([-4 0 -5 5]) w=[0.1, 0.5, 1, 1.41, 5, 10] p=horner(g2,w*%i) plot(real(p),imag(p),'o') [w;p] k1=5/6; l=num2+k1*den2; r=roots(l); h=pdiv(l,s^2+2);
Añadir nuevo comentario