Marcillo Salcedo Odegil4 Paralelo 701

UNIVERSIDAD DE GUAYAQUIL
FACULTAD DE INGENIERÍA QUÍMICA

CARRERA DE INGENIERÍA QUÍMICA
SIMULACION DE PROCESOS II
ESTUDIANTE:
Marcillo Tigua Gabriela Noemi

Salcedo Chango Andrea
Elizabeth
DOCENTE:
Ing. Gonzalo Villa Manosalvas
TEMA:
Taller 5
FECHA DE ELABORACION:
06/02/2021
FECHA DE ENTREGA:
08/02/2021
PARALELO:
7-1
PERIODO:
2020 –2021 CII
SOLUCIONES NUMÉRICAS DE ECUACIONES DIFERENCIALES
EMPLEO DEL CODIGO ODEGIL4
El método más común de resolución numérica de ecuaciones diferenciales es el método

de Runge-Kutta de cuarto orden, cuya aplicación se resume a continuación.
 Sean las ecuaciones diferenciales tales como:
dy
=f (t , y 1 , y 2 ,… , y n )
dx 1
dy 2
=f 2 (t , y 1 , y 2 , … , y n )
dx
dy n
=f n(t , y 1 , y 2 , … , y n )
dx
 Su integración se realiza evaluando los siguientes parámetros:
pk =h∗f ∗(t k , y k )
h q
(
q k =h∗f t k +
2
yk + k
2 )
sk =h∗f ( t k + hy k +r )
 Las ecuaciones integradas numéricamente se obtienen mediante:
pk +2 q k + 2r k + s k
y k , i+1= y k , i+
6
Para integración de ecuaciones diferenciales por el método de Runge-Kuta de cuatro

parámetros se elaboró la función odegil4.m de paso de integración fijo.
RESOLVER EL SIGUIENTE EDO
d2 y dy
2
+3 +6 y=0
dt dt
Para su resolución numérica se realiza la siguiente transformación: z= y I
Con lo que la ecuación se transforma en dos ecuaciones diferenciales de primer grado
z ' =−3 z−6 y
y ' =z
Se pretende integrar desde –π hasta π, para lo cual es necesario conocer los valores de y (-π) y
de dy (-π), que se obtiene según el scrip.
SCRIPT 1: odegil4.m
SCRIPT 2: ecudif.m
SCRIPT 3: EcudifP1.m
COMANDOS UTILIZADOS
COMANDO: spline
Interpolación de datos de spline cúbico.
Sintaxis:
s = spline(x,y,xq)
pp = spline(x,y)
Descripción:
s = spline(x,y,xq) devuelve un vector de valores interpolados s correspondiente a los

puntos de consulta de xq. Los valores de s se determinan mediante la interpolación de
spline cúbico de x y y.
pp = spline(x,y) devuelve una estructura polinómica tramos para su uso por ppval y la
utilidad spline unmkpp.
Ejemplo 7
x = [0 1 2.5 3.6 5 7 8.1 10];

y = sin(x);
xx = 0:.25:10;
yy = spline (x, y, xx);
plot(x,y,'o',xx,yy)
COMANDO: reshape
Reformular un array.
Sintaxis:
B = reshape(A,sz)
B = reshape(A,sz1,...,szN)
Descripción:
B = reshape(A,sz) reformula A utilizando el vector de tamaño, sz, para definir size(B).

Por ejemplo, reshape(A,[2,3]) reformula A como una matriz de 2 por 3. sz debe
contener al menos 2 elementos y prod(sz) debe ser igual que numel(A).
B = reshape(A,sz1,...,szN) reformula A como un array de sz1 por ... por szN,

donde sz1,...,szN indica el tamaño de cada dimensión. Puede especificar un único
tamaño de dimensión de [] para que el tamaño de la dimensión se calcule
automáticamente, de forma que el número de elementos de B coincida con el número de
elementos de A. Por ejemplo, si A es una matriz de 10 por 10, reshape (A,2,2,
[]) reformula los 100 elementos de A en un array de 2 por 2 por 25.
Ejemplo 8:
Reformular un vector como matriz.
Reformular un vector de 1 por 10 como una matriz de 5 por 2.
A = 1:10;
B = reshape (A, [5,2])
B = 5×2
1 6
2 7
3 8
4 9
5 10
COMANDO: feval
Evaluar la función.
Sintaxis: [y1,...,yN] = feval(fun,x1,...,xM)
Descripción: [ y1,...,yN ] = feval(fun,x1,...,xM ) evalúa una función usando su nombre

o su manejador, y utilizando los argumentos de entrada x1,...,xM .
La función feval sigue las mismas reglas de alcance y precedencia que llamar
directamente a un identificador de función.
Ejemplo 9:
Evaluar función con nombre de función como vector de caracteres
Redondee el valor de pi al entero más cercano utilizando el nombre de la función.
diversión = 'ronda';
x1 = pi;
y = feval (diversión, x1)
y=3
Redondee el valor de pi a dos dígitos a la derecha del punto decimal.
x2 = 2;
y = feval (diversión, x1, x2)
y = 3.1400
COMANDO: tic
Iniciar cronómetro.
Sintaxis:
tic
timerVal = tic
Descripción:
Tic inicia un cronómetro para medir el rendimiento. La función registra el tiempo

interno en la ejecución del comando tic. Mostrar el tiempo transcurrido con la función
de toc.
timerVal = tic devuelve el valor del temporizador interno en la ejecución del

comando tic, por lo que puede registrar tiempo para palmos del tiempo simultáneo.
Ejemplo 10:
Mida cómo varía el tiempo requerido para resolver un sistema lineal con el orden de una
matriz:
t = zeros (1,100);
for n = 1: 100
A = rand (n, n);
b = rand (n, 1);
tic;
x = A \ b;
t(n) = toc;
end
COMANDO: toc
Leer el tiempo transcurrido desde el cronómetro.
Sintaxis:
toc
elapsedTime = toc
toc(timerVal)
elapsedTime = toc(timerVal)
Descripción:
Toc lee el tiempo transcurrido desde el cronómetro iniciado por el tic función. La
función lee el tiempo interno en la ejecución del comando y muestra el tiempo
transcurrido desde la última llamada a la función que no tuvo salida, en segundos.toc tic
elapsedTime = toc Devuelve el tiempo transcurrido en una variable.
toc(timerVal)muestra el tiempo transcurrido desde el ticcomando correspondiente a

timerVal.
elapsedTime = toc(timerVal)devuelve el tiempo transcurrido desde el ticcomando

correspondiente a timerVal.
Ejemplo 11:
Mida cómo varía el tiempo requerido para resolver un sistema lineal con el orden de una
matriz:
t = ceros (1,100);
for n = 1: 100
A = rand (n, n);
b = rand (n, 1);
tic;
x = A \ b;
t(n) = toc;
end
trama (t)
CONCLUSIONES.
Se utilizó la función odegil4, con el fin de obtener resultados óptimos en el

Command Windows, ya que esta función da una mayor precisión en los cálculos.
Es una potente herramienta que dispone para realizar cálculos con capacidades que
gradualmente podemos ir incluyendo a nuestro conocimiento y aplicaciones del
quehacer laboral e intelectual.
La resolución de ecuaciones diferenciales de forma numérica es una herramienta útil
para encontrar la solución de un sistema dinámico representado por ese tipo de
expresiones matemáticas.
RECOMENDACIONES.
Es muy importante conocer las condiciones iniciales del ejercicio, para evitar así
cualquier tipo de error al momento de resolver el ejercicio.
Tener en consideración los comandos su sintaxis, esta información se la puede
encontrar en la ayuda que nos ofrece Matlab.
La función odegil4, es de suma importancia, dado a que guarda los parámetros del
problema a resolver.
Conocer que hacen los comandos tic y toc (dan el tiempo que transcurre en correr un
ejercicio).

Marcillo Salcedo Odegil4 Paralelo 701

Cargado por

Copyright:

Formatos disponibles

Marcillo Salcedo Odegil4 Paralelo 701

Cargado por

Información del documento

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

Marcillo Salcedo Odegil4 Paralelo 701

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE GUAYAQUIL

FACULTAD DE INGENIERÍA QUÍMICA

Marcillo Tigua Gabriela Noemi

Ing. Gonzalo Villa Manosalvas

2020 –2021 CII

SOLUCIONES NUMÉRICAS DE ECUACIONES DIFERENCIALES

EMPLEO DEL CODIGO ODEGIL4

El método más común de resolución numérica de ecuaciones diferenciales es el método

 Sean las ecuaciones diferenciales tales como:

 Su integración se realiza evaluando los siguientes parámetros:

 Las ecuaciones integradas numéricamente se obtienen mediante:

Para integración de ecuaciones diferenciales por el método de Runge-Kuta de cuatro

RESOLVER EL SIGUIENTE EDO

Para su resolución numérica se realiza la siguiente transformación: z= y I

Con lo que la ecuación se transforma en dos ecuaciones diferenciales de primer grado

z ' =−3 z−6 y

Interpolación de datos de spline cúbico.

s = spline(x,y,xq) devuelve un vector de valores interpolados s correspondiente a los

x = [0 1 2.5 3.6 5 7 8.1 10];

B = reshape(A,sz) reformula A utilizando el vector de tamaño, sz, para definir size(B).

B = reshape(A,sz1,...,szN) reformula A como un array de sz1 por ... por szN,

Reformular un vector como matriz.

Reformular un vector de 1 por 10 como una matriz de 5 por 2.

B = reshape (A, [5,2])

Sintaxis: [y1,...,yN] = feval(fun,x1,...,xM)

Descripción: [ y1,...,yN ] = feval(fun,x1,...,xM ) evalúa una función usando su nombre

Evaluar función con nombre de función como vector de caracteres

Redondee el valor de pi al entero más cercano utilizando el nombre de la función.

y = feval (diversión, x1)

Redondee el valor de pi a dos dígitos a la derecha del punto decimal.

y = feval (diversión, x1, x2)

Tic inicia un cronómetro para medir el rendimiento. La función registra el tiempo

timerVal = tic devuelve el valor del temporizador interno en la ejecución del

A = rand (n, n);

b = rand (n, 1);

Leer el tiempo transcurrido desde el cronómetro.

elapsedTime = toc Devuelve el tiempo transcurrido en una variable.

toc(timerVal)muestra el tiempo transcurrido desde el ticcomando correspondiente a

elapsedTime = toc(timerVal)devuelve el tiempo transcurrido desde el ticcomando

A = rand (n, n);

b = rand (n, 1);

Se utilizó la función odegil4, con el fin de obtener resultados óptimos en el

También podría gustarte