Universidad Tecnológica de Santiago (UTESA)
Taller de Controles de Sistemas Automáticos
Guía de la Práctica No. 04
Matemática Simbólica
Nombre: Matrícula:
Docente: José M. Solís Clave de la Asignatura: _IEL – 905 – 003 _ Ciclo: 01 – 2025 _
Fecha de inicio: / / Fecha de término: / / Duración: 90 minutos
1. Objetivo
• Utilizar la matemática simbólica para obtener las funciones de
transferencias
• Aplicar la matemática simbólica para resolver ecuaciones simultáneas
• Utilizar el concepto de polinomio en operaciones básicas de matemáticas.
Matemática simbólica
Matlab permite la realización de cálculos de manera simbólica si previamente se
definen símbolos las variables y funciones con las que se va a operar. Todos los
comandos se encuentran en la toolbox de matemática simbólica, incorporando un
nuevo tipo de dato en Matlab: el objeto simbólico (sym object). Las variables se
han de definir como símbolos mediante syms. Las funciones se definen
directamente con el comando sym.
Con variables y funciones simbólicas se puede realizar todo tipo de cálculo:
diferencial, integral y matricial; también es posible la resolución de ecuaciones
simbólicas y diferenciales, y la obtención de distintas transformadas como Laplace
o transformada Z.
2. Equipos
Equipo Característica Cantidad
PC Instalado con programa MATLAB 1 por estudiante
3.Indicaciones
• Los estudiantes trabajan individualmente.
• Cada estudiante presentará un informe de la práctica de laboratorio en
la próxima sesión.
Operación con variables simbólicas
Dadas las funciones de transferencia
Obtener:
a) Suma de G1 y G2
>> s=sym('s');
>> G1=1/(s^2+9*s+20);
>> G2=(s+2)/(s^3+15*s^2+71*s+105);
>> sum_G=G1+G2
sum_G =
1/(s^2 + 9*s + 20) + (s + 2)/(s^3 + 15*s^2 + 71*s + 105)
>> pretty(sum_G)
1 s + 2
------------- + -----------------------
2 3 2
s + 9 s + 20 s + 15 s + 71 s + 105
b) Producto de G1 y G2
>> prod_G=(G1*G2)
prod_G =
(s + 2)/((s^2 + 9*s + 20)*(s^3 + 15*s^2 + 71*s + 105))
>> pretty(prod_G)
s + 2
-----------------------------------------
2 3 2
(s + 9 s + 20) (s + 15 s + 71 s + 105)
c) solve(G1)
>> solve(G1)
ans =
Empty sym: 0-by-1
d) solve(G2)
>> solve(G2)
ans =
-2
2. Utilizar el comando pzmap para graficar los polos y ceros en el
plano s de las expresiones siguientes:
>> p=[1 4];
>> d=[1 6 17 13];
>> roots(p)
ans =
-4
>> roots(d)
ans =
-2.4344 + 2.3593i
-2.4344 - 2.3593i
-1.1312 + 0.0000i
>> pzmap(p,d)
>> num=10*[1 2];
>> den=conv([1 0 0], conv([1 1], [1 10]));
>> roots(num)
ans =
-2
>> roots(den)
ans =
0
0
-10
-1
>> pzmap(num,den)
>> num=10*[1 2];
>> den=conv([1 0], [1 25 2]);
>> roots(num)
ans =
-2
>> roots(den)
ans =
0
-24.9197
-0.0803
>> pzmap(num,den)
>> num=10*conv([1 0], [1 1]);
>> den=conv([1 2], [1 35 2]);
>> roots(num)
ans =
0
-1
>> roots(den)
ans =
-34.9428
-2.0000
-0.0572
>> pzmap(num,den)
3. Realice la conversión de los polinomios siguientes a polinomios
numéricos utilizando el comando sym2poly
>> syms s
>> G=10*(s+2)/(s*(s^2+2*s+2);
>> G=10*(s+2)/(s*(s^2+2*s+2));
>> [n,d]=numden(G)
n =
10*s + 20
d =
s*(s^2 + 2*s + 2)
>> num_1=sym2poly(n)
num_1 =
10 20
>> den_1=sym2poly(d)
den_1 =
1 2 2 0
>> syms s
>> G=(10*s*(s+1))/((s+2)*(s^2+3*s+2));
>> [n,d]=numden(G)
n =
10*s
d =
(s + 2)^2
>> num_1=sym2poly(n)
num_1 =
10 0
>> den_1=sym2poly(d)
den_1 =
1 4 4
𝑛𝑢𝑚 = [1 0 −20.6];
4. Dado los polinomios
𝑑𝑒𝑛 = [1 19.6 151.2];
4.1 𝒄 = 𝒄𝒐𝒏𝒗(𝒂, 𝒃)
Obtenga:
>> num= [1 0 -20.6];
>> den=[1 19.6 151.2];
>> c=conv(num,den)
c =
1.0e+03 *
0.0010 0.0196 0.1306 -0.4038 -3.1147
4.2 N1=poly2sym(num,'s')
>> syms s;
>> N1= poly2sym(num,s)
N1 =
s^2 - 103/5
4.3 D1=poly2sym(den,'s')
>> syms s;
>> D1=poly2sym(den,s)
D1 =
s^2 + (98*s)/5 + 756/5
4.4 G=N1/D1
>> G=N1/D1
G =
(s^2 - 103/5)/(s^2 + (98*s)/5 + 756/5)
>> pretty(G)
2 103
s - ---
5
---------------
2 98 s 756
s + ---- + ---
5 5
