2
Universidade Federal do Ceará – Centro de Tecnologia Data: 09/09/2014 Curso de Engenharia de Energias e Meio Ambiente – Métodos Numéricos para EEMA Atividade Computacional – Raízes de Polinômios – 2014.2 Prof. Tarcisio Ferreira Maciel, Dr.-Ing. Exercício 1. Considere um polinômio p n (x) de grau n com n +1 coeficientes reais da forma p n (x)= a n x n + a n-1 x n-1 + a n-2 x n-2 + ... + a 2 x 2 + a 1 x + a 0 , = n X i=0 a i x i , com a n 6=0, (1) e n raízes reais x 1 , x 2 , x 3 , ..., x n . Um algoritmo eficiente para o cálculo das raízes de p n (x) é o algoritmo de Briot-Ruffini-Horner, também chamado de algoritmo de deflação polinomial, que é obtido reescrevendo p n (x) como p n (x)=(... ( b 0 z }| { b 2 z }| { b n-1 z }| { a n |{z} bn x + a n-1 ) x + a n-2 ) | {z } b n-2 x + ... + a 2 ) x + a 1 ) | {z } b 1 x + a 0 . (2) De acordo com a equação acima, podemos definir b n = a n , b n-1 = a n x + a n-1 = b n x + a n-1 , b n-2 =(a n x + a n-1 )x + a n-2 = b n-1 x + a n-1 , ... b 2 =(... (a n x + a n-1 )x + a n-2 )x + ... + a 2 = b 3 x + a 2 , b 1 =(... (a n x + a n-1 )x + a n-2 )x + ... + a 2 )x + a 1 )= b 2 x + a 1 , b 0 =(... (a n x + a n-1 )x + a n-2 )x + ... + a 2 )x + a 1 )x + a 0 = b 1 x + a 0 = p n (x), (3) Na formulação acima, ¯ x é uma raiz de p n (x) p n x)=0 b 0 =0. Aplicando a b i para i =1, 2,...,n,o mesmo algoritmo aplicado a a i para i =0, 1, 2,...,n, podemos gerar c n = b n , c n-1 = c n x + b n-1 , c n-2 = c n-1 x + b n-2 , ... c 2 = c 3 x + b 2 , c 1 = c 2 x + b 1 . (4) É possível mostrar que c 1 = p 0 n (x). Como p n (x) é um polinômio, o método de Newton para o mesmo pode ser expresso como x i+1 = x i - b 0 (x i ) c 1 (x i ) . Se ¯ x é uma raiz de p n (x), então os coeficientes b n ,b n-1 ,...,b 2 ,b 1 , são os coeficientes do polinômio q n-1 (x) de grau n - 1 onde p n (x)=(x - ¯ x)q n-1 (x), (5) e todas as raízes de q n-1 (x) são também raízes de p n (x). Assim, o processo acima pode ser repetido para q n-1 (x) gerando um novo polinômio q n-2 (x) de grau n - 2, e assim sucessivamente, permitindo determinar todas as raízes de p(x).

AtividadeRaizesPolinomios

Embed Size (px)

Citation preview

  • Universidade Federal do Cear Centro de Tecnologia Data: 09/09/2014Curso de Engenharia de Energias e Meio Ambiente Mtodos Numricos para EEMAAtividade Computacional Razes de Polinmios 2014.2Prof. Tarcisio Ferreira Maciel, Dr.-Ing.

    Exerccio 1. Considere um polinmio pn(x) de grau n com n + 1 coeficientes reais da forma

    pn(x) = anxn + an1xn1 + an2xn2 + . . . + a2x2 + a1x + a0,

    =

    ni=0

    aixi, com an 6= 0,

    (1)

    e n razes reais x1, x2, x3, . . ., xn. Um algoritmo eficiente para o clculo das razes de pn(x) o algoritmode Briot-Ruffini-Horner, tambm chamado de algoritmo de deflao polinomial, que obtido reescrevendopn(x) como

    pn(x) = (. . . (

    b0 b2

    bn1 anbn

    x + an1)x + an2)

    bn2

    x + . . . + a2)x + a1)

    b1

    x + a0 . (2)

    De acordo com a equao acima, podemos definir

    bn = an,

    bn1 = anx + an1 = bnx + an1,bn2 = (anx + an1)x + an2 = bn1x + an1,. . .

    b2 = (. . . (anx + an1)x + an2)x + . . . + a2 = b3x + a2,b1 = (. . . (anx + an1)x + an2)x + . . . + a2)x + a1) = b2x + a1,b0 = (. . . (anx + an1)x + an2)x + . . . + a2)x + a1)x + a0 = b1x + a0 = pn(x),

    (3)

    Na formulao acima, x uma raiz de pn(x) pn(x) = 0 b0 = 0. Aplicando a bi para i = 1, 2, . . . , n, omesmo algoritmo aplicado a ai para i = 0, 1, 2, . . . , n, podemos gerar

    cn = bn,

    cn1 = cnx + bn1,cn2 = cn1x + bn2,. . .

    c2 = c3x + b2,

    c1 = c2x + b1.

    (4)

    possvel mostrar que c1 = pn(x). Como pn(x) um polinmio, o mtodo de Newton para o mesmo pode ser

    expresso como xi+1 = xi b0(xi)c1(xi)

    . Se x uma raiz de pn(x), ento os coeficientes bn, bn1, . . . , b2, b1, so

    os coeficientes do polinmio qn1(x) de grau n 1 ondepn(x) = (x x)qn1(x), (5)

    e todas as razes de qn1(x) so tambm razes de pn(x). Assim, o processo acima pode ser repetido paraqn1(x) gerando um novo polinmio qn2(x) de grau n 2, e assim sucessivamente, permitindo determinartodas as razes de p(x).

  • 1. Com base no exposto acima e nas referncias bibliogrficas do curso, implemente em Octave (ou MATLAB)o algoritmo de Briot-Ruffini-Horner para o clculo das razes de um polinmio.

    2. Apresente um relatrio sucinto contendo uma introduo sobre o tema "clculo de razes de polinmios",uma descrio matemtica do mtodo Briot-Ruffini-Horner apresentando seus principais aspectos, umaimplementao do algoritmo em Octave (ou MATLAB) e um exemplo de aplicao comparando o resultadodo seu cdigo (e.g., em termos de tempo de execuo e preciso dos resultados) ao das funes proprietriasdo Octave (ou MATLAB).

    3. A atividade dever ser desenvolvida em grupo de no mais que 03 alunos. Todos os membros do grupodevem ter domnio completo sobre todo o contedo da atividade. A atividade deve ser entregue na forma derelatrio obedecendo o Guia de Normalizao de Trabalhos Acadmicos da UFC.