1. Considere los dos polinomios:
p(s) = s2+2s+1
q(s) = s+1
Calcule lo siguiente:
a). p(s).q(s)
b). Los polos y ceros de:
c). p(-1)
2. Considere el sistema realimentado mostrado en la figura:
a). Calcule la función de transferencia de lazo cerrado.
b). Obtenga la respuesta al impulso y la respuesta al escalón del sistema.
Grafique dichas respuestas y halle el valor final en el caso de la respuesta al escalón.
3. Considere la ecuación diferencial:
Donde 
y u(t) es un escalón unitario. Determine la solución y(t) solucionando por transformada inversa de Laplace y por el comando step. Compare los resultados, dibujando simultáneamente las dos respuestas.
4. Un sistema de control de altitud de eje sencillo de una nave espacial puede ser representado por el diagrama de bloques de la siguiente figura.
Para ver el gráfico seleccione la opción "Descargar" del menú superior
Las variables k, a y b son parámetros del controlador, y J es el momento de inercia de la nave espacial. Suponga que el momento de inercia es J=10.8×108 (slug-ft2), y los parámetros de control son k=10.8×108, a=1, b=8.
a). Encuentre la función de transferencia de lazo cerrado 
b). Encuentre la respuesta al escalón del sistema cuando la entrada corresponde a un escalón de 10º.
c). El momento de inercia exacto es generalmente desconocido y puede variar ligeramente con el tiempo. Compare el desempeño de la respuesta al escalón de la nave espacial cuando J se reduce en un 20% y un 50%. Utilice los mismos parámetros de control anteriores y la misma entrada escalón de 10º. Analice los resultados.
5. Considere el diagrama de bloques mostrado en la siguiente figura.
Para ver el gráfico seleccione la opción "Descargar" del menú superior
a). Reduzca el diagrama de bloques y encuentre la función de transferencia de lazo cerrado.
b). Genere el diagrama de polos y ceros de la función de transferencia de lazo cerrado.
c). Determine explícitamente los polos y ceros de la función de transferencia de lazo cerrado utilizando "roots". Compare con los resultados obtenidos en el numeral anterior.
6. Un sistema tiene función de transferencia:
Dibuje la respuesta del sistema cuando r(t) equivale a un impulso unitario para el parámetro z = 3, 6, 9, 12.
Programa que desarrolla el taller
close all
%EJERCICIO 1
%p(s).q(s)
a1=conv([1 2 1],[1 1]);
%Los polos y ceros de G(s)=q(s)/p(s)
g=tf([1 1],[1 2 1]);
b1=zpk(g);
p1=b1.p{:}; %Polos de G
z1=b1.z{:}; %Ceros de G
%p(-1)
c1=polyval([1 2 1],(-1));
%EJERCICIO 2
%Función de transferencia de lazo cerrado
a2=feedback((series((tf(1,[1 1])),(tf([1 2],[1 3])))),1,'+');
%Respuesta al impulso
b2i=impulse(a2,[1:11]);
figure
plot([1:11],b2i), title('Respuesta al impulso'), grid, pause;
b2e=step(a2,[9:15]);
figure
plot([9:15],b2e), title('Respuesta al escalón'), grid, pause;
%EJERCICIO 3
%Por transformada inversa de Laplace
syms ys s
ys=1/(s*(s+1)^2);
yt=ilaplace(ys);
pretty(ys)
pretty(yt), pause;
figure
ezplot(yt,[0:0.1:6]), grid, pause;
%Por el comando step
ys1=tf(1,[1 2 1]);
yt1=step(ys1,[0:0.1:6]);
figure
plot([0:0.1:6],yt1), title('Respuesta al escalón'), grid, pause;
%EJERCICIO 4
k=10.8e8; a=1; b=8; j=10.8e8;
a4=tf(feedback((series((k*zpk(-a,-b,1)),((1/j)*zpk([],[0 0],1)))),1,'+'))
disp('Función de transferencia de lazo cerrado'), pause;
%Respuesta al escalón de 10°
b4s=a4*10;
b4=step(b4s,[0:0.1:6]);
figure
plot([0:0.1:6],b4), title('Respuesta al escalón de 10°'), grid, pause;
%Cuando j se reduce en un 20%
a4j20=tf(feedback((series((k*zpk(-a,-b,1)),((1/(j*0.8))*zpk([],[0 0],1)))),1,'+'));
b4sj20=a4j20*10;
b4j20=step(b4sj20,[0:0.1:6]);
figure
plot([0:0.1:6],b4j20), title('Respuesta al escalón de 10° con 0.8j'), grid, pause;
a4j50=tf(feedback((series((k*zpk(-a,-b,1)),((1/(j*0.5))*zpk([],[0 0],1)))),1,'+'));
b4sj50=a4j50*10;
b4j50=step(b4sj50,[0:0.1:6]);
figure
plot([0:0.1:6],b4j50), title('Respuesta al escalón de 10° con 0.5j'), grid, pause;
%EJERCICIO 5
m5=tf(1,[1 0 0]);
a5=tf([1 0 2],[1 0 0 14]);
u5=tf([4 2],[1 2 1]);
r5=tf([1 0],[1 0 2]);
o5=tf(1,[1 1]);
t5=series(4,(feedback((series((feedback(series(o5,r5),u5,'+')),(feedback(m5,50)))),a5,'+')));
[num5,den5]=tfdata(t5,'v');
%Diagrama de polos y ceros
figure
zplane(num5,den5), title('Diagrama de polos y ceros'), grid, pause;
%Polos y ceros con el comando roots
z5=roots(num5) %Ceros de la función de transferencia
disp('Ceros de la función de transferencia'), pause;
p5=roots(den5) %Polos de la función de transferencia
disp('Polos de la función de transferencia'), pause;
%EJERCICIO 6
for y=3:3:12
t6=tf([(15/y),15],[1 3 15])
t6imp=impulse(t6,[0:0.1:6]);
figure
plot([0:0.1:6],t6imp), title('Respuesta al impulso unitario'), grid, pause;
end
Resultados obtenidos en MATLAB
1. a). Transfer function P(s).Q(s):
s^3 + 3 s^2 + 3 s + 1
b).
Polos de G(s).
Ceros de G(s).
c). P (-1) = 0
2. a). Función de transferencia de lazo cerrado:
b). Respuesta al impulso:
Respuesta al escalón:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
3. Análisis:
y(t) = 1+(-t-1)*exp(-t)
Por transformada inversa de Laplace:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
Por el comando Step:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
4. a). Función de transferencia T(s):
b). Respuesta al escalón de 10°:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
c). Respuesta al escalón de 10° con 0.8J:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
Respuesta al escalón de 10° con 0.5J:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
5. a). Función de transferencia de lazo cerrado:
b). Diagrama de polos y ceros:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
c). Polos de la función de transferencia:
Ceros de la función de transferencia:
6. Respuesta del sistema cuando Z=3:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
Respuesta del sistema cuando Z=6:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
Respuesta del sistema cuando Z=9:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
Respuesta del sistema cuando Z=12:
Para ver el gráfico seleccione la opción "Descargar" del menú superior
Mauro Baquero
Rigoberto Hernando Olarte
ING Mecatronico. BUCARAMANGA – SANTANDER – COLOMBIA
2005