Fundamento Matemático de los Sistemas de Control (I)

1. Mediante el MATLAB, CC5, VISSIM, realizar los mismos pasos del desarrollo para la solución de la siguiente función de transferencia

- Para el programa MATLAB tenemos:

clc
%Lab.2
%Name: Chávez Choque Alvaro
%Adress: Pedro Ferrari #4, Toledo, Madrid
%Thelephone: 5264639
%I.N. 5775394 or.
%e-mail: alvaro.chavez_choque@hotmail.com
%Desarrollo:
'Introducción de la función de transferencia'
num=[1 4 3]% numeradores de la función de transferencia
den1=[1 1.5]
den2=[1 0.3]
den3=[1 3 3]
den4=conv(den1,den2)
den=conv(den3,den4)%denominador de la función de transferencia
%conversión de la función de transferencia
G=tf(num,den)
display (G)
'a)polos y ceros'
[ceros,polos,ganancia]=tf2zp(num,den)
'b)Expansión de fracciones parciales de la función de transferencia'
[residue,polos,directo]=RESIDUE(num,den)
'c)la transformada de la salida, si la entrada es un escalón de amplitud 10'
%introducción de la entrad r(s)
numr=10
denr=[1 0]
R=tf(numr,denr)
'la salida es'
Y=G*R
'd)expansión en fracciones parciales de la función Y(s)'
'conversión de Y(s)a vectores'
[numY,denY]=tfdata(Y,'v')
[residueY,polosY,directoY]=RESIDUE(numY,denY)
'e)la inversa de transformada de Laplace'
syms g y r s
g=(s^2+4*s+3)/((s+1.5)*(s+0.3)*(s^2+3*s+3))
pretty(g)
r=10/s
pretty(r)
y=g*r
pretty(y)
pretty(simplify(ilaplace(y)))
'f) La inversa en el tiempo'
%modificación de la amplitud
G1=10*G
subplot(1,2,1)
step(G1)
%modification
subplot(1,2,2)
step(G1,10)
grid
title('respuesta del sistema')
'g)El modelo en variables de estado'
[A,B,C,D]=tf2ss(num,den)
'h)A partir de variables de estados, encontrar la función de transferencia'
[numc,denc]=ss2tf(A,B,C,D)
Gc=tf(numc,denc)
display(Gc)