Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 127 vistas

    Método de Bairstow

    Cargado por

    Isaías Cocherán Montenegro
    Este documento presenta dos métodos numéricos para encontrar las raíces de un polinomio: el método de Bairstow y el método de Muller. El método de Bairstow es un algoritmo iterativo que encuentra las raíces de un polinomio real de grado arbitrario dividiéndolo en un polinomio cuadrático y un residuo. El método de Muller aproxima la gráfica de la función por una parábola que pasa por tres puntos iniciales para obtener una nueva aproximación a la raíz. También se proporcion

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd

    Método de Bairstow

    Cargado por

    Isaías Cocherán Montenegro
    127 vistas19 páginas
    Este documento presenta dos métodos numéricos para encontrar las raíces de un polinomio: el método de Bairstow y el método de Muller. El método de Bairstow es un algoritmo iterativo que encuentra las raíces de un polinomio real de grado arbitrario dividiéndolo en un polinomio cuadrático y un residuo. El método de Muller aproxima la gráfica de la función por una parábola que pasa por tres puntos iniciales para obtener una nueva aproximación a la raíz. También se proporcion

    Descripción original:

    métodos numéricos

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta dos métodos numéricos para encontrar las raíces de un polinomio: el método de Bairstow y el método de Muller. El método de Bairstow es un algoritmo iterativo que encuentra las raíces de un polinomio real de grado arbitrario dividiéndolo en un polinomio cuadrático y un residuo. El método de Muller aproxima la gráfica de la función por una parábola que pasa por tres puntos iniciales para obtener una nueva aproximación a la raíz. También se proporcion

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    127 vistas19 páginas

    Método de Bairstow

    Cargado por

    Isaías Cocherán Montenegro
    Este documento presenta dos métodos numéricos para encontrar las raíces de un polinomio: el método de Bairstow y el método de Muller. El método de Bairstow es un algoritmo iterativo que encuentra las raíces de un polinomio real de grado arbitrario dividiéndolo en un polinomio cuadrático y un residuo. El método de Muller aproxima la gráfica de la función por una parábola que pasa por tres puntos iniciales para obtener una nueva aproximación a la raíz. También se proporcion

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    de 19

    Mtodos Numricos

    Investigacin

    Profesor
    Grupo: 2IL121

    Estudiantes:

    Fecha de entrega:

    Martes 3 de septiembre de 2013

    Mtodo de Bairstow

    En anlisis numrico, el mtodo de Bairstow es un algoritmo eficiente


    de bsqueda de las races de un polinomio real de grado arbitrario. Es
    un mtodo iterativo, basado en el mtodo de Mller y de Newton
    Raphson. Dado un polinomio
    se encuentran dos factores, un
    polinomio cuadrtico
    y
    El procedimiento general para el mtodo de Bairstow es:
    Dado

    1.-Utilizando el mtodo de Newton Raphson


    calculamos
    y
    , tal que, el residuo
    de
    sea igual a cero.
    2.-Se determinan la races

    , utilizando la formula general.

    3.-Se calcula
    4.-Hacemos
    5.-Si el grado del polinomio es mayor que tres regresamos al paso 2;
    en caso contrario, terminamos.

    La principal diferencia de este mtodo, respecto a otros, es que


    permite calcular todas las races de un polinomio (reales e
    imaginarias).
    Para calcular la divisin de polinomios, hacemos uso de la divisin
    sinttica. As dado

    Al dividir entre
    siguiente polinomio

    con un residuo
    son.

    , tenemos como resultado el

    , el residuo ser cero solo si

    lo

    Los trminos b, los calculamos utilizamos divisin sinttica, la cual


    puede resolverse utilizando la siguiente relacin de recurrencia

    Una manera de determinar los valores de r y s que hacen cero el


    residuo es utilizar el Mtodo de Newton-Raphson. Para ello
    necesitamos una aproximacin lineal de
    respecto a r y s la cual
    calculamos utilizando la serie de Taylor

    donde los valores de r y s estn dados y calculamos los incrementos

    dr y ds que hacen a
    y
    igual a cero.
    El sistema de ecuaciones que tenemos que resolver es:

    Bairtow muestra que las derivadas parciales pueden obtener haciendo


    un procedimiento similar a la divisin sinttica, as

    donde

    Sustituyendo trmino
    ( Tomado de
    http://lc.fie.umich.mx/~calderon/programacion/Mnumericos/Bairstow.ht
    ml )
    Dado el polinomio f

    5
    (x) = x5- 3.5x4+ 2.75x3 + 2.125x2- 3.875x + 1.25, determinar los
    valores de r y s que hacen el resido igual a cero. Considere r0 = -1 y
    s0 = 2.
    Solucin.
    Iteracin 1.
    La divisin sinttica con el polinomio f2(x) = x2-x + 2.0 da como
    resultadof3(x) = x3- 4.5x2 + 9.25x -16.125 Residuo = {30.75, -61.75}
    Aplicando el mtodo de Newton tenemos
    -43.875

    16.75

    dr -30.75

    108.125

    -43.875 ds 61.75

    de donde
    r1 = -1.0 + 2.7636812508572213 =1.7636812508572213
    s1 = 2.0 + 5.403374022767796 =7.403374022767796
    Iteracin 2.
    La divisin sinttica con el polinomio f2(x) = x2-1.7636812508572213x
    - 7.403374022767796
    da como resultado
    f3(x) = x3- 1.7363187491427787x2 + 7.091061199392814x 1.776754563401905
    Residuo = {51.75640698828836, 105.68578319650365}
    Aplicando el mtodo de Newton tenemos
    27.628006

    14.542693 dr -51.75640

    208.148405

    27.62800

    de donde

    ds -105.68578

    r2 = 1.7636812508572213 - 0.04728019113442016 =
    1.7164010597228012
    s2 = 7.403374022767796 - 3.469106187802152 =
    3.934267834965644
    Iteracin 3.
    La divisin sinttica con el polinomio f2(x)= x2- 1.7164010597228012x
    - 3.934267834965644
    da como resultado
    f3(x) = x3- 1.7835989402771988x2 + 3.622896723753395x +
    1.3261878347051992
    Residuo = {12.654716254544885, 28.1881465309956}
    Aplicando el mtodo de Newton tenemos
    13.83497

    7.44182

    65.679212 13.83497

    dr -12.65471
    ds -28.18814

    de donde r3 = 1.7164010597228012 - 0.11666951305731528 =


    1.599731546665486
    s3 = 3.934267834965644 - 1.4835870659929915 =
    2.4506807689726524
    En resumen
    k rs

    Residuo

    0 -1 2 30.75 -61.75
    1 1.76368 7.403374 51.756406 105.68578
    2 1.71640 3.93426 12.65471 28.18814
    3 1.599731 2.450680 2.89958 8.15467
    4 1.33354 2.18666 0.760122 2.522228

    5 1.11826 2.11302 0.271940 0.607688


    6 1.02705 2.02317 0.04313 0.11185
    7 1.00165 2.00153 0.00277 0.00634
    8 1.00000 2.00000 1.13930E-5 2.67534E-5
    La solucin es:
    f3(x) = x3- 2.53x2 + 2.25x - 0.625 y f2(x) = x2- x 2
    Las races de f2(x) = x2- x - 2, son x1 = 2

    x2 = -1

    Si repetimos el ejemplo pero ahora considerando el polinomio f3(x) =


    x3- 2.53x2 + 2.25x - 0.625 ,
    Podemos calcular el total de las races del polinomio original.
    Ejemplo realizado por el estudiante
    Use el mtodo de Baiestow para determinar las races de
    f ( x) 2 x 6,2 x 4 x 2 0,7 x 3

    con r=1 y s=-1 y un porcentaje de error de

    1%
    Demostracin de la primera iteracin y encontrar los valores de b y c
    b0 2 (1)( 2.2) (1)( 3.3) 3.5
    b1 6.2 (1)( 3.3) (0.7)( 1) 2.2
    b2 4 (1)(0,7) 3.3
    b3 0.7
    c1 2.2 (1)( 2.6) ( 1)(0.7) 1.1
    c 2 3.3 (1)(0.7) 2.6
    c 3 0 .7

    Por lo tanto se debe resolver

    2.6r 0.7 s 2.2


    1.2r 2.6s 3.5
    2.2 2.6r
    0.7
    2.2 2.6r
    1.1r 2.6(
    ) 3.5
    0.7
    r 1.0846
    s

    Remplazando
    s

    2.2 2.6(1.0846)
    0.8856
    0.7

    Corrigiendo los valores iniciales


    r 1 1.0846 2.0846
    s 1 0.8856 0.1144

    El porcentaje de error aproximado


    ea ,r

    1.0847
    100 52%
    2.0846

    ea ,s

    0.8856
    100 774%
    0.1144

    Otras iteraciones

    Al llegar a la convergencia deseada, para determinar los coeficientes r


    y s, procedemos a calcular dos races de la funcin original a travs de

    r r 2 4s
    2
    r 2.436
    s 0.871
    x

    2.436 (2.436) 2 4(0.871)


    2
    x1 2.006
    x

    x 2 0.4354

    Para calcular la raz faltante:


    A0
    0.7

    A1
    -2.295

    A2
    -0.0002

    (2.295) (2.295) 2 4(0.7)( 0.0002)


    2(0.7)
    x3 3.2787
    x

    El otro valor de la raz se desprecia por ser tan pequea solo se


    buscan races grandes.

    Mtodo de Muller

    Este mtodo utilizado para encontrar races de ecuaciones con races


    mltiples, y consiste en obtener los coeficientes de la parbola que
    pasa por tres puntos elegidos. Dichos coeficientes son sustituidos en
    la formula cuadrtica para obtener el valor donde la parbola
    intersecta al eje X; es decir, la raz estimada. La aproximacin se
    puede facilitar, si se escribe la ecuacin de la parbola en una forma
    conveniente.
    Una de las mayores ventajas de este mtodo, es que al trabajar con la
    formula cuadrtica es posible localizar tanto races reales, como races
    complejas.
    Una de las mayores ventajas de este mtodo, es que al trabajar con la
    formula cuadrtica es posible localizar tanto races reales, como races
    complejas.
    Caractersticas del mtodo Muller.
    Converge cuadrticamente en un intervalo cercano a la raz.
    No requiere evaluar la primera derivada.
    Se obtiene races reales y complejas (cuando las races sean
    repetidas.
    Requiere valores inciales
    Extensin del mtodo de la secante; aprox. la grfica de la funcin
    f(x) por una lnea recta que pasa por x los puntos (xi -1,f(x1)) su
    interseccin con el eje x da una nueva aproximacin (Xi -1).
    n=2 segundo
    Se tom 3 valores inciales X0, X1, X2 y est ene. Polinomio P(x) de
    segundo grado que pasa por los puntos (X0,f(X0)), (X1,f(X1)) y
    (X2,f(X2)).
    Se toman sus races de p(x), la ms cercana a X 2, como la
    siguiente aproximacin a X3

    Frmula
    Los tres valores iniciales necesitados son denotados como xk, xk-1 yxk-2.
    La parbola pasa a travs de los puntos: (xk, f(xk)), (xk-1, f(xk-1)) y (xk2, f(xk-2)), si se escribe en la forma de Newton, entonces:
    donde f[xk, xk-1] y f[xk, xk-1, xk-2] denotan restas divididas. Esto puede
    ser escrito como:
    Donde:
    La prxima iteracin esta dada por la raz que brinda la ecuacin y =
    0.

    Determina la siguiente aproximacin x3 encontrando la


    interseccin con el eje x de la parbola definida por los puntos
    (x0,f(x0)), (x1,f(x1)), (x2,f(x2)).

    Procedimiento

    Se determina un X0, X1 y un X2.


    Segundo paso :
    h0 = X1 X0
    h1 = X2 X1
    Tercer paso:

    0 = F (X1) - F (X0)
    h0
    1 = F (X2) - F (X1)
    h1
    Cuarto paso:
    Se obtienen:
    a = 1 0
    h1 + h0
    b = a * h1 + 0
    c = F (X2)
    Quinto paso:
    x3 x 2

    2c
    b b 2 4ac

    Sexto paso:
    b b 2 4ac b b 2 4ac

    Si

    Se escoge b+ o si no b Calculo de error


    ea

    x3 x 2
    100%
    x3

    METODO DE MULLER

    E para:
    Parte a para F(x) = x^3 + x^2 4x 4 :

    X0 = 1
    X1 = 1.5

    X2 = 1.75

    Iteraciones

    X3

    Ea (%)

    1.75

    ---------------

    2.0112

    12.9863

    1.999882423

    0.5648

    1.99999997

    0.0059

    1.3686 * 10 ^ 6

    Parte b para: F(x) = x^3 0.5x^2 + 4x - 2

    X0 = 0.4

    X1 = 0.6
    Iteraciones

    X3

    Ea (%)

    0.8

    ---------------

    0.5007

    59.7750

    0.49999

    0.141817

    0.500000

    0.00100

    X2 = 0.8

    Ejemplo realizado por el estudiante


    Se ah F(x) = x^3 13x -12 encuentre las races con un porcentaje de
    error de 1%

    X0 = 4.5
    X1 = 5.5

    Iteraciones

    X3

    Ea (%)

    ---------------

    3.9765

    25.7391

    4.0011

    0.6139

    4.0000

    0.0262

    4.0000

    1.7631 * 10 ^ - 5

    MTODO DE HORNER
    Los Pasos a seguir para encontrar las Races son los siguientes:
    1. Determinar por medio de la Regla de los Signos de
    Descartes y hasta donde sea posible el nmero de Races
    Positivas, Negativas e Imaginarias que la Ecuacin Pueda tener.
    2. Encontrar por medio de la Relacin P en Q ( p/q ) ( P son
    factores del termino independiente y Q los factores del
    coeficiente principal respectivamente) las posibles Races
    Racionales.
    3. Probar las posibles Races por medio de la Divisin Sinttica:

    4. Las Races sern aquellas que al calcularlas por medio de


    la Divisin Sinttica el Residuo sea CERO.
    5. Cada vez que obtengas una Raz debes de Separarla, y
    utilizars la Ecuacin Reducida ( tercera lnea ) para buscar la
    siguiente.
    6. Continuars este proceso hasta que la Ecuacin Reducida sea
    de Grado 2 (cuadrtica ) la cual se podr resolver por
    Factorizacin por la Formula General.
    Ejemplo: Encontrar las Races de la siguiente Ecuacin.

    Regla de los Signos de Descartes


    +

    Calcular la Relacin

    Probar cada una de estas posibles Races por medio de


    la Divisin Sinttica sea calcular +1, +2, +4, -1, -2 y 4:

    Como te puedes dar cuenta si ya las calculaste 2 y 1 son


    Races, puesto que el Residuo en ambas es Cero, Calcula
    primero el 2 despus utiliza la ecuacin Reducida ( tercera
    lnea ) que es de Grado 3 para calcular el nmero 1, y
    obtendrs la siguiente Ecuacin Reducida que es de Grado 2.
    Las Races faltantes las calculars por medio de la Formula
    General y son

    Bibliografa

    http://mat156.wordpress.com/category/muller/
    http://es.scribd.com/doc/50377833/METODO-DE-HORNER
    https://lc.fie.umich.mx/~calderon/programacion/Mnumericos/Muller.html
    http://es.wikipedia.org/wiki/M%C3%A9todo_de_Bairstow

    https://lc.fie.umich.mx/~calderon/programacion/Mnumericos/Bairstow.html

    También podría gustarte