Węzły Czebyszewa

Węzły Czebyszewa są równoważne współrzędnym n {\displaystyle n} na osi x {\displaystyle x} równo rozmieszczonych punktów na półokręgu jednostki (tutaj dla n = 10 {\displaystyle n=10} )[1].

W analizie numerycznej węzły Czebyszewa są specyficznymi rzeczywistymi liczbami algebraicznymi, mianowicie pierwiastkami wielomianów Czebyszewa pierwszego rodzaju. Są często używane jako węzły w interpolacji wielomianowej, ponieważ wynikowy wielomian interpolacyjny minimalizuje efekt Rungego, czyli duże oscylacje wielomianu interpolacyjnego przy krańcach przedziału[2]. Fakt, że miejsca zerowe wielomianów Czebyszewa zagęszczają się ku krańcom przedziału, pozwala lepiej związać wielomian zapobiegając naturalnym dla wielomianów wysokiego rzędu oscylacjom.

Definicja

Miejsca zerowe pierwszych 50 wielomianów Czebyszewa pierwszego rodzaju

Dla danej liczby całkowitej n {\displaystyle n} węzły Czebyszewa w przedziale ( 1 , 1 ) {\displaystyle (-1,1)} to

x k = cos ( 2 k 1 2 n π ) , k = 1 , , n . {\displaystyle x_{k}=\cos \left({\frac {2k-1}{2n}}\pi \right),\quad k=1,\dots ,n.}

Są to pierwiastki wielomianu Czebyszewa pierwszego rodzaju stopnia n . {\displaystyle n.} Dla węzłów w dowolnym przedziale [ a , b ] {\displaystyle [a,b]} można zastosować przekształcenie afiniczne:

x k = 1 2 ( a + b ) + 1 2 ( b a ) cos ( 2 k 1 2 n π ) , k = 1 , , n . {\displaystyle x_{k}={\frac {1}{2}}(a+b)+{\frac {1}{2}}(b-a)\cos \left({\frac {2k-1}{2n}}\pi \right),\quad k=1,\dots ,n.}

Aproksymacja

Węzły Czebyszewa są ważne w teorii aproksymacji, ponieważ tworzą szczególnie dobry zestaw węzłów do interpolacji wielomianowej. Biorąc pod uwagę funkcję f {\displaystyle f} na przedziale [ 1 , + 1 ] {\displaystyle [-1,+1]} i n punktów x 1 , x 2 , , x n , {\displaystyle x_{1},x_{2},\dots ,x_{n},} w tym przedziale, wielomian interpolacyjny jest unikalnym wielomianem P n 1 {\displaystyle P_{n-1}} stopnia co najwyżej n 1 , {\displaystyle n-1,} który ma wartość f ( x i ) {\displaystyle f(x_{i})} w każdym punkcie x i . {\displaystyle x_{i}.} Błąd interpolacji w x {\displaystyle x} jest równy

f ( x ) P n 1 ( x ) = f ( n ) ( ξ ) n ! i = 1 n ( x x i ) {\displaystyle f(x)-P_{n-1}(x)={\frac {f^{(n)}(\xi )}{n!}}\prod _{i=1}^{n}(x-x_{i})}

dla pewnego ξ {\displaystyle \xi } [3] Więc logiczne jest, aby próbować zminimalizować

max x [ 1 , 1 ] | i = 1 n ( x x i ) | . {\displaystyle \max _{x\in [-1,1]}\left|\prod _{i=1}^{n}(x-x_{i})\right|.}

Ten produkt jest unormowanym wielomianem stopnia n . {\displaystyle n.} Można wykazać, że maksymalna wartość bezwzględna (maksymalna norma) dowolnego takiego wielomianu jest ograniczona od dołu przez 2 1 n . {\displaystyle 2^{1-n}.} Ta granica jest osiągana przez skalowane wielomiany Czebyszewa 2 1 n T n , {\displaystyle 2^{1-n}T_{n},} które również są moniczne. Zauważmy, że | T n ( x ) | 1 {\displaystyle |T_{n}(x)|\leqslant 1} dla x [ 1 , 1 ] {\displaystyle x\in [-1,1]} [4]. Dlatego też, gdy węzły interpolacji x i {\displaystyle x_{i}} są pierwiastkami T n , {\displaystyle T_{n},} błąd spełnia:

| f ( x ) P n 1 ( x ) | 1 2 n 1 n ! max ξ [ 1 , 1 ] | f ( n ) ( ξ ) | . {\displaystyle \left|f(x)-P_{n-1}(x)\right|\leqslant {\frac {1}{2^{n-1}n!}}\max _{\xi \in [-1,1]}\left|f^{(n)}(\xi )\right|.}

Dla dowolnego przedziału [ a , b ] {\displaystyle [a,b]} zmiana zmiennej pokazuje, że:

| f ( x ) P n 1 ( x ) | 1 2 n 1 n ! ( b a 2 ) n max ξ [ a , b ] | f ( n ) ( ξ ) | . {\displaystyle \left|f(x)-P_{n-1}(x)\right|\leqslant {\frac {1}{2^{n-1}n!}}\left({\frac {b-a}{2}}\right)^{n}\max _{\xi \in [a,b]}\left|f^{(n)}(\xi )\right|.}

Zobacz też

  • Richard L.R.L. Burden Richard L.R.L., J. DouglasJ.D. Faires J. DouglasJ.D., Numerical analysis, wyd. 8, s. 503–512, ISBN 0-534-39200-8 .

Przypisy

  1. Lloyd N. Trefethen: Approximation Theory and Approximation Practice. SIAM, 2012.
  2. Kurtis D. Fink, John H. Mathews: Numerical Methods using MATLAB. Wyd. 3. Upper Saddle River, NJ: Prentice Hall, 1999, s. 236–238.
  3. Afternotes on Numerical Analysis. A series of lectures on elementary numerical analysis presented at the University of Maryland at College Park and recorded after the fact. 1996. ISBN 978-0-89871-362-6. (20.3).
  4. Afternotes on Numerical Analysis. A series of lectures on elementary numerical analysis presented at the University of Maryland at College Park and recorded after the fact. 1996. ISBN 978-0-89871-362-6., Lecture 20, § 14.