Info - Lab1 - Candia Zapata Estefany Klenny

UNIVERSIDAD NACIONAL DE SAN ANTONIO ABAD DEL
CUSCO
FACULTAD DE INGENIERIA DE PROCESOS
ESCUELA PROFESIONAL DE INGENIERIA QUIMICA
CONTROL DE PROCESOS
LABORATORIO N° 1:
“SOLUCION DE EDO’S CON MATLAB Y LABVIEW”
Alumna:
- Candia Zapata Estefany Klenny / 124695
Ing. Franklin A. Salas Camacho

Docente
22/06/2017
1) Resumen
En la presente práctica de laboratorio se desarrolló la solución analítica de ecuaciones
diferenciales ordinarias lineales con el uso de Matlab y Labview. En el software Matlab se
desarrollaron dos métodos de solución para los ejemplos dados: dsolve y ODE45 Runge
Kutta de 4to Orden, mientras que en el Software Labview solo se desarrolló un método
Runge Kutta de 4to Orden.
Con el fin de reforzar más el conocimiento adquirido se desarrolló un ejercicio para resolver
una ecuación diferencial ordinaria por el método dsolve en el Software Matlab.
2) Introducción
2.1) Objetivo
 Aprender a resolver numéricamente las ecuaciones diferenciales haciendo uso de los
software Matlab y Labview
2.2) Ecuaciones Diferenciales
2.2.1) Definición
Una ecuación diferencial es una igualdad que contiene derivadas de una función
desconocida con respecto a una o más variables independientes.
Existen diferentes tipos de ecuaciones diferenciales, por lo que se hace necesario
realizar una clasificación de ellas.
2.2.1) Clasificación de ecuaciones diferenciales
2.2.2.1) Clasificación por tipo
a) Ecuación diferencial Ordinaria (EDO)
Si una ecuación contiene solo derivadas de una o más variables
dependientes respecto a una sola variable independiente se dice que
es una ecuación diferencial ordinaria
b) Ecuación diferencial Parcial (EDP)
Una ecuación que involucra derivadas parciales de una o más
variables dependientes de dos o más variables independientes se
llama ecuación diferencial parcial
2.2.2.3) Clasificación por orden

El orden de una ecuación diferencial (ya sea EDO o EDP) es el orden de la
mayor derivada en la ecuación
3.2.2.3) Clasificación por linealidad

a) Ecuación diferencial Lineal
Una ecuación diferencial es lineal, si se puede escribir de la forma:
𝑎𝑛 (𝑥)𝑦 (𝑛) + 𝑎𝑛−1 (𝑥)𝑦 (𝑛−1) + ⋯ + 𝑎2 (𝑥)𝑦 ′′ + 𝑎1 (𝑥)𝑦 ′ + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥)
Esto implica que debe cumplir con las siguientes condiciones:

(i) La función desconocida y sus derivadas son a lo sumo de
primer grado, es decir, de potencia 1.
(ii) Los coeficientes de la función desconocida y sus derivadas
dependen solo de la variable independiente.
b) Ecuación diferencial No-Lineal
En caso de que no se cumpla alguna de las condiciones anteriores, se
dice que la Ecuación diferencial es no lineal. Por ejemplo:
𝑦 ′′ − (𝑦 2 + 1)𝑦 ′ = 𝑥
Es no lineal, ya que el coeficiente 𝑦 ′ depende de 𝑦
𝑑4 𝑦 𝑑𝑦
+𝑥 − 𝑦2 = 0
𝑑𝑥 4 𝑑𝑥
No es lineal, ya que el termino 𝑦 no es de primer grado
2.3) Metodo de Runge Kutta de 4to orden
Este metodo trata de resolver ecuaciones diferenciales ordinarias, se utiliza para
calcular la aproximacion de 𝑦𝑛+1 (que es la variable dependiente). El metodo de
cuarto orden mas habital es el determinado por las siguientes formulas:
ℎ
𝑦𝑛+1 = 𝑦𝑛 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 )
6
𝑘1 = 𝑓(𝑥𝑛 ; 𝑦𝑛 )
ℎ ℎ
𝑘2 = 𝑓 (𝑥𝑛 + ; 𝑦𝑛 + 𝑘1 )
2 2
ℎ ℎ
𝑘3 = 𝑓 (𝑥𝑛 + ; 𝑦𝑛 + 𝑘2 )
2 2
𝑘2 = 𝑓(𝑥𝑛 + ℎ; 𝑦𝑛 + ℎ 𝑘3 )
3) Materiales y métodos
3.1) Materiales
 Laptop
 Software Matlab
 Software Labview
3.2) Métodos
3.2.1) Solución de ecuaciones diferenciales ordinarias en Matlab
3.2.1.1) Con dsolve solucionador de EDO’s
Ejemplo 1
Se procederá a hallar la solución analítica de la ecuación (1), con la
siguiente condición inicial: 𝑦(0) = 1
5𝑦 ′ + 4𝑦 − 2 = 0 (1)
Ahora procedemos a realizar la solución en Matlab. Escribimos los códigos
mostrados en la Figura 1 en la ventana de comandos del software Matlab,
Figura 1. Secuencia de códigos dsolve para Matlab

Al escribir el código ezplot (A), Matlab graficará la solución de la ecuación
(1), lo cual se observa en la Figura 2
Figura 2. Grafica resultante del código para la solución de la ecuación (1)
3.2.1.2) Con ODE45 Runge de 4to Orden

Ejemplo 2
El perfil de las concentraciones de una reacción en serie efectuada en un
reactor discontinuo e isotérmico, cuyas ecuaciones resultantes del modelado
son (2), (3) y (4):
𝑑[𝑎]
= −𝑘1 [𝑎] (2)
𝑑𝑡
𝑑[𝑟]
𝑑𝑡
= 𝑘1 [𝑎] − 𝑘2 [𝑟] (3)
𝑑[𝑠]
𝑑𝑡
= 𝑘2 [𝑟] (4)
Con las siguientes condiciones iniciales:
[𝑎] = 1 y [𝑟] = [𝑠] = 0
−1
𝑘1 = 0.064074 𝑠 y 𝑘2 = 0.04344 𝑠 −1
Remplazando estos valores en las ecuaciones anteriores, se obtienen las
ecuaciones (5), (6) y (7):
𝑑[𝑎]
𝑑𝑡
= −0.064074[𝑎] (5)
𝑑[𝑟]
𝑑𝑡
= 0.064074[𝑎] − 0.04344[𝑟] (6)
𝑑[𝑠]
= 0.04344[𝑟] (7)
𝑑𝑡
Ahora procedemos a realizar la solución en Matlab. En la ventana de

comandos del software escribiremos los códigos que se muestran en la Figura
3
Figura 3. Código para la solución de un sistema de EDO’s

Al ejecutar el código plot (t,x), obtenemos la Figura 4, donde se representa el
perfil de las concentraciones de una reacción en serie en función del tiempo.
Figura 4. Grafica resultante del código para la solución de un sistema de EDO’s
3.2.2) Solución de ecuaciones diferenciales ordinarias con Labview

Iniciaremos un NewVI en Labview, luego hacemos clic derecho en la
ventana de Block Diagram, buscamos la opción Mathematics/Differrential
Equations/Ordinary Differential Equations/RK 4th Order, como lo podemos
observar en la Figura 5
Figura 5. Secuencia de despliegue de los bloques para la solución de EDO’s en Labview
Posteriormente colocamos el bloque RK 4th Order en la ventana de Block

Diagram, se procede a agregar los controles para el bloque como se muestra
en la Figura 6,y los indicadores como se muestra en la Figura 7, además de
agregar Waveform Graph para realizar la gráfica de las respuestas.
Figura 6. Bloque RK 4th Order con sus respectivos controles, indicadores y Waveform Graph
En el Front Panel se acomodan los elementos de entradas y salidas del programa, para
luego insertar la información disponible de las ecuaciones diferenciales así como los
valores de tiempos y nombres de las constantes, y por último se da a correr el
programa para poder observar los resultados, como se muestra en la Figura 7
Figura 7. Captura del Front Panel con el resultado de la simulación del sistema de EDO’s
4) Resultados
Hallar la solución analítica de la siguiente Ecuación Diferencial Ordinaria, por el método de
dsolve en el software Matlab, con las condiciones iniciales 𝑦(0) = 1 ; 𝐷𝑦(0) = 0:
𝑑2 𝑦
+ 16𝑦 = 0 (8)
𝑑𝑥 2
Procederemos a realizar la solución en Matlab. Escribimos los códigos mostrados en la

Figura 8 en la ventana de comandos del software Matlab,
Figura 8. Secuencia de códigos dsolve para Matlab
Al escribir el código ezplot (A, [0 10]), Matlab graficará la solución de la

ecuación (8), lo cual se observa en la Figura 9
Figura 9. Grafica resultante del código para la solución de la ecuación (8)
5) Conclusiones
Luego de haber desarrollado este laboratorio que tuvo como objetivo aprender a usar los
softwares Matlab y Labview para resolver Ecuaciones Diferenciales Ordinarias, se llegó a la
siguiente conclusión, estos dos softwares son complejos debido a que cada uno tiene su
propio lenguaje de programación, a mi parecer el más complejo es el software Labview, pero
a pesar de su complejidad son de gran utilidad para el control de procesos en ingenieria
6) Cuestionario
a) ¿Qué es una Ecuación Diferencial Ordinaria?
Una ecuación diferencial es una ecuación que involucra derivadas de una función desconocida
de una o más variables. Si la función desconocida depende sólo de una variable (de tal modo
que las derivadas son derivadas ordinarias) la ecuación se llama una ecuación diferencial
ordinaria.
7) Bibliografía
Spiegel, M. (1983). Ecuaciones Diferenciales Aplicadas (3ra ed.). México: Prentice Hall
Hispanoamerica.
Zill, D. G. (2009). Ecuaciones Diferenciales con aplicaciones de Modelado (9na ed.). México
D.F.: Cengage Learning Latinoamerica.
8) Páginas Webs
Campus USAL. (01 de 10 de 2009). Recuperado el 22 de 06 de 2017, de
http://campus.usal.es/~mpg/Personales/PersonalMAGL/Docencia/MetNumTema4Teo(0
9-10).pdf
Castillo U., C. (02 de 2010). Udomatematica. Recuperado el 22 de 06 de 2017, de
https://udomatematica.files.wordpress.com/2010/02/ecuaciones-diferenciales-
ordinarias1.pdf
MTY ITESM. (01 de 06 de 2008). Recuperado el 22 de 06 de 2017, de
http://www.mty.itesm.mx/dtie/deptos/cb/cb00854-1/Apuntes/HMA/MN06fED.pdf

Info - Lab1 - Candia Zapata Estefany Klenny

Cargado por

Copyright:

Formatos disponibles

Info - Lab1 - Candia Zapata Estefany Klenny

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Info - Lab1 - Candia Zapata Estefany Klenny

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE SAN ANTONIO ABAD DEL

FACULTAD DE INGENIERIA DE PROCESOS

ESCUELA PROFESIONAL DE INGENIERIA QUIMICA

“SOLUCION DE EDO’S CON MATLAB Y LABVIEW”

- Candia Zapata Estefany Klenny / 124695

Ing. Franklin A. Salas Camacho

2.2.2.3) Clasificación por orden

3.2.2.3) Clasificación por linealidad

𝑎𝑛 (𝑥)𝑦 (𝑛) + 𝑎𝑛−1 (𝑥)𝑦 (𝑛−1) + ⋯ + 𝑎2 (𝑥)𝑦 ′′ + 𝑎1 (𝑥)𝑦 ′ + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥)

Esto implica que debe cumplir con las siguientes condiciones:

Figura 1. Secuencia de códigos dsolve para Matlab

Figura 2. Grafica resultante del código para la solución de la ecuación (1)

3.2.1.2) Con ODE45 Runge de 4to Orden

Ahora procedemos a realizar la solución en Matlab. En la ventana de

Figura 3. Código para la solución de un sistema de EDO’s

Figura 4. Grafica resultante del código para la solución de un sistema de EDO’s

3.2.2) Solución de ecuaciones diferenciales ordinarias con Labview

Figura 5. Secuencia de despliegue de los bloques para la solución de EDO’s en Labview

Posteriormente colocamos el bloque RK 4th Order en la ventana de Block

Procederemos a realizar la solución en Matlab. Escribimos los códigos mostrados en la

Figura 8. Secuencia de códigos dsolve para Matlab

Al escribir el código ezplot (A, [0 10]), Matlab graficará la solución de la

Figura 9. Grafica resultante del código para la solución de la ecuación (8)

También podría gustarte