Polinomios de Chebyshev

De testwiki
Revisión feita o 19 de xullo de 2024 ás 17:08 por imported>Breogan2008
(dif) ← Revisión máis antiga | Revisión actual (dif) | Revisión máis nova → (dif)
Saltar á navegación Saltar á procura

Os polinomios de Chebyshev son dúas secuencias de polinomios relacionadas coas funcións coseno e seno, escritas como Tn(x) e Un(x). Pódense definir de varias formas equivalentes, unha das cales comeza con funcións trigonométricas:

Os polinomios de Chebyshev do primeiro tipo Tn están definidos por

Tn(cosθ)=cos(nθ).

Os polinomios de Chebyshev do segundo tipo Un están definidos por

Un(cosθ)sinθ=sin((n+1)θ).

Que estas expresións definen polinomios en cosθ pode non ser obvio a primeira vista, mais conséguese escribindo cos(nθ) e sin((n+1)θ) coa fórmula de de Moivre ou mediante as fórmulas de suma de ángulos para cos e sin repetidamente. Por exemplo, coas fórmulas de ángulo duplo, obtemos

T2(cosθ)=cos(2θ)=2cos2θ1,
U1(cosθ)sinθ=sin(2θ)=2cosθsinθ,

que son respectivamente un polinomio en cosθ e un polinomio en cosθ multiplicado por sinθ. De aí

T2(x)=2x21 e
U1(x)=2x.

Unha propiedade importante e conveniente dos Modelo:Math é que son ortogonais con respecto ao produto interno

f,g=11f(x)g(x)dx1x2,

e os Modelo:Math son ortogonais con respecto a outro produto interno análogo, que se indica máis abaixo.

Os polinomios de Chebyshev Modelo:Math son polinomios co maior coeficiente posible cuxo valor absoluto no intervalo [−1, 1] está limitado por 1. Tamén son os polinomios "extremos" de moitas outras propiedades.[1]

En 1952, Cornelius Lanczos demostrou que os polinomios de Chebyshev son importantes na teoría da aproximación para a solución de sistemas lineares;[2] as raíces de Modelo:Math, que tamén se denominan nós de Chebyshev, úsanse como puntos de correspondencia para optimizar a interpolación polinómica. O polinomio de interpolación resultante minimiza o problema do fenómeno de Runge e proporciona unha aproximación que se aproxima á mellor aproximación polinómica a unha función continua baixo a norma máxima, tamén chamada criterio minimax. Esta aproximación leva directamente ao método da cuadratura de Clenshaw-Curtis.

Estes polinomios recibiron o nome de Pafnuty Chebyshev [3] A letra Modelo:Mvar úsase polas transliteracións alternativas do nome Chebyshev como Modelo:Lang, Modelo:Lang(francés) ou Modelo:Lang (alemán).

Definicións

Definición de recorrencia

Gráfico dos primeiros cinco polinomios Modelo:Mvar de Chebyshev (primeiro tipo)

Os polinomios de Chebyshev do primeiro tipo obtéñense da relación de recorrencia:T0(x)=1T1(x)=xTn+1(x)=2xTn(x)Tn1(x). A recorrencia permite tamén unha representación explícita como o determinante dunha matriz tridiagonal de tamaño k×k : Tk(x)=det[x10012x1012x010012x]

A función xeradora ordinaria ven dada por n=0Tn(x)tn=1tx12tx+t2. A función xeradora exponencial é: n=0Tn(x)tnn!=12(et(xx21)+et(x+x21))=etxcosh(tx21).

Gráfico dos primeiros cinco polinomios Modelo:Mvar de Chebyshev (segundo tipo)

Os polinomios de Chebyshev do segundo tipo están definidos por unha relación de recorrencia:U0(x)=1U1(x)=2xUn+1(x)=2xUn(x)Un1(x). Observemos que os dous conxuntos de relacións de recorrencia son idénticos, agás T1(x)=x versus Modelo:Nowrap

A función xeradora ordinaria para Modelo:Mvar é: n=0Un(x)tn=112tx+t2, e a función xeradora exponencial é: n=0Un(x)tnn!=etx(cosh(tx21)+xx21sinh(tx21)).

Definición trigonométrica

Como se describe na introdución, os polinomios de Chebyshev do primeiro tipo poden definirse como os polinomios únicos que satisfán: Tn(cosθ)=cos(nθ) ou, noutras palabras, como os únicos polinomios que satisfán: Tn(x)={cos(narccosx) if |x|1cosh(narcoshx) if x1(1)ncosh(narcosh(x)) if x1

Os polinomios do segundo tipo satisfán: Un1(cosθ)sinθ=sin(nθ), que ven sendo Un(cosθ)=sin((n+1)θ)sinθ,para Modelo:Math.

Os polinomios de Chebyshev pódense definir desta forma cando se estudan polinomios trigonométricos.[4]. Dita forma ten unha estrutura bastante semellante ao núcleo de Dirichlet Dn(x)=sin((2n+1)x2)sinx2=U2n(cosx2). (O núcleo de Dirichlet, de feito, coincide co que agora se coñece como o polinomio de Chebyshev do cuarto tipo).

Unha forma equivalente de ver isto é mediante a exponenciación dun número complexo, dado un número complexo Modelo:Math con valor absoluto de 1: zn=Tn(a)+ibUn1(a).

Que Modelo:Math é un polinomio de grao Modelo:Mvar en Modelo:Math pódese ver observando que Modelo:Math é a parte real do lado dereito da fórmula de Moivre:cosnθ+isinnθ=(cosθ+isinθ)n.A parte real é un polinomio en Modelo:Math e Modelo:Math, no que todas as potencias de Modelo:Math son pares e así substituíbles mediante a identidade Modelo:Math. Polo mesmo razoamento, Modelo:Math é a parte imaxinaria do polinomio, na que todas as potencias de Modelo:Math son impares e, polo tanto, se se factoriza un factor de Modelo:Math, os restantes factores pódense substituír para crear un polinomio de grao Modelo:Math en Modelo:Math.

Definición mediante a ecuación de Pell

Podemos definir os polinomios de Chebyshev como as solucións da ecuación de Pell : Tn(x)2(x21)Un1(x)2=1 nun anel Modelo:Math[5]. Así, pódense xerar mediante a técnica estándar para as ecuacións de Pell de toma de potencias dunha solución fundamental:Tn(x)+Un1(x)x21=(x+x21)n.

Relacións entre os dous tipos de polinomios de Chebyshev

Os polinomios de Chebyshev do primeiro e do segundo tipo corresponden a un par complementario de secuencias de Lucas Modelo:Math e Modelo:Math cos parámetros Modelo:Math e Modelo:Math: U~n(2x,1)=Un1(x),V~n(2x,1)=2Tn(x). Así conseguimos tamén un par de ecuacións de recorrencia mutua: [6] Tn+1(x)=xTn(x)(1x2)Un1(x),Un+1(x)=xUn(x)+Tn+1(x).

Podemos reorganizar o segundo usando a definición de recorrencia para os polinomios de Chebyshev do segundo tipo para dar: Tn(x)=12(Un(x)Un2(x)).

Usando iterativamente esta ecuación dá a fórmula da suma:Un(x)={2 odd jnTj(x) for odd n.2 even jnTj(x)+1 for even n,aquí substituímos Un(x) e Un2(x) utilizando a fórmula da derivada para Tn(x) dá a relación de recorrencia para a derivada de Tn: 2Tn(x)=1n+1ddxTn+1(x)1n1ddxTn1(x),n=2,3, Os polinomios de Chebyshev cumpren as desigualdades de Turán coa seguinte expresión: [7] Tn(x)2Tn1(x)Tn+1(x)=1x2>0 for 1<x<1 and Un(x)2Un1(x)Un+1(x)=1>0. As relacións con integrais son Modelo:R Modelo:Sfn 11Tn(y)yxdy1y2=πUn1(x),11Un1(y)yx1y2dy=πTn(x) onde as integrais son consideradas como valor principal.

Expresións explícitas

Diferentes enfoques para definir os polinomios de Chebyshev conducen a diferentes expresións explícitas. A definición trigonométrica dá unha fórmula explícita como segue:Tn(x)={cos(narccosx) para 1x1cosh(narcoshx) para 1x(1)ncosh(narcosh(x)) para x1A partir desta forma trigonométrica, podemos calcular a recorrencias con condicións iniciais:T0(cosθ)=cos(0θ)=1T1(cosθ)=cosθ,e tendo en conta que a identidade trigonométrica de produto a suma cumpre:2cosnθcosθ=cos[(n+1)θ]+cos[(n1)θ].Usando a definición de exponenciación de números complexos dos polinomios de Chebyshev, pódese obter as seguintes expresións equivalentes:Tn(x)=12((xx21)n+(x+x21)n) para x Tn(x)=12((xx21)n+(xx21)n) para x As dúas son equivalentes porque (x+x21)(xx21)=1.

Tamén temos unha expresión con monomios Modelo:Math que se obtén da fórmula de Moivre:Tn(cos(θ))=Re(cosnθ+isinnθ)=Re((cosθ+isinθ)n),onde Modelo:Math denota a parte real. Se expandimos a fórmula, obtemos: (cosθ+isinθ)n=j=0n(nj)ijsinjθcosnjθ. A parte real obtense cos sumandos dos índices pares. Decatándonos de que i2j=(1)j e sin2jθ=(1cos2θ)j, obtemos a fórmula explícita: cosnθ=j=0n/2(n2j)(cos2θ1)jcosn2jθ, que representa: Tn(x)=j=0n/2(n2j)(x21)jxn2j. E á súa vez pode ser escrita como unha función hiperxeométrica Modelo:Math: Tn(x)=k=0n2(n2k)(x21)kxn2k=xnk=0n2(n2k)(1x2)k=n2k=0n2(1)k(nk1)!k!(n2k)!(2x)n2k para n>0=nk=0n(2)k(n+k1)!(nk)!(2k)!(1x)k para n>0=2F1(n,n;12;12(1x)) e os seus inversos:[8][9]

xn=21nj=0jn(mod2)n(nnj2)Tj(x), onde a vírgual da suma indica que para Modelo:Math toma a metade do valor.

Relacionado con isto Modelo:Math como suma de monomios con coeficientes binomiais e potencias de 2 é Tn(x)=m=0n2(1)m((nmm)+(nm1n2m))2n2m1xn2m.

Modelo:Math tamén pode expresarse como unha función hiperxeométrica: Un(x)=(x+x21)n+1(xx21)n+12x21=k=0n/2(n+12k+1)(x21)kxn2k=xnk=0n/2(n+12k+1)(1x2)k=k=0n/2(2k(n+1)k)(2x)n2k para n>0=k=0n/2(1)k(nkk)(2x)n2k para n>0=k=0n(2)k(n+k+1)!(nk)!(2k+1)!(1x)k para n>0=(n+1) 2F1(n,n+2;32;12(1x)).

Propiedades

Simetría

Tn(x)=(1)nTn(x)={Tn(x) para n parTn(x) para n imparUn(x)=(1)nUn(x)={Un(x) para n parUn(x) para n impar

Raíces e extremos

Un polinomio de Chebyshev de calquera tipo con grao Modelo:Mvar ten Modelo:Mvar raíces simples diferentes, chamadas raíces de Chebyshev, no intervalo [ −1, 1]. Usando a definición trigonométrica e o feito de que:cos((2k+1)π2)=0pódese demostrar que as raíces de Modelo:Mvar son: xk=cos(π(k+1/2)n),k=0,,n1. As raíces de Modelo:Mvar son: xk=cos(kn+1π),k=1,,n. os máximos e mínimos de Modelo:Mvar no intervalo Modelo:Math son: xk=cos(knπ),k=0,,n. Unha propiedade única dos polinomios de Chebyshev do primeiro tipo é que no intervalo Modelo:Math todos os extremos teñen valores que son −1 ou 1. Tanto o primeiro como o segundo tipo de polinomio de Chebyshev teñen extremos nos puntos finais, dados por:Tn(1)=1Tn(1)=(1)nUn(1)=n+1Un(1)=(1)n(n+1).Os extremos de Tn(x) no intervalo 1x1 están situados en n+1 valores de x. Son ±1, ou cos(2πkd) onde d>2, d|2n, 0<k<d/2 e (k,d)=1, é dicir, k e d son números coprimos.[10][11]

Diferenciación e integración

Ao diferenciar os polinomios nas súas formas trigonométricas, pódese demostrar que:dTndx=nUn1dUndx=(n+1)Tn+1xUnx21d2Tndx2=nnTnxUn1x21=n(n+1)TnUnx21.Para a derivada Modelo:Mvar do polinomio Modelo:Mvar nos puntos ±1 temosdpTndxp|x=±1=(±1)n+pk=0p1n2k22k+1,que é de gran utilidade na resolución numérica de problemas de valores propios.

En canto á integración, a derivada primeira de Modelo:Mvar implica que:Undx=Tn+1n+1a maiores a relación de recorrencia para os polinomios de primeiro tipo que inclúen derivadas establece que para Modelo:Math:Tndx=12(Tn+1n+1Tn1n1)=nTn+1n21xTnn1.Podemos manipular a última fórmula para expresar a integral de Modelo:Mvar só en función dos polinomios de Chebyshev do primeiro tipo:Tndx=nn21Tn+11n1T1Tn=nn21Tn+112(n1)(Tn+1+Tn1)=12(n+1)Tn+112(n1)Tn1.Ademais, temos:11Tn(x)dx={(1)n+11n2 if n10 if n=1.

Propiedades de composición e divisibilidade

As definicións trigonométricas de Modelo:Math e Modelo:Math implican a composición ou propiedades de aniñamento: Tmn(x)=Tm(Tn(x)),Umn1(x)=Um1(Tn(x))Un1(x).Para Modelo:Math a orde de composición pódese inverter, e conseguimos que as funcións polinómicas Modelo:Math sexan un semigrupo conmutativo baixo composición.

Dado que Modelo:Math é divisible por Modelo:Mvar se Modelo:Mvar é impar, deducimos que Modelo:Math é divisible por Modelo:Math se Modelo:Mvar é impar. Alén diso, Modelo:Math é divisible por Modelo:Math, e para o caso de que Modelo:Mvar sexa par, é divisible por Modelo:Math.

Ortogonalidade

Tanto Modelo:Mvar como Modelo:Mvar forman unha secuencia de polinomios ortogonais. Os polinomios do primeiro tipo Modelo:Mvar son ortogonais con respecto ao peso:11x2,no intervalo [ −1, 1 ], é dicir, temos:11Tn(x)Tm(x)dx1x2={0 se nm,π se n=m=0,π2 se n=m0.Isto pódese probar con Modelo:Math e usando a identidade da definición Modelo:Math .

Do mesmo xeito, os polinomios Modelo:Mvar son ortogonais con respecto ao peso:1x2no intervalo [ −1, 1 ], é dicir, temos:11Un(x)Um(x)1x2dx={0 se nm,π2 se n=m.(Con unha constante normalizadora Modelo:Math é a distribución do semicírculo de Wigner).

Estas propiedades de ortogonalidade danse porque os polinomios de Chebyshev resolven as ecuacións diferenciais de Chebyshev:(1x2)TnxTn+n2Tn=0,(1x2)Un3xUn+n(n+2)Un=0,que son ecuacións diferenciais de Sturm-Liouville. É unha característica xeral destas ecuacións diferenciais que hai un conxunto significativo de solucións ortonormais.

O Modelo:Mvar tamén satisfai unha condición de ortogonalidade discreta:k=0N1Ti(xk)Tj(xk)={0 se ij,N se i=j=0,N2 se i=j0,onde o número Modelo:Mvar é calquera enteiro maior que Modelo:Math,Modelo:Sfn e Modelo:Math son os Modelo:Mvar nós de Chebyshev (ver arriba) de Modelo:Math :xk=cos(π2k+12N) para k=0,1,,N1.Para os polinomios Modelo:Mvar e calquera número enteiro Modelo:Math cos mesmos nós de Chebyshev Modelo:Math, hai sumas similares:k=0N1Ui(xk)Uj(xk)(1xk2)={0 se ij,N2 se i=j,e sen a función de peso:k=0N1Ui(xk)Uj(xk)={0 se i≢j(mod2),N(1+min{i,j}) se ij(mod2).Para todo número enteiro Modelo:Math, baseado nos Modelo:Mvar ceros de Modelo:Math: yk=cos(πk+1N+1) parak=0,1,,N1,pódese obter a suma:k=0N1Ui(yk)Uj(yk)(1yk2)={0 if ij,N+12 se i=j,e de novo sen a función de peso:k=0N1Ui(yk)Uj(yk)={0 se i≢j(mod2),(min{i,j}+1)(Nmax{i,j}) se ij(mod2).

Exemplos

Primeiro tipo

Os primeiros polinomios de Chebyshev do primeiro tipo son Modelo:OEIST0(x)=1T1(x)=xT2(x)=2x21T3(x)=4x33xT4(x)=8x48x2+1T5(x)=16x520x3+5xT6(x)=32x648x4+18x21T7(x)=64x7112x5+56x37xT8(x)=128x8256x6+160x432x2+1T9(x)=256x9576x7+432x5120x3+9xT10(x)=512x101280x8+1120x6400x4+50x21T11(x)=1024x112816x9+2816x71232x5+220x311x

Segundo tipo

Os primeiros polinomios de Chebyshev do segundo tipo son Modelo:OEISU0(x)=1U1(x)=2xU2(x)=4x21U3(x)=8x34xU4(x)=16x412x2+1U5(x)=32x532x3+6xU6(x)=64x680x4+24x21U7(x)=128x7192x5+80x38xU8(x)=256x8448x6+240x440x2+1U9(x)=512x91024x7+672x5160x3+10x

Como conxunto de bases

No espazo de Sobolev apropiado, o conxunto de polinomios de Chebyshev forma unha base ortonormal, de xeito que unha función no mesmo espazo pode, en Modelo:Math, expresarse mediante a expansión:[12] f(x)=n=0anTn(x).Os polinomios de Chebyshev forman unha base ortogonal que implica que os coeficientes Modelo:Math poden determinarse facilmente mediante a aplicación dun produto interno. Esta suma chámase unha serie de Chebyshev ou unha expansión de Chebyshev .

Unha serie de Chebyshev está relacionada cunha serie coseno de Fourier mediante un cambio de variables. Así todos os teoremas, identidades, etc. que se aplican ás series de Fourier teñen unha contraparte de Chebyshev. [12] Estes atributos inclúen:

A abundancia de teoremas e identidades herdados das series de Fourier fan que os polinomios de Chebyshev sexan ferramentas importantes na análise numérica; por exemplo, son as funcións básicas de propósito xeral máis populares utilizadas no método espectral,[12] moitas veces máis favorables que as series trigonométricas debido á converxencia xeralmente máis rápida para funcións continuas (o fenómeno de Gibbs segue sendo un problema).

Exemplo 1

Vemos a expansión de Chebyshev de Modelo:Math: log(1+x)=n=0anTn(x).Podemos atopar os coeficientes Modelo:Math mediante a aplicación dun produto interno ou pola condición de ortogonalidade discreta. Para o produto interno:1+1Tm(x)log(1+x)1x2dx=n=0an1+1Tm(x)Tn(x)1x2dx,que dá:an={log2 para n=0,2(1)nn para n>0.Como alternativa, cando non se pode avaliar o produto interno da función que se está a aproximar, a condición de ortogonalidade discreta dá un resultado frecuentemente útil para obter os coeficientes aproximados:an2δ0nNk=0N1Tn(xk)log(1+xk),onde o Modelo:Mvar é a función delta de Kronecker e os Modelo:Mvar son os Modelo:Mvar ceros de Gauss–Chebyshev de Modelo:Math :xk=cos(π(k+12)N).Así podemos calcular os coeficientes aproximados Modelo:Mvar de forma moi eficiente mediante a transformada do coseno discreta:an2δ0nNk=0N1cos(nπ(k+12)N)log(1+xk).

Exemplo 2

Para poñer outro exemplo:(1x2)α=1πΓ(12+α)Γ(α+1)+212αn=0(1)n(2ααn)T2n(x)=22αn=0(1)n(2α+1αn)U2n(x).

Sumas parciais

As sumas parciais de:f(x)=n=0anTn(x)son moi útiles na aproximación de varias funcións e na resolución de ecuacións diferenciais (ver método espectral). Dous métodos comúns para determinar os coeficientes Modelo:Mvar son o uso do produto interno como no método de Galerkin e o uso da colocación que está relacionada coa interpolación.

Como interpolante, os Modelo:Mvar coeficientes das Modelo:Math primeiras sumas parciales adoitan obterse nos puntos Chebyshev–Gauss–Lobatto[13] (ou retícula de Lobatto), o que resulta nun erro mínimo e evita o fenómeno de Runge asociado a unha retícula uniforme. Esta colección de puntos ven dada por: xk=cos(kπN1);k=0,1,,N1.

Familias de polinomios relacionados cos polinomios de Chebyshev

Moitas veces son usados os polinomios denotados como Cn(x) e Sn(x) relacionados cos de Chebyshev do seguinte modo: Cn(x)=2Tn(x2),Sn(x)=Un(x2)Modelo:OEIS e satisfai:Cn(x)=Sn(x)Sn2(x).AF Horadam chamou aos polinomios Cn(x) Polinomios de Vieta–Lucas con nomenclatura vn(x) e chamou aos Sn(x) Polinomios de Vieta–Fibonacci con nomenclatura Modelo:Nowrap As listas dos dous conxuntos de polinomios aparecen na Opera Mathematica de Viète, capítulo IX, teoremas VI e VII.[14] Os polinomios de Vieta–Lucas e Vieta–Fibonacci do argumento real son, ata unha potencia de i e un desprazamento de índice no caso deste último, igual aos polinomios de Lucas e Fibonacci Modelo:Math e Modelo:Math de argumento imaxinario.

Notas

Modelo:Reflist

Véxase tamén

Bibliografía

Outros artigos

Ligazóns externas

Modelo:Control de autoridades