Mensaje de error

Deprecated function: implode(): Passing glue string after array is deprecated. Swap the parameters en drupal_get_feeds() (línea 394 de /home1/montes/public_html/drupal/includes/common.inc).

plot

Undefined

Respuesta Transitoria

 

Español

Programa 7.5 OGATA 4edicion pag455

Vamos a dibujar la respuesta a un escalon para los valores de $K=0.4490$ y $K=1.4130$ mediante el Scilab del sistema en lazo cerrado:

 

Funcion de transferencia en lazo cerrado, transformada de Laplace


Programa en Scilab:
s=%s;

k1=0.4490;

g1=20/(s*(s+1)*(s+4)+20*(1+k1*s));

k2=1.4130;

g2=20/(s*(s+1)*(s+4)+20*(1+k2*s));

t=0:0.1:10;

y1=csim('step',t,g1);

y2=csim('step',t,g2);

clf;

plot(t,y1,'k');

plot(t,y2,'b');

xgrid;

xtitle('Respuesta a un escalon unitario de un sistema paralelo','t','salida');

legend(['K=0.4490';'K=1.4130'],style=4);
Respuesta del sistema a un escalon para dos valores de ganancia con Scilab
Español

Ejemplo 7.3 OGATA 4edicion pag444

Vamos a calcular el lugar de las raices mediante Scilab de un sistema en lazo cerrado con realimentación unitaria cuyo funcion de transferencia en lazo abierto es:Vamos a solucionar el problema del libro de otra forma. La funcion de transferencia en lazo abierto es:

 

Funcion de transferencia del sistema en lazo abierto


 

  1. Calculamos el polo dominante
    Tenemos el factor amortiguamient de 0.5 y la frecuencia natural no amortiguada 5. Por lo tanto el polo dominante tiene que estar en:

     

    polo dominante


     

  2. Calculamos el angulo a corregir:

     

    Valor de la funcion de transferencia para el polo dominante

     

     

    angulo a corregir

     

     

    180-angulo a corregir


     

     

    representacion de los angulos

 

 

 

  1. Calculamos $T_{1}$ y $\gamma$

    compensador para corregir el angulo


    Vamos a posicionar el cero en -2.5 .Alineado con el polo dominante.

     

    posicion de los polos y ceros del compensador

     

     

    distancia entre polos


     

    distancia entre el polo y el cero


     

    funcion de transferencia del compensador de adelanto


     

    Tenemos un cero en -2.5 y un polo en -8.63. Por lo tanto el compensador adelanto queda:

     

    \begin{displaymath}K_{c}\cdot \frac{s+\frac{1}{T_{1}}}{s+\frac{\gamma}{T_{1}}}=K_{c}\cdot \frac{s+2.5}{s+8.63} \end{displaymath}


     

    Por lo que podemos calcular el valor de $\gamma$

     

    ecuaciones para calcular la gamma

     

    valor de la gamma


     

  2. Vamos a calcular el valor de $K_{c}$

     

    ecuacion del valor de la ganancia del compensador de adelanto

     

     

    valor de la ganancia del compensador de adelanto


     

    Con lo que el compensador de adelanto queda:

     

    Funcion de transferencia del compensador de adelanto


     

  3. Vamos a calcular el $\beta$ del compensador de atraso. Sabiendo que $Kv=80$

     

    funcion de transferencia del compensador de atraso

     

    ganancia del compensador de atraso

     

     

    valor de la beta

     

  4. Escogiendo $T_{2}=5$ como en el libro vamos a comprobar que verifica las condiciones el compensador

     

    funcion de transferencia obtenida para el compensador de atraso


     

    valor del compensador de atraso para el polo dominante

     

    angulo del compensador de atraso para el polo dominante

     

     

    \begin{displaymath}-5º<-1.42<0º\end{displaymath}


     

    Con lo que verifica las condiciones. El compensador adelanto-atraso nos quedaria:

     

    Funcion de transferencia para el compensador adelanto atraso


     

Vamos a sacar las graficas de respuesta del sistema(el que calculamos,el del libro y el del ejemplo 7.4 del libro) a una entrada escalon . Tambien se mostrara la programacion de todos los calculos obtenidos.

Programa en Scilab:
s=%s;
g=4/(s*(s+0.5));
s1=-2.5+5*sqrt(1-0.5^2)*%i;
gs=syslin('c',g);
gs1=horner(gs,s1);
angulo=180-360*atan(abs(imag(gs1))/abs(real(gs1)))/(2*%pi);
angulocorregir=180-angulo;
l=imag(s1)*tan(2*%pi*angulocorregir/360);
p1=-2.5-l;
z1=-2.5;
gc=(s-z1)/(s-p1);
gma=p1/z1;
gt=gc*g;
aux1=(abs(horner(gt,s1)));
kc=1/aux1;
gct=kc*gc;
gt2=kc*gt;
aux3=s*gt2;
aux4=horner(aux3,0);
b=80/aux4;
gc2=(s+(1/5))/(s+(1/(5*b)));
aux5=horner(gc2,s1);
aux6=abs(aux5);
angulo2=-360*atan(abs(imag(aux5))/abs(real(aux5)))/(2*%pi);
gt3=gc2*gt2;
gct2=6.26*((s+0.5)/(s+5.02))*((s+0.2)/(s+0.01247));
gt4=g*gct2;
gct3=10*((s+2.38)/(s+8.34))*((s+0.1)/(s+0.0285));
gt5=g*gct3;
t=0:0.01:5;
glc=g /. 1;
glc1=gt3 /. 1;
glc2=gt4 /. 1;
glc3=gt5 /. 1;
y=csim('step',t,glc);
y1=csim('step',t,glc1);
y2=csim('step',t,glc2);
y3=csim('step',t,glc3);
clf;
//negro sistema sin compensar
plot(t,y,'k');
//verde, sistema compensado que se ha calculado
plot(t,y1,'g');
//azul,sistema compensado del libro ejemplo 7.3
plot(t,y2,'b');
//cyan sistema compensado del libro ejemplo 7.4
plot(t,y3,'c');

legend(['sin compensar';'compensado';'compensado libro 7.3';'compensado 
libro 7.4']);

xtitle('Respuesta a un escalon unitario de un sistema con compensacion de 
adelanto-atraso','t','salida');

xgrid;
respuesta del sistema compensado y no compensado a un impulso con Scilab

Ahora vamos a dibujar la respuesta a una rampa.

 

Añadimo al programa anterior en Scilab el siguiente codigo:
y=csim(t,t,glc);
y1=csim(t,t,glc1);
y2=csim(t,t,glc2);
y3=csim(t,t,glc3);
clf;
plot(t,t,'r');
//negro sistema sin compensar
plot(t,y,'k');
//verde, sistema compensado que se ha calculado
plot(t,y1,'g');
//azul,sistema compensado del libro ejemplo 7.3
plot(t,y2,'b');
//cyan sistema compensado del libro ejemplo 7.4
plot(t,y3,'c');

legend(['rampa';'sin compensar';'compensado';'compensado libro 7.3';'compensado
 libro 7.4'],style=4);

xtitle('Respuesta a una rampa unitaria de un sistema con compensacion de 
adelanto-atraso','t','salida');

xgrid;

respuesta del sistema compensado y no compensado a una rampa con Scilab
Español

Ejercicio B.6.15

[adsense:336x280:9156825571]
Vamos a dibujar el lugar de las raíces del sistema en lazo abierto y calcular las raíces en lazo cerrado cuando $K=2$

 

Funcion de transferencia en lazo abierto, Tranformada de Laplace

 

Programa en Scilab:
//Borramos la ventana grafica

clf;
//Declaramos la variable s y la funciones 'g','h' y 'g*h'

s=%s;

g=(s+1)/(s*(s^2+2*s+6));

h=1/(s+1);

gh=g*h;

//Declaramos la funcion g*h lineal y continua

ghs=syslin('c',gh)

//Declaramos ghs2 como la funcion anterior
 pero con K=2

ghs2=2*ghs;
[adsense:336x280:9156825571]
//Calculamos las raices de (1+ghs2)*(s+1), este ultimo es el polo que se anulo
 con el cero

glc=(1+ghs2)*(s+1);

r=roots(numer(glc))

//Dibujamos el lugar de las raices de la funcion ghs

evans(ghs);
mtlb_axis([-5 5 -5  5])

//Dibujamos el polo y el cero que se anularon:

plot(-1,0,'o');

plot(-1,0,'x');

//Dibujamos las raices para K=2;

plot(real(r),imag(r),'+');

xgrid;

xtitle('Lugar de las raices, los polos en lazo cerrado -> + para K=2','Eje real'
,'Eje imaginario');


Resultados
-->r=roots(numer(glc))
 r  =
 
  - 0.3706383              
  - 1.                     
  - 0.8146808 + 2.175406i  
  - 0.8146808 - 2.175406i



Lugar de las raices en lazo cerrado con Scilab
Español

Programa 7.3 OGATA 4edicion pag437

Vamos a dibujar la respuesta a un rampa mediante el Scilab de un sistema compensado y no compensado. El sistema en lazo abierto a compensar es el siguiente:

 

Funcion de transferencia del sistema en lazo abierto


 

El sistema de compensacion (compensador de atraso) obtenido es:

 

Funcion de transferencia del compensador de atraso

 

Programa en Scilab:
clf;

s=%s;

g=1.06/(s*(s+1)*(s+2));

gc=0.9656*(s+0.05)/(s+0.005);

gt=gc*g;

gc=g /. 1;

gct=gt /. 1;

t=0:0.1:50;

gs=syslin('c',gc);

gcs=syslin('c',gct);

y=csim(t,t,gs);

yt=csim(t,t,gcs);

plot(t,y);

plot(t,yt,'g');

plot(t,t,'r');

xgrid;

xtitle('Respuesta a una rampa unitaria de un sistema con compensacion de atraso','t'
,'salida');

legend(['no compensado';'compensado';'rampa'],style=4);

Respuesta en el tiempo del sistema a una entrada rampa unitaria con Scilab
Español

Programa 7.1 OGATA 4edicion pag428

Vamos a dibujar mediante el Scilab un sistema compensado y no compensado a una entrada escalon. El sistema en lazo abierto a compensar es el siguinte:

 

Funcion de transferencia en lazo abierto, Transformada de Laplace


 

El sistema de compensacion (compensador de adelanto) obtenido es:

 

Compensador de adelanto del sistema


Programa en Scilab:
clf;
 
s=%s;
 
g=4/(s*(s+2));
 
gs=4.68*(s+2.9)/(s+5.4);
 
gc=g/. 1;
 
gt=gs*g;
 
gtc=gt/. 1;
 
t=0:0.05:5;
 
gcs=syslin('c',gc);
 
gtcs=syslin('c',gtc);
 
y=csim('step',t,gcs);
 
yt=csim('step',t,gtcs);
 
plot(t,y);
 
plot(t,yt,'g');
 
xgrid;
 
legend(['no compensado';'compensado'],style=4);


 
Salida del sistema compensado y no compensado a una entrada escalon con Scilab
Español

Problema A6.15 OGATA 4edicion pag407

Vamos a dibujar el lugar de las raices mediante Scilab de un sistema en lazo cerrado con realimentación unitaria cuyo funcion de transferencia en lazo abierto es:

 

Funcion de transferencia, Transformada de Laplace


 

Programa en Scilab

s=%s;

num=(s^2+25)*s;

den=(s^4+404*s^2+1600);

g=syslin('c',num/den);

evans(g)

v=[-20 20 -21 21];

mtlb_axis(v)
Lugar de las raices del sistema con Scilab

El resultado que obtendremos no nos valdra la grafica esta deformada tendremos que programarlo de la siguiente manera

 
Programa en Scilab

s=%s;

num=(s^2+25)*s;

den=(s^4+404*s^2+1600);

g=syslin('c',num/den);

for k1=0.2:0.2:20,

  gl=1+k1*g;

  numl=numer(gl);

  r=roots(numl);

  x=real(r);

  y=imag(r);

  plot(x,y,'o');

end;

for k2=20.02:0.2:30,

  g2=1+k2*g;

  num2=numer(g2);

  r2=roots(num2);

  x2=real(r2);

  y2=imag(r2);

  plot(x2,y2,'o');

end;

for k3=35:5:1000,

  g3=1+k3*g;

  num3=numer(g3);

  r3=roots(num3);

  x3=real(r3);

  y3=imag(r3);

  plot(x3,y3,'o');

end;

v=[-20 20 -21 21];

mtlb_axis(v)

xtitle('lugar de las raices','Eje Real','Eje imaginario')

xgrid;

 

Lugar de las raices con puntos del sistema con Scilab

 

Español

Programa 6.3 OGATA 4edicion pag362 (Lugar de las Raices)

Vamos a calcular el lugar de las raices mediante Scilab de un sistema en lazo cerrado con realimentación unitaria cuyo funcion de transferencia en lazo abierto es:

 

Funcion de tranferencia del sistema. Laplace


Programa en Scilab

num=poly([1 0 0 0 0],'s','coeff');

den=poly([0 5 10.3 1.1 1],'s','coeff');

g=syslin('c',num/den);

for k1=0.2:0.2:20,
  gl=1+k1*g;
  numl=numer(gl);
  r=roots(numl);
  x=real(r);
  y=imag(r);
  plot(x,y,'o');
end;

for k2=20.02:0.2:30,
  g2=1+k2*g;
  num2=numer(g2);
  r2=roots(num2);
  x2=real(r2);
  y2=imag(r2);
  plot(x2,y2,'o');
end;

for k3=35:5:1000,
  g3=1+k3*g;
  num3=numer(g3);
  r3=roots(num3);
  x3=real(r3);
  y3=imag(r3);
  plot(x3,y3,'o');
end;

v=[-20 20 -25 25];

mtlb_axis(v)

xgrid;

Lugar de las raices por puntos de la trasformada de Laplace con Scilab

Español

Páginas

Pin It
Subscribe to RSS - plot