torres Valero Ricardo, Corredor Alipio Gustavo
Estudiantes séptimos Semestre Ing. Electrónica.
>> numg=20*[1 5] EN ESTA LINEA PRESENTAMOS LOS COEFICIENTES DEL NUMERARDOR
numg =
20 100
>> deng=poly([0 -1 -4]) EN ESTA LINEA PRESENTAMOS LOS COEFICIENTES DEL POLINOMIO DEL DENOMINADOR
deng =
1 5 4 0
>> G=tf(numg, deng) EN ESTA LINEA HALLAMOS LA FUNCION DE TRANFERENCIA
Transfer function:
20 s + 100
-----------------
s^3 + 5 s^2 + 4 s
>> pos=input('type desired setting time%'); EN ESTA LINEA HALLAMOS EL COEFICIENTE DE AMORTIGUAMIENTO DEL SISTEMA
type desired setting time%95
>> ts=input('tipe value'); EN ESTA LINEA SE HALLA EL TIEMPO ESTABLE DEL SISTEMA
tipe value0.75
>> z=(-log(pos/100))/(sqrt(pi^2+log(pos/100)^2));
>> wn=4/(z*ts) EN ESTA LINEA SE COLOCA LA FRECUENCIA ANGULAR
wn =
326.6976
>> [num, den]=ord2(wn, z)
num =
1
den =
1.0e+005 *
0.0000 0.0001 1.0673
>> r=roots(den) SE HALLA LAS RAICES GENERALES O LOS POLOS DE LA FUNCION DE TRANSFERENCIA
r =
1.0e+002 *
-0.0533 + 3.2665i
-0.0533 - 3.2665i
>> [Ac, Bc, Cc, Dc]=tf2ss(numg, deng) SE HALLA LA MATRIZ Ac, Bc, Cc, Dc
Ac =
-5 -4 0
1 0 0
0 1 0
Bc =
1
0
0
Cc =
0 20 100
Dc =
0
>> p=[0 0 1; 0 1 0; 1 0 0] SE INGRESA UNA MATRIZ
p =
0 0 1
0 1 0
1 0 0
>> Ap=inv(p)*Ac*p
Ap =
0 1 0
0 0 1
0 -4 -5
LO ANTERIOSR FUE PARA HALLAR LA FUNCION DE TRANSFERENCIA EL LUGAR DE POLOS Y OBSERVAR LA ESTABILIDAD DE UN SISTEMA.
2) HALLAR LAS RAICEZ DE LA FUNCION DE TRANSFERENCIA:
numf=20*[1 5]
numf =
20 100
>> denf=poly([0 -1 -4 ])
denf =
1 5 4 0
>> denf=poly([0 -1 -4 ]);
>> denf=poly(0 -1 -4 )
denf =
1 5
>> den=poly(0 -1 -4 )
den =
1 5
>> ca= conv([1 0 ], [1 1 ]) SE HACE LA CONVOLUCION
ca =
1 1 0
>> cb=conv(ca, [1 4]) SE HACE OTRA CONVOLUCION
cb =
1 5 4 0
>> res=(cb+ [0 0 20 100]) SE HALLA LOS COEFICIENTES DE LA FUNCION DE TRANSFERENCIA
res =
1 5 24 100
>> m=roots(res) EN ESTA LINEA SE REPRESENTA LAS RAICEZ O SE HALLA LAS RAICEZ
m =
-4.5529
-0.2236 + 4.6813i
-0.2236 - 4.6813i
LO ANTERIOR ES PARA HALLAR EL LUGAR DE LAS RAICEZ DE LA FUNCION DE TRANSFERENCIA
