Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Prepara tus exámenes
Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Prepara tus exámenes con los documentos que comparten otros estudiantes como tú en Docsity
Encuentra los documentos específicos para los exámenes de tu universidad
Estudia con lecciones y exámenes resueltos basados en los programas académicos de las mejores universidades
Responde a preguntas de exámenes reales y pon a prueba tu preparación
Consigue puntos base para descargar
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Comunidad
Pide ayuda a la comunidad y resuelve tus dudas de estudio
Ebooks gratuitos
Descarga nuestras guías gratuitas sobre técnicas de estudio, métodos para controlar la ansiedad y consejos para la tesis preparadas por los tutores de Docsity
ejercicios de la tarea 3 y ene este esta el desarrollo y tambien sus diferentes aportes que se realisan mas el del aplicativo
Tipo: Guías, Proyectos, Investigaciones
1 / 14
Esta página no es visible en la vista previa
¡No te pierdas las partes importantes!
En oferta
2
Con esto vamos a empezar por diseñar la ecuación para poder empezar a estabilizar el sistema con la metodología del compensador de ganancia, lo primero será hacer el diagrama de bode clc clear all Mos=20; %%compensador de ganancia nump=5; a0=1; a1=3; a2=7; Gp= tf((nump),[a0 a1 a2]); figure(1) margin(Gp) grid on Vamos a encontrar lo sigiente El máximo de obershut el cual es 0,095 por que es de -9.5 porciento zeta=sqrt(log(Mos/100)^2/(pi^2+log(Mos/100)^2)); El coeficiente de amortiguamiento pm=atand(2zeta/(sqrt(-2zeta^2+sqrt(1+4*zeta^4))));
Con esto vamos a empezar por diseñar la ecuación para poder empezar a estabilizar el sistema con la metodología del compensador de ganancia, lo primero será hacer el diagrama de bode Mos=15; %%compensador de ganancia s=tf('s'); nump=3; denp=s((7s)+3); Gp= nump/denp; figure(1) margin(Gp) grid on Vamos a encontrar lo sigiente El máximo de obershut el cual es 0,095 por que es de -9.5 porciento zeta=sqrt(log(Mos/100)^2/(pi^2+log(Mos/100)^2)); El coeficiente de amortiguamiento pm=atand(2zeta/(sqrt(-2zeta^2+sqrt(1+4*zeta^4)))); El margen de fase p=-180+pm; Vamos a buscar en el diagrama de bode donde la fase Recordemos que la fase es de -180 mas la margen de fase el cual nos daría
[mag,phase,wout]=bode(Gp); wm=interp1(squeeze(phase),wout,p); Am= interp1(squeeze(phase), squeeze(mag),p); Kc=1/Am; ganancia=20log10(Am); opt=stepDataOptions('StepAmplitude',10); figure(2); step(feedback(KcGp,1),opt)
clc clear all Mos=15; Kc=39.9910; %%compensador de ganancia s=tf('s'); nump=3; denp=s((7s)+3); Gp= nump/denp; figure(1) margin(KcGp) grid on %% Vamos a encontrar lo sigiente El máximo de obershut el cual es 0,095 por que es de -9.5 porciento zeta=sqrt(log(Mos/100)^2/(pi^2+log(Mos/100)^2)); El coeficiente de amortiguamiento pm=atand(2zeta/(sqrt(-2zeta^2+sqrt(1+4*zeta^4)))); El margen de fase p=-180+pm; Vamos a buscar en el diagrama de bode donde la fase Recordemos que la fase es de -180 mas la margen de fase el cual nos daría
p=-180+pm; [mag,phase,wout]=bode(KcGp); wm=interp1(squeeze(phase),wout,p); Am= interp1(squeeze(phase), squeeze(mag),p); alpha=10^(Am/20); wz=0.1wm; wp=wz/alpha; k=Kcwp/wz; Gc=k(s+wz)/(s+wp); Gol=series (Gp,Gc); opt=stepDataOptions('StepAmplitude',10); figure(2); step(feedback(Gol,1),opt) Como podemos ver en este no se estabilisa por lo tanto ahora utilizaremos la ultima opción
clc clear all Mos=15; tp=0.1 ; Kc=39.9910; %%compensador de ganancia s=tf('s'); nump=3; denp=s((7*s)+3);
Gp= nump/denp; figure(1) margin(KcGp) grid on zeta=sqrt(log(Mos/100)^2/(pi^2+log(Mos/100)^2)); wBW =Wbw(zeta,tp,'tp'); pm=atand(2zeta/(sqrt(-2zeta^2+sqrt(1+4zeta^4)))); [MG,MF,Wcg,Wcp]=margin(KcGp); pmax=pm+10-MF; beta = (1-sind(pmax))/(1+sind(pmax)); m=1/sqrt(beta); [mag,phase,wout]=bode(Gp); wmax=interp1(squeeze(mag),wout,p); T=1/(wmaxsqrt(beta)); Gc=(s+1/T)/(s+1/(betaT)); k=1/dcgain(Gc); Gol=series (KcGp,k*Gc); opt=stepDataOptions('StepAmplitude',10); figure(2); step(feedback(Gol,1),opt)
Mos=20; tp=0.1 ; kv=40; %%compensador de ganancia s=tf('s'); nump=100;
