Metodo de Steffensen

1.
- Introduccin
En la resolucin de ecuaciones no lineales nos aparecen muchos problemas en forma natural, con la necesidad de calcular el valor de x donde una funcin f se anula, es decir, una raz de f. En general, con las herramientas analticas que se usan para estudiar y graficar funciones suaves (derivables) slo podemos analizar si hay un intervalo a, b! donde el gr"fico de f cruza el e#e x. En muchas ocasiones, slo tiene sentido encontrar una solucin aproximada. $ veces, el c"lculo exacto no es posible ya sea porque se trata de una raz irracional o porque la funcin viene dada por coeficientes cuyos valores se conocen slo en forma aproximada. %o importante al utilizar m&todos que estimen el valor deseado es poder controlar el error que se comete al utilizar un valor aproximado en lugar del exacto.
2.- Fundamento Terico

'ara una sucesin (xn), obtenida por el m&todo del punto fi#o xn*+ , f(xn), partimos de tres puntosy., f(x.) z., f(y.) donde x. es el punto inicial. /bteniendo asx+ , x. 0 (y. 0x.)+12 z. 0 23y. 0 x. En forma general4n*+ , xn 0 (yn 0 xn)+12 zn 0 23 yn 0 xn 5onde si 6xn*+ 0 xn6 , error 7 8ol entonces se satisface la tolerancia.
3.- Planteamiento del Problema

Obtener una raz de la funcin f(x) = x - cos x por la iteracin del punto fijo y aplicando el mtodo de Steffensen, si tomamos g(x) = cos x = 0 y error<0 0!
=0
= g(0) = !
4.- Desarrollo del mtodo

El mtodo de Steffensen consiste en aplicar un mtodo de Aitken modificado a una sucesin convergente obtenida de la iteracin de punto fijo. Para ello partiremos de tres puntos obtenidos mediante la iteracin en el punto fijo, que cumplen:
a!ora podemos aplicar Aitken a estos tres trminos
Para la obtencin de los siguientes trminos de la sucesin basta considerar los siguientes valores . .
. . . en la siguiente formula genrica
5.- Diseo del modelo matemtico

Deflexin de vigas en voladizo
"a ecuacin diferencial que determina el despla#amiento de una viga puede determinarse a partir de la relacin e$istente entre el momento flector %&$', el radio de curvatura k&$' de la curva (&$' que define la forma de la viga, el mdulo de elasticidad E &que depende del material' ( el momento de inercia dic!a relacin es la siguiente:
)eniendo en cuenta que el radio de curvatura depende de las dos primeras derivadas de (&$', puede deducirse la siguiente ecuacin no lineal:
donde " es la longitud de la viga. *uando la fle$in de la viga es peque+a, se tiene que (,&$' -. ,. Por tanto, asumiendo esta !iptesis de peque+as deformaciones, obtenemos:
.- !eri"icacin de los modelos matemticos

Si consideramos el caso de una viga en voladi#o, debemos imponer las siguientes condiciones iniciales: (&,' - ,, (,&,' - ,. "a condicin (&,' - , indica que el e$tremo i#quierdo de la viga est/ situado en el origen de coordenadas, mientras que la condicin (,&,' - , nos dice que dic!o e$tremo permanece fijo. 0e esta forma !emos e$presado (&$' como solucin de un problema de valores iniciales. Para resolverlo numricamente ser/ necesario escribirlo en forma de sistema. "a ecuacin diferencial obtenida depende de la forma del momento flector, que a su ve# depende de la distribucin de la carga que soporta la viga. 1amos a tener en cuenta dos casos: *arga uniformemente distribuida P &figura 2.3,':
*arga puntual P en el e$tremo de la viga en voladi#o &figura 2.33':
Figura 5.10: 0efle$in de una viga en voladi#o con carga uniformemente distribuida: esquema.
Figura 5.11: 0efle$in de una viga en voladi#o con carga puntual en un e$tremo: esquema.
*onsideremos una viga de longitud " - 4 m ( tomemos los valores E - 4 $3,33 5.m4 ( 6 - 7,82 $ 3,99 m:. Se representan los resultados obtenidos al aplicar el mtodo ;<: con 5 - 4, particiones ( P - 32,. En la figura 2.34 se !a supuesto que la carga est/ uniformemente distribuida a lo largo de la viga, representa el caso de una carga puntual situada en el e$tremo derec!o de la viga.
*omo puede observarse, en este =ltimo caso la !iptesis de peque+as deformaciones es menos plausible.
#.- Pro$ramacin del mtodo a%licado

function [x,i,tolf]=steffensen(x0,f,tolx,nmax) err=tolx+1; x=x0; phi=0; while(i<nmax & err>tolx) xx=x; fx =fe!al(f,x); tolf=tolx"a#s(phi); if a#s(fx )<=tolf #rea en$ fx %=fe!al(f,x+fx ); phi=(fx %&fx )'fx ; x=xx&fx 'phi; err=a#s(x&xx); i=i+1; en$
&.- '%licacin a un %roblema de su es%ecialidad

(inclu$e (inclu$e (inclu$e (inclu$e <cmath> <iomanip> <iostream> <conio%)h>
usin* namespace st$; ''prec+ error en la rai, ''fprec+ !alor maximo $e la funcion en el !alor acepta$o para la rai, ''numiter+ numero maximo $e iteraciones const const const const $ou#le prec=1)e&-,fprec=1)e&-; int nprec=10,nwi$th=1.; int numiter=100; char" ta#=/0t/;
'' f+funcion, f$+ $eri!a$a, fs+se*un$a $eri!a$a, *+ funcion para meto$o ''iterati!o $e 1 punto $ou#le $ou#le $ou#le $ou#le f($ou#le x); f$($ou#le x); fs($ou#le x); *($ou#le x);
int steffensen($ou#le , $ou#le);
'' 1efinicion $e la funcion 2 $eri!a$as $ou#le f($ou#le x)3 ''funcion return x"x&exp(&x"x"x); 4
$ou#le f$($ou#le x)3 ''primera $eri!a$a return %"x&exp(&x); 4 $ou#le fs($ou#le x)3 '' secon$ $eri!ati!e return %+exp(&x); 4 $ou#le *($ou#le x)3 '' * function x=*(x) return exp(&x"x"x)'x; 4 int main()3 '" 5ro*rama para comparar $i!ersos m6to$os $e calculo $e raices"' '' 5oner los puntos $e hor7uilla$o; x0 es el punto $e parti$a en meto$os $e 1 punto ''steffensen(x0,rai,); int steffensen( $ou#le x1, $ou#le rai,)3 $ou#le f($ou#le x); int n=0; $ou#le x%=0,$; $ou#le f1=f(x1); while(a#s(f1)>=fprec )3 if(n>=numiter)3 cout<</ numero maximo $e iteraciones exce$i$o en steffensen/<<en$l; return 1; 4 $= f(x1+f1)&f1; if ($==0)3cout<</$i!ision por 0 en steffensen/<<en$l; return 1;4 '' cout<</f1/<<f1<</$=/<<$<<en$l; x%=x1&f1"f1'$; x1=x%; f1=f(x1); '' cout<<x%<</f1/<<f1<<en$l; n++; 4 rai,=x%; cout<</steffensen /<<n&1<</ iteraciones /<< /rai,= /<<rai,<</ funcion= /<<f(rai,)<<en$l; return 0; 4
(.- )onclusiones
:in duda el proyecto y su programa podra solucionar problemas de tiempo y efectividad en el desarrollo de ecuaciones convirti&ndose en una gran venta#a para los estudiantes de ingenieria en general. El dise;o de este proyecto de investigacin formativa fue de gran importancia para demostrar que lo ense;ado durante el ciclo universitario fue muy bien aprendido y fue y ser" aplicado en mi futuro profesional.
1*.- +iblio$ra",a - ane.o

"nalisis #umrico $urden, % & ' (aires, ) *rupo +ditorial ,beroamericana, ! -./ "nalysis of #umerical 0et1ods ,saacson, + ' $is1op 2eller, 3 )o1n 4iley and Sons, ! -55 ,terati6e Solution Of #onlinear +7uations ,n Se6eral 8ariables ) 0 Orte9a ' 4 : %1einboldt "cademic ;ress ,nc

Metodo de Steffensen

Cargado por

Copyright:

Formatos disponibles

Metodo de Steffensen

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

Metodo de Steffensen

Cargado por

Copyright:

Formatos disponibles

1.

2.- Fundamento Terico

3.- Planteamiento del Problema

4.- Desarrollo del mtodo

a!ora podemos aplicar Aitken a estos tres trminos

. . . en la siguiente formula genrica

5.- Diseo del modelo matemtico

.- !eri"icacin de los modelos matemticos

*arga puntual P en el e$tremo de la viga en voladi#o &figura 2.33':

#.- Pro$ramacin del mtodo a%licado

&.- '%licacin a un %roblema de su es%ecialidad

int steffensen($ou#le , $ou#le);

1*.- +iblio$ra",a - ane.o

También podría gustarte