Metodos Numericos - 2018
Metodos Numericos - 2018
Metodos Numericos - 2018
MÉTODOS
NUMÉRICOS
INGENIERIA DE SISTEMAS
Métodos Numéricos
sfs
© Corporación Universitaria
Remington
Tercera edición
2018
MÉTODOS NUMÉRICOS
Elkin Ceballos Gómez.
FACULTAD DE CIENCIAS BÁSICAS E INGENIERÍA
Editorial Uniremington
Medellín, Colombia
Derechos Reservados ©2011
Responsables
Jorge Mauricio Sepúlveda Castaño
Decano de la Facultad de Ciencias Básicas e Ingeniería
[email protected]
Edición y Montaje
Vicerrectoría de Educación a Distancia y Virtual
Equipo de diseño gráfico
www.uniremington.edu.co
[email protected]
3
TABLA DE CONTENIDO
Pág.
4
2.5.3 EJERCICIOS DE APRENDIZAJE 121
2.5.4 TALLER DE ENTRENAMIENTO 123
2.5.5 TALLER DE APRENDIZAJE 125
2.5.6 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO 129
2.5.7 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO. 130
2.5.8 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO 132
2.5.9 TALLER DE ENTRENAMIENTO 132
5 GLOSARIO 164
6 BIBLIOGRAFÍA 165
Métodos Numéricos
sfs
PROPÓSITO GENERAL
MÉTODOS NUMÉRICOS
MÉTODOS NUMÉRICOS
OBJETIVO GENERAL
Aplicar los métodos de manera eficiente en la solución de problemas que
involucran modelos matemáticos, procurando que la solución obtenida
mediante la aplicación de los diferentes algoritmos sea óptima, precisa y exacta.
OBJETIVOS ESPECÍFICOS
Utilizar algunos métodos para la solución numérica de ecuaciones no
lineales.
Solución Numérica De
Cálculo infinitesimal y
Ecuaciones No E Interpolación
Métodos numéricos
Lineales
Métodos Numéricos
sfs
7
1 UNIDAD 1 SOLUCIÓN NUMÉRICA DE ECUACIONES NO
LINEALES
Video: Método de Newton para ecuaciones
”Ver Video”: https://www.youtube.com/watch?v=5boNNInwkhk
8
1.2.1.1 MANEJO DE DECIMALES
Cuando se utilizan métodos numéricos, los cálculos realizados, por lo general contienen decimales
infinitos, estos decimales se pueden tratar de la siguiente manera:
TRUNCAMIENTO:
REDONDEO:
9
1.2.2 EJERCICIOS DE APRENDIZAJE
1. Redondee el número 3,469823... con cuatro cifras decimales.
Solución
La quinta cifra es 2, por lo tanto, la cuarta cifra, que es 8, se deja igual, entonces:
3,469823... 3,4698
2. Redondee el número 0,324686... con cuatro cifras decimales.
Solución
La quinta cifra decimal es ocho, por lo tanto, se debe sumar uno a la cuarta cifra, que es 6, entonces:
0,324686... 0,3247
3. Redondee el número 2,356342 con una cifra decimal.
Solución
La segunda cifra decimal es cinco, se debe sumar uno a la primera cifra decimal, que es 3, entonces:
2,356342... 2,4
4. Trunque el número 63,45678...con dos cifras decimales.
Solución
Como el número se trunca, simplemente se escribe con las dos primeras cifras decimales:
63,45678... 63,45
1.2.2.1 CONVERGENCIA
Se dice que un método converge cuando se llega al valor buscado, esto se
puede identificar por qué se puede observar que cada cálculo efectuado lleva a
un valor cada vez más parecido al anterior.
Métodos Numéricos
sfs
10
1.2.2.2 DIVERGENCIA
UN MÉTODO DIVERGE CUANDO NO SE LLEGA AL VALOR BUSCADO, ESTO SE IDENTIFICA,
PORQUE CADA CÁLCULO EFECTUADO LLEVA A UN VALOR CADA VEZ DIFERENTE.
E Vverdadero Vcalculado
Error absoluto: Es el valor absoluto de la diferencia entre el valor verdadero y el
valor calculado o medido. Posee las mismas unidades que la magnitud que se mide.
E= Vverdadero Vcalculado
Otra manera de calcular este error es restando los dos últimos cálculos.
𝑬𝒂 = 𝑬𝒗 − 𝑬𝑪
𝑬𝒂 = 𝒙𝒏 − 𝒙𝒏−𝟏
Error relativo: Es el cociente entre el error absoluto y el valor absoluto del valor
verdadero, es un número a dimensional y en algunos casos se expresa en tanto por
ciento.
|𝑽𝒗𝒆𝒓𝒅𝒂𝒅𝒆𝒓𝒐 − 𝑽𝒄𝒂𝒍𝒄𝒖𝒍𝒂𝒅𝒐 |
𝑬𝒓 =
|𝑽𝒗𝒆𝒓𝒅𝒂𝒅𝒆𝒓𝒐 |
|𝒙𝒏 − 𝒙𝒏−𝟏 |
𝑬𝒓 =
𝒙𝒏
Métodos Numéricos
sfs
11
1.2.3 TALLER DE ENTRENAMIENTO
1) De acuerdo a lo trabajado en el desarrollo del tema anterior, realice los siguientes
ejercicios:
12
1.3.1 CONTINUIDAD Y DISCONTINUIDAD
Video: Tipos de discontinuidad AINTE Mat 1 Bach
”Ver Video”: https://www.youtube.com/watch?v=RU4v-7u4blI
Se dice que una función es continua en todo su dominio, cuando se puede recorrer toda la gráfica
sin tener que levantar la mano, cuando no hay huecos o espacios entre sus gráficas, si algo de
esto se llega a presentar se dice que la función es discontinua.
La función y = f (x) mostrada en la figura 1 es continua, porque se puede recorrer toda su gráfica
sin necesidad de levantar la mano
13
14
Función Polinómica.
El dominio de las funciones polinómicas son todos los números reales, esto se escribe así:
Dom. x lR
Las funciones polinómicas son siempre continuas
Función racional.
𝒈(𝒙)
𝒚 = 𝒇(𝒙) = 𝒉(𝒙), con 𝒉(𝒙) ≠ 𝟎
Su dominio está formado por todos los números reales menos las asíntotas verticales de la función.
Video: Problema sobre el cálculo del dominio y asíntotas de una función que es
un cociente de polinomios
”Ver Video”: https://www.youtube.com/watch?v=gSDOXzMx_h8
Para determinar las asíntotas verticales de una función racional (si tiene) se procede de la siguiente
manera:
15
SOLUCIÓN:
DOMINIO:
𝟕 𝟕 𝟕
𝒙 ≠ → 𝑫𝒇 = (−∞, ) ∪ ( , +∞)
𝟐 𝟐 𝟐
CONTINUIDAD:
𝟓𝒙 𝟕
La función: 𝒚 = 𝒘(𝒙) = es discontinua en 𝒙 =
𝟐𝒙−𝟕 𝟐
𝟒𝒙 − 𝟕
𝟐. 𝒚 = 𝒘(𝒙) =
𝟐𝒙𝟐 − 𝒙 − 𝟔
Solución:
Dominio:
Métodos Numéricos
sfs
16
Se iguala el denominador a cero:
𝟐𝒙𝟐 − 𝒙 − 𝟔 = 𝟎
Factorizando el polinomio se tiene:
𝟐(𝒙 − 𝟐)(𝟐𝒙 + 𝟑)
= 𝟎 → (𝒙 − 𝟐)(𝟐𝒙 + 𝟑) = 𝟎 →
𝟐
(𝒙 − 𝟐) = 𝟎 → 𝒙 = 𝟐
𝟑
𝟐𝒙 + 𝟑 = 𝟎 → 𝒙 = −
𝟐
Por lo tanto, el Dominio de la función está dado por:
𝟑
𝑫𝑭 = ℝ − {− , 𝟐}
𝟐
También se puede escribir con el siguiente intervalo:
𝟑 𝟑
𝑫𝑭 = (−∞, − ) ∪ (− , 𝟐) ∪ (𝟐, +∞)
𝟐 𝟐
Continuidad:
𝟒𝒙−𝟕
𝑒𝑠 𝒅𝒊𝒔𝒄𝒐𝒏𝒕𝒊𝒏𝒖𝒂 𝑒𝑛 𝑙𝑜𝑠 𝑝𝑢𝑛𝑡𝑜𝑠: 𝒙 = − 𝟐 𝒚 𝒙 = 𝟐
𝟑
La función: y w (x)
𝟐𝒙𝟐 −𝒙−𝟔
𝟑
𝒚 = 𝒇(𝒙) =
𝒙𝟐 − 𝟏
Métodos Numéricos
sfs
17
Dominio: Se toma el denominador y se hace igual a cero:
𝒙𝟐 − 𝟏 = 𝟎 → (𝒙 + 𝟏)(𝒙 − 𝟏) = 𝟎 →
𝒙 + 𝟏 = 𝟎 → 𝒙 = −𝟏
𝒙−𝟏=𝟎→𝒙=𝟏
El dominio de la función está dado por:
𝑫𝑭 = ℝ − {−𝟏, 𝟏}
También se puede escribir con el siguiente intervalo:
𝟑
La función: 𝒚 = 𝒇(𝒙) = :
𝒙𝟐 −𝟏
𝒏
𝒚 = 𝒇(𝒙) = √𝒈(𝒙)
DOMINIO:
Caso #1:
18
Caso #2:
Si la raíz es par, se debe garantizar que los valores que se le asignen a “x”
hagan positivo todo dentro de la raíz.
Ejercicios de Aprendizaje
Métodos Numéricos
sfs
19
1. Determine el dominio y la continuidad de:
𝒚 = 𝒇(𝒙) = √𝒙𝟐 + 𝟐𝒙 − 𝟏𝟓
SOLUCIÓN:
PASOS:
1. Todo lo que está dentro de la raíz debe ser mayor o igual que cero, por lo tanto, se debe
plantear la inecuación:
𝒙𝟐 + 𝟐𝒙 − 𝟏𝟓 ≥ 𝟎
1. Encuentre los números ceros de la expresión resultante. Esto es, iguale a cero y
resuelva la ecuación resultante, los valores obtenidos son los ceros de la expresión.
Un cero es un número donde posiblemente hay cambio de signo
en la expresión.
Para el ejemplo:
x 5 x 3
Estos son los ceros o números críticos.
2. Evalúe el signo de la expresión obtenida en el paso uno. Para ello se toma un número que
se encuentre a la izquierda del primer cero, se toma un número que se encuentre entre
ambos ceros y se toma un número que se encuentre a la derecha del segundo cero. Estos
números se reemplazan en la expresión obtenida en el paso uno y el signo del resultado
se escribe en la recta numérica. Véase la figura 3.
Métodos Numéricos
sfs
20
3. La respuesta o solución de la inecuación, se da tomando los intervalos que cumplan con
el sentido de la desigualdad. Para ello observe el sentido de la desigualdad de la expresión
obtenida en el paso dos.
𝒚 = 𝒇(𝒙) = √𝟑𝒙 − 𝟓
Métodos Numéricos
sfs
21
SOLUCIÓN
𝟓
3x 5 0 3x 5 0 3x 5 x
𝟑
Función Exponencial.
Es una función de la forma:
𝒚 = 𝒇(𝒙) = 𝒃𝒈(𝒙)
Métodos Numéricos
sfs
22
Con: 𝒃 >𝟎𝒚𝒃≠𝟏
DOMINIO:
CONTINUIDAD:
Ejercicios de Aprendizaje
Determine el dominio y la continuidad de las siguientes funciones:
1. 𝒚 = 𝒇(𝒙) = 𝟓𝒙
Es una función exponencial que tiene en el exponente un polinomio, su dominio
corresponde a todos los números reales.
𝒙
𝑫𝒇 = { ∈ ℝ}
𝒙
La función es siempre continua.
2. y f (x) e2x1
SOLUCIÓN:
Dominio: Como en el exponente hay una función lineal, entonces su dominio está dado por todos
𝑫𝒇 = {𝒙/𝒙 ∈ ℝ}
Continuidad: La función es continua en todo su dominio.
𝟓𝒙
3. y 𝒇(𝒙) =
𝒆𝒙−𝟒
Métodos Numéricos
sfs
23
SOLUCIÓN
Como en el exponente hay una función racional, el dominio será todos los reales menos las asíntotas
verticales.
x40x4
Dominio:
𝑫𝒇 = ℝ − {𝟒} →= (−∞, 𝟒) ∪ (𝟒, +∞)
𝒈(𝒙) > 𝟎
Continuidad:
Las funciones logarítmicas son continuas en todo su dominio.
Métodos Numéricos
sfs
24
Taller de Entrenamiento
SOLUCIÓN
DOMINIO:
𝟑
2x 3 0 2x 3 0 2x 3 𝒙 =
𝟐
𝟑
Continuidad: La función es continua en el mismo intervalo del dominio, (𝟐 , +∞)
𝟐. 𝒚 = 𝒇(𝒙) = 𝐥𝐨𝐠(𝒙𝟐 − 𝟒)
SOLUCIÓN:
Métodos Numéricos
sfs
25
DOMINIO: Para calcular el dominio se hace: x2 40
Factorizando se tiene:
x2 4 0 x 2x 2 0 x 2 0 x 2 0
x 2 x 2
Funciones trigonométricas.
Se analizará el dominio y la continuidad de las funciones trigonométricas para el
caso en que el ángulo es x.
26
𝒚 = 𝒇(𝒙) = 𝒄𝒐𝒔𝒙
Y su gráfica está determinada por:
El dominio de esta función es todos los reales, por lo tanto, esta función
es siempre continua.
Métodos Numéricos
sfs
27
1.3.3 TALLER DE ENTRENAMIENTO
Para cada una de las siguientes funciones determine: Dominio, Continuidad o discontinuidad
(intervalos donde es discontinua) y su respectiva gráfica:
Nota: Recuerde que estas funciones se pueden expresar en función de sen 𝑥 𝑦 cos 𝑥.
𝑦 = sec 𝑥
𝑦 = csc 𝑥
𝑦 = 𝑡𝑎𝑛𝑥
𝑦 = cot 𝑥
1.3.3.1 SEMILLA
Es una estimación inicial de la raíz. Es el valor inicial con el cual se empieza a operar el
método numérico que se va a aplicar. La estimación de la semilla se puede hacer de dos
formas diferentes:
Haciendo un bosquejo de la gráfica y el (los) punto(s) donde la gráfica corte el eje x es
(son) la(s) raíz (o raíces) de la ecuación, como el ojo no nos puede dar buena exactitud,
este valor sólo sirve como semilla para iniciar los cálculos.
Otra forma puede ser evaluando la función en diferentes valores de x, entre los puntos
donde se observe cambio de signo se encuentra una raíz, por lo tanto, se escoge como
semilla un valor de x que se encuentre entre dichos puntos.
3x3 9x2 3x 2 0
Su gráfica se observa en la figura 8.
Métodos Numéricos
sfs
28
La ecuación tiene raíces en los siguientes intervalos:
1, 0, 0, 1 2, 3,
Esto es cierto, ya que en estos intervalos la gráfica de la función corta
el eje x.
x -1 0 1 2 3
y -13 2 -1 -4 11
En resumen:
La ecuación 3x3 9x2 3x 2 0 tiene raíces en:
1, 0 Porque:
Métodos Numéricos
sfs
29
f (1) f (0) 0 (o tienen signo diferente) y
−𝟏 + 𝟎
𝒙𝟏 = = −𝟎, 𝟓
𝟐
0,1 Porque f (0) f (1) 0 (o tienen signo diferente) yf (x) es continua en este intervalo.
𝟏+𝟎
𝒙𝟏 = = 𝟎, 𝟓
𝟐
2, 3 Porque
f (2) f (3) 0 (o tienen signo diferente) y f (x) es continua en este intervalo
𝟐+𝟑
𝒙𝟐 = = 𝟐, 𝟓
𝟐
Métodos Numéricos
sfs
30
Si una función f (x) es continua en un intervalo a, b y además, f a y f b tienen:
diferente signo (esto se simboliza como f af b 0 ) entonces, la ecuación f (x) 0 tiene
una raíz en a, b ; la semilla para determinar esta raíz se calcula como:
𝒂+𝒃
𝒙=
𝟐
2. Determine en que intervalos se encuentran las raíces de la ecuación:
4x4 9x 7 0.
Determine la semilla para cada raíz.
SOLUCIÓN
Por ser
Métodos Numéricos
sfs
31
f (x) 4x4 9x 7 una función polinómica, es continua para cualquier valor de x:
x 2 y 4 2 9(2) 7 y 75
4
x 1 y 41 9(1) 7 y 6
4
x 0 y 40 9(0) 7 y 7
4
x 2 y 42 9(2) 7 y 39
4
x -3 -2 -1 0 1 2 3
Se concluye que:
−𝟏+𝟎
𝒙𝟏 = = −𝟎, 𝟓
𝟐
Métodos Numéricos
sfs
32
𝟏+𝟐
𝒙𝟐 = = 𝟏, 𝟓
𝟐
3. Determine en que intervalos se encuentran las raíces de la ecuación 5 ln x 3x 4 0
Determine la semilla para cada raíz.
SOLUCIÓN:
x 1 y 5ln1 3(1) 4 1
33
Se llena la siguiente tabla.
x -0,01 0,5 1 2 3 4
−𝟎. , 𝟏 + 𝟎, 𝟓
𝒙𝟏 = = 𝟎, 𝟐𝟒𝟓
𝟐
Raíz número 2 está en el intervalo [3, 4, la semilla es:
𝟑+𝟒
𝒙𝟐 = = 𝟑. 𝟓
𝟐
1. 5x3 10x2 8x 3 0
2. 2 x 3x 2 5 0
3. x2 3 x 6 0
34
1.4 TEMA 3 MÉTODOS NUMÉRICOS PARA SOLUCIONAR
ECUACIONES NO LINEALES
En los cinco métodos que se explica a continuación, se parte del mismo problema:
f (x) 0
Método Bisección
𝒂 𝟏 + 𝒃𝟏
𝒙𝟏 =
𝟐
Cerciórese que se cumpla que:
35
PASO TRES: Calcule f x1 . Pueden resultar tres posibilidades:
POSIBILIDAD CARACTERÍSTICA
f (x1) f a1 0 Esto quiere decir que la raíz se encuentra en el intervalo a1 , x1 .
o Para este caso haga a2 =a1 b2 x1
f (x1 ) f b1 0
Esto quiere decir que la raíz se encuentra en el intervalo x1 ,
b1 .
𝒂𝟐 + 𝒃𝟐
𝒙𝟐 =
𝟐
Regrese al paso 3.
Encuentre una raíz real de la ecuación: f (x) 5x3 16x2 6x 4 0 con una precisión de
0.01,
36
x -3 -2 -1 0 1 2 3
y -13 8 1 -4 23 112 293
Al observar la tabla se puede ver que la ecuación tiene tres aíces, así: Raíz número 1
Raíz número 2
Está en el intervalo 1, 0esto se puede afirmar, ya que f (x) 5x3 16x2 6x 4
continua en el intervalo 1, 0 y porque:
Raíz número 3
f 0 f 1 0
Procedimiento para hallar la raíz número 3 con una precisión de 0,01 o 10-2.
𝟎+𝟏
𝒙𝟏 = = 𝟎, 𝟓
𝟐
Paso tres. Se calcula f (x1 )
Métodos Numéricos
sfs
37
Como f 0,5 f 10
𝟎 + 𝟎, 𝟓
𝒙𝟏 = = 𝟎, 𝟐𝟓
𝟐
Comparando x2 0,25 con x1 0,5 , no se repite ninguna cifra decimal, por lo tanto, se debe
continuar el proceso en el paso tres:
Cálculo de f x2
[0,25, 0,5]
Cálculo de x3
𝟎, 𝟐𝟓 + 𝟎, 𝟓
𝒙𝟑 = = 𝟎, 𝟑𝟕𝟓
𝟐
Observando 𝑥3 ∧ 𝑥2 no se repite decimal
Como
Cálculo de
𝟎, 𝟐𝟓 + 𝟎, 𝟑𝟕𝟓
𝒙𝟒 = = 𝟎, 𝟑𝟏𝟐𝟓
𝟐
Métodos Numéricos
sfs
38
Observando x4 x3 se repite un decimal (deben ser dos). Se continúa el proceso.
Se calcula: f x4
Cálculo de x5
𝟎, 𝟑𝟏𝟐𝟓 + 𝟎, 𝟑𝟕𝟓
𝒙𝟓 = = 𝟎, 𝟑𝟒𝟑𝟕𝟓
𝟐
𝟎, 𝟑𝟏𝟐𝟓 + 𝟎, 𝟑𝟒𝟑𝟕𝟓
𝒙𝟔 = = 𝟎, 𝟑𝟐𝟖𝟏𝟐𝟓
𝟐
Observando x6 x5 se repite un decimal (deben ser dos). Se continúa el proceso.
39
Observando x7 x6 se repite un decimal (deben ser dos). Se continúa el proceso.
f 0.3359375
50.3359375 160.3359375 60.3359375 4 0.0108485218..
3 2
La raíz es 𝒙𝟖 0.33203125
Queda como ejercicio hallar las otras dos raíces utilizando el mismo error.
2. Con una precisión de 10-3 encuentre una raíz de real de la ecuación ln x cos x 0
Solución:
40
cada función que la conforman, es decir, Dom. x 0,
Para determinar en qué intervalo hay una raíz, se debe dar valore a x desde cero
(sin incluirlo) hasta que se obtenga un cambio de signo.
Antes de dar valores a x, no olvide que la calculadora debe estar en modo radianes.
Cálculo de la semilla
𝟏 + 𝟏, 𝟓
𝒙𝟏 = = 𝟏, 𝟐𝟓
𝟐
Cálculo de f x1 f x1 ln1,25 cos1,25 -0. 09217818108..
Cálculo de x2
𝟏, 𝟐𝟓 + 𝟏, 𝟓
𝒙𝟐 = = 𝟏, 𝟑𝟕𝟓
𝟐
Métodos Numéricos
sfs
41
Comparando x2 1,375con x1 1,25, no se repite ninguna cifra decimal, por lo tanto, se debe
continuar el proceso.
Cálculo de x3
𝟏, 𝟐𝟓 + 𝟏, 𝟑𝟕𝟓
𝒙𝟑 = = 𝟏, 𝟑𝟏𝟐𝟓
𝟐
Cálculo de 𝒙𝟒
𝟏, 𝟐𝟓 + 𝟏, 𝟑𝟏𝟐𝟓
𝒙𝟒 = = 𝟏, 𝟐𝟖𝟏𝟐𝟓
𝟐
42
𝟏, 𝟐𝟖𝟏𝟐𝟓 + 𝟏, 𝟑𝟏𝟐𝟓
𝒙𝟓 = = 𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓
𝟐
Cálculo de x6
𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓 + 𝟏, 𝟑𝟏𝟐𝟓
𝒙𝟔 = = 𝟏, 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝟐
Observando x6 x5 no se repite decimal (deben ser tres). Se continúa el proceso.
Cálculo de x7
𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓 + 𝟏, 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝒙𝟕 = = 𝟏, 𝟑𝟎𝟎𝟕𝟖𝟏𝟐𝟓
𝟐
43
-0. 003780291734..
(-) (+)
Cálculo de x8
𝟏. 𝟑𝟎𝟎𝟕𝟖𝟏𝟐𝟓 + 𝟏. 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝒙𝟖 = = 𝟏, 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓
𝟐
(-) (+)
Cálculo de x9
𝟏, 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓 + 𝟏. 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝒙𝟗 = = 𝟏. 𝟑𝟎𝟑𝟕𝟏𝟎𝟗𝟑 𝟕
𝟐
44
Como f 1.302734375 f 1.303710937 0 , la raíz está en el intervalo:
(-) (+)
Cálculo de x10
𝟏, 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓 + 𝟏. 𝟑𝟎𝟑𝟕𝟏𝟎𝟗𝟑 𝟕
𝒙𝟏𝟎 = = 𝟏, 𝟑𝟎𝟑𝟐𝟐𝟐𝟔𝟓𝟓
𝟐
Observando x9 x10
𝒙𝟏𝟎 == 𝟏, 𝟑𝟎𝟑𝟐𝟐𝟐𝟔𝟓𝟓
Si la función 𝑓(𝑥) es continua en algún intervalo que contenga la raíz buscada, este método
siempre converge.
Sea f (x) una función continua en el intervalo a, bcon f af b 0, sea
r una raíz de f x 0 en el intervalo a, b.
Véase la figura 9
Métodos Numéricos
sfs
45
El método consiste en encontrar la ecuación de la recta que pasa por los puntos:
46
((𝒂𝟏, 𝒇 (𝒂𝟏 ) ∧ (𝒃𝟏,𝒇( 𝒃𝟏 ))
𝑓(𝑏)1 −𝑓(𝑎1 )
La pendiente de la recta está dada por: 𝑚 = 𝑏1 −𝑎1
La ecuación es:
𝒃𝟏 −𝒂𝟏
Cuando:
𝒙 = 𝒙𝟏 , 𝒚 = 𝟎
La ecuación queda:
𝑓 (𝑏1 ) − 𝑓 (𝑎1 )
0 − 𝑓(𝑎1 ) (𝑥1 − 𝑎1 )
𝑏1 − 𝑎1
Despejando x1
𝑏1 − 𝑎1
𝑥1 = 𝑎1 − 𝑓(𝑎1 ) [ ]
𝑓 (𝑏1 ) − 𝑓 (𝑎1 )
𝑎1 𝑓 (𝑏1 ) − b1 𝑓 (𝑎1 )
𝑥1 =
𝑓 (𝑏1 ) − 𝑓 (𝑎1 )
PASO TRES: Halle
Si f a1 f x1 0 , quiere decir que la raíz se encuentra en a1 , x1 en este caso,
haga:
a2 a1 b2 x1
Métodos Numéricos
sfs
47
De lo contrario tome a2 x1 b2 b1
PASO CUATRO: Se aplica nuevamente el proceso anterior al intervalo a2 , b2 , es decir, haga.
𝑎2 𝑓 (𝑏2 ) − b2 𝑓 (𝑎2 )
𝑥=
𝑓 (𝑏2 ) − 𝑓 (𝑎2 )
PASO CINCO: Después de varias iteraciones se tiene que:
𝑎n 𝑓 (𝑏n ) − bn 𝑓 (𝑎n )
𝑥=
𝑓 (𝑏n ) − 𝑓 (𝑎n )
El proceso se termina cuando:
1. Utilizando el método “Regula Falsi” con un error de 10-2 obtenga una de las
raíces reales de la ecuación:
3x4 2x 9 0
PROCEDIMIENTO
La función polinómica. f x 3x4 2x 9 es continua para cualquier valor de x, por ser una
función polinómica.
x -2 -1 0 1 2
y 35 -8 -9 -4 43
NOTA:
Métodos Numéricos
sfs
48
La ecuación no tiene más raíces reales, ya que, es posible comprobar, dando
más valores a x, que a la izquierda de x = - 2 la función asignará a y solo valores
positivos (no hay más cambio de signos hacia la izquierda) y a la derecha de x =
2 tampoco habrá más cambio de signos. Esto se puede asegurar, porque el signo
que predomina en ambos costados es el de 3x 4 que es positivo.
Raíz número 1: Esta en [−2 − 2.] ya que f 2 f 1 0 f x es continua en 2, 1
+ −
Raíz número 2: Esta en [1; 2.] ya que f 1 f 2 0 f x es continua en 1, 2
− +
x 1 32 22 9 2 31 21 9 / 32 22 9 31 21 9
4 4 4 4
𝒙𝟏 1.08510638
1
Cálculo de f (x1 )
[1.085106382 … , 2]
− +
Cálculo de x2
Métodos Numéricos
sfs
49
x 1 .085106328... 32 22 9 2 31.085106328... 21.085106328... 9
4 4
2
/ 32 22 9 31.085106328... 21.085106328... 9
4 4
x2 1.138604649...
Cálculo de f x2
[1.138604649 … , 2]
− +
Cálculo de x3
x 1 .138604694... 32 22 9 2 31.138604694... 21.295463593... 9
4 4
3
/ 32 22 9 31.138604694... 21.138604694... 9
4 4
x3 1.17100619...
50
1.1 1.1
𝑓(𝑥3 ) = 3 ( ) + 2( ) − 9 = −1.01696214
1.00619 … 1.00619 …
1.1
La raíz está en: [ … 2]
1.006199…
− +
Cálculo de x4
𝑥4 = 1 .17100619... 32 22 9 2 31.17100619... 21.17100619... 9
4 4
4
/ 32 22 9 31.17100619... 21.17100619... 9
4 4
𝒙𝟒 = 𝟏. 𝟏𝟗𝟎𝟏𝟓𝟗𝟏𝟕𝟎 …
Comparando 𝑥4 ∧ 𝑥3 Se repite un decimal, el proceso continúa.
Cálculo de 𝑓(𝑥4 ):
𝑓(𝑥4 ) = 3(1.190159170)4 + 2(1.190159170) − 9 = −0.600444630
− +
Cálculo de x5
x 1 .190159107... 32 22 9 2 31.190159107... 21.190159107... 9
4 4
5
/ 32 22 9 31.190159107... 21.190159107... 9
4 4
Métodos Numéricos
sfs
51
𝑥5 1.20131194...
Comparando 𝑥5 ∧ 𝑥4 Se repite un decimal, el proceso continúa.
Cálculo de 𝑓(𝑥5 ):
𝑓(𝑥5 ) = 3(1.20131194 … ) + 2(1.20131194 … ) − 9 = 1.20131194 …
− +
Cálculo de x6
x 1 .201311941... 32 22 9 2 31.201311941... 21.201311941... 9
4 4
6
/ 32 22 9 31.201311941... 21.201311941... 9
4 4
x6 1.20774808…
Comparando 𝑥6 ∧ 𝑥5 Se repite dos decimales, el proceso terminó.
Con una precisión de 10-2, la raíz de la ecuación 3x4 2x 9 0 es x6 1.20774808..
Queda como ejercicio encontrar la raíz número 1 utilizando el mismo método y la misma
precisión.
Si f (x) es una función continua en un intervalo a, b, existe un número x en a, b tal
que:
Métodos Numéricos
sfs
52
𝑓(𝑏) − 𝑓(𝑎)
𝑓 (𝑥 ) =
𝑏−𝑎
Utilizando el método de la regla falsa, con una precisión de 10-5, encuentre el número
( 2
´( )
5𝑥 𝑥 − 4) − 2𝑥(5𝑥 − 3)
𝑓 𝑥 =
(𝑥 2 − 4)2
Realizando las operaciones indicadas y reduciendo términos semejantes, se tiene:
2
´( )
−5𝑥 + 6𝑥 − 20
𝑓 𝑥 =
(𝑥 2 − 4)2
Se sabe que:
𝑓(𝑏) − 𝑓(𝑎) 83
=−
𝑏−𝑎 160
Métodos Numéricos
sfs
53
se debe plantear la ecuación:
−5𝑥 2 + 6𝑥 − 20 83
= −
(𝑥 2 − 4)2 160
Ejercicio de Entrenamiento:
54
Video: MÉTODO DE NEWTON
”Ver Video”: https://www.youtube.com/watch?v=_qNllgwlXwc
Se tiene que 𝒙𝟏 es una aproximación inicial para la raíz 𝒓 como se ve en la figura 10.
Observe que la línea tangente a la curva en
(𝒙𝟏 , f (𝒙𝟏 )) corta el eje x en el punto (𝒙𝟐 , 𝟎) ,0) y que 𝒙𝟐 es una mejor aproximación a 𝒓
que 𝒙𝟏 .
La ecuación de esta línea tangente es: y f (𝒙𝟏 ) f '(𝒙𝟏 )(x 𝒙𝟏 ) Se puede encontrar
𝒙𝟐 utilizando la ecuación de la línea tangente.
Como 𝒙𝟐 , 𝟎 está en la línea tangente, sus coordenadas deben satisfacer la
ecuación anterior, por lo tanto:
𝒇 (𝒙𝟏 )
𝟎 − f (𝒙𝟏 )= f '(𝒙𝟏 )( 𝒙𝟐 𝒙𝟏 )− =𝒙𝟐 − 𝒙𝟏 ⟹
𝒇 ′(𝒙𝟏 )
𝒇 (𝒙𝟏 )
𝒙𝟐 = 𝒙𝟏 −
𝒇 ′(𝒙𝟏 )
𝒇 (𝒙𝟐 )
𝒙𝟑 = 𝒙𝟐 −
𝒇 ′(𝒙𝟐 )
Generalizando el proceso se puede decir que:
𝒇 (𝒙𝒏 )
𝒙𝒏+𝟏 = 𝒙𝒏 − , donde 𝒏 = 𝟏, 𝟐, 𝟑, …
𝒇 ′(𝒙𝒏 )
55
Figura 10. Determinación del método de Newton (Autor Elkin Ceballos Gómez)
Para determinar intervalos para las raíces, se da valores a x. De esta manera se llena la siguiente tabla:
x 0 1 2 3 4 5
y f 0 20 f 1 16.8 f 2 12.9 f (3) 8.1 f 4 3.2 f 5 -0.9
Métodos Numéricos
sfs
56
4, 5
Se observa el cambio de signo en el intervalo: [ ]
+ −
4, 5
La ecuación 𝑓 (𝑥 ) = 𝑠𝑒𝑛𝑥 − 4𝑥 + 20 = 0 tiene una raíz en el intervalo[ ]
+ −
ya que:
4, 5
𝑓(4). 𝑓(5) < 0 𝑦 𝑓 (𝑥 )𝑒𝑠 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑎 𝑒𝑛 [ ]
+ −
𝑓 1 (𝑥 ) = 𝑐𝑜𝑠𝑥 − 4
Se realiza el cálculo de la semilla:
𝟒+𝟓 𝟗
𝒙𝟏 = = = 𝟒. 𝟓
𝟐 𝟐
Cálculo de 𝒙𝟐 :
𝒇 (𝒙𝟏 ) 𝒔𝒆𝒏(𝟒. 𝟓) − 𝟒(𝟒. 𝟓) + 𝟐𝟎
𝒙𝟐 = 𝒙𝟏 − = 𝟒. 𝟓 − = 𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 …
𝒇𝟏 (𝒙𝟏 ) 𝒄𝒐𝒔(𝟒. 𝟓) − 𝟒
57
𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … − (𝒔𝒆𝒏(𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … ) − 𝟒(𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … ) + 𝟐𝟎
𝑥3 =
(𝒄𝒐𝒔(𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … ) − 𝟒)
𝑥3 = 4.75017167 …
𝑥4 = 4.7501784845
𝑥5 = 4.75017848.5…
Utilizando el método de Newton, con una precisión de 10-8 Encuentre una raíz real de la ecuación:
3x2 e x 6 0
Por tener una combinación de función polinómica y función exponencial (con polinomio en el
exponente la función f (x) 3x2 ex 6 es continua en todo su dominio.
Es fácil comprobar que la función tiene una raíz en el intervalo 2, 3, ya que:
Métodos Numéricos
sfs
58
f 2 32 e 2 6 -1.38905698
2
Determinación de la derivada:
Determinación de la semilla:
2+3 5
𝑥1 = = = 2.5
2 2
Nota: Para escribir la expresión e 2.5 en la calculadora, debe digitar las
teclas:
SHIFT ln 2.5
Cálculo de 𝑥2
𝒇 (𝒙𝟏 ) 𝟑(𝟐. 𝟓)𝟐 − 𝒆𝟐.𝟓 − 𝟔
= 𝒙𝟏 − = 𝟐. 𝟓 −
𝒇𝟏 (𝒙𝟏 ) 𝟔(𝟐. 𝟓) − 𝒆𝟐.𝟓
𝒙𝟐 = 𝟐. 𝟐𝟗𝟖𝟓𝟕𝟖𝟓𝟗 …
x 2.5 32.5 SHIFTln2.5 6 /62.5 SHIFTln2.5
2
2
x2 2.29857859...
Métodos Numéricos
sfs
59
Cálculo de 𝑥3 :
𝑥3 = =
2.29857859... 32.29857859... SHIFTln2.29857859... 6
2
/ 62.29857859... SHIFTln2.29857859...
𝑥3 = 2.3271903...
Cálculo de 𝑥4 :
𝑥4 = =
2.3271903... 32.3271903...) SHIFTln2.3271903…) 6
2
/ 62.3271903… SHIFTln2.3271903…)
𝑥4 = 2.327637268...
Cálculo de 𝑥5 :
𝑥5 = 2.327637268... 32.327637268...) 2
SHIFTln2.327637268... ) 6
/ 62.327637268... SHIFTln2.327637268... )
𝑥5 = 2.3276373282..
Cálculo de 𝑥6 :
2
𝑥6 = 2.327637282... 32.327637282...) SHIFTln2.327637282...) 6
/ 62.327637282... SHIFTln2.327637282... )
𝑥6 = 2.327637382..
Métodos Numéricos
sfs
60
Observando los dos últimos cálculos, se repiten todos los decimales.
3. Problemas de cajas:
Se desea construir una caja sin tapa y de forma rectangular con una lámina de cartón d e 50 cm *
35 cm. Para ello se cortarán cuadrados idénticos en las cuatro esquinas y se doblarán los lados hacia
arriba
Si la caja será hecha para contener un volumen 𝟑𝟓𝟎𝟎𝒄𝒎𝟑 con una precisión de 8 decimales
determine.
a. Las dimensiones de la caja.
b. La cantidad de material utilizado.
PROCEDIMIENTO:
x: Altura.
50 2x: Largo de la caja
Se plantea la ecuación:
61
Solucionando la ecuación presentada se tiene que las raíces se encuentran en 2, 3, 11,
12 27, 28
Queda como ejercicio de entrenamiento encontrar las raíces con el método de
Newton.
El valor de x 27.9829 queda descartado, ya que daría longitudes negativas para el largo y
para el ancho.
Con x 2.631
x: Altura. 2.631
Con x 11.8864
x: Altura. 11.886
La función de ingreso en cierta empresa es: rq 150q 0.2q2 $ y la función de costo
total es: cq q3 80q 600 $ donde q es el número de docenas de cierto artículo
Métodos Numéricos
sfs
62
producidas y vendidas. Determine el punto de equilibrio con una precisión de 4
decimales.
rq cq
Esto es:
63
1.4.6 MÉTODO DE LA SECANTE
1 𝐿𝑖𝑚 𝑓 (𝑥 ) − 𝑓((𝑥𝑛+1 )
𝑓 (𝑥𝑛+1 ) =
𝑥 → (𝑥𝑛+1 ) 𝑥 − 𝑥𝑛+1
Haciendo x xn y evaluando el límite, tenemos que:
𝒇(𝒙𝒏 ) − 𝒇((𝒙𝒏+𝟏 )
𝒇𝟏 (𝒙𝒏+𝟏 ) ≈
𝒙𝒏 − 𝒙
𝒇(𝒙𝒏+𝟏 )
𝒙𝒏+𝟐 = 𝒙𝒏+𝟏 − 𝟏
𝒇 (𝒙𝒏+𝟏 )
𝟏
Reemplazando 𝒇 (𝒙𝒏+𝟏 ) Se tiene que:
𝒇(𝒙𝒏+𝟏 )
𝒙𝒏+𝟐 = 𝒙𝒏+𝟏 −
𝒇(𝒙𝒏 ) − 𝒇((𝒙𝒏+𝟏 )
𝒙𝒏 − 𝒙𝒏+𝟏
Métodos Numéricos
sfs
64
Utilizando el método de la secante, con una precisión de 10−6 encuentre una raíz real para la
ecuación:
32x 35x2 3 0
Procedimiento:
x = -1 y x = 0,
Métodos Numéricos
sfs
65
x=0yx=1,y
x = 2 y x = 3.
]
−1, 0
Raíz número 1: Está en [
− +
0, 1
Raíz número 2: Está en [ ]
+ −
]
2, 3
Raíz número 3: Está en [
− +
𝒙𝒐 2 𝒙𝟏 3
Cálculo de 𝒙𝟐
𝟑. 𝒇(𝟐) − 𝟐. 𝒇(𝟑)
𝒙𝟐 =
𝒇(𝟐) − 𝒇(𝟑)
𝒙𝟐 2.118393243...
Actividad: Verifique este resultado calculando 𝒇(𝟐)⋀𝒇(𝟑)
Cálculo de 𝒙𝟑
Métodos Numéricos
sfs
66
(𝒙𝟐 ). 𝒇(𝒙𝟏 ) − 𝒙𝟏 . 𝒇(𝒙𝟐 )
𝒙𝟑 =
𝒇(𝒙𝟏 ) − 𝒇(𝒙𝟐 )
Reemplazando los valores correspondientes, se obtiene que:
𝒙𝟑 = 𝟐. 𝟐𝟏𝟏𝟎𝟗𝟒𝟑𝟕𝟏. ..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.
Cálculo de 𝒙𝟒
𝒙𝟓 = 𝟐. 𝟑𝟓𝟓𝟓𝟑𝟖𝟐. 𝟗..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.
Métodos Numéricos
sfs
67
Comparando 𝒙𝟓 𝒙𝟒 no se repiten decimales, el proceso continua.
Cálculo de 𝒙𝟔
𝒙𝟔 = 𝟐. 𝟑𝟗𝟖𝟑𝟒𝟎𝟖. 𝟓..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.
𝒙𝟖 =2.4132248.9..
Métodos Numéricos
sfs
68
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.
𝒙𝟗 = 𝟐. 𝟒𝟏𝟑𝟐𝟔𝟐𝟔. 𝟐..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.
𝒙𝟏𝟎 =2.4132627.2..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.
69
Procedimiento
f (x) x cosx
Al hacer dar valores a x se tiene que hay cambio de signo en x = 0.5 y x = 0.9, se pueden tomar:
𝒙𝟎 = 0.5 𝒙𝟏 = / 4 0.78539816
Para 𝒏 = 𝟎 se tiene:
𝒙𝟐 0.7363841
n xn
0 x0 0,5
1 x1 0,78539816
2 x2 0,73638414
3 x3 0,73905814
4 x4 0,73908515
5 x5 0,73908513
Métodos Numéricos
sfs
70
6 x6 0,73908513
7 x7 0,73908513
Se puede observar que después de la quinta iteración se obtiene la respuesta con la exactitud
deseada, por lo tanto, la raíz de esta ecuación es: x5 0.7390851.
y f x 3x 4
71
x g(x) consiste en encontrar una raíz r que no cambie mediante la función g. Al número r se
le llama punto fijo.
𝑎+𝑏
2. Halle la semilla para la primera raíz como: 𝑥1 = 2
3. Lleve la ecuación 𝒇(𝒙) = 𝟎 𝑎 𝑙𝑎 𝑓𝑜𝑟𝑚𝑎 𝒙 = 𝒈(𝒙), es decir despeje 𝒙. Tenga presente que
pueden existir muchas formas de obtener 𝒙 = 𝒈(𝒙).
Dando valores a x en la ecuación original, es decir, en ecuación tiene una raíz en el intervalo 0, 1, ya
que:
Cálculo de la semilla:
Métodos Numéricos
sfs
72
𝑎+𝑏 1+0 1
𝑥1 = → 𝑥1 = → 𝑥1 = → 𝑥1 = 0.5
2 2 2
Aplicando el algoritmo, resulta que:
73
Luego se lleva la ecuación a la forma x g(x)
Una posible forma es x 2 cos x pero al calcular los siguientes valores de x se observa una
oscilación, lo que indica que no hay convergencia, por lo tanto con esta forma de:
x gx 2cosx, no se obtiene una solución a la ecuación:
x 2 cos x 0.
Se debe escribir la ecuación de otra forma:
𝑥+2𝑐𝑜𝑠𝑥
Puede ser: 𝑔(𝑥) = 2
1+2𝑐𝑜𝑠1
x2 g(x1 ) x2𝑔(1) = = 1.040302306
2
1.040302306 +2𝑐𝑜𝑠1.040302306
x3 g(x2 ) x3= = 1.026110667
2
1.026110667+2𝑐𝑜𝑠1.026110667
x4 g(x3 ) x4= =
2
1.031204594+2𝑐𝑜𝑠1.031204594
x5 g(x4 ) x5= =
2
1.029388072+2𝑐𝑜𝑠1.029388072
x6 g(x5 ) x6= =
2
1.03003739+2𝑐𝑜𝑠1.03003739
x7 g(x6 ) x7= =
2
1.029805485+2𝑐𝑜𝑠1
x8 g(x7 ) x8= = 1.02988833.5
2
raíz x8 1.02988833.5
CONVERGENCIA DEL MÉTODO
Prueba de convergencia.
La g(x) que usted encontró se puede utilizar si y sólo si g'(x) 1 para todo número en a, b
Actividad: comprobar la convergencia en los dos ejemplos anteriores.
Taller de Entrenamiento
1. Utilizando al menos dos métodos diferentes, con una precisión de 10-4, encuentre todas
Métodos Numéricos
sfs
74
las raíces de las siguientes ecuaciones:
a. 5cosx 7x 2 4 0
b. 4x 4 10x3 5x 2 9x 15 0
c. ln(3x) 5x2 4 0
3. ANUALIDADES
Al comprar a crédito, en muchos casos el pago se efectúa mediante anualidades; Las
ANUALIDADES, son cuotas de dinero periódicas e iguales que se entregan o se reciben al
comienzo o al final de cada período de tiempo. Para calcular el valor presente de las
anualidades se utiliza la fórmula:
(𝟏 + 𝒊)𝒏 − 𝟏
𝒑 = 𝑨[ ]
𝒊(𝟏 + 𝒊)𝒏
Donde:
p : El valor presente, es el dinero que se tendría que pagar de contado, es decir, sin financiación.
i : Es la tasa de interés.
Con una precisión de 10-4 y utilizando el método que desee solucione: Un automóvil de
contado cuesta $ 25’000 000, pero se puede pagar de la siguiente manera: Una cuota
75
Encuentre el número 𝒙 garantizado por el teorema del valor medio para la función
𝒙
𝒇(𝒙) = 𝟒 𝒆𝒏 𝒆𝒍 𝒊𝒏𝒕𝒆𝒓𝒗𝒂𝒍𝒐 [𝟎, 𝟏]
𝒙 +𝟏
7. A partir de una lámina de cartón de 15 cm por 10 cm se desea construir una caja, para ello
se debe cortar cuadrados idénticos en las cuatro esquinas y se doblan los lados hacia
arriba. La caja debe ser hecha para contener un volumen de 60 cm3, determine las
dimensiones de la caja. Seleccione aquellas dimensiones que permitan obtener la mayor
altura.
Métodos Numéricos
76
2 UNIDAD 2 REGRESIÓN E INTERPOLACIÓN
Definición de Determinante
77
Si A es una matriz cuadrada de orden n, su determinante se puede simbolizar como
det A o |𝐀|
Cálculo de determinantes
𝑎 𝑏
𝐴=[ ]
𝑐 𝑑
Recuerde que el determinante se calcula realizando la diferencia entre el producto de los
elementos de la diagonal principal y el producto de los elementos de la diagonal secundaria, esto
es:
𝑑𝑒𝑡 𝐴 = 𝑎. 𝑑 − 𝑏. 𝑐
Ejercicios:
6 −2
1. 𝑑𝑒𝑡 = [ ] = (6). (1) − (−2). (−3) = 6 − 6 = 0
−3 1
2 −3
2. 𝑑𝑒𝑡 = [ ] = (2). (−5) − (−3). (4) = −10 + 12 = 2
4 −5
78
Sea la matriz A de orden 3 × 3: Se tiene que para calcular su determinante se repiten las dos
primeras columnas (o las dos primeras filas) y luego se realiza la suma del producto de las
diagonales menos la suma del producto de las diagonales secundarias.
79
(𝒂𝟏𝟏 . 𝒂𝟐𝟐 . 𝒂𝟑𝟑 + 𝒂𝟏𝟐 . 𝒂𝟐𝟑 . 𝒂𝟑𝟏 + 𝒂𝟏𝟑 . 𝒂𝟐𝟏 . 𝒂𝟑𝟐 ) − (𝒂𝟏𝟑 . 𝒂𝟐𝟐 . 𝒂𝟑𝟏 + 𝒂𝟏𝟏 . 𝒂𝟐𝟑 . 𝒂𝟑𝟐 +
𝒂𝟏𝟐 . 𝒂𝟐𝟏 . 𝒂𝟑𝟑 )
Ejercicio de Aprendizaje
2 −3 1
𝐴 = [4 0 −2]
3 −1 −3
2 −3 1 2 −3
𝑑𝑒𝑡𝐴 = [3 0 −2] 4 0
3 −1 −3 3 −1
𝑑𝑒𝑡𝐴 = 14 − 40 = −26
Métodos Numéricos
sfs
80
Actividad: calcula el determinante para la siguiente matriz, siguiendo el proceso anterior:
7 −2 5
𝐴 = [9 3 −2]
8 7 5
2.2.1.3 CASO 3: DETERMINANTE DE MATRICES 𝒏 × 𝒏. DESARROLLO POR
COFACTORES
Dada una matriz cuadrada 𝐴 de orden 𝒏 × 𝒏:
7 −2 5
𝐴 = [9 3 −2]
8 7 5
Obtenga la matriz menor para cada una de las filas:
7 −2 5
𝐴11 = [9 3 −2]
8 7 5
Quedaría:
3 −2
𝐴11 = [ ]
7 5
Para la segunda fila: Se elimina su fila y su respectiva columna:
Métodos Numéricos
sfs
81
7 −2 5
𝐴21 = [9 3 −2]
8 7 5
7 5
𝐴21 = [ ]
8 5
Para la tercera fila: Se elimina su fila y su respectiva columna:
7 −2 5
𝐴31 = [9 3 −2]
8 7 5
−2 5
𝐴31 = [ ]
3 −2
Nota: Para calcular el determinante por menores, se realiza el siguiente proceso:
Los signos de los productos se colocan alternos empezando por el signo más.
En el ejemplo, se tiene:
3 −2 7 5 −2 5
7[ ] − 9[ ]+8 [ ]=
7 5 8 5 3 −2
𝟕(𝟑 × 𝟓 − (−𝟐) × 𝟕) − 𝟗(𝟕 × 𝟓 − 𝟓 × 𝟖) + 𝟖(−𝟐 × (−𝟐) − 𝟓 × 𝟑) =
82
2.2.2 EJERCICIO DE ENTRENAMIENTO
1. Dada la siguiente matriz, encuentra el determinante por menores de los elementos de
cada una de las columnas y compara el resultado con el ejercicio anterior.
7 −2 5
𝐴 = [9 3 −2]
8 7 5
2. ¿Cómo se calcularía el determinante de una matriz de orden 4 × 4?Realiza un ejemplo.
COFACTOR: Cij
𝟑 𝟓 𝟑
𝑨 = [𝟑 𝟎 𝟖]
𝟐 𝟐 𝟑
Métodos Numéricos
sfs
83
Procedimiento:
𝟑 𝟓
𝑪𝟐𝟑 = (−𝟏)𝟐+𝟑 [ ] = (−𝟏)𝟓 [𝟑(𝟐) − 𝟐(𝟓)] =
𝟐 𝟐
−𝟏 (𝟔 − 𝟏𝟎) = −𝟏(−𝟒) = 𝟒
---------------------------------------------------------
𝟎 𝟖
𝑪𝟏𝟏 = (−𝟏)𝟏+𝟏 [ ] = (−𝟏)𝟐 [𝟎(𝟑) − 𝟐(𝟖)] =
𝟐 𝟑
𝟔 −𝟑 𝟓
𝑨 = [−𝟒 −𝟗 𝟖]
𝟓 𝟑 𝟕
Procedimiento:
𝟔 −𝟑
𝑪𝟐𝟑 = (−𝟏)𝟐+𝟑 [ ] = (−𝟏)𝟓 [𝟔(𝟑) − 𝟓(−𝟑)] =
𝟓 𝟑
84
Para calcular determinantes de matrices n x n se utiliza el método: CÁLCULO DE DETERMINANTES
POR EXPANSIÓN POR COFACTORES que consiste en:
𝟐 −𝟏 𝟑
𝑨 = [𝟑 𝟎 −𝟓]
𝟐 𝟏 𝟏
Procedimiento
𝟐 −𝟏 𝟑 𝟐 −𝟏
[𝟑 𝟎 −𝟓] 𝟑 𝟎
𝟐 𝟏 𝟏 𝟐 𝟏
𝒅𝒆𝒕 𝑨 = 𝟐(𝟎)(𝟏) + (−𝟏)(−𝟓)(𝟐) + 𝟑(𝟑)(𝟏) −2031 52 131
𝒅𝒆𝒕 𝑨 = 𝟎 + 𝟏𝟎 + 𝟗 − 𝟎 + 𝟏𝟎 + 𝟑 → 𝒅𝒆𝒕 𝑨 = 𝟑𝟐
𝒅𝒆𝒕𝑨 = 𝒂𝟏𝟏 (−𝟏)𝟏+𝟏 𝑴𝟏𝟏 + 𝒂𝟏𝟐 (−𝟏)𝟏+𝟐 𝑴𝟏𝟐 + 𝒂𝟏𝟑 (−𝟏)𝟏+𝟑 𝑴𝟏𝟑
Métodos Numéricos
sfs
85
𝟎 −𝟓 𝟑 −𝟓 𝟑 𝟎
𝒅𝒆𝒕𝑨 = 𝒂𝟏𝟏 (−𝟏)𝟏+𝟏 [ ] + 𝒂𝟏𝟐 (−𝟏)𝟏+𝟐 [ ] + 𝒂𝟏𝟑 (−𝟏)𝟏+𝟑 [ ]
𝟏 𝟏 𝟐 𝟏 𝟐 𝟏
𝒅𝒆𝒕𝑨 = 𝟐(−𝟏)𝟐 × [(𝟎)(𝟏) − (𝟏)(−𝟓)] + (−𝟏)(−𝟏)𝟑 × [(𝟑)(𝟏) − (𝟐)𝟖 − 𝟓)] + 𝟑(−𝟏)𝟒 × [(𝟑)(𝟏) − (𝟐)(𝟎)]
𝒅𝒆𝒕𝑨 = 𝟏𝟎 + 𝟏𝟑 + 𝟗 = 𝟑𝟐
𝒅𝒆𝒕𝑨 = 𝒂𝟏𝟐 (−𝟏)𝟏+𝟏 𝑴𝟏𝟐 + 𝒂𝟐𝟐 (−𝟏)𝟏+𝟐 𝑴𝟐𝟐 + 𝒂𝟑𝟐 (−𝟏)𝟏+𝟑 𝑴𝟑𝟐
𝟑 −𝟓 𝟐 𝟑 𝟐 𝟑
𝒅𝒆𝒕𝑨 = (−𝟏)(−𝟏)𝟑 [ ] + 𝟎(−𝟏)𝟒 [ ] + 𝟏(−𝟏)𝟑+𝟐 [ ]
𝟐 𝟏 𝟐 𝟏 𝟑 −𝟓
𝒅𝒆𝒕𝑨 = 𝟏𝟑 + 𝟎 + 𝟏𝟗 = 𝟑𝟐
2. Utilizando cofactores a lo largo del segundo renglón (fila) obtenga el determinante de la matriz:
12 −1 3
𝐴=[2 3 2]
0 −1 3
Métodos Numéricos
sfs
86
Det(A) 3(1) 21 1(3) 2*3 1(1) 22 12(3) (10)3 1123 12* 2 (10)
(1)
det(𝐴) = −1
3. Utilizando cofactores, calcule el determinante de la matriz:
0 1 1
𝐴 = [2 3 2]
0 −1 3
Nota: Como en la primera columna hay dos entradas iguales a cero, vamos a calcular el determinante
utilizando cofactores en la primera columna
1 0 3 0 0 1
1+1 4
det 𝐵 = 2(−1) × 𝑑𝑒𝑡 [0 1 2] + 1(−1) 𝑑𝑒𝑡 [1 0 3]
2 3 0 0 1 2
det(𝐵) = 2 × [−(3 × 1 × 2 + 2 × 3 × 1]
87
5 4 −3
𝐴=[ 2 𝑥 0]
−8 6 3
A det A 5x3 40 8 326 8x 3 605 324 105
−45
−9𝑥 = −45 → 𝑥 = →𝑥=5
−9
88
Si el det A, matriz de coeficientes, es diferente de cero, el sistema tiene única solución que se
obtiene como:
𝐷1 𝐷2 𝐷𝑛
𝑥1 = , 𝑥2 = , … 𝑥𝑛 =
𝐷 𝐷 𝐷
Dónde:
2x 5 y 4z 13
3x 2 y 3z 30
4x 6 y 5z 4
Procedimiento
89
2 −5 4
𝐷 = [3 2 3]
4 6 −5
b. Se repiten las dos primeras columnas:
2 −5 4 2 −5
𝐷 = [3 2 3 ]3 2
4 6 −5 4 6
D 151
c. Para calcular el determinante D1 se cambia las entradas de la primera columna en la matriz de
coeficientes por las entradas del término independiente de cada ecuación:
13 −5 4
𝐷1 = [ 30 2 3]
−4 6 −5
Recuerde: Para realizar el determinante se repiten las dos primeras columnas
D1 302
2 13 4
𝑫𝟐 = [3 30 3]
4 −4 −5
Recuerde: Para realizar el determinante se repiten las dos primeras columnas
Métodos Numéricos
sfs
90
𝑫𝟐 =230 5 1334 43 4 4304 432 133 5
D2 453
2 −5 13
𝐷 = [3 2 30 ]
4 6 −4
Recuerde: Para realizar el determinante se repiten las dos primeras columnas
x 2, y 3 z 6
g. Prueba:
91
2.3 TALLER DE ENTRENAMIENTO
1. Siguiendo el procedimiento presentado en el ejercicio de aprendizaje, resuelva el
siguiente sistema, dándole la respectiva prueba para verificar los resultados obtenidos:
2x y z 0
4x 3 y 2z 2
2x y 3z 0
Respuesta: x=-1/2, y=2, z=-1
2. Determine el valor de x, de tal manera que el determinante de la matriz sea igual a -1
−6 3 𝑥
𝐴 = [ 7 2 −5] 𝑅𝑒𝑠𝑝𝑢𝑒𝑠𝑡𝑎. 𝑥 = 4
5 9 −𝑥
3. Determine el valor de x, de tal manera que el determinante de la matriz sea igual a -141
𝑥 6 5
𝐴 = [4 𝑥2 2] 𝑅𝑒𝑠𝑝𝑢𝑒𝑠𝑡𝑎. 𝑥 = 3
6 𝑥 5
𝑥−3 9 3
𝐴=[ 7 8 𝑥 − 2] 𝑅𝑒𝑠𝑝𝑢𝑒𝑠𝑡𝑎. 𝑥 = 8
5 0 4
Métodos Numéricos
sfs
92
2.4 TEMA 2 REGRESIÓN POLINÓMICA
Regresión polinómica.
Según Asmar (1999), suponga que existe una relación funcional y = f (x) entre dos cantidades
x e y, con f desconocida y se conoce los distintos valores x e y.
Suponga que se conocen los datos x0, y0 , x1, y1 ...,xn, yn distintos y se desea encontrar un
polinomio:
El grado m del polinomio 𝑷𝒎 (x) puede ser el que cada quien deseé, teniendo en cuenta que
no puede superar a n 1 que es el número de parejas x e y .
Video: MinimosCuadrados
”Ver Video”: https://www.youtube.com/watch?v=TiLeQ7rDvSc
Métodos Numéricos
sfs
93
Video: REGRESION LINEAL
”Ver Video”: https://www.youtube.com/watch?v=rF66WNFYAoM
Consiste en encontrar una fórmula o ecuación lineal que muestre como se relacionan las variables
𝒚 = 𝒃𝟎 + 𝒃𝟏 𝒙
Donde:
𝒃𝟎 𝒚 𝒃𝟏 𝒔𝒐𝒏 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆𝒔
𝒃𝟏 : 𝑷𝒆𝒏𝒅𝒊𝒆𝒏𝒕𝒆 𝒅𝒆 𝒍𝒂 𝒍í𝒏𝒆𝒂 𝒓𝒆𝒄𝒕𝒂
𝒃𝟎 : 𝑰𝒏𝒕𝒆𝒓𝒔𝒆𝒄𝒄𝒊ó𝒏 𝒄𝒐𝒏 𝒆𝒍 𝒆𝒋𝒆 𝒚.
Para obtener los valores de b0 y b1 de esta ecuación lineal, y utilizando el método de los mínimos
cuadrados, se tiene que:
𝒆 = 𝒚 − 𝒃 𝟎 − 𝒃𝟏 𝒙
Calculando el error para cada valor, se tiene que:
e0 y0 b0 b1 x0 ,
e1 y1 b0 b1 x1 ,
e2 y2 b0 b1 x2 ,...,
en yn b0 b1 xn
Elevando al cuadrado cada lado de la ecuación y sumando, se tiene que:
Métodos Numéricos
sfs
94
𝑛 𝑛
∑(𝑒𝑖 )2 = ∑(𝑦𝑖 − 𝑏0 − 𝑏1 𝑥𝑖 )2
𝑖=1 𝑖=1
Lo que es lo mismo:
𝑺 = ∑(𝒚𝒊 − 𝒃𝟎 − 𝒃𝟏 𝒙𝒊 )𝟐
𝝏𝑺
= −𝟐 ∑(𝑦𝑖 − 𝑏0 − 𝑏1 𝑥𝑖 )
𝝏𝒃𝟎
𝝏𝑺
= −𝟐 ∑(𝑦𝑖 − 𝑏0 − 𝑏1 𝑥𝑖 )𝒙
𝝏𝒃𝟏
Igualando estas derivadas a cero, se encuentra el S mínimo:
→ ∑(𝑦𝑖 − ∑ 𝑏0 − ∑ 𝑏1 𝑥𝑖 ) = 𝟎
→ ∑ 𝑦𝑖 𝑥𝑖 − ∑ 𝑏0 𝑥𝑖 − ∑ 𝑏1 𝑥𝑖 𝟐 = 𝟎
∑ 𝑦 = 𝑏0 𝑛 + 𝑏1 ∑ 𝑥
∧
Métodos Numéricos
sfs
95
∑ 𝑥𝑦 = 𝑏0 ∑ 𝑥 + 𝑏1 ∑ 𝑥 2
Cuya solución es:
𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2
𝑏0 = 𝑦̅ − 𝑏1 𝑥̅
EJERCICIOS DE APRENDIZAJE:
1. Este es parte de un ejemplo planteado por Asmar (1999). Para la tabla:
x y
0 -1
2 0
3 2
5 1
x y xy
x2
0 -1 0 0
2 0 4 0
3 2 9 6
5 1 25 5
x 10 y 2 x 2
38 x y 11
Con n = 4, Ya que hay cuatro parejas x, y
Cálculo de 𝑥̅ :
Métodos Numéricos
sfs
96
∑ 𝑥 10 5
𝑥̅ = = =
𝑛 4 2
Cálculo de 𝑦̅:
∑𝑦 2 1
𝑦̅ = = =
𝑛 4 2
Cálculo de 𝑏1 :
𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2
𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2
𝟒 ∗ 𝟏𝟏 − 𝟏𝟎 ∗ 𝟐 𝟐𝟒 𝟔
𝒃𝟏 = = =
𝟒 ∗ 𝟑𝟖 − (𝟏𝟎)𝟐 𝟓𝟐 𝟏𝟑
Cálculo de:
̅ − 𝒃𝟏 𝒙
𝒃𝟎 = 𝒚 ̅
𝟏 𝟔 𝟓 𝟏 𝟑𝟎 𝟏𝟑 − 𝟑𝟎
𝒃𝟎 = − × = − =
𝟐 𝟏𝟑 𝟐 𝟐 𝟐𝟔 𝟐𝟔
𝟏𝟕
𝒃𝟎 = − 𝟐𝟔
La recta regresión queda entonces:
17 6
𝑦= + 𝑥
26 13
Métodos Numéricos
sfs
97
2. Para los datos de la siguiente tabla determine la recta de regresión
x y
-5 15
-2 9
1 3
4 -3
6 -7
PROCEDIMIENTO:
X Y X2 X.Y
-5 15 25 -75
-2 9 4 -18
1 3 1 3
4 -3 16 -12
6 -7 36 -42
x 4 y 17 x 2
82 x y 144
Con n = 5, Ya que hay cinco parejas x, y
Cálculo de 𝑥̅ :
∑𝑥 4
𝑥̅ = =
𝑛 5
Cálculo de 𝑦̅:
∑ 𝑦 17
𝑦̅ = =
𝑛 5
Cálculo de 𝑏1 :
𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2
Métodos Numéricos
sfs
98
𝟓 ∗ (−𝟏𝟒𝟒) − 𝟒 ∗ 𝟏𝟕 𝟕𝟖𝟖
𝒃𝟏 = = − = −𝟐
𝟓 ∗ 𝟖𝟐 − (𝟒)𝟐 𝟑𝟗𝟒
Cálculo de:
̅ − 𝒃𝟏 𝒙
𝒃𝟎 = 𝒚 ̅
𝟏𝟕 𝟒 𝟏𝟕 𝟖 𝟐𝟓
𝒃𝟎 = − (−𝟐) × = + = =𝟓
𝟓 𝟓 𝟓 𝟓 𝟓
𝒃𝟎 = 𝟓
La recta regresión queda entonces:
𝑦 = 5 − 2𝑥
REGRESIÓN CUADRÁTICA O REGRESIÓN PARABÓLICA
𝒚 𝒃𝟎 𝒃𝟏 x 𝒃𝟐 x 2
Para encontrar los valores de 𝒃𝟎, 𝒃𝟏 , 𝒃𝟐 el método de los mínimos cuadrados plantealas siguientes
ecuaciones:
2
y 𝒃𝟎 n 𝒃𝟏 x 𝒃𝟐 x
2 3
xy 𝒃𝟎 x 𝒃𝟏 x 𝒃𝟐 x
0 1 2
2 2 3 4
x y 𝒃𝟎 x 𝒃𝟏 x 𝒃𝟐 x
Al solucionar este sistema 3X3 se obtienen las constantes b0 , b1 b2
Métodos Numéricos
sfs
99
2.4.2 EJERCICIOS DE APRENDIZAJE
1. Para los datos de la tabla:
x y
-5 113
-2 23
1 5
4 59
6 135
10 383
X Y X2 X.Y x2.y x3 x4
-5 113 25 -565 2825 -125 625
-2 23 4 -46 92 -8 16
1 5 1 5 5 1 1
4 59 16 236 944 64 256
6 135 36 810 4860 216 1296
10 383 100 3830 38300 1000 10000
2 2 3 4
x 14 y 718 x 182 x y 4270 x y 47026 x 1148 x 12194
Con n6
Se debe plantear y resolver el siguiente sistema:
6 14 182
𝐷 = [ 14 182 1148 ] = 2.840.040
182 1148 −12194
Métodos Numéricos
sfs
100
718 14 182
𝐷1 = [ 4270 182 1148 ] = 8.520.120
47026 1148 −12194
6 718 182
𝐷2 = [ 14 4270 1148 ] = −5.680.080
182 47026 −12194
6 14 718
𝐷3 = [ 14 182 4270 ] = 11.360.160
182 1148 47026
Recuerde que:
𝐷
𝑏𝑛 =
𝐷𝑛
Reemplazando los valores obtenidos, se tiene que:
𝐷1 8.520.120
𝑏0 = = =3
𝐷 2.840.040
𝐷2 −5.680.080
𝑏1 = = = −2
𝐷 2.840.040
𝐷3 11.360.160
𝑏2 = = =4
𝐷 2.840.040
La parábola de regresión está dada por:
𝒚 = 𝟑 − 𝟐𝒙 − 𝟒𝒙𝟐
101
Procedimiento:
X Y X2 X.Y x2y x3 x4
-4 10 16 -40 160 -64 256
-1 -2 1 2 -2 -1 1
3 10 9 30 90 27 81
5 28 25 140 700 125 625
8 70 64 560 4480 512 4096
y 11 6 x x y 69 2 x y 542 8 x 599 x
2 3 4
115 5059
2
x 11
5 11 115
𝐷 = [ 11 115 599 ] = 497.376
115 599 5059
116 11 115
𝐷1 = [ 692 115 599 ] = −994.752
5428 599 5059
5 116 115
𝐷2 = [ 11 692 599 ] = 497.376
115 5428 5059
5 11 116
𝐷3 = [ 11 115 692 ] = 497.376
115 599 5248
Métodos Numéricos
sfs
102
Recuerde que:
𝐷
𝑏𝑛 =
𝐷𝑛
Reemplazando los valores obtenidos, se tiene que:
𝐷1 −994.752
𝑏0 = = = −2
𝐷 497.376
𝐷2 497.376
𝑏1 = = =1
𝐷 497.376
𝐷3 497.376
𝑏2 = = =1
𝐷 497.376
La parábola de regresión es:
y 2 x x2
Polinomio de regresión de grado 3 o polinomio cúbico de regresión
El polinomio de grado 3 o polinomio cúbico que se debe obtener tiene la siguiente forma:
𝒚 𝒃𝟎 𝒃𝟏 x 𝒃𝟐 x2+𝒃𝟑 x3
Para encontrar las constantes b0 , b1 , b2 b3 el método de los mínimos cuadrados plantea las
siguientes ecuaciones:
2 3
y 𝒃𝟎 n 𝒃𝟏 x 𝒃𝟐 x 𝒃𝟑 x
0 1 2 3
2 3 4
xy 𝒃𝟎 x 𝒃𝟏 x 𝒃𝟐 x 𝒃𝟑 x
0 1 2 3
2 2 3 4 5
x y 𝒃𝟎 x 𝒃𝟏 x 𝒃𝟐 x 𝒃𝟑 x
0 1 2 3
3 3 4 5 6
x y 𝒃𝟎 x 𝒃𝟏 x 𝒃𝟐 x 𝒃𝟑 x
0 1 2 3
Métodos Numéricos
sfs
103
Al solucionar este sistema 4X4 se obtienen las constantes b0 , b1 , b 2 b3 .
Polinomio de regresión de grado n.
El polinomio de grado 3 o polinomio cúbico que se debe obtener tiene la siguiente forma:
𝒚 𝒃𝟎 𝒃𝟏 x 𝒃𝟐 x2+𝒃𝟑 x3+ ⋯ + 𝒃𝒏 xn
Para encontrar las constantes b0 , b1, b2 , b3 .... bn el método de los mínimos cuadrados plantea
las siguientes ecuaciones:
3 n
y b0 n b1 x b2 x b3 x ... bn x
0 1 2 3 n
2 3 4 n1
xy b0 x b1 x b2 x b3 x ... b x
0 1 2 3 n
2 2 3 4 5 n2
x y b0 x b1 x b2 x b3 x ... bn x
0 1 2 3 n
3 3 5 6 n3
x y b0 x b1 x b x b3 x ... bn x
0 1 2 3 n
En forma general:
x -2 3 4 9
y 7 4 -2 -3
Métodos Numéricos
sfs
104
2. Encuentre la parábola de regresión para el siguiente conjunto de datos.
x -3 -1 1 2 4 5
y 9 2 1 3 7 9
x -9 -6 -4 -1 1 2 3 5 7 10
y -921 -303 -101 7 19 25 39 115 295 865
Este tema consiste en dado un conjunto de valores x e y encontrar una función polinómica
que represente estos valores, que pase exactamente por cada uno de ellos y que permita
estimar valores no observados (es decir, que permita hacer predicciones o interpolaciones o
aproximaciones).
También permite aproximar mediante una función polinómica funciones que no son polinómicas. Se
utiliza la aproximación polinomial ya que los polinomios son continuos en todo su dominio, con esto
se garantizan funciones continuas, que son las que interesa estudiar en el análisis numérico y además
la derivada y la integral de polinomios son muy fáciles de determinar y también son polinómicas.
Video: Interpolación 1b
”Ver Video”: https://www.youtube.com/watch?v=S3QEyUN3sn0
Métodos Numéricos
sfs
105
Polinomio interpolante de LaGrange
Video: Interpolación 2
”Ver Video”: https://www.youtube.com/watch?v=h5s6lZ9EDbk
Video: Calculo y métodos Numéricos. Interpolación de Lagrange
”Ver Video”: https://www.youtube.com/watch?v=S4KdYhTzaxE
La idea es encontrar un polinomio de primer grado que pasa por dos puntos distintos:
x0 , y0 x1 , y1
O lo que es lo mismo aproximar una función f, para la cual:
f (x0 ) y0 y f (x1 ) y1
Por medio de un polinomio de primer grado que interpole los valores de f en los puntos dados o que
coincida con ellos.
P1 (x) a0 a1 x
Se tiene que:
𝑃1 (𝑥0 ) = 𝑎0 + 𝑎1 𝑥0 = 𝑦0
𝑃1 (𝑥1 ) = 𝑎0 + 𝑎1 𝑥1 = 𝑦1
Se tiene un sistema 2X2, donde las incógnitas son a0 a1
Solución del sistema:
1 𝑥0
𝐷=[ ] = 𝑥1 − 𝑥0
1 𝑥1
𝑦0 𝑥0
𝐷1 = [𝑦 𝑥1 ] = 𝑥1 𝑦0 − 𝑥0 𝑦1
1
Métodos Numéricos
sfs
106
1 𝑦0
𝐷2 = [ ] = 𝑦1 − 𝑦0
1 𝑦1
𝐷1 𝑥1 𝑦0 −𝑥0 𝑦1 𝐷 𝑦 −𝑦
𝑎0 = = , 𝑎1 = 𝐷2 = 𝑥1−𝑥0
𝐷 𝑥1 −𝑥0 1 0
𝑥1 𝑦0 −𝑥0 𝑦1 𝑦1 −𝑦0
P1 (x) +( )𝑥
𝑥1 −𝑥0 𝑥1 −𝑥0
Factorizando:
𝒙−𝒙𝟏 𝒙−𝒙𝟎
P1 (x)= ( ) 𝒚𝟎 + ( ) 𝒚𝟏
𝒙𝟎 −𝒙𝟏 𝒙𝟏 −𝒙𝟎
𝒙−𝒙𝟏 𝒙−𝒙𝟎
L0 (x)= ( ) ∧ 𝑳𝟏 (𝒙) ( )
𝒙𝟎 −𝒙𝟏 𝒙𝟏 −𝒙𝟎
Y se define:
Se puede afirmar que P(x) es la única función lineal que pasa por los puntos x0 , y0 x1 , y1 .
107
Procedimiento
𝑥 8 𝟒 𝒙
L0 (x)=− + → L0 (x)= −
6 6 𝟑 𝟔
Cálculo de:
𝒙 − 𝒙𝟎 𝒙−𝟐 𝑥 2
𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = − →
𝒙𝟏 − 𝒙𝟎 𝟖−𝟐 6 6
𝒙 𝟏
L1 (x)=𝟔 − 𝟑
Reemplazando en: P1 (x) L0 (x) f (x0 ) L1 (x) f (x1 )
𝒙 𝟏𝟕
𝟒 𝒙 𝒙 𝟏
P1 (x) ( − ) 𝟓 + ( − ) 𝟑 →P1 (x)
𝟑 𝟔 𝟔 𝟑
− +
𝟑 𝟑
2. Encuentre el polinomio lineal que pasa por los puntos (1,7) (6,27)
3. Procedimiento
108
Cálculo de:
𝒙 − 𝒙𝟎 𝒙−𝟏 𝒙−𝟏
𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = ( )
𝒙𝟏 − 𝒙𝟎 𝟔−𝟏 𝟓
𝒙 𝟏
L1 (x)=𝟔 − 𝟑
Reemplazando en: P1 (x) L0 (x) f (x0 ) L1 (x) f (x1 )
𝒙 𝟏𝟕
P1 (x) (
−𝒙+𝟔
𝟓
𝒙−𝟏
) 𝟕 + ( 𝟓 ) 𝟐𝟕 →P1(x) − +
𝟑 𝟑
−7𝑥+42 27𝑥−27
P1 (x)
5
+
5
→P1 (x)=−7𝑥+42+27𝑥−27
5
=
20𝑥+15
5
20𝑥 15
P1 (x)= + = 4𝑥 + 3
5 5
Los puntos por los cuales debe pasar el polinomio son: x0 , y0 , x1 , y1 , x2 , y2
Métodos Numéricos
sfs
109
Recuerde que x0 , x1 , x2 deben ser diferentes.
Se desea encontrar el polinomio interpolante de grado menor o igual que dos, el polinomio se simboliza
como P2 (x) y tiene la forma:
𝑃2 (𝑥 ) = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2
Se debe encontrar el valor de las constantes: a0 , a1 a2
Deducción del polinomio de LaGrange de grado 2 (Esta demostración puede ser omitida)
Como los puntos x0 , y0 , x1 , y1 , x2 , y2 pertenecen al polinomio P2 (x) , se debe cumplir que:
𝒚𝟎 = 𝒂𝟎 + 𝒂𝟏 𝒙𝟎 + 𝒂𝟐 𝒙𝟎 𝟐
𝒚𝟏 = 𝒂𝟎 + 𝒂𝟏 𝒙𝟏 + 𝒂𝟐 𝒙𝟏 𝟐
𝒚𝟐 = 𝒂𝟎 + 𝒂𝟏 𝒙𝟐 + 𝒂𝟐 𝒙𝟐 𝟐
1 𝑥0 𝑥0 2 𝑎0 𝑦0
𝐴 = [1 𝑥1 𝑥1 2 ] , 𝑥 = [𝑎1 ] ∧ 𝐵 = [𝑦1 ]
1 𝑥2 𝑥2 2 𝑎2 𝑦2
Sea D el det(A)
Para simplificar el cálculo de del determinante D se aplican operaciones elementales en la matriz A:
Métodos Numéricos
sfs
110
R1 R1 , R2 R2 R1 (1) R3 R3 R1 (1)
1 𝑥0 𝑥0 2
[0 𝑥1 − 𝑥0 𝑥1 2 −𝑥0 2 ]
0 𝑥2 − 𝑥0 𝑥2 2 −𝑥0 2
Se obtiene D por expansión por cofactores a lo largo de la primera columna:
𝑥1 − 𝑥0 𝑥1 2 −𝑥0 2 2 2
𝐷=[ 2 2 ]2 = (𝑥1 − 𝑥0 ) (𝑥2 −𝑥0 ) − (𝑥2 −
𝑥2 − 𝑥0 𝑥2 −𝑥0
𝑥0 )(𝑥1 2 −𝑥0 2 )
𝑦0 𝑥0 𝑥0 2
𝑫𝟏 = [𝑦1 𝑥1 𝑥1 2 ]
𝑦2 𝑥2 𝑥2 2
Se obtiene este determinante por expansión a lo largo de la primera columna:
1+1 𝑥1 𝑥1 2 2+1 𝑥0 𝑥0
2
𝑫𝟏 = 𝑦0 (−1) [ ] + 𝑦1 (−1) [ ]
𝑥2 𝑥2 2 𝑥2 𝑥2 2
2
3+1 𝑥0 𝑥0
+ 𝑦2 (−1) [ ]
𝑥1 𝑥1 2
𝑫𝟏 = 𝑦0 (𝑥1 𝑥2 2 − 𝑥1 2 𝑥2 ) − 𝑦1 (𝑥0 𝑥2 2 − 𝑥0 2 𝑥2 ) + 𝑦2 (𝑥0 𝑥1 2 − 𝑥0 2 𝑥1 )
111
𝑦0 𝑥0 2
𝟏
𝑫𝟐 = [𝟏 𝑦1 𝑥1 2 ]
𝟏 𝑦2 𝑥2 2
Se soluciona este determinante por expansión a lo largo de la segunda columna:
1+2 1 𝑥1 2 2+2 1 𝑥0
2
3+2 1 𝑥0
2
𝑫𝟐 = 𝑦0 (−1) [ ] + 𝑦1 (−1) [ ] + 𝑦2 (−1) [ ]
1 𝑥2 2 1 𝑥2 2 1 𝑥1 2
𝑫𝟐 = −𝑦0 (𝑥2 2 − 𝑥1 2 ) + 𝑦1 (𝑥2 2 − 𝑥0 2 ) − 𝑦2 (𝑥1 2 − 𝑥0 2 )
𝟏 𝑥0 𝑦0
𝑫𝟑 = [𝟏 𝑥1 𝑦1 ]
𝟏 𝑥2 𝑦2
Se soluciona este determinante por expansión a lo largo de la segunda columna:
Se tiene que:
𝐷1 𝐷2 𝐷3
𝑎0 = ; 𝑎1 = ; 𝑎2 =
𝐷 𝐷 𝐷
Reemplazando estas expresiones en el polinomio 𝑷𝟐 (𝒙) = 𝒂𝟎 + 𝒂𝟏 𝒙 + 𝒂𝟐 𝒙𝟐 queda:
Métodos Numéricos
sfs
112
𝑦 𝑥1 𝑥2 (𝑥2 − 𝑥1 ) − 𝑦1 𝑥0 𝑥2 (𝑥2 − 𝑥0 ) + 𝑦2 𝑥0 𝑥1 (𝑥1 − 𝑥0 )
𝑷𝟐 (𝒙) = 0 +
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )
Factorizando, se obtiene:
Se tiene que:
113
−𝑦1 (𝑥2 − 𝑥0 )[𝑥(𝑥 − 𝑥0 ) − 𝑥2 (𝑥 − 𝑥0 )] = −𝑦1 (𝑥2 − 𝑥0 )(𝑥 − 𝑥0 )(𝑥 − 𝑥2 )
Simplificando, se tiene:
𝒚𝟎 (𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝑷𝟐 (𝒙) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )
𝒚𝟏 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
+
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
𝒚𝟐 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
+
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )
Métodos Numéricos
sfs
114
Haciendo el mismo análisis que para el polinomio interpolante lineal, se tiene que:
(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝑳𝟎 (𝒙) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝑳𝟏 (𝒙) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝑳𝟐 (𝒙) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )
1 2
Y el polinomio interpolante de grado dos sería:
Cálculo de:
(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝐿0 (𝑥 ) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )
Métodos Numéricos
sfs
115
(𝒙 − 𝟓)(𝒙 − 𝟕) 𝒙𝟐 − 𝟕𝒙 − 𝟓𝒙 + 𝟑𝟓
𝐿0 (𝑥 ) = = =
(−𝟐 − 𝟓)(−𝟐 − 𝟕) (−𝟕)(−𝟗)
𝒙𝟐 − 𝟏𝟐𝒙 + 𝟑𝟓
𝐿0 (𝑥 ) =
𝟔𝟑
Cálculo de:
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝐿1 (𝑥 ) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
(𝒙 − (−𝟐))(𝒙 − 𝟕) (𝒙 + 𝟐))(𝒙 − 𝟕)
𝐿1 (𝑥 ) = = =
(𝟓 − (−𝟐))(𝟓 − 𝟕) (𝟓 + 𝟐)(𝟓 − 𝟕)
𝒙𝟐 − 𝟓𝒙 − 𝟏𝟒
𝐿1 (𝑥 ) =
−𝟏𝟒
Cálculo de:
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝐿2 (𝑥 ) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )
𝒙𝟐 − 𝟑𝒙 − 𝟏𝟎
𝐿2 (𝑥 ) =
𝟏𝟖
Cálculo de 𝑃2 (𝑥): Recuerde que:
116
Simplificando:
𝟏 𝟐 𝟏 −𝟏𝟐 𝟏𝟎 𝟑 𝟑𝟓 𝟐𝟖 𝟏𝟎
P2 (x)= ( − + ) 𝒙𝟐 − ( + − )𝒙 + ( + − )
𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗
𝟖 𝟏𝟏 𝟒𝟏
P2 (x)=− 𝒙𝟐 + 𝒙+
𝟔𝟑 𝟐𝟏 𝟗
𝒙𝟐 − 𝟐𝟐𝒙 + 𝟏𝟐𝟎
𝐿0 (𝑥 ) =
𝟒𝟖
Métodos Numéricos
sfs
117
Cálculo de:
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝐿1 (𝑥) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
(𝒙 − (−𝟐))(𝒙 − 𝟕) (𝒙 + 𝟐))(𝒙 − 𝟕)
𝐿1 (𝑥 ) = = =
(𝟓 − (−𝟐))(𝟓 − 𝟕) (𝟓 + 𝟐)(𝟓 − 𝟕)
𝒙𝟐 − 𝟓𝒙 − 𝟏𝟒
𝐿1 (𝑥 ) =
−𝟏𝟒
Cálculo de:
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝐿2 (𝑥 ) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )
𝒙𝟐 − 𝟑𝒙 − 𝟏𝟎
𝐿2 (𝑥 ) =
𝟏𝟖
Cálculo de 𝑃2 (𝑥): Recuerde que:
Simplificando:
118
𝟏 𝟐 𝟏 −𝟏𝟐 𝟏𝟎 𝟑 𝟑𝟓 𝟐𝟖 𝟏𝟎
P2 (x)= ( − + ) 𝒙𝟐 − ( + − )𝒙 + ( + − )
𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗
𝟖 𝟏𝟏 𝟒𝟏
P2 (x)=− 𝒙𝟐 + 𝒙+
𝟔𝟑 𝟐𝟏 𝟗
1. Encuentre el polinomio de grado 2 o parábola de interpolación que pase por los puntos:
𝒙𝟐 − 𝟐𝟐𝒙 + 𝟏𝟐𝟎
𝐿0 (𝑥 ) =
𝟒𝟖
Cálculo de:
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝐿1 (𝑥 ) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
Métodos Numéricos
sfs
119
(𝒙 − 𝟒)(𝒙 − 𝟏𝟐) 𝒙𝟐 − 𝟏𝟔𝒙 + 𝟒𝟖
𝐿1 (𝑥 ) = = =
(𝟏𝟎 − 𝟒)(𝟏𝟎 − 𝟏𝟐) −𝟏𝟐
𝒙𝟐 − 𝟏𝟔𝒙 + 𝟒𝟖
𝐿1 (𝑥 ) =
−𝟏𝟐
Cálculo de:
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝐿2 (𝑥 ) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )
𝒙𝟐 − 𝟏𝟒 + 𝟒𝟎
𝐿2 (𝑥 ) =
𝟏𝟔
Cálculo: P2 (x) y0 L0 (x) y1L1 (x) y2 L2 (x)
P2 (x) = 𝟓𝒙𝟐 − 𝒙 + 𝟗
120
CASO 3: Polinomio de grado máximo tres que pasa por cuatro puntos, donde los valores de x son
todos diferentes.
(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )
𝑳𝟎 (𝒙) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )(𝒙𝟎 − 𝒙𝟑 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )
𝑳𝟏 (𝒙) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )(𝒙𝟏 − 𝒙𝟑 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟑 )
𝑳𝟐 (𝒙) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )(𝒙𝟐 − 𝒙𝟑 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝑳𝟑 (𝒙) =
(𝒙𝟑 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )(𝒙𝟑 − 𝒙𝟐 )
Caso 4: Polinomio de grado máximo cuatro que pasa por cinco puntos con la condición que todos
los valores de x son diferentes.
121
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )(𝒙 − 𝒙𝟒 )
𝑳𝟏 (𝒙) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )(𝒙𝟏 − 𝒙𝟑 )(𝒙𝟏 − 𝒙𝟒 )
𝒇(𝟐, 𝟓) 𝒔𝒊:
𝒇(𝟐) = 𝟎, 𝟓𝟏𝟎𝟑𝟕𝟓
𝒇(𝟐, 𝟐) = 𝟎, 𝟓𝟐𝟎𝟕𝟖𝟒
𝒇(𝟐, 𝟒) = 𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
𝒇(𝟐, 𝟔) = 𝟎, 𝟒𝟖𝟏𝟑𝟑𝟎
𝒇(𝟐, 𝟖) = 𝟎, 𝟒𝟑𝟓𝟗𝟏𝟔
Métodos Numéricos
sfs
122
Procedimiento:
(𝟐, 𝟓 − 𝟐, 𝟔) (𝟐, 𝟓 − 𝟐, 𝟒)
𝑷𝟏 (𝟐, 𝟓) = 𝟎, 𝟓𝟏𝟒𝟏𝟒𝟕 + 𝟎, 𝟒𝟖𝟏𝟑𝟑𝟎𝟔
(𝟐, 𝟒 − 𝟐, 𝟔) (𝟐, 𝟔 − 𝟐, 𝟒)
(𝟐, 𝟓 − 𝟐, 𝟒)(𝟐, 𝟓 − 𝟐, 𝟔)
𝑷𝟐 (𝟐, 𝟓) = 𝟎, 𝟓𝟐𝟎𝟕𝟖𝟒
(𝟐, 𝟐 − 𝟐, 𝟒)(𝟐, 𝟐 − 𝟐, 𝟔)
(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟔)
+𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
(𝟐, 𝟒 − 𝟐, 𝟐)(𝟐, 𝟒 − 𝟐, 𝟔)
(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟒)
+𝟎, 𝟒𝟖𝟏𝟑𝟑𝟎𝟔
(𝟐, 𝟔 − 𝟐, 𝟐)(𝟐, 𝟔 − 𝟐, 𝟒)
=
𝟎, 𝟒𝟗𝟖𝟐𝟏𝟐𝟎 ≈ 𝒇(𝟐, 𝟓)
El otro polinomio interpolante de grado dos se obtiene tomando los nodos:
Métodos Numéricos
sfs
123
𝒙𝟎 = 𝟐, 𝟒, 𝒙𝟏 = 𝟐, 𝟔 ∧ 𝒙𝟐 = 𝟐, 𝟖
(𝟐, 𝟓 − 𝟐, 𝟔)(𝟐, 𝟓 − 𝟐, 𝟖)
𝑷𝟐 (𝟐, 𝟓) = 𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
(𝟐, 𝟐 − 𝟐, 𝟒)(𝟐, 𝟐 − 𝟐, 𝟔)
(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟔)
+𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
(𝟐, 𝟒 − 𝟐, 𝟐)(𝟐, 𝟒 − 𝟐, 𝟔)
(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟒)
+𝟎, 𝟒𝟑𝟓𝟗𝟏𝟔𝟎
(𝟐, 𝟔 − 𝟐, 𝟐)(𝟐, 𝟔 − 𝟐, 𝟒)
==
𝟎, 𝟒𝟗𝟕𝟗𝟏𝟒𝟎 ≈ 𝒇(𝟐, 𝟓)
𝒙𝟎 = 𝟐, 𝟐 𝒙𝟏 = 𝟐, 𝟒 𝒙𝟐 = 𝟐, 𝟔 ∧ 𝒙𝟑 = 𝟐, 𝟖
Ya que 𝟐, 𝟓 ∈ [𝟐. 𝟐, 𝟐. 𝟖] y 𝟐, 𝟐; 𝟐, 𝟒; 𝟐, 𝟔 ∧ 𝟐, 𝟖 son los nodos más cercanos a 2,5
Nota: Resuelve el problema teniendo como modelo los ejercicios realizados, debes tener
presente que es para el grado tres el polinomio interpolante.
124
Dados n + 1 puntos x0 , y0 , x1 , y1 ,...,xn , yn con 𝒙𝟎 , 𝒙𝟏 , . . , 𝒙𝒏 números distintos.
Para obtener los coeficientes de los polinomios de interpolación utilizando el método de Newton,
es necesario efectuar cálculos con resultados anteriores, para facilitar esto, se utilizarán las
diferencias divididas:
Para dos puntos (𝒙𝟎 , 𝒇(𝒙𝟎 )) ∧ (𝒙𝟏 , 𝒇(𝒙𝟏 )) se tiene que las primeras y únicas diferencias
divididas son:
𝒇(𝒙𝟏 ) − 𝒇(𝒙)𝟎
𝒇[𝒙𝟏 , 𝒙𝟎 ] =
𝒙𝟏 − 𝒙𝟎
Para tres puntos (𝒙𝟎 , 𝒇(𝒙𝟎 )), (𝒙𝟏 , 𝒇(𝒙𝟏 )) ∧ (𝒙𝟐 , 𝒇(𝒙𝟐 )) se tiene que las primeras
diferencias divididas son:
𝒇(𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟏 , 𝒙𝟎 )
𝒇[𝒙𝟐 , 𝒙𝟏, 𝒙𝟎 ] =
𝒙𝟏 − 𝒙𝟎
Para cuatro puntos (𝒙𝟎 , 𝒇(𝒙𝟎 )), (𝒙𝟏 , 𝒇(𝒙𝟏 )), (𝒙𝟐 , 𝒇(𝒙𝟐 )) ∧ (𝒙𝟑 , 𝒇(𝒙𝟑 )) se tiene que las
primeras diferencias
divididas son:
𝒇(𝒙𝟑 ) − 𝒇(𝒙𝟐 )
𝒇[ 𝒙 𝟑 , 𝒙 𝟐 ] = ,
𝒙𝟑 − 𝒙𝟐
𝒇(𝒙𝟐 ) − 𝒇(𝒙𝟏 )
𝒇[ 𝒙 𝟐 , 𝒙 𝟏 ] = ∧
𝒙𝟐 − 𝒙𝟏
Métodos Numéricos
sfs
125
𝒇(𝒙𝟏 ) − 𝒇(𝒙𝟎 )
[ 𝒙𝟏 , 𝒙𝟎 ] =
𝒙𝟏 − 𝒙𝟎
Las segundas diferencias divididas son:
𝒇(𝒙𝟑 , 𝒙𝟐 ) − 𝒇(𝒙𝟐 , 𝒙𝟏 )
𝒇[𝒙𝟑 , 𝒙𝟐, 𝒙𝟏 ] =
𝒙𝟑 − 𝒙𝟏
∧
𝒇(𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟏 , 𝒙𝟎 )
𝒇[𝒙𝟐 , 𝒙𝟏, 𝒙𝟎 ] =
𝒙𝟐 − 𝒙𝟎
Las Terceras diferencias divididas son:
𝒇(𝒙𝟑 , 𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 )
𝒇[𝒙𝟑 , 𝒙𝟐, 𝒙𝟏 , 𝒙𝟎 ] =
𝒙𝟑 − 𝒙𝟎
Y así sucesivamente.
Procedimiento:
𝒙𝟎 = 𝟏, 𝒇(𝒙𝟎 ) = 𝟎
𝒙𝟏 = 𝟒, 𝒇(𝒙𝟏 ) = 𝟏, 𝟑𝟖𝟔𝟑
𝒙𝟐 = 𝟖, 𝒇(𝒙𝟐 ) = 𝟏, 𝟕𝟗𝟏𝟖
126
𝟐−𝟕 𝟓
𝒇[𝒙𝟐 , 𝒙𝟏 ] = =−
𝟖−𝟒 𝟒
𝟕−𝟑 𝟒
[𝒙𝟏 , 𝒙𝟎 ] = =
𝟒−𝟏 𝟑
Las segundas diferencias divididas son:
𝟓 𝟒 𝟑𝟏
𝒇(𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟏 , 𝒙𝟎 ) − 𝟒 − 𝟑 − 𝟏𝟐 𝟑𝟏
𝒇[𝒙𝟐 , 𝒙𝟏, 𝒙𝟎 ] = = = =−
𝒙𝟐 − 𝒙𝟎 𝟖−𝟏 𝟕 𝟖𝟒
(𝒙𝟎 , 𝒇(𝒙𝟎 )) ∧ (𝒙𝟏 , 𝒇(𝒙𝟏 )), (𝒙𝟐 , 𝒇(𝒙𝟐 )), … , (𝒙𝒏 , 𝒇(𝒙𝒏 ))
Con:
𝒙𝟎 ≠ 𝒙𝟏 ≠ 𝒙𝟐 ≠ ⋯ ≠ 𝒙𝒏
El polinomio que pasa exactamente por todos los puntos es:
(𝒙𝟎 , 𝒇(𝒙𝟎 )) ∧ (𝒙𝟏 , 𝒇(𝒙𝟏 )), (𝒙𝟐 , 𝒇(𝒙𝟐 )), … , (𝒙𝒏 , 𝒇(𝒙𝒏 ))
Con:
𝒙𝟎 ≠ 𝒙𝟏 ≠ 𝒙𝟐 ≠ ⋯ ≠ 𝒙𝒏
El polinomio de Newton que pasa exactamente por todos los puntos es:
Con:
𝒃𝒐 = 𝒇(𝒙𝟎 )
𝒃 𝟏 = 𝒇[ 𝒙 𝟏 , 𝒙 𝟎 ]
Métodos Numéricos
sfs
127
𝒃 𝟐 = 𝒇[ 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝟑 = 𝒇[ 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝟒 = 𝒇[ 𝒙 𝟒 , 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝟓 = 𝒇[ 𝒙 𝟓 , 𝒙 𝟒 , 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝒏 = 𝒇[ 𝒙 𝒏 … , 𝒙 𝟓 , 𝒙 𝟒 , 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
EJERCICIO DE APRENDIZAJE:
1. Para 3, 40, 1, 8 5, 40, encuentre el polinomio de Newton de grado máximo que
pasa por todos los puntos.
Procedimiento
Datos:
𝒙𝟎 = −𝟑 ∧ 𝒇(𝒙𝟎 ) = 𝟒𝟎
𝒙𝟏 = 𝟏 ∧ 𝒇(𝒙𝟏 ) = −𝟖
𝒙𝟐 = 𝟓 ∧ 𝒇(𝒙𝟐 ) = 𝟒𝟎
El polinomio de Newton de grado máximo que pasa por todos los puntos es:
P2 (x) b0 b1 x x0 b2 x x0 x x1
Con:
𝒃𝟎 = 𝒇(𝒙𝟎 )
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ]
𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ]
Determinación de 𝒃𝟎 = 𝒇(𝒙𝟎 ):
𝒃𝟎 = 𝒇(𝒙𝟎 ) = 𝟒𝟎
Métodos Numéricos
sfs
128
Cálculo de todas las diferencias divididas:
−𝟖 − 𝟒𝟎 −𝟒𝟖
𝒇[𝒙𝟏 , 𝒙𝟎 ] = = = −𝟏𝟐
𝟓−𝟏 𝟒
Las segundas diferencias divididas son:
Determinación de:
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ] = −𝟏𝟐
Determinación de:
𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ] = 𝟑
Determinación de:
P2 (x) b0 b1 x x0 b2 x x0 x x1
𝑷𝟐 (𝒙) = 𝟑𝒙𝟐 − 𝟔𝒙 − 𝟓
Métodos Numéricos
sfs
129
2.5.6 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO
2. Para 10, 4, 3, 8 5, 10, encuentre el polinomio de Newton de grado máximo que
pasa por todos los puntos.
Procedimiento
Datos:
𝒙𝟎 = −𝟏𝟎 ∧ 𝒇(𝒙𝟎 ) = 𝟒
𝒙𝟏 = 𝟑 ∧ 𝒇(𝒙𝟏 ) = 𝟖
𝒙𝟐 = 𝟓 ∧ 𝒇(𝒙𝟐 ) = 𝟏𝟎
El polinomio de Newton de grado máximo que pasa por todos los puntos es:
P2 (x) b0 b1 x x0 b2 x x0 x x1
Con:
𝒃𝟎 = 𝒇(𝒙𝟎 )
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ]
𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ]
Determinación de 𝒃𝟎 = 𝒇(𝒙𝟎 ):
𝒃𝟎 = 𝒇(𝒙𝟎 ) = 𝟒
Cálculo de todas las diferencias divididas:
𝟖−𝟒 𝟒
𝒇[𝒙𝟏 , 𝒙𝟎 ] = = = −𝟏𝟐
𝟑 − (−𝟏𝟎) 𝟏𝟑
130
𝟒 𝟗
𝒇[𝒙𝟐 , 𝒙𝟏 ] − 𝒇[𝒙𝟏 , 𝒙𝟎 ] 𝟏−
𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ] =] = 𝟏𝟑 = 𝟏𝟑 = 𝟑
𝒙𝟐 − 𝒙𝟎 𝟓 − (−𝟏𝟎) 𝟏𝟓 𝟔𝟓
Determinación de:
𝟒
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ] =
𝟏𝟑
Determinación de:
𝟑
𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ] =
𝟔𝟓
Determinación de:
P2 (x) b0 b1 x x0 b2 x x0 x x1
𝟒 𝟑
𝑷𝟐 (𝒙) = 𝟒 + (𝒙 − (−𝟏𝟎)) + (𝒙 − (−𝟏𝟎))(𝒙 − 𝟑)
𝟏𝟑 𝟔𝟓
𝟒 𝟑
𝑷𝟐 (𝒙) = 𝟒 + (𝒙 + 𝟏𝟎) + (𝒙 + 𝟏𝟎)(𝒙 − 𝟑)
𝟏𝟑 𝟔𝟓
𝟒 𝟒𝟎 𝟑 𝟐
𝑷(𝒙) = 𝟒 + 𝒙+ + (𝒙 + 𝟕𝒙 − 𝟑𝟎)
𝟏𝟑 𝟏𝟑 𝟔𝟓
𝟒 𝟒𝟎 𝟑 𝟐 𝟐𝟏 𝟗𝟎
𝑷𝟐 (𝒙) = 𝟒 − 𝒙+ + 𝒙 + 𝒙−
𝟏𝟑 𝟏𝟑 𝟔𝟓 𝟔𝟓 𝟔𝟓
𝟑 𝟒𝟏 𝟕𝟒
𝑷𝟐 (𝒙) = 𝒙𝟐 + 𝒙+
𝟔𝟓 𝟔𝟓 𝟏𝟑
Procedimiento
Métodos Numéricos
sfs
131
Datos:
𝒙𝟎 = 𝟒 ∧ 𝒇(𝒙𝟎 ) = 𝟑
𝒙𝟏 = 𝟔 ∧ 𝒇(𝒙𝟏 ) = 𝟏𝟏
El polinomio de Newton de grado máximo que pasa por todos los puntos es:
P1(x) b0 b1 x x0
Con:
𝒃𝟎 = 𝒇(𝒙𝟎 )
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ]
Determinación de 𝒃𝟎 = 𝒇(𝒙𝟎 ):
𝒃𝟎 = 𝒇(𝒙𝟎 ) = 𝟑
Determinación de:
𝟏𝟏 − 𝟑 𝟖
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ] = = =𝟒
𝟔−𝟒 𝟐
Determinación de:
P1(x) b0 b1 x x0
𝑷𝟏 (𝒙) = 𝟑 + 𝟒(𝒙 − 𝟒)
𝑷𝟏 (𝒙) = 𝟑 + 𝟒𝒙 − 𝟏𝟔
𝑷𝟏 (𝒙) = −𝟏𝟑 + 𝟒𝒙
Métodos Numéricos
sfs
132
2.5.8 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO
2.5.9 TALLER DE ENTRENAMIENTO
Teniendo presente los conceptos vistos y los ejercicios de aprendizaje, resuelva el siguiente taller:
2. x -6 4 10
y -1 10 3
x 10 20
y -3 15
𝒙 -4 -1 4 5 8 9 11 13 20 25
𝒚 10 6 3 4 7 10 9 5 1 -4
133
3 UNIDAD 3 CÁLCULO INFINITESIMAL Y MÉTODOS
NUMÉRICOS
Video: Resolución Ecuaciones Diferenciales 1er orden por Euler, Euler Modificado
y Runge-Kutta 2da Parte
”Ver Video”: https://www.youtube.com/watch?v=aypLi3gIepU
Sea la función y f x, la primera derivada de esta función se obtiene aplicando la siguiente fórmula:
Métodos Numéricos
sfs
134
1(
𝑓 (𝑥 + ℎ) − 𝑓(𝑥)
𝑓 𝑥 ) = lim
ℎ→0 ℎ
Como la derivada se va a hallar numéricamente:
1(
𝑓(𝑥 + ℎ) − 𝑓(𝑥)
𝑓 𝑥 ) = lim −𝐸
ℎ→0 ℎ
Si se quiere evaluar la derivada en un valor 𝑥0 , la derivada se obtiene como:
1(
𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 )
𝑓 𝑥0 ) = lim −𝐸
ℎ→0 ℎ
Donde el error se obtiene como 𝐸 :
ℎ "
𝐸 = 𝑓 (𝑐)
2
Con: 𝑥0 < 𝑐 < 𝑥0 + ℎ
Se toma el valor de c que produce el mayor error, esto se llama encontrar una cota para el error.
Pero este error se puede calcular de manera exacta sólo si se conoce la función f x
Sean 𝒇(𝒙) = 𝐥𝐧 𝒙 ∧ 𝒙𝟎 = 𝟏, 𝟖
a. Como 𝒇(𝒙) = 𝐥𝐧 𝒙 ∧ 𝒙𝟎 = 𝟏, 𝟖
Métodos Numéricos
sfs
135
1 𝟏
𝑓 1 (𝑥) = ⟹ 𝒇(𝟏, 𝟖) = = 𝟎, 𝟓𝟓𝟓𝟓𝟓𝟓𝟓 …
𝑥 𝟏, 𝟖
b. Para h 0.1
1 𝒇(𝟏, 𝟖 + 𝟎, 𝟏) − 𝒇(𝟏, 𝟖) 𝐥𝐧 𝟏, 𝟗 − 𝐥𝐧 𝟏, 𝟖
𝑓 (1,8) = = = 𝟎, 𝟓𝟒𝟎𝟔𝟕𝟐𝟐𝟏 …
𝟎, 𝟏 𝟎, 𝟏
1 1
𝑓 1 (𝑥 ) = ⇒ 𝑓 " (𝑥 ) = − 2
𝑥 𝑥
0,1 " 0,1 1 0,1 1
𝐸𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑓 (𝑐) = × |− 2 | ≤ × ≤ 0,015432098
2 2 𝑐 2 1,82
Nota: El error se calculó con c 1.8 y no con c 1.9, ya que con c 1.8 se obtiene
un error mayor y de esta manera se puede garantizar que el error verdadero sea
menor que este error calculado
a. Para h 0.01
136
1(
1 "( )
1
𝑓 𝑥) = ⇒𝑓 𝑥 =− 2
𝑥 𝑥
b. Para h 0.01
1(
1 "
1
𝑓 𝑥 ) = ⇒ 𝑓 (𝑥 ) = − 2
𝑥 𝑥
d. Para h 0.001
137
El Error real es:
1(
1 "( )
1
𝑓 𝑥) = ⇒𝑓 𝑥 =− 2
𝑥 𝑥
___________________________________________
Dada:
𝑓(3+0,1)−𝑓(3)
𝒇(𝒙) = , Halle:
0,1
c) Error real.
d) Una cota para el error cometido.
Procedimiento:
138
40,55 − 38 2,55
𝑓(3) = = 𝑓(3) = = 25,5
0,1 0,1
c. Error Real:
𝑓 1 (𝑥 ) = 10𝑥 − 5 ⟹ 𝑓"(𝑥) = 10
𝑓"(3) = 10 ∧ 𝑓"(3,1) = 10
0,1 0,1
𝐸𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑓"(𝑐) = × |10|
2 2
𝟏
𝒇𝟏 (𝒙𝟎 ) = [−𝟑𝒇(𝒙𝟎 ) + 𝟒𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 + 𝟐𝒉)]
𝟐𝒉
Que se conoce como Fórmula de Tres Puntos.
𝒙𝟎
𝒙𝟎 + 𝒉
Métodos Numéricos
sfs
139
𝒙𝟎 + 𝟐𝒉
El Error está dado por:
𝒉𝟐 𝟏𝟏𝟏
𝑬= 𝒇 (𝒄)
𝟑
Con:
𝒙𝟎 < 𝒄 < 𝒙𝟎 + 𝟐𝒉
El valor de c se toma en este intervalo de tal manera que, con este se, obtenga el mayor error.
𝟏(
𝟏
𝒇 𝒙𝟎 ) = [𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 − 𝒉)]
𝟐𝒉
También se conoce como fórmula de. Los tres puntos son: x0 h, x0 x0 h aunque
x0 no aparece en la fórmula.
El error es:
𝒉𝟐 𝟏𝟏𝟏
𝑬= 𝒇 (𝒄)
𝟔
𝒙𝟎 − 𝒉 < 𝒄 < 𝒙𝟎 + 𝒉
El valor de c se toma en este intervalo de tal manera que, con este, se obtenga el mayor error.
𝟏
𝒇" (𝒙𝟎 ) = [𝒇(𝒙𝟎 − 𝒉) − 𝟐𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟎 + 𝒉)
𝒉𝟐
Métodos Numéricos
sfs
140
El Error es:
𝒉𝟐 (𝟒)
𝑬= 𝒇 (𝒄)
𝟏𝟐
Con 𝒙𝟎 − 𝒉 < 𝒄 < 𝒙𝟎 + 𝒉
x 1 1.5 2 2.5 3
y f x 6 9.5 15 22.5 32
Procedimiento:
𝟏(
𝟏
𝒇 𝒙𝟎 )−= [𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 − 𝒉)]
𝟐𝒉
𝟏
𝒇𝟏 (𝟐) = [𝒇(𝟐 + 𝟎, 𝟓) − 𝒇(𝟐 − 𝟎, 𝟓)] = 𝟐𝟐, 𝟓 − 𝟗, 𝟓
𝟐(𝟎, 𝟓)
𝒇𝟏 (𝟐) = 𝟏𝟑
141
𝟏
𝒇" (𝒙𝟎 ) = [𝒇(𝒙𝟎 − 𝒉) − 𝟐𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟎 + 𝒉)
𝒉𝟐
Reemplazando los valores conocidos, se tiene que:
𝟏
𝒇" (𝟐) = [𝒇(𝟐 − 𝟎, 𝟓) − 𝟐𝒇(𝟐) + 𝒇(𝟐 + 𝟎, 𝟓)
(𝟎, 𝟓)𝟐
𝟏
𝒇" (𝟐) = [𝒇(𝟏, 𝟓) − 𝟐𝒇(𝟐) + 𝒇(𝟐, 𝟓)
𝟎, 𝟐𝟓
𝒇" (𝟐) = 𝟒[𝟗, 𝟓 − 𝟐(𝟏𝟓) + 𝟐𝟐, 𝟓] = 𝟖
𝒙=𝟑
La curva corresponde a los datos que se muestran en la tabla.
Procedimiento:
𝟏
𝒎 = 𝒇𝟏 (𝟑) = [𝒇(𝟑, 𝟓) − 𝒇(𝟐, 𝟓)] = 𝟏𝟑
𝟐(𝟎, 𝟓)
𝒚 − 𝒚𝟎 = 𝒎(𝒙 − 𝒙𝟎 )
Métodos Numéricos
sfs
142
Donde:
𝒎 = 𝟏𝟑, 𝒙𝟎 = 𝟑 ∧ 𝒚𝟎 = 𝟏𝟑
Reemplazando en la ecuación de la recta se tiene:
Procedimiento
𝟏
𝒗 = 𝒇𝟏 (𝒙𝟎 )−= [𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 − 𝒉)]
𝟐𝒉
Con: 𝒙𝟎 = 𝟏, 𝟒 ∧ 𝒉 = 𝟎, 𝟐
Cálculo de la velocidad:
𝟏 𝟏
𝒗 = 𝒇𝟏 (𝟏, 𝟒)−= [𝒇(𝟏, 𝟒 + 𝟎, 𝟐) − 𝒇(𝟏, 𝟒 − 𝟎, 𝟐)] = [𝒇(𝟏, 𝟔) − 𝒇(𝟏, 𝟐)]
𝟐(𝟎, 𝟐) 𝟎, 𝟒
143
𝑚
𝑣 = 54
𝑠
La aceleración se obtiene como la segunda derivada de la posición:
𝟏
𝒂 = 𝒇" (𝒙𝟎 ) = [𝒇(𝒙𝟎 − 𝒉) − 𝟐𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟎 + 𝒉)
𝒉𝟐
Cálculo de la aceleración:
𝟏
𝒂 = 𝒇" (𝟏, 𝟒) = [𝒇(𝟏, 𝟒 − 𝟎, 𝟐) − 𝟐𝒇(𝟏, 𝟒) + 𝒇(𝟏, 𝟒 + 𝟎, 𝟐)
𝒉𝟐
𝟏
𝒂 = 𝒇" (𝟏, 𝟒) = [𝒇(𝟏, 𝟐) − 𝟐𝒇(𝟏, 𝟒) + 𝒇(𝟏, 𝟔)]
𝒉𝟐
144
3.3 TEMA 2 INTEGRACIÓN NUMÉRICA
En algunos casos al intentar resolver una integral por los métodos de integración conocidos es
una tarea difícil y en algunas ocasiones no es posible resolver dicha integral, adicionalmente si se
tiene una serie de datos y es necesario integrarlos es muy práctico aplicar algún método numérico
para resolverla.
𝒃
∫ 𝒇(𝒙)𝒅𝒙
𝒂
Deducción de la fórmula de la regla del trapecio para integrar.
Para ello se divide el área en áreas más pequeñas, y observe, en la figura 16, que las áreas más
pequeñas corresponden a trapecios.
Métodos Numéricos
sfs
145
Figura 15. Gráfica de una función y= f(x) (Autor Elkin Ceballos Gómez abril 18 de 2011)
Figura 16. División del área en n trapecios iguales. (Autor Elkin Ceballos Gómez abril 18 de 2011)
Métodos Numéricos
sfs
146
Se pude observar que el intervalo a, b está dividido en n sub intervalos de igual longitud. Como
la longitud de a, bes 𝑏 − 𝑎 la longitud de cada sub intervalo es:
𝒃−𝒂
El cual se llamará ℎ.
𝒏
Observe que:
𝒙𝟎 = 𝒂
𝒙𝟏 = 𝒂 + 𝒉
𝒙𝟐 = 𝒂 + 𝟐𝒉
𝒙𝟑 = 𝒂 + 𝟑𝒉
⋮
𝒙𝒏 = 𝒂 + 𝒏𝒉 = 𝒃
También se puede ver que en cada uno de los sub intervalos se forma un trapecio, el área A de la
b
𝟏
𝑨𝟏 = 𝒉[𝒇(𝒂) + 𝒇(𝒂 + 𝒉)]
𝟐
El segundo trapecio tiene área igual a:
Métodos Numéricos
sfs
147
𝟏
𝑨𝟐 = 𝒉[𝒇(𝒂 + 𝒉) + 𝒇(𝒂 + 𝟐𝒉)]
𝟐
El tercer trapecio tiene área igual a:
𝟏
𝑨𝟑 = 𝒉[𝒇(𝒂 + 𝟐𝒉) + 𝒇(𝒂 + 𝟑𝒉)]
𝟐
El área total es:
𝒃 𝟏
∫𝒂 𝒇(𝒙)𝒅𝒙 ≈ 𝟐h f (a) 2 f a h 2 f (a 2h) 2 f (a 3h) ... 2 f (a (n 1) h) f
(b)
Con:
𝒃−𝒂
𝒉= 𝒏
El error cometido en la regla trapezoidal está dado por:
ℎ3
𝐸𝑎 = | 𝑁𝑓 " (𝑐)| , 𝑐 ∈ [𝑎, 𝑏]
12
Se escoge el valor de c que produzca el máximo error; esto es lo que se llama: ENCONTRAR UNA COTA
PARA EL ERROR. Para utilizar esta fórmula se debe conocer la función a integrar.
1 1
Use la regla del trapecio para estimar el valor de ∫0 1+𝑥 2
Métodos Numéricos
sfs
148
Con: 𝒂 = 𝟎, 𝒃 = 𝟏 ∧ 𝑵 = 𝟓
PROCEDIMIENTO
𝟏−𝟎
𝒉=
𝟓
𝑓(𝑎) = 𝑓(0) = 𝟏
1
2𝑓(𝑎 + ℎ) = 2 × = 𝟏, 𝟗𝟐𝟑𝟏
1 + 0,22
1
2𝑓(𝑎 + 2ℎ) = 2 × = 𝟏, 𝟕𝟐𝟒𝟏
1 + 0,42
1
2𝑓(𝑎 + 3ℎ) = 2 × = 𝟏, 𝟒𝟕𝟎𝟔
1 + 0,62
1
2𝑓(𝑎 + 4ℎ) = 2 × = 𝟏, 𝟐𝟏𝟗𝟓
1 + 0,82
2𝑓(𝑏) = 𝑓(1) = 𝟎, 𝟓
1
1 1
∫ 2
≈ × 0,2 × 0,78373 = 0,784
0 1 + 𝑥 2
Error cometido:
Métodos Numéricos
sfs
149
𝟐(𝟑𝒙𝟐 − 𝟏)
𝒇"(𝒙) =
(𝟏 + 𝒙𝟐 )𝟑
Reemplazando los valores respectivos, se tiene:
𝟐(𝟑(𝟎)𝟐 − 𝟏)
𝒇"(𝟎) = = −𝟐
(𝟏 + (𝟎)𝟐 )𝟑
𝟐(𝟑(𝟏)𝟐 − 𝟏) 𝟒 𝟏
𝒇"(𝟏) = = =
(𝟏 + (𝟏)𝟐 )𝟑 𝟖 𝟐
(1 − 0)3 " 𝟏𝟑
𝐸𝑎 = | 𝑓 (0)| < × 𝟐 = 𝟔, 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟒 × 𝟏𝟎−𝟑 = 𝟎, 𝟎𝟎𝟔𝟔𝟔𝟔 …
12 × 52 12 × 52
ℎ3 (𝟎, 𝟐)𝟑
𝐸𝑎 = | 𝑁𝑓 " (𝑐)| < × 𝟓 × 𝟐 = 𝟎, 𝟎𝟎𝟔𝟔𝟔𝟔𝟔𝟔𝟔
12 12 ×
𝟏
𝟐
∫ 𝒆−𝒙 𝒅𝒙
𝟎
Procedimiento
Datos:
Métodos Numéricos
sfs
150
𝟐
𝒇(𝒙) = 𝒆−𝒙 , 𝒂 = 𝟎, 𝒃 = 𝟏, 𝒏 = 𝟔
1−0 1
ℎ= = = 0,1667
6 6
𝟐
𝒇(𝒂) = 𝒇(𝟎) = 𝒆−(𝟎) = 𝟏
𝟐
𝟐 × 𝒇(𝒂 + 𝒉) = 𝒇(𝟎, 𝟏𝟔𝟔𝟕) = 𝟐 × 𝒆−(𝟎,𝟏𝟔𝟔𝟕) = 𝟏, 𝟗𝟒𝟓𝟐
𝟐
𝟐 × 𝒇(𝒂 + 𝟐𝒉) = 𝒇(𝟎, 𝟑𝟑𝟑𝟑) = 𝟐 × 𝒆−(𝟎,𝟑𝟑𝟑𝟑) = 𝟏, 𝟕𝟖𝟗𝟕
𝟐
𝟐 × 𝒇(𝒂 + 𝟑𝒉) = 𝒇(𝟎, 𝟓) = 𝟐 × 𝒆−(𝟎,𝟓) = 𝟏, 𝟓𝟓𝟕𝟔
𝟐
𝟐 × 𝒇(𝒂 + 𝟒𝒉) = 𝒇(𝟎, 𝟔𝟔𝟔𝟕) = 𝟐 × 𝒆−(𝟎,𝟔𝟔𝟔𝟕) = 𝟏, 𝟐𝟖𝟐𝟒
𝟐
𝟐 × 𝒇(𝒂 + 𝟓𝒉) = 𝒇(𝟎, 𝟖𝟑𝟑𝟑) = 𝟐 × 𝒆−𝟎,𝟖𝟑𝟑𝟑) = 𝟎, 𝟗𝟗𝟖𝟖
𝟐
𝒇(𝒃) = 𝒇(𝟏) = 𝒆−(𝟏) = 𝟎, 𝟑𝟔𝟕𝟗
La suma de estos valores es: 𝟖, 𝟗𝟒𝟏𝟔
𝟏
−𝒙𝟐
𝟏 𝟏
∫ 𝒆 𝒅𝒙 = × × 𝟖, 𝟗𝟒𝟏𝟔 = 𝟎, 𝟕𝟒𝟓𝟏𝟑𝟑𝟑𝟑 ≈ 𝟎, 𝟕𝟒𝟓𝟏
𝟎 𝟐 𝟔
El error cometido en la regla trapezoidal está dado por:
ℎ3
𝐸𝑎 = | 𝑁𝑓 " (𝑐)| , 𝑐 ∈ [𝑎, 𝑏]
12
Estime el máximo error posible para el ejemplo anterior
𝟐 𝟐 𝟐
𝒇(𝒙) = 𝒆−𝒙 , 𝒇𝟏 (𝒙) = −𝟐𝒙𝒆−𝒙 ∧ 𝒇" (𝒙) = 𝒆−𝒙 (𝟒𝒙𝟐 − 𝟐)
Se tiene que:
Métodos Numéricos
sfs
151
" −𝒙 𝟐
|𝒇 (𝒙)| = |𝒆 (𝟒𝒙𝟐 − 𝟐)|
𝟐
|𝒇" (𝟎)| = |𝒆−𝟎 (𝟒(𝟎)𝟐 − 𝟐)| ≤ |𝟏(𝟎) − 𝟐| = |−𝟐| = 𝟐
𝟐
|𝒇" (𝟏)| = |𝒆−𝟏 (𝟒(𝟏)𝟐 − 𝟐)| = |𝟎, 𝟑𝟔𝟕𝟖𝟕𝟗𝟒𝟒𝟏(𝟎 − 𝟐)| = |𝟎, 𝟑𝟔𝟕𝟖𝟕𝟗𝟒𝟒𝟏(−𝟐)|
= 𝟎, 𝟕𝟑𝟓𝟕𝟖𝟖𝟖𝟐
Entonces:
𝟏 𝟑
( 𝟔)
|𝑬| ≤ × 𝟔 × 𝟐 ⇒ |𝑬| = 𝟎, 𝟎𝟎𝟒𝟔𝟐𝟗𝟔𝟗 …
𝟏𝟐
𝟏 𝟐 − 𝟏 𝒛𝟐
∫ 𝒆 𝟐 𝒅𝒛
𝟐𝝅 −𝟐
Procedimiento
𝟏 𝟐 − 𝟏 𝒛𝟐
∫ 𝒆 𝟐 𝒅𝒛
𝟐𝝅 −𝟐
𝟐 − (−𝟐) 𝟒
𝒉= = = 𝟎, 𝟓
𝟖 𝟖
𝒙𝟎 = −𝟐
𝒙𝟏 = −𝟐 + 𝟎, 𝟓 = −𝟏, 𝟓
𝒙𝟐 = −𝟏, 𝟓 + 𝟎, 𝟓 = −𝟏
𝒙𝟑 = −𝟏 + 𝟎, 𝟓 = −𝟎, 𝟓
𝒙𝟒 = −𝟎, 𝟓 + 𝟎, 𝟓 = 𝟎
Métodos Numéricos
sfs
152
𝒙𝟓 = 𝟎 + 𝟎, 𝟓 = 𝟎, 𝟓
𝒙𝟔 = 𝟎, 𝟓 + 𝟎, 𝟓 = 𝟏
𝒙𝟕 = 𝟏 + 𝟎, 𝟓 = 𝟏, 𝟓
𝒙𝟑 = 𝟏, 𝟓 + 𝟎, 𝟓 = 𝟐
Métodos Numéricos
sfs
153
FIGURA 17
𝟓 𝟑𝒙
2. Resuelva ∫𝟎 𝒙𝟐 +𝟒
a. De manera exacta, aplicando el método de integración por sustitución.
154
El problema a resolver en esta sección es que:
Dada una ecuación diferencial de la forma:
𝒅𝒚
= 𝒇(𝒙, 𝒚)
𝒅𝒙
Con la condición inicial:
𝒚(𝒙𝟎 ) = 𝒚𝟎
Halle el valor que da la función f(x) para otro valor dado de x, es decir, halle:
𝑦(𝑥𝑑𝑎𝑑𝑜 )
Para solucionar este problema que se llama P.V.I. (Problema de valor inicial) se analizarán dos
métodos.
𝒅𝒚
= 𝒇(𝒙, 𝒚)
𝒅𝒙
Con la condición inicial:
𝒚(𝒙𝟎 ) = 𝒚𝟎
Se tienen los siguientes datos:
En
𝒙 = 𝒙𝟎 , 𝒚𝟎 = 𝒚(𝒙𝟎 )
Como:
Métodos Numéricos
sfs
155
𝒅𝒚
= 𝒇(𝒙𝟎 , 𝒚𝟎 ) en 𝒙 = 𝒙𝟎 ∧ 𝒚 = 𝒚𝟎
𝒅𝒙
𝒅𝒚 𝒚𝟏 − 𝒚𝟎
𝒎= = 𝒇(𝒙𝟎 , 𝒚𝟎 ) =
𝒅𝒙 𝒙𝟏 − 𝒙𝟎
Tomando la igualdad:
𝒚𝟏 − 𝒚𝟎
𝒇(𝒙𝟎 , 𝒚𝟎 ) =
𝒙𝟏 − 𝒙𝟎
Despejando 𝒚𝟏 , quedaría:
𝒚𝟏 = 𝒚𝟎 + 𝒇(𝒙𝟎 , 𝒚𝟎 )(𝒙𝟏 − 𝒙𝟎 )
𝒚𝟐 = 𝒚𝟏 + 𝒇(𝒙𝟏 , 𝒚𝟏 )(𝒙𝟐 − 𝒙𝟏 )
En forma general:
𝒚𝒊+𝟏 = 𝒚𝒊 + 𝒇(𝒙𝒊 , 𝒚𝒊 )𝒉
Para: 𝒊 = 𝟏, 𝟐, 𝟑, … , 𝒏 − 𝟏
Llamada Fórmula de Euler.
Métodos Numéricos
sfs
156
Enlaces método de Euler.
𝒚∗𝑰+𝟏 = 𝒚𝒊 + 𝒇(𝒙𝒊 , 𝒚𝒊 )𝒉
Este valor se ingresa al corrector, que es la fórmula de EULER MEJORADO.
𝒅𝒚 𝟏 𝟐
Para el P.V.I: = (𝒚 + 𝒚) con condición inicial:
𝒅𝒙 𝒙
Procedimiento
Métodos Numéricos
sfs
157
Con 𝒉 = 𝟎, 𝟐, 𝒙𝟎 = 𝟏 ∧ 𝒚𝟎 = −𝟐
Se tiene que:
𝒙𝟏 = 𝟏, 𝟐; 𝒙𝟐 = 𝟏, 𝟒; 𝒙𝟑 = 𝟏, 𝟔
𝟏
𝒚𝟏 = 𝒚𝟎 + 𝒇(𝒙𝟎 , 𝒚𝟎 )𝒉 ⇒ 𝒚𝟏 = −𝟐 + [ ((−𝟐)𝟐 + (−𝟐))] 𝟎, 𝟐 = 𝟏, 𝟔
𝟏
𝒇(𝒙𝟎 , 𝒚𝟎 ) + 𝒇(𝒙𝟏 , 𝒚𝟏 )
𝒚𝟏 = 𝒚𝟎 + 𝒉⇒
𝟐
𝟏 𝟐 𝟏 𝟐
((−𝟐) + (−𝟐)) + (−𝟏, 𝟔))
𝟏 𝟏, 𝟐 ((−𝟏, 𝟔) +
𝒚𝟏 = −𝟐 + × 𝟎, 𝟐 ⇒
𝟐
𝒚𝟏 = −𝟏, 𝟕𝟐 = 𝒚(𝟏, 𝟐)
𝟏
𝒚∗𝟐 = 𝒚𝟏 + 𝒇(𝒙𝟏 , 𝒚𝟏 )𝒉 ⇒ 𝒚∗𝟏 = −𝟏, 𝟕𝟐 + [ ((−𝟏, 𝟕𝟐)𝟐 + (−𝟏, 𝟕𝟐))] 𝟎, 𝟐 = −𝟏, 𝟓𝟏𝟑𝟔
𝟏, 𝟐
158
𝒅𝒚
1. Para el P.V.I: = 𝟐𝒙𝒚 + 𝒚 𝒄𝒐𝒏 𝒚(𝟐) = 𝟐
𝒅𝒙
Use la fórmula de Euler mejorado con tamaños de paso de h 0.2 para determinar un valor
aproximado de 𝒚(𝟐, 𝟒)
Procedimiento:
h 0,2; 𝒙𝟎 2 𝒚𝟎 2
𝒙𝟏 = 𝟐 + 𝟎, 𝟐 = 𝟐, 𝟐
𝒙𝟐 = 𝟐, 𝟐 + 𝟎, 𝟐 = 𝟐, 𝟒
Cálculo del Predictor:
𝒚∗𝟏 = 𝒚𝟎 + 𝒇(𝒙𝟎 , 𝒚𝟎 )𝒉
𝒇(𝒙𝟎 , 𝒚𝟎 ) + 𝒇(𝒙𝟏 , 𝒚𝟏 )
𝒚𝟏 = 𝒚𝟎 + 𝒉
𝟐
[𝟐(𝟐)(𝟐) + 𝟐] + [𝟐(𝟐, 𝟐)(𝟒) + 𝟒]
𝒚𝟏 = 𝟐 + 𝟎, 𝟐
𝟐
𝒚𝟏 = 𝟓, 𝟏𝟔 = 𝒚(𝟐, 𝟐)
Cálculo del Predictor:
𝒚∗𝟐 = 𝒚𝟏 + 𝒇(𝒙𝟏 , 𝒚𝟏 )𝒉
159
Cálculo del Corrector:
𝒇(𝒙𝟏 , 𝒚𝟏 ) + 𝒇(𝒙𝟐 , 𝒚𝟐 )
𝒚𝟐 = 𝒚𝟏 + 𝒉
𝟐
[𝟐(𝟐, 𝟐)(𝟓, 𝟏𝟔) + 𝟓, 𝟏𝟔] + [𝟐(𝟐, 𝟒)(𝟏𝟎, 𝟕𝟑𝟐𝟖) + 𝟏𝟎, 𝟕𝟑𝟐𝟖]
𝒚𝟐 = 𝟓, 𝟏𝟔 + 𝟎, 𝟐
𝟐
160
𝑑𝑦 𝑥
= 𝑐𝑜𝑛 𝑦(2) = 3
𝑑𝑥 𝑦
𝑑𝑦 𝑥
= 𝑐𝑜𝑛 𝑦(2) = 4
𝑑𝑥 𝑦
𝟔
𝟐𝒙
∫ 𝟐 𝒅𝒙
𝟎 𝒙 + 𝟏
3. Aproximando a cuatro cifras decimales, el resultado exacto de la integral es:
a. 3.5835
b. 3.6109
c. 3.4174
d. 6.8348
4. Aplicando la regla del trapecio con 𝑛 = 6 y dando la respuesta final aproximada con
cuatro decimales, el resultado obtenido es:
a. 3.5835
b. 2.6333
c. 1.7918
d. 3.4174
161
X Y
(Tiempo) (distancia)
2 -46
2.5 -1,875
3 62
3.5 149,375
a. y (2) 1.36218
b. y (2) 1.13520
c. y (2) 1.25996
d. y (2) 1.08250
e. y (2) 1.12604
Métodos Numéricos
162
4 PISTAS DE APRENDIZAJE
No olvide Que para encontrar la raíz de una ecuación no lineal f x 0 , Se sabe que hay una
raíz en un intervalo a, bporque: La función f x es continua en a, by porque f a f b 0
Traer a la memoria La división entre cero no está definida en ningún campo numérico. Cuando
en el numerador hay un número diferente de cero y en el denominador está el cero se dice que
el resultado no existe; si en el numerador y en el denominador está el cero, se dice que el
resultado es indefinido.
El grado máximo del polinomio de interpolación debe ser menor que el número de puntos.
Tener en cuenta Que cuando se realiza la interpolación por el método de LaGrange y por el
método de Newton, el polinomio obtenido por ambos métodos debe ser el mismo.
Tenga Presente Que, al utilizar el método del Punto Fijo para solucionar una ecuación no lineal,
para que la
función gx, encontrada por usted, funcione, se debe cumplir que |𝑔1 (𝑥)| < 1 en el intervalo
donde encuentra la raíz
Traer a la memoria Que en física si f xes la función de posición de un objeto, la velocidad instantánea del
objeto se halla como v f 'x y función de aceleración instantánea se obtiene como a f ''x.
Métodos Numéricos
sfs
163
Traer a la memoria Para encontrar la ecuación de una recta conocidos; un punto de coordenadas
x0 , y0 y la pendiente m de la recta, se utiliza la ecuación punto pendiente de la línea recta
cuya forma es:
y y0 mx x0
m f ' x0
Tenga presente El orden en que se efectúa operaciones es: Primero potencias o raíces, luego
multiplicaciones o divisiones y por último sumas y restas.
Tenga presente Al utilizar la calculadora para resolver una ecuación lineal que involucra funciones
trigonométricas, la calculadora se debe llevar al modo radianes
Métodos Numéricos
164
5 GLOSARIO
Métodos numéricos: Son técnicas que permiten solucionar problemas utilizando operaciones
repetitivas de suma, resta, multiplicación y división.
Interpolación: Dado un conjunto de valores 𝑥, 𝑦 encontrar una función que represente estos
valores, que pase exactamente por cada uno de ellos y que permita estimar valores no observados
(es decir, que permita hacer predicciones o interpolaciones o aproximaciones).
Semilla: Es una estimación inicial de la raíz de una ecuación de la forma f x 0. Es el valor
inicial con el cual se empieza a operar el método numérico que se va a aplicar.
Punto fijo: Dada una función y=f xun punto fijo es un punto en el cual 𝑦 = 𝑥
Convergencia: Se dice que un método converge cuando se llega al valor buscado, esto se puede
identificar por qué se puede observar que cada cálculo efectuado lleva a un valor cada vez más
parecido al anterior.
Divergencia: Un método diverge cuando no se llega al valor buscado, esto se identifica, porque
cada cálculo efectuado lleva a un valor cada vez diferente.
Asíntota vertical: Es un valor de x en el cual el denominador de una función se hace igual a cero.
Métodos Numéricos
165
6 BIBLIOGRAFÍA
Asmar, I. (1999). Métodos numéricos un
primer curso. Medellín. http://books.google.com.co/books?id=65J
SdxLG4M4C&printsec=frontcover&dq=m%
Burden, & Faires, D. (2003). Análisis C3%A9todos
numérico. México: International
+num%C3%A9ricos&source=bl&ots=oe_U
Thomson Editores, S.A. 1iXJDU&sig=rfk15RFvg6S3O1Vx8wvgZXFC
8sE&hl=es&ei=
Chapra, S., & Canale, R. (2003). Métodos a_qLTeuuKsaCtgfahsm1DQ&sa=X&oi=boo
numéricos para ingenieros. . México D.F.: k_result&ct=result&resnum=10&ved=0CG
Mc Graw-Hill Interamericana. AQ6AEwCQ#v
=onepage&q&f=false
Haeussler, E., & Richard, P. S. (1997).
Mendoza, F. Matemáticas para ingenieros.
Matemáticas para Administración,
Recuperado el 24 de marzo de 2011 de,
Economía, Ciencias sociales y de la vida. http://matematicas5.galeon.com/index.htm
México: Prentice hall. l
Métodos numéricos. Recuperado
Purcell, E., & Varberg, D. (1993). Cálculo
el 24 de
con geometría analítica. México: Prentice
marzo de 2011,
Hall.
de
Asmar,I. (2000). Métodos numéricos. http://www.mitecnologico.com/Main/Met
Recuperado el 24 odosNumericos
de marzo de 2011,
Quintana, P, Villalobos, E “y” Cornejo, M.
de
http://www.unalmed.edu.co/~ifasmar/int (2005). Métodos numéricos con
ro.html aplicaciones en Excel. Recuperado el 24 de
marzo de 2011 de,
Díaz,W. (1998). Métodos numéricos.
http://books.google.com.co/books?id=wX
Recuperado el 24
de marzo de 2011, P3VzSa-
de http://www.uv.es/diaz/mn/fmn.html HIC&printsec=frontcover&dq=m%C3%A9to
dos+num%C3%A9ricos&source=bl&ots=j9J
4hx8T_Q&si
Huerta, A, Sarrate, J “y” Rodríguez, A. g=BnEOwJ1CgzurjEudVbSP2qAdTJk&hl=es
(2001). Métodos numéricos introducción,
&ei=a_qLTeuuKsaCtgfahsm1DQ&sa=X&oi
aplicaciones y programación.
Recuperado el =book_result&
24 de ct=result&resnum=9&ved=0CF0Q6AEwCA#
marzo de v=onepage&q&f=false
2011 de,