NIT.

890500622 - 6

www.ufps.edu.co
OPERACIONES MATEMATICAS UTILIZANDO
MATLAB
1. Objetivos.
 Aprender a utilizar los comandos de Matlab para la resolver operaciones matemáticas,
que permitan la solución de problemas orientados al control.
 Aprender a utilizar la herramienta de Matlab para el análisis matemático que soportarán
el diseño de sistemas de control.

2. Elementos necesarios
 Computador de ultimas especificaciones
 Cuaderno de notas
 MATLAB 2017

3. Introducción
La simulación es una disciplina que provee la estructura para resolver problemas desafiantes
de la mayoría de las áreas de aplicación en casos donde la experimentación y la experiencia
son necesarias. Ambas características se pueden ofrecer bajo todos los parámetros reales y
extremos. La simulación es una experimentación dirigida a objetivos con modelos de
sistemas dinámicos. Con la habilidad para proveer experiencia, la simulación juega un papel
importante en la capacitación al proporcionar posibilidades para desarrollar y/o mejorar la
competencia en uno de los tres tipos de habilidades; Habilidades motrices, la toma de
decisiones y las habilidades de comunicación, y las habilidades operativas.
Por tal motivo el modelamiento y la simulación computacional son aspectos muy importantes
de la ingeniería moderna aplicada, ya que son vitales en todas las áreas de la ciencia y la
ingeniería basadas en la simulación.[1]

4. Procedimiento

En esta práctica se bosquejaran algunas operaciones a fines con el análisis de los sistemas
de control utilizando la herramienta Matlab.

4.1 Matrices:
Una matriz es un arreglo rectangular de números:

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co

1 3 5 2 4 6
𝐴=[ ] 𝐵=[ ]
7 9 11 8 10 12
2
𝐶 = [1 3 5] 𝐹 = [4]
6
Operaciones con Matrices

Sintaxis Definición

D = A+B Suma de matrices.

D = A-B Resta de matrices.

D = A*B Multiplicación de matrices

Retorna la multiplicación de
D = k*A un escalar k por cada elemento
de la matiz A.
Devuelve un vector fila que
Size (A) corresponde con la longitud de
la dimensión de A
Retorna una matriz identidad
I = eye (n) nxn con unos en la diagonal
principal.
Intercambia las filas por las
columnas. Si A es compleja
B = A.’
devuelve la traspuesta no
conjugada de A.
Devuelve el determinante de la
D = det (A)
matriz cuadrada A.
Calcula la inversa de la
D = inv (A)
matriz cuadrada A.
Multiplica los arreglos A y B
D = A.*B
elemento por elemento.
Devuelve el rango de la matriz
D = rank (A)
A.

 Suma y Resta: Esta operación entre matrices solo se puede realizarse, si A = B.

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co
D = A+B ó E=A – B

1+2 3+4 5+6 3 7 11

𝐷 =𝐴+𝐵 = [ ]=[ ]
7+8 9 + 10 11 + 12 15 19 23
1−2 3−4 5−6 −1 −1 −1
𝐸 =𝐴−𝐵 =[ ]=[ ]
7 − 8 9 − 10 11 − 12 −1 −1 −1

MATLAB

%%%%%% Suma y Resta %%%%%%%%%

clc D=
Clear 3 7 11
A= [1 3 5; 7 9 11]; 15 19 23
B= [2 4 6; 8 10 12]; E=
D = A+B % Suma -1 -1 -1
E = A-B % Resta -1 -1 -1

 Multiplicación:

Multiplicación por un escalar. Si k es un escalar, entonces la multiplicación de la matriz

A por un escalar, es la multiplicación de k por cada uno de los elementos de A.

Multiplicación de Matrices: Se dice que dos matrices son conformes para la multiplicación
A*B, cuando el número de columnas de la matriz A es igual al número de filas de la matriz
B. Esto es, el producto de A*B (no B*A), es conformable para la multiplicación únicamente
si A es una matriz de m x p y la matriz B es p x n, entonces el producto será una matriz
m x n. Así la multiplicación es fila por columna.[2]

Muestra que A y B son conformables Indica la dimensión

para la multiplicación del producto B*A

Indica la dimensión del producto A*B

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co
2
𝐶 ∗ 𝐹 = [1 3 5 4] = [1 ∗ 2 + 3 ∗ 4 + 5 ∗ 6] = 44
] [
6
𝐷𝑖𝑚𝑒𝑛𝑠𝑖ó𝑛 𝑑𝑒𝑙 𝑎𝑟𝑟𝑒𝑔𝑙𝑜 𝑓𝑖𝑛𝑎𝑙 1𝑥3 3𝑥1 = 1𝑥1

2 2∗1 2∗3 2∗5 2 6 10

𝐹 ∗ 𝐶 = [4] [1 ]
3 5 = [4 ∗ 1 4 ∗ 3 4 ∗ 5] = [4 12 20]
6 6∗1 6∗3 6∗5 6 18 30
𝐷𝑖𝑚𝑒𝑛𝑠𝑖ó𝑛 𝑑𝑒𝑙 𝑎𝑟𝑟𝑒𝑔𝑙𝑜 𝑓𝑖𝑛𝑎𝑙 3𝑥1 1𝑥3 = 3𝑥3

MATLAB
%%%%%% Multiplicación %%%%%%%%% D=
10 15 30
k=5; 35 45 55
A = [2 3 6; 7 9 11; 13 15 18]% Matriz A 65 75 90
C = [1 3 5]; E=
F = [2; 4; 6]; 44
D = k*A % Multiplicación por constante H=
E = C*F % Multiplicación C por F 2 6 10
H = F*C % Multiplicación F por C 4 12 20
6 18 30

 Matriz escalar: Una matriz con un escalar se llama matriz diagonal. Si ese escalar
es uno (1), se llama matriz identidad.

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co

MATLAB
L=
%%%%%% Matriz Identidad %%%%%%%%% 1 0
0 1
P= [1 3 5; 7 9 11; 13 15 18]; M=
L = eye (2) % Matriz identidad de 2x2 1 0 0 0
0 1 0 0
M = eye (4) % Matriz identidad de 4x4
0 0 1 0
N = eye (size (P)) % matriz identidad 0 0 0 1
del tamaño de P
N=
1 0 0
0 1 0
0 0 1

Matriz traspuesta: Es la matriz que se obtiene de intercambiar las filas por las columnas de
una matriz.

MATLAB
%%%%%% Matriz Traspuesta %%%%%%%%% B=
2 4 6
8 10 12
B= [2 4 6; 8 10 12] % matriz B
Q=
Q=B' % Traspuesta de la matriz B 2 8
4 10
6 12

Determinantes: Los determinantes permiten obtener la solución de ecuaciones simultáneas.

En nuestro caso particular nos permite establecer si una matriz es invertible o no. Una matriz

MATLAB
%%%%%% Determinante %%%%%%%%%

A= [1 3 5; 7 9 11; 13 15 18]; % R=
Matriz A
R = det(A) % determinante de A -12.0000

es invertible si el determinante de ella es ≠ 0. Para matrices mayores o iguales a 3 x 3 se

resuelve por el método de los Cofactores. El determinante de una matriz cuadrada A de
orden n se denomina como el determinante de orden n.

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co
Matriz inversa (𝐴−1 ): Una matriz cuadrada A es invertible si su determinante es ≠ 0. Y se
obtiene así.

𝑎𝑑𝑗 𝐴
𝐴−1 =
| 𝐴|

La matriz inversa nos permite solucionar un sistema de ecuaciones simultáneas de la forma

AX = B en 𝑿 = 𝑨−𝟏 ∗ 𝑩 así:

2𝑥1 + 3𝑥2 + 𝑥3 = 9 2 3 1 𝑥1 9
𝑥
Dado el sistema: [ 1 + 2𝑥2 + 3𝑥3 = 6] = [1 2 𝑥
3] [ 2 ] = [6]
3𝑥1 + 𝑥2 + 2𝑥3 = 8 3 1 2 𝑥3 8

MATLAB
%%%%%% Matriz Inversa %%%%%%%%% x= w= z=
1.9444 1.9444 1.9444
A = [2 3 1; 1 2 3; 3 1 2] % Matriz A 1.6111 1.6111 1.6111
B = [9; 6; 8] % Matriz B 0.2778 0.2778 0.2778
x = inv (A)*B % solución del sistema
w = A^-1*B %Otra forma - Multiplicando
z = A\B % Dividiendo B entre A

Matriz adjunta (𝑎𝑑𝑗 𝐴): La matriz adjunta de una matriz cuadrada A, es la traspuesta de
la matriz de los cofactores de A.

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co
Por tanto 𝑎𝑑𝑗 𝐴 = (𝐴−1 ∗ det 𝐴)

MATLAB
%%%%%% Matriz adjunta %%%%%%%%% K=
-2
b = [1 2 3; 1 3 4; 1 4 3] % Matriz b inv_B =
a = [4 3; 2 2] % Matriz a
3.5000 -3.0000 0.5000
-0.5000 0 0.5000
K=det (b) % determinante de la Matriz b -0.5000 1.0000 -0.5000
inv_B=inv (b) % Matriz inversa de b
adj_B = (K*inv_B) % adjunta de b adj_B =

Otra forma: -7 1 1
6 0 -2
-1 -1 1
inversa_B2= b^-1

Multiplicación de dos vectores:

MATLAB
%%% Multiplicación de Vectores %%

A = [1 3 5 7]; % Vector A mul_AB =

B = [-3 2 6 -9]; % Vector B
mul_AB = A.*B
-3 6 30 -63

Rango de una Matriz

El rango de una matriz es el número máximo de filas y columnas que son linealmente
independientes de una matriz. Si una matriz A tiene un número máximo de columnas
linealmente independientes 𝑟𝑎 , este valor hace referencia al rango de columnas de A. De
igual manera se puede establecer un rango de filas 𝑟𝑏 . Cuando estos dos valores son iguales
𝑟𝑎 = 𝑟𝑏 entonces se refieren al valor del rango de la matriz A. Para matrices de dimensiones
grandes es muy difícil calcularlo. En Matlab es cuestión del comando Rank (A).

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co

MATLAB
A=
%%% Rango de una matriz %%
2 3 1

A = [2 3 1; 1 2 3; 3 1 2] % Matriz A 1 2 3
rango= rank(A) % rango de A 3 1 2

rango = 3

4.2 Polinomios:[3]
Es frecuente tener que modificar los polinomios para su tratamiento y posterior solución, a
continuación veremos algunos parámetros que se pueden obtener de manera rápida con
Matlab:

Sintaxis Definición
Devuelve las raíces del
R = roots (p)
polinomio (p).
Devuelve los coeficientes del
P = poly (r) polinomio P, con las raíces
que conforman el vector r.
Crea el polinomio simbólico P
del vector de coeficientes b.
P = poly2sym (b)
la variable del polinomio es
x.
Usa la variable var como
P = poly2sym(b,var)
variable del polinomio.
Presenta todos los factores
F = factor (x) irreductibles de x en el
vector F.
Reescribe f en formato de
P = collect (f,var) potencias de términos
semejantes de la variable var.
Expande multiplicación de
P = expand (f) polinomios y productos
notables.
Trabajo en matemática
Syms var simbólica, con la variable
especificada.

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co
Sintaxis Definición
Multiplica los coeficientes
P = conv (A,B) del polinomio A con los
coeficientes del polinomio B.
Devuelve un vector Q que
contiene los coeficientes del
polinomio de la división. Y un
[Q,R] = deconv(A,B)
vector R que contiene los
coeficientes del polinomio
del residuo de la división.
Los polinomios en Matlab son expresados como un vector fila de sus coeficientes:

𝑝1 (𝑥) = 𝑥 5 − 25𝑥 4 + 230𝑥 3 − 950𝑥 2 + 169𝑥 − 945

𝑝2 (𝑥) = 𝑥 3 + 3𝑥 2 + 9𝑥 − 13

𝐷𝑒𝑠𝑐𝑟𝑖𝑝𝑐𝑖ó𝑛 𝑒𝑛 𝑀𝑎𝑡𝑙𝑎𝑏

𝑝1 = [1 −25 230 −950 169 − 945]

𝑝2 = [1 3 9 − 13]

Raíces de un polinomio

Las raíces son los valores de x que hacen que el polinomio tome un valor de cero.

MATLAB
%%%%%% Raíces de un polinomio %%%%%%% raices_p1 = raices_p2 =

p1 = [1 -25 230 -950 1689 -945] 9.0000 1.0000 + 0.0000i

% Polinomio 7.0000 -2.0000 + 3.0000i
raices_p1= roots (p1) % Raíces 5.0000 -2.0000 – 3.0000i
p2 = [1 3 9 -13] % Polinomio 3.0000
raices_p2= roots (p2) % Raices 1.0000

𝑅𝑎𝑖𝑐𝑒𝑠 𝑑𝑒𝑙 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑜

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co
𝑝1 (𝑥) = (𝑥 − 9)(𝑥 − 7)(𝑥 − 5)(𝑥 − 3)(𝑥 − 1)

𝑝2 (𝑥) = (𝑥 − 1) ∗ (𝑥 + 2 − 3𝑗) ∗ (𝑥 + 2 + 3𝑗)

Reconstrucción de un polinomio:

Se puede calcular el polinomio a partir de sus raíces

MATLAB
%%%% Reconstrucción a partir de %%% polin_r1 =
%%%% las Raíces de un polinomio %%%
3 9 -13
r1 = [1 -2+3j -2-3j] % Raíces de p1
polin_r1 = poly (r1) % polinomio p1 psimb =
psimb=poly2sym(polin_r1) % convierte el
polinomio vectorial en simbólico. x^3 + 3*x^2 + 9*x – 13
syms s
psimb1=poly2sym(polin_r1,s) % Cambia la psimb 1=
variable de x a s.
s^3 + 3*s^2 + 9*s - 13

Factores de un polinomio

En este proceso el polinomio se descompone en el producto de otros polinomios.

MATLAB
%%%% Factores de un polinomios %%%%% D2 =

syms s; % se declara "s" simbólica, (s + 1)*(s^2 + 4*s + 8)

D2=factor (s^3+5*s^2+12*s+8) % Los
factores no incluye el cuadrático

Multiplicación de polinomios

Convolucionar dos vectores es equivalente a multiplicarlos, el resultado es un vector cuya

longitud es igual a la Longitud de P1 + Longitud de P2 - 1

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co

MATLAB
%%% Multiplicación de dos polinomios %%%%

p1 = [2 -4 8 12] % Polinomio p1
p2 = [1 -3 5 7] % Polinomio p2
p1p2 = conv (p1, p2) % Multiplicación

%%% Multiplicación forma polinomial %%%%

syms s % convierte a s en variable simbol

simb_p1p2 = collect ((2*s^3-4*s^2+8*s+12)*(s^3-3*s^2+5*s+7)) %
forma para multiplicar dos o más funciones simbólicas

p1p2 =
2 -10 30 -18 -24 116 84
simb_p1p2 =
2*s^6 - 10*s^5 + 30*s^4 - 18*s^3 - 24*s^2 + 116*s + 84

División de dos polinomios

El resultado de la división de P1 por P2 se almacena en la variable Q, y el resto se almacena

en la variable R.

MATLAB
%%%%%%%% División de dos polinomios %%%%%%%%
p1 = [2 -4 8 12] % 2*s^3-4*s^2+8*s+12 Polinomio p1
p2 = [1 -3 5 7] % s^3-3*s^2+5*s+7 Polinomio p2
[Q, R] = deconv (p1, p2) % División de p1/p2
p1 = conv (Q, p2) + R % comprobar el resultado
syms s % Establece a s como simbólica
psimb=poly2sym(p1) % convierte el polinómio simbólico

Q= R=
2 0 2 -2 -2
p1 = psimb=
2 -4 8 12 2*x^3 - 4*x^2 + 8*x + 12

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co
3 2
𝑝1 = 2𝑥 − 4𝑥 + 8𝑥 + 12

𝑝2 = 𝑥 3 − 3𝑥 2 + 5𝑥 + 7

𝑝1 2𝑥 2 − 2𝑥 − 2
=2+ 3
𝑝2 𝑥 − 3𝑥 2 + 5𝑥 + 7

5. ACTIVIDADES:

1. Se quiere analizar un sistema eléctrico y para ello se aplica la segunda ley de

Kirchhoff Σ 𝑣 = 0 . El análisis arroja el siguiente sistema de ecuaciones lineales
cuyas variables son el 𝑣1 , 𝑣2 , 𝑣3 .

−𝑣3 + 3𝑣1 − 2𝑣4 = 0

6𝑣1 + 80 − 𝑣3 − 2𝑣4 = 0
−16𝑣4 + 6𝑣1 − 5𝑣3 = 40

2. Exprese el sistema de ecuaciones en forma matricial.

3. Si el sistema es de la forma AX = B. cuáles serían los valores de los coeficientes

de A. Cuáles serían los valores de los coeficientes de B.

4. Determine el valor de los coeficientes de la matriz adjunta de A. Para ello realice los
procedimientos analíticos en el cuaderno y utilizando la herramienta de Matlab
compruebe si corresponden a los obtenidos en el proceso anterior.

5. Calcule de manera analítica en el cuaderno el determinante de A. Utilice la

herramienta de Matlab para comprobar su valor.

6. Utilizando los determinantes, use la regla de Cramer’s para resolver el sistema

A X = B.

7. Encuentre el rango de la Matriz A.

8. Determine la matriz traspuesta de A.

9. Para el circuito de la figura aplique la primera ley de Kirchhoff Σv = 0 y encuentre

el sistema de ecuaciones lineales que relacionan las corrientes de malla. Responda
para este caso las preguntas (2, 3, 4, 5, 6, 7,8).

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co

10. De los polinomios 𝑝1 = 𝑥 3 − 𝑥 2 + 11𝑥 − 51 ; 𝑝2 = 𝑥 2 − 5𝑥 + 45 encuentre las

raíces de los polinomios.

11. Encuentre los factores de cada polinomio si existen.

12. Multiplique el polinomio p1 y p2.

13. Divida el polinomio p1 entre p2.

6. BIBLIOGRAFÍA.

[1] C. S. Chin, Computer Aided Control Systems Design. Boca Raton London
New York: Taylor & Francis Group, an informa business, 2013.
[2] S. T. Karris, Numerical Analysis Using MATLAB and Excel, 3th ed., vol. 1, no.
7. Fremont, California, USA, 2007.
[3] S. S. Haykin and B. Van Veen, Signals and systems. 1998.
[4] K. Ogata, Ingeniería de Control Moderna, vol. 1, no. 9. Madrid, 2010.
[5] White Robert: Computational Mathematics: Models, Methods, and Analysis
with MATLAB and MPI. Chapman and Hall/CRC; 1 edition, 2003.
[6] Hanselman, D.; Littlefield, B. MATLAB 5: Versão do Estudante, Guia do
Usuário, Makron Books, 1999.
Paginas de Internet.

 https://www.mathworks.com/

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia
 NIT. 890500622 - 6

www.ufps.edu.co

Av. Gran Colombia No. 12E-96 Colsag

Teléfono: 5776655
Cúcuta - Colombia