Aplicaciones de Matlab y Simulink
Aplicaciones de Matlab y Simulink
Aplicaciones de Matlab y Simulink
DE INGENIERIA INDUSTRIAL
TEMA
APLICACIONES DE MATLAB Y
SIMULINK
CURSO: SISTEMAS DINÁMICOS
DOCENTE:
ALUMNO:
CICLO: IX
2022
1
i
ÍNDICE
REFERENCIAS BIBLIOGRÁFICAS..................................................................................... 16
iii
LISTA DE ILUSTRACIONES
Matlab ofrece varios algoritmos numéricos para resolver una extensa variedad de ecuaciones
diferenciales. Esta demostración enseña la formulación y solución para dos tipos distintos de
La ecuación de Van Der Pol tiene una larga historia de utilización tanto en física como en
sismología para modelizar las dos placas en una falla geológica. En la dinámica, el oscilador
escalar que indica la fuerza de la amortiguación. Se puede demostrar a través del Teorema de
Lienard que existe un ciclo límite para el oscilador de Van Der Pol no forzado, por lo que es
Dos regímenes de interés para las características del oscilador no forzado son:
1
Esta es una forma del oscilador armónico simple, y por tanto se conserva la energía
mecánica.
Cuando µ> 0, el sistema entrará en un ciclo límite, donde la energía sigue siendo
conservada, pero parte de la energía mecánica se perderá en forma de calor. Cerca del
𝑑𝑥
origen 𝑥 = = 0 el sistema es inestable, y lejos del origen el sistema se
𝑑𝑡
amortigua.
El oscilador de Van Der Pol forzado o dirigido, toma la función original, y añade una
función de corrección 𝐴. 𝑠𝑒𝑛 (𝜔𝑡) para dar una ecuación diferencial de la forma:
de dos funciones ODE de primer orden. Estas son evaluadas para distintos valores del
parámetro Mu. Para una integración más rápida, elegimos un método de solución basado en
Para Mu = 1, cualquiera de los métodos de solución ODE Matlab puede resolver la ecuación
2
• ODE: es un método de solución que emplea Matlab para resolver ecuaciones
diferenciales ordinarias.
• [t,y]: Nos dice el método empleado por Matlab para la resolución del problema, y las
% Gráfico de la solución
• 𝑝𝑙𝑜𝑡: Comando que ordena a Matlab representar gráficamente la función que le damos.
𝑥𝑙𝑎𝑏𝑒𝑙(′𝑡 ′ )
𝑦𝑙𝑎𝑏𝑒𝑙(′𝑆𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑦 ′ )
3
Para magnitudes más grandes de Mu, el problema se vuelve más rígido, es decir, mas dificil
de resolver numericamente. Para una integración rapida son necesarios métodos numéricos
especiales como ODE15S, ODE23S, ODE23T, y ODE23TB, que pueden resolver problemas
rígidos eficientemente. Aquí hay una solución a la ecuación de Van Der Pol para Mu = 1000
𝑢𝑠𝑎𝑛𝑑𝑜 𝑂𝐷𝐸15𝑆.
𝑀𝑢 = 1000;
4
función de ejemplo TWOODE tiene una ecuación diferencial escrita como un sistema de dos
Antes de utilizar BVP4C, tenemos que proporcionar una suposición para la solución que
que refina la solución. BVPINIT ensambla la suposición inicial en la forma que el método de
solución BVP4C tendrá. Para una malla inicial de [0 1 2 3 4] y una suposición constante de
Con esta suposición inicial, podemos resolver el problema con BVP4C. La solución (abajo)
sol: Nos dice el método empleado por Matlab para la resolución del problema, y las variables
➢ plot: Comando que ordena a Matlab representar gráficamente la función que le damos.
5
➢ plot(xint,yint(1,:),’r’): representa gráficamente la función en color rojo.
Este problema particular del valor límite tiene exactamente dos soluciones. La otra solución
es obtenida para una suposición inicial de: 𝑦(𝑥) = −1, 𝑦′(𝑥) = 0 , y representando
𝑥𝑖𝑛𝑡 = 𝑙𝑖𝑛𝑠𝑝𝑎𝑐𝑒(0,4,50);
6
Ilustración 4. Grafica utilizando la función de BVP4C
Además, podemos añadir algunos comandos más para modificar las gráficas:
Es una representación gráfica de las funciones que lleva a cabo cada componente y el flujo de
señales. Tales diagramas muestran las relaciones existentes entre los diversos componentes.
tiene la ventaja de indicar de forma más realista el flujo de las señales del sistema real. En un
diagrama de bloques todas las variables del sistema se enlazan unas con otras mediante
bloques funcionales.
7
Definición de bloque
Es un símbolo para representar la operación matemática que sobre la señal de entrada hace el
bloque para producir la salida. Las funciones de transferencia de los componentes por lo
general se introducen en los bloques correspondientes, que se conectan mediante flechas para
realimenta al punto de suma para compararse con la entrada, es necesario convertir la forma
La salida C(s) se realimenta al punto de suma, donde se compara con la entrada de referencia
R(s). La naturaleza en lazo cerrado del sistema se indica con claridad en la figura. La salida
del bloque, C(s) en este caso, se obtiene multiplicando la función de transferencia G(s) por la
8
Función de transferencia
lineal, cada una de ellas puede tratarse de forma independiente; y las salidas correspondientes
a cada entrada pueden sumarse para obtener la salida completa. La forma en que se introduce
cada entrada en el sistema se muestra en el punto de suma mediante un signo más o un signo
menos.
9
Fuente: (Roldán Guzmán, 2016).
1. Expresar lar reglas básicas del álgebra de bloques para que Matlab realice los cálculos
2. Expresar las conexiones presentes en el diagrama mediante los comandos: series, 𝑝𝑎𝑟𝑎𝑙𝑙𝑒𝑙
y 𝑓𝑒𝑒𝑑𝑏𝑎𝑐𝑘.
𝑐𝑜𝑛𝑛𝑒𝑐𝑡.
Ejemplo. Reducir el diagrama de bloques de la ilustración 10. Para las siguientes funciones
de transferencia:
…………. (𝒆𝒄𝒖𝒂𝒄𝒊ó𝒏 𝟏. 𝟐)
>> 𝐺3 = 2;
comando 𝑚𝑖𝑛𝑟𝑒𝑎𝑙.
CASO PRÁCTICO
determinando a partir de las ecuaciones físicas un modelo no lineal que lo describe con
linealización por Taylor, dos modelos lineales que lo aproximan en torno a sendos puntos de
equilibrio. Tanto el modelo no lineal como los dos modelos linealizados se simularán
11
Asimismo, se realizará una animación del péndulo que permitirá corroborar visualmente el
El sistema se compone de una bola de masa m situada en el extremo de una barra de masa
despreciable con una longitud l. Además, se sabe que el momento de inercia del péndulo
El ángulo girado θ, que será la variable de salida y, se toma según indica la figura adjunta.
Modelado matemático
El ángulo θ queda determinado por la ecuación (1). El par T aplicado sobre el péndulo se
del péndulo.
12
Ecuaciones del espacio de estados
espacio de estados. En muchos sistemas físicos una elección adecuada consiste en tomar
como variables de estado la salida y sus derivadas. En este caso tomaremos el ángulo girado θ
la figura adjunta que define la variable de salida θ ante una entrada de par T.
𝑙 = 1; % 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑒𝑛 𝑚𝑒𝑡𝑟𝑜𝑠
13
𝑔 = 9.8; % 𝑎𝑐𝑒𝑙𝑒𝑟𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑎 𝑔𝑟𝑎𝑣𝑒𝑑𝑎𝑑 𝑚. 𝑠^2
𝑚 = 3; % 𝑚𝑎𝑠𝑎 𝑒𝑛 𝑘𝑔
implementación tiene el aspecto mostrado en la figura. Hay que tener en cuenta que las
condiciones iniciales quedan determinadas por los valores iniciales de los integradores.
14
ANÁLISIS CRÍTICO
Los temas resumidos son de base fundamental que se debe conocer para poder realizar
Simulink, así como también el sistema masa resorte, la cual de estos sistemas sus valores se
puede colocar en el mismo Simulink para poder simular o estos valores se les puede iniciarlas
con ayuda de un script, ya que de este modo las variables se almacenan en el workspace,
aunque también se puede realizar códigos y guardarlos en scripts para poder simularlos en
Matlab. También debemos tener en consideración que las ecuaciones diferenciales son de
gran importancia ya que gracias a ello podemos representar estos sistemas dinámicos
mediante ecuaciones matemáticas, aunque por lo general estas ecuaciones matemáticas son
son de gran ayuda para poder simplificar sistemas que contengan diagrama de bloques
15
REFERENCIAS BIBLIOGRÁFICAS
➢ García Reyes, L., & Rincón Pasaye, J. J. (s.f.). Reducción de Diagramas de Bloques
SISTEMAS FÍSICOS .
http://mathgene.usc.es/matlab/trabajos/ecuaciones_diferenciales.pdf
diagramas-de-bloques.html
https://la.mathworks.com/help/matlab/math/differential-equations.htm l
de 2022, de
http://isa.uniovi.es/wiki/isa/index.php/Modelado_y_Simulaci%C3%B3n_de_un_P%C
3%A9ndulo
16