Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
128 6. POLINOMIOS ORTOGONALES Y APROXIMACIÓN POR CUADRADOS MÍNIMOS<br />
Terminamos el capítulo dando una forma más eficiente de encontrar los polinomios ortogonales<br />
asociados a un producto interno. En efecto, el siguiente teorema muestra que cada polinomio qn<br />
se escribe en función de los dos anteriores y por lo tanto la sucesión de polinomios ortogonales<br />
mónicos puede obtenerse por recurrencia.<br />
Teorema 6.20. Si un producto interno en C[a, b] satisface 〈xf, g〉 = 〈f, xg〉 entonces los polinomios<br />
ortogonales mónicos qn satisfacen la relación de recurrencia<br />
donde an y bn están dados por<br />
qn(x) = (x − an)qn−1(x) − bnqn−2(x) , ∀n ≥ 2 (6.5)<br />
an = 〈xqn−1, qn−1〉<br />
〈qn−1, qn−1〉<br />
y bn = 〈qn−1, qn−1〉<br />
〈qn−2, qn−2〉 .<br />
Demostración. Sea n ≥ 2. Como cero es raíz del polinomio qn(x) − qn(0) podemos escribir<br />
qn(x) − qn(0) = xrn−1<br />
donde rn−1 es un polinomio de grado menor o igual que n − 1. Además, como qn es mónico, rn−1<br />
también lo es. Tenemos entonces,<br />
qn(x) = xrn−1(x) + qn(0) = xqn−1(x) + x(rn−1(x) − qn−1(x)) + qn(0). (6.6)<br />
Pero como rn−1 y qn−1 son mónicos su diferencia resulta un polinomio de grado menor o igual<br />
que n − 2 y por lo tanto, como q0, . . . , qn−1 forman una base de Pn−1, existen coeficientes βj<br />
tales que<br />
n−1 <br />
x(rn−1(x) − qn−1(x)) + qn(0) = βjqj(x)<br />
y reemplazando en (6.6) obtenemos<br />
Ahora, para i < n − 2, tenemos<br />
j=0<br />
n−1 <br />
qn(x) = xqn−1(x) + βjqj(x). (6.7)<br />
n−2 <br />
0 = 〈qn, qi〉 = 〈(x + βn−1)qn−1 + 〈βjqj, qi〉 = 〈xqn−1, qi〉 + βi〈qi, qi〉<br />
j=0<br />
donde en el último paso hemos usado la ortogonalidad de los qj. Pero, como xqi es un polinomio<br />
de grado menor que n − 1, resulta<br />
j=0<br />
〈xqn−1, qi〉 = 〈qn−1, xqi〉 = 0<br />
y en consecuencia βi = 0 para todo i < n−2. Por lo tanto, definiendo an = −βn−1 y bn = −βn−2,<br />
(6.5) se obtiene de (6.7).<br />
Finalmente, usando (6.5) y la ortogonalidad de los qj tenemos,<br />
0 = 〈qn, qn−1〉 = 〈xqn−1, qn−1〉 − an〈qn−1, qn−1〉<br />
de donde se obtiene la expresión para an. Análogamente,<br />
0 = 〈qn, qn−2〉 = 〈xqn−1, qn−2〉 − bn〈qn−2, qn−2〉