UNIVERSIDAD NACIONAL DE JAÉN
FACULTAD DE INGENIERÍA MECÁNICA
ELÉCTRICA
CARRERA PROFESIONAL DE INGENIERÍA MECÁNICA Y ELÉCTRICA
SEMANA 08
Asignatura: MÉTODOS NUMÉRICOS
Docente: DR. LIC. LENIN QUIÑONES HUATANGARI
Estudiantes:
INTEGRANTES EJERCICIOS DESEMPEÑO
Angel Huaman Villalobos 1-16 Presentó
Benjamín Samaniego Adrianzén 17-32 Presentó
Luis Campos Guerrero 33-48 No presentó
Alexander Vega Ocupa 49-64 Presentó
Semestre: 2023-I
Ciclo: V
Jaén, 24 de Mayo del 2023
ACTIVIDADES Y EVALUACIÓN
ACTIVIDAD 1: REVISAR LA GUÍA DE APRENDIZAJE DE LA SEMANA
Leer detenidamente, además realizar los ejercicios en el software Octave y el manual de la
Universidad Politécnica de Madrid.
EJEMPLO 01
Escribir en OCTAVE:
EJEMPLO 02
Escribir en OCTAVE:
EJEMPLO 03
Escribir en OCTAVE:
EJEMPLO 04
Escribir en OCTAVE:
Se considera z1 = 3 + 4i , z2 = 4 + 3i . Determinar, con OCTAVE, los siguientes apartados:
1. Realizar las siguientes operaciones:
2. Calcular el módulo y el argumento de z1 .
3. Escribir la forma trigonométrica y exponencial de z1 . 4. Calcular sin(z1 ), cos(z1 )
EJEMPLO 05
Escribir en OCTAVE:
EJEMPLO 06
Escribir en OCTAVE:
EJEMPLO 07
Escribir en OCTAVE:
EJEMPLO 08
Escribir en OCTAVE:
EJEMPLO 09
Escribir en OCTAVE:
Dados los vectores:
Calcular, con OCTAVE, los siguientes apartados:
1. El producto escalar de ~a y ~b.
2. El producto vectorial de ~a y ~b.
3. El vector unitario perpendicular al plano formado por los dos vectores.
EJEMPLO 10
Escribir en OCTAVE:
El tensor de la figura, se ajusta hasta que la tensión del cable AB es de 2, 5kN. Determinar, con
OCTAVE, el momento respecto al punto O de la tensión del cable ,que actúa en el punto A, y la
magnitud de ese momento.
EJEMPLO 11
Escribir en OCTAVE:
Resolver el sistema lineal de ecuaciones:
EJEMPLO 12
Escribir en OCTAVE:
Una compañía petrolera dispone de tres refinerías de petróleo. Estas se denominan de la siguiente
forma: Refinería 1, Refinería 2 y Refinería 3. Cada refinería produce tres productos basados en el
crudo: Alquitrán, Gasóleo y Gasolina. Supongamos que, de un barril de petróleo, se sabe que:
• la primera refinería produce 4 litros de alquitrán, 2 de gasóleo, y 1 de gasolina.
• la segunda refinería produce 2 litros de alquitrán, 5 de gasóleo y 2.5 de gasolina.
• y la tercera refinería produce 2 litros de alquitrán, 2 de gasóleo y 5 de gasolina.
Supongamos que hay una demanda de estos productos de la siguiente manera:
• 600 litros de alquitrán.
• 800 litros de gasóleo.
• 1000 litros de gasolina.
¿Cuántos barriles de crudo necesitará cada refinería para satisfacer la demanda? El enunciado se
puede representar de la siguiente forma:
EJEMPLO 13
Escribir en OCTAVE:
Dibuja la gráfica de la función:
EJEMPLO 14
Escribir en OCTAVE:
Dibuja la gráfica de la función:
EJEMPLO 15
Escribir en OCTAVE:
Dibuja la gráfica de la función:
EJEMPLO 16
Escribir en OCTAVE:
Dibuja la gráfica de la función:
EJEMPLO 17
Modificaciones de las graficas
≫ 𝑒𝑧𝑝𝑙𝑜𝑡(´𝑠𝑞𝑟𝑡(1 − 𝑥 ∧ 2)´
≫ 𝑎𝑥𝑖𝑠 ([−2,2, −1,2])
EJEMPLO 18
Se puede añadir un titulo y etiquetas a los ejes :
≫ 𝑡𝑖𝑡𝑙𝑒(′𝐹𝑖𝑔𝑢𝑟𝑎 𝑝𝑟𝑢𝑒𝑏𝑎 ′ )
≫ 𝑥𝑙𝑎𝑏𝑒𝑙(′𝐸𝑗𝑒 𝑑𝑒 𝑙𝑎𝑠 𝑥´)
≫ 𝑦𝑙𝑎𝑏𝑒𝑙(´𝐿𝑜 𝑞𝑢𝑒 𝑞𝑢𝑖𝑒𝑟𝑎𝑠´)
EJEMPLO 19
𝑥
≫ 𝑒𝑧𝑝𝑙𝑜𝑡(´ sin(𝑥 ∧ 2) ∗ 2 ´ , [−4,4]
≫ ℎ𝑜𝑙𝑑 𝑜𝑛
≫ 𝑒𝑧𝑝𝑙𝑜𝑡(′𝑥 − 2′ , [−4,4])
1′
≫ 𝑒𝑧𝑝𝑙𝑜𝑡(′𝑥 ∧ 2 − 𝑥 , [−4,4])
≫ ℎ𝑜𝑙𝑑 𝑜𝑓𝑓
EJEMPLO 20
𝑥′
≫ 𝑒𝑧𝑝𝑙𝑜𝑡 (′ sin(𝑥 ∧ 2) ∗ 2 , [−4,4]
≫ ℎ𝑜𝑙𝑑 𝑜𝑛
≫ 𝑒𝑧𝑝𝑙𝑜𝑡(′𝑥 − 2′ , [−4,4])
1′
≫ 𝑒𝑧𝑝𝑙𝑜𝑡 (′𝑥 ∧ 2 − 𝑥 , [−4,4])
≫ ℎ𝑜𝑙𝑑 𝑜𝑓𝑓
EJEMPLO 21
≫ 𝑠𝑢𝑏𝑝𝑙𝑜𝑡(2,2,1)
′𝑥 ′
≫ 𝑒𝑧𝑝𝑙𝑜𝑡( 2 , [0,4])
≫ 𝑠𝑢𝑏𝑝𝑙𝑜𝑡(2,2,2)
≫ 𝑒𝑧𝑝𝑙𝑜𝑡(′ sin(3 ∗ 𝑥)′ )
≫ subplot(2,2,3)
≫ ezplot(’x^2’, [0,4])
≫ subplot(2,2,4)
≫ ezplot(’cos(x/2)’, [−1,4])
EJEMPLO 22
≫ x = 0: 0.01: 3;
≫ y = abs(exp(−0.5 ∗ x) .∗
sin(5 ∗ x));
≫ subplot(221);
≫ plot(x, y)
≫ title(’lineal’)
≫ hold on
≫ subplot(222)
≫ loglog(x, y)
≫ title(’logaritmica’)
≫ subplot(223)
≫ semilogx(x, y)
≫ title(’semilogaritmico en eje x’)
≫ subplot(224)
≫ semilogy(x, y)
>> title(’semilogaritmico en eje y’)
EJEMPLO 23
E>> t = 0: 900; A = 1000;
≫ a = 0.005; b = 0.005;
≫ z1 = A ∗ exp(−a ∗ t) ;
≫ z2 = sin(b ∗ t) ;
≫ [haxes, hline1, hline2] = plotyy(t, z1, t, z2, ’ semilogy’, ’plot’);
≫ axes(haxes(1)), ylabel(’Semilogaritmico’)
≫ axes(haxes(2)), ylabel(’Lineal’)
>> set(hline2, ’LineStyle’, ’ − −’)
EJEMPLO 24
>> ezplot3(’3 ∗ cos(t)’, ’t ∗ sin(t^2)’, ’sqrt(t)’)
EJEMPLO 25
≫ ezmesh(f)
≫ ezmesh(f, [a, b])
>> ezmesh(f, [a, b, c, d])
donde, f es una expresión de dos variables. Representa la superficie z = f (x, y) para (x, y)
considerando los intervalos indicados en la función..
>> ezmesh(’x* exp(−x^2 − y^2)’)
EJEMPLO 26
≫ ezcontour(f)
≫ ezcontour(f, [a, b])
>> ezcontour(f, [a, b, c, d]
EJEMPLO 27
>> ezmeshc(f)
Representa simultáneamente las líneas de nivel y la superficie
>> ezmeshc(’sin(u/2) ∗ sin(v/2)’)
EJEMPLO 28
>> ezsurf(f)
Representa una superficie coloreada z = f (x, y)
>> ezsurf(’sin(sqrt(x^2 + y^2))/sqrt(x^2 + y^2)’)
EJEMPLO 29
>> ezsurfc(f)
Representa simultáneamente las curvas de nivel y la superficie.
>> ezsurfc(’sin(sqrt(x^2 + y^2))/sqrt(x^2 + y^2)’)
EJEMPLO 30
Calcular, con OCTAVE, una raíz de la ecuación:
𝑥
x + log ( ) = 0
3
>> ezplot(’x + log(x/3)’)
Vemos, a simple vista, que la raíz está cerca de x=1.
x
≫ fun = @(x)x + log (3)
>> z = fzero(fun, 1)
Solución OCTAVE: z = 1.0499
EJEMPLO 31
Calcular, con OCTAVE, las raíces de la ecuación:
𝑥
sin ( ) + cos(√𝑥) = 0 𝑒𝑛 [𝜋, 3𝜋]
2
>> ezplot(’sin(x/2) ∗ cos(sqrt(x))’, [−pi, 3 ∗ pi])
A “simple vista"se observa que tiene 3 raíces: una “cerca"de x=0, otra “cerca"de x=2 y otra
“cerca"de x=6
x
≫ fun = @(x) sin (2) .∗ cos(sqrt(x))
≫ z1 = fzero(fun, 0)
≫ z2 = fzero(fun, 2)
>> z3 = fzero(fun, 6)
Solución OCTAVE: z1 = 0.0 ; z2 = 2.4674 ; z3 = 6.2832
EJEMPLO 32
Un cable de una línea de transmisión de energía eléctrica, cuyos anclajes están a la misma
altura tiene una cuerda de 240 [m]. El cable pesa 83 [N/m] y la flecha en su punto medio es
de 60 [m]. Determinar, con OCTAVE, la tensión del cable en el punto medio. La ecuación
cartesiana de la catenaria se expresa de la siguiente forma :
La tensión del cable en el punto medio puede determinarse utilizando la ecuación de la curva
catenaria, que representa la forma que adopta el cable bajo la acción exclusiva de su peso.
x
≫ fun = @(x) sin ( ) .∗ cos(sqrt(x))
2
≫ w = 83; ≫ L = 240; ≫ h = 60;
L2
≫ T0 = w ∗ 8∗ h ;
T w∗ L
≫ fun = @(To ) (( wo ) ∗ (cosh (2∗ T ) − 1)) − h;
o
>> T_o = fzero(fun, T_0)
Solución OCTAVE: Tensión del cable en su punto medio T = 10700.23 N
[EJERCICIOS 33-48]
EJEMPLO 49 (LEY DE ENFRIAMIENTO DE NEWTON)
Insertamos el código para la resolución del ejercicio:
Y nos da la siguiente gráfica:
EJEMPLO 50 (TEOREMA DE TORRICELLI)
La ecuación diferencial para expresar la altura del agua en cualquier instante t es:
dz
A( z ) Cd A0 2 gz
dt
Si el depósito es esférico, 𝐴(𝑧) = 𝜋 ∗ (2𝑅𝑧 − 𝑧 2 ), y la ecuación a integrar queda de la siguiente
forma: 2
dz r
Cd 2 gz
dt (2 Rz z 2 )
z (0) 1,5
Ahora ingreso el código: (en el intervalo de [1,100])
El cual nos da el siguiente gráfico:
EJEMPLO 51 (TANQUE DE MEZCLA)
Ecuación que modela el proceso:
dMs (t ) Ms (t )
Qe Cse (t ) Qs
dt V
0, t 0
Cse (t ) 4, 0 t 10
2, t 10
Ms (t 0) 300
Ahora ingreso el código:
EJEMPLO 52 (OSCILADOR DE VAN DER POL)
El problema es equivalente a :
z1 z2
z2 (1 z12 ) z2 z1
z1 (0) y00
z2 (0) y10
Ahora ingreso el código:
El cual nos da los siguientes gráficos:
EJEMPLO 53 (OSCILACIONES LIBRES AMORTIGUADAS)
Ingresamos el código:
El cual me da el siguiente gráfico:
EJEMPLO 54 (VARIACIÓN DE CONDICIONES INICIALES)
Ingreso el código:
El cual me da el siguiente gráfico:
EJEMPLO 55 (VARIACIÓN DE UN PARÁMETRO)
Ingresamos el código:
El cual nos da la siguiente función:
EJEMPLO 56 (MODELO DE LOTKA-VOLTERRA)
Ingresamos el código:
El cual nos da los siguientes gráficos:
EJEMPLO 57 (MAXIMIZACIÓN DE LA FUNCIÓN OBJETIVO)
Ahora ingresamos el código:
EJEMPLO 58 (MINIMIZACIÓN DE LA FUNCIÓN OBJETIVO)
Ahora ingresamos el código:
EJEMPLO 59 (MODELO DE TRANSPORTE)
Ahora ingreso el código:
EJEMPLO 60 (MODELO DE MEZCLA DE PRODUCTOS)
Ahora ingreso el código:
EJEMPLO 61 (OPTIMIZACIÓN DE PARÁMETROS EDO. CINÉTICA QUÍMICA)
Ahora ingreso el código:
EJEMPLO 62 (MINIMIZACIÓN DE LA FUNCIÓN OBJETIVO)
Ahora ingreso el código:
EJEMPLO 63 (UBICACIÓN DE CENTRO LOGÍSTICO)
EJEMPLO 64 (DISEÑO DE UNA TORRE TUBULAR DE UN AEROGENERADOR)
Ahora ingreso el código:
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
Scanned by TapScanner
