| Autor |
Mensaje |
leandrob_90
Nivel 9
Edad: 34
Registrado: 17 Ago 2009
Mensajes: 1586
Ubicación: Mundo de los Ryuo Shin
Carrera: Mecánica
|
|
Buenas, estaba viendo algunos parciales viejos de numérico y me topé con el siguiente problema:
Dado un numero positivo ![[tex]a>=2[/tex]](images/latex/45d7d9b6d2c01d41bd1ee12d09dc9e1a20a2bea7_0.png "[tex]a>=2[/tex]") y la función ![[tex]f(x)=e^x-a[/tex]](images/latex/ac42c094c78dc8dec0c744c8acf5cf2db1a46c61_0.png "[tex]f(x)=e^x-a[/tex]")
Me pide que determine el intervalo de convergencia y la cantidad de iteraciones para que el error absoluto sea ![[tex]10^{-4}[/tex]](images/latex/d2bcca4c951d1c8369012777d5f2d97df4270366_0.png "[tex]10^{-4}[/tex]") para el método NR
Cómo se determinan esas dos cosas? hay alguna forma? algún método? porque por lo que estuve viendo, el intervalo se saca a pulmón, y la cantidad de iteraciones la única formula que vi fue para el método de bisección (dudo que valga acá también)
|
|
|
|
_________________
leandrob_90
Revivamos el Chat-FIUBA
¿Qué te pasó foro? Antes eras chévere.
Por un ping-pong libre, popular y soberano.
|
|
   |
     |
|
Sebastian Santisi
Administrador Técnico
Edad: 42
Registrado: 23 Ago 2005
Mensajes: 17451
|
|
| ¿A qué te referís con lo de "a pulmón"?, digamos, NR te propone resolver el problema como un problema de punto fijo, las hipótesis que tenés que cumplir son las de PF y podés sacarlas analíticamente.
|
|
|
|
_________________
 ![[tex] ${. \ \ \ \ \ \ \ \ \ .}$ [/tex]](images/latex/d678cf273c99d2546c259db3422b3d968b184721_0.png "[tex] ${. \ \ \ \ \ \ \ \ \ .}$ [/tex]") ![[tex] ${\Large Usá \LaTeX, no seas foro...}$ [/tex]](images/latex/1e6580e1ed9e054ddefd65c0551c670c44f57f38_0.png "[tex] ${\Large Usá \LaTeX, no seas foro...}$ [/tex]")
|
|
 |
 |
|
4WD
Administrador
Edad: 39
Registrado: 07 Sep 2006
Mensajes: 2430
Ubicación: Ingeniero
Carrera: Mecánica
|
|
Para el intervalo, sí, tenés que encontrar "a mano" un intervalo que cumpla las condiciones del teorema de punto fijo, o hay algunas condiciones más restrictivas pero fáciles para NR.
Para estimar la cantidad de iteraciones, podés estimarlas sabiendo que NR converge cuadráticamente (casi siempre), en base al primer error que obtengas con iterar usando NR. Para esto necesitás saber la definición de orden de convergencia, que relaciona los errores en iteraciones sucesivas.
|
|
|
|
_________________
 
 
 
|
|
  |
|
|
leandrob_90
Nivel 9
Edad: 34
Registrado: 17 Ago 2009
Mensajes: 1586
Ubicación: Mundo de los Ryuo Shin
Carrera: Mecánica
|
|
Sí, con "a pulmón" me refería a "a mano"... Entonces, no hay forma de hallar el intervalo que no sea probando hasta que los valores den.
El orden de convergencia, para n "grandes" lo calculo como ![[tex]\alpha \approx \frac{ \log \left( |e_{n+2}|/|e_{n+1}| \right) }{ \log \left( |e_{n+1}|/|e_{n}| \right) }[/tex]](images/latex/ee6611ff1f04b1e451a9f8f7c7872284425a3645_0.png "[tex]\alpha \approx \frac{ \log \left( |e_{n+2}|/|e_{n+1}| \right) }{ \log \left( |e_{n+1}|/|e_{n}| \right) }[/tex]") al menos así lo vi para bisección y regula falsi, acá aplica también?
No entiendo qué tengo que hacer ahora... 
|
|
|
|
_________________
leandrob_90
Revivamos el Chat-FIUBA
¿Qué te pasó foro? Antes eras chévere.
Por un ping-pong libre, popular y soberano.
|
|
| |
|
|
leandrob_90
Nivel 9
Edad: 34
Registrado: 17 Ago 2009
Mensajes: 1586
Ubicación: Mundo de los Ryuo Shin
Carrera: Mecánica
|
|
Acá encontré algo, con la fórmula para acotar el error:
![[tex]|x_n-p| \leq \frac{ K^n }{|1-K|} |x_1-x_0|<\varepsilon[/tex]](images/latex/c31d3b1883bf0e3c7102ae2326e6ad0ebe583b37_0.png "[tex]|x_n-p| \leq \frac{ K^n }{|1-K|} |x_1-x_0|<\varepsilon[/tex]")
usando el término del medio y la derecha despejo la cantidad de iteraciones n, siendo K la cota de la derivada de g en el intervalo
|
|
|
|
_________________
leandrob_90
Revivamos el Chat-FIUBA
¿Qué te pasó foro? Antes eras chévere.
Por un ping-pong libre, popular y soberano.
|
|
| |
|
|
4WD
Administrador
Edad: 39
Registrado: 07 Sep 2006
Mensajes: 2430
Ubicación: Ingeniero
Carrera: Mecánica
|
|
Correcto!
Lo primero que pusiste no te sirve mucho: ya sabés que NR converge con orden 2 (excepto casos particulares), no hace falta estimarlo.
Yo había pensado en algo más burdo y aproximado:
![[tex]\left|x_{n+1} - r\right| \leq c \cdot \left|x_n - r\right|^2[/tex]](images/latex/ca34ebe80fe1452686a1cb4046fec1f45f321bd5_0.png "[tex]\left|x_{n+1} - r\right| \leq c \cdot \left|x_n - r\right|^2[/tex]")
Es decir, cada iteración te eleva al cuadrado el error (que es menor que 1).
Pero la expresión que pusiste es mejor, y la correcta para hacer lo que te piden.
Recordar que ![[tex]g[/tex]](images/latex/ec5d8a9a3f42318688bd76c52018a51b90659195_0.png "[tex]g[/tex]") es la función de punto fijo (la función de NR en este caso) y no la función a resolver.
|
|
|
|
_________________
|
|
| |
|
|
|
|
|
Ver tema siguiente
Ver tema anterior
Podés publicar nuevos temas en este foro
No podés responder a temas en este foro
No podés editar tus mensajes en este foro
No podés borrar tus mensajes en este foro
No podés votar en encuestas en este foro
No Podéspostear archivos en este foro
No Podés bajar archivos de este foro
|
Todas las horas son ART, ARST (GMT - 3, GMT - 2 Horas)
Protected by CBACK CrackerTracker365 Attacks blocked.
|
|
|
[ Tiempo: 0.3589s ][ Pedidos: 20 (0.2657s) ] |
