Trabajos1

Matlab ejemplo.

'(ch4p1) Example 4.6'
ans =

(ch4p1) Example 4.6
% desplega etiqueta de titulo%

p1=[1 3+7*i];
p2=[1 3-7*i];
% definicion de los polinomios con el 1er y 2do polo respectivamente%

deng=conv(p1,p2);
%multiplicacion de los polinomios para determinar el polinomio de 2do
%grado%

omegan= sqrt(deng(3)/deng(1))

omegan =

7.6158
%calculo de la frecuencianatural no amortiguada, sqrt(c/a)%

zeta=(deng(2)/deng(1))/(2*omegan)

zeta =

0.3939
%factor de amortiguamiento relativo, ((b/a)/2*wn).%

Ts=4/(zeta*omegan)

Ts =

1.3333
%calculo tiempo de asentamiento, (4/z*wn)%

Tp=pi/(omegan*sqrt(1-zeta^2))

Tp =

0.
%calculo tiempo pico, pi/wn*sqrt(1-z^2).%

pos=100*exp(-zeta*pi/sqrt(1-zeta^2))

pos =26.0176
%calculo sobrepaso en porcentaje (100*e^(z*pi/sqrt(1-z^2))

pause



'(ch4p2) Ejemplo 4.8'

ans =

(ch4p2) Ejemplo 4.8


'corrida de prueba'

ans =

corrida de prueba
%despliegan etiqueta de titulo%
clf
% borrar grafica%

numt1=[24.542];
%definir numerador de T1%

dent1=[1 4 24.542];
%def. denominador de T1%

'T1(s)'

ans =

T1(s)%desplegar etiqueta%

T1=tf(numt1,dent1)

Transfer function:
24.54
-----------------
s^2 + 4 s + 24.54
%crear y desplegar T1(s)%


step(T1)
%correr demostracion de la grafica a la respuesta escalon%

title('corrida de prueba de T1(s)')
%agregar titulo a la grafica%

pause

'corrrida completa'

ans =

corrrida completa
% desplegar etiqueta%[y1,t1]=step(T1);
%correr la respuesta escalon de T1%

numt2=[245.42];
%definicion numerador T2%
p1=[1 10];
%definir (s+10) en denominador de T2%

p2=[1 4 24.542];
%definir (s^2+4s+24.542)en denom. T2%

dent2=conv(p1,p2);
% multiplicar (s+10)(s^2+4s+24.542) po el denominador T2%

'T2(s)'

ans =

T2(s)

T2=tf(numt2,dent2)

Transfer function:
245.4------------------------------
s^3 + 14 s^2 + 64.54 s + 245.4
%crear y desplegar funcion transf. T2%

[y2,t2]=step(T2);
%correr la respuesta escalon de t1%
numt3=[73.626];
%definir numerador T3%

p3=[1 3];
%definir (s+3) en el denominador p3%

dent3=conv(p3,p2);
%mult. (s+3)(s^2+4s+24.542) por el denominador t3%

'T3(s)'

ans =

T3(s)

T3=tf(numt3,dent3)

Transfer function:73.63
-----------------------------
s^3 + 7 s^2 + 36.54 s + 73.63
%crear y desplegar ft de T3%

[y3,t3]=step(T3);
%correr resp. escalon T3%

clf
plot(t1,y1,t2,y2,t3,y3)
title('Respuestas escalón de T1(s),T2(s), y T3(s)')
xlabel('Tiempo(segundos)')
ylabel('Respuesta normalizada')
text(0.7,0.7,'c3(t)')
text(0.7,1.1,'c2(t)')
text(0.5,1.3,'c1(t)')
%graficar los puntos de las 3graficas sobre los mismos ejes, adicional
%titulo a la grafica, adicionar etiqueta del tiempo y del eje de respuesta,
%etiquetar las respuestas ecalon respectivamente%

pause
step(T1,T2,T3)
title('Respuesta escalon de T1(s),T2(s), y T3(s)')
%adicionar titulo a la grafica%

'(ch4p3)'

ans =

(ch4p3)

clf
%borrar grafica%

A=[0 1 0;0 0 1;-24 -26 -9];
B=[0;0;1];
C=[2 7 1];D=0;
%generar la matriz A, vector B, vector C y D respectivamente%

T=ss(A,B,C,D)

a =
x1 x2 x3
x1 0 1 0
x2 0 0 1
x3 -24 -26 -9

b =
u1
x1 0
x2 0
x3 1

c =
x1 x2 x3
y1 2 7 1

d =
u1
y1 0


Continuous-time model.
%generar el objeto LTI,T, en el espacio de estado ydesplegarlo%

t=0:0.1:10;
%intervalo de tiempo para la grafica%

grid on
%dibujar reticula sobre la grafica%

step(T,t)
%grafica resp. escalon para el intervalo de tiempo dado%

pause
'(ch4p4) Estudio de caso del control de la antena'

ans =

(ch4p4) Estudio de caso del control de la antena

clf
numg=20.83;
deng=[1 101.71 171];
%definir numerador y denominador G(s)...