[adsense:responsibe:9545213979]
s=%s; g=1/(s*(s^2+s+4)); gs=syslin('c',g); ang=40*2*%pi/360; //Vamos a calcular wc aux=tan(ang); aux1=aux*s^2+s-4*aux; aux2=roots(aux1) wc=aux2(1) //Valor de K aux3=horner(gs,wc*%i); aux4=abs(aux3); k=1/aux4 //Vamos a calcular wcf aux5=tan(%pi/2); aux6=aux5*s^2+s-4*aux5; aux7=roots(aux6) wcf=aux7(1) //Vamos a calcular el margen de ganancia para K aux8=horner(k*gs,wcf*%i); aux9=abs(aux8); aux10=20*log10(aux9) //comprobaciones aux11=horner(k*gs,wc*%i) aux12=abs(aux11); aux13=atan(imag(aux11),real(aux11)); aux14=aux13*360/(2*%pi) s1=s/(2*%pi) gb=1/((s1+0.0000000000001)*(s1^2+s1+4)); gbs=syslin('c',gb); clf; bode(gbs) show_margins(gbs) [mp,frp]=p_margin(gbs) [mg,frg]=g_margin(gbs)
En este problema pasa una cosa curiosa debido al pico de resonancia de aumenta la frecuencia de cruce de ganancia y se reduce el margen de fase al calculado.
Añadir nuevo comentario