Stała Lebesgue'a

W matematyce stałe Lebesgue'a ( zależne od zbioru węzłów i jego wielkości) dają wyobrażenie o tym, jak dobry jest interpolant funkcji (w danych węzłach) w porównaniu z najlepszym wielomianowym przybliżeniem funkcji (stopień wielomiany są ustalone). Stała Lebesgue'a dla wielomianów stopnia co najwyżej n i dla zbioru n + 1 węzłów T jest ogólnie oznaczana przez Λ n ( T ) . Stałe te zostały nazwane na cześć Henriego Lebesgue'a .

Definicja

Naprawiamy węzły interpolacji i interwał zawierający wszystkie węzły interpolacji Proces na . To definiuje odwzorowanie przestrzeni C [ a , b ]) wszystkich funkcji ciągłych na [ a , b ] do siebie. Mapa X jest liniowa i jest rzutem na podprzestrzeń Π n wielomianów stopnia n lub mniejszego.

Stała norma X . _ _ Ta definicja wymaga od nas określenia normy na C ([ a , b ]). Jednolita norma jest zwykle najwygodniejsza.

Nieruchomości

Stała Lebesgue'a ogranicza błąd interpolacji: niech p oznacza najlepsze przybliżenie f spośród wielomianów stopnia n lub mniejszego. Innymi słowy, p minimalizuje || p fa || wśród wszystkich p w Π n . Następnie

Udowodnimy tutaj to stwierdzenie z maksymalną normą.

przez nierówność trójkąta . Ale X jest rzutem na Π n , więc

p - X ( fa ) = X ( p ) - X ( fa ) = X ( p - fa ) .

To kończy dowód, ponieważ . Zauważ, że ta relacja jest również szczególnym przypadkiem lematu Lebesgue'a .

Innymi słowy, wielomian interpolacyjny jest co najwyżej o czynnik Λ n ( T ) + 1 gorszy od najlepszego możliwego przybliżenia. Sugeruje to, że szukamy zbioru węzłów interpolacji z małą stałą Lebesgue'a.

Stałą Lebesgue'a można wyrazić za pomocą wielomianów bazowych Lagrange'a :

W rzeczywistości mamy funkcję Lebesgue'a

a stała Lebesgue'a (lub liczba Lebesgue'a) dla siatki jest jej wartością maksymalną

Niemniej jednak nie jest łatwo znaleźć jednoznaczne wyrażenie dla Λ n ( T ) .

Minimalne stałe Lebesgue'a

W przypadku równoodległych węzłów stała Lebesgue'a rośnie wykładniczo . Dokładniej, mamy następujące oszacowanie asymptotyczne

Z drugiej strony stała Lebesgue'a rośnie tylko logarytmicznie, jeśli używane są węzły Czebyszewa , ponieważ mamy

Ponownie dochodzimy do wniosku, że węzły Czebyszewa są bardzo dobrym wyborem do interpolacji wielomianowej. Istnieje jednak łatwa (liniowa) transformacja węzłów Czebyszewa, która daje lepszą stałą Lebesgue'a. Niech t i oznacza i -ty węzeł Czebyszewa. Następnie zdefiniuj

Dla takich węzłów:

Węzły te nie są jednak optymalne (tj. nie minimalizują stałych Lebesgue'a), a poszukiwanie optymalnego zbioru węzłów (który już przy pewnych założeniach okazał się unikalny) jest nadal intrygującym tematem współczesnej matematyki. Jednak ten zestaw węzłów jest optymalny do interpolacji po zbiorze n razy różniczkowalnych funkcji, których n -te pochodne do są ograniczone w wartościach bezwzględnych przez stałą M jak pokazuje NS Hoang. Za pomocą komputera można przybliżyć wartości minimalnych stałych Lebesgue'a, tutaj dla przedziału kanonicznego [−1, 1] :

N 1 2 3 4 5 6 7 8 9
Λ n ( T ) 1.0000 1.2500 1.4229 1,5595 1,6722 1.7681 1.8516 1,9255 1.9917

Istnieje niezliczona liczba nieskończenie wielu zbiorów węzłów w [-1,1], które minimalizują, dla ustalonego n > 1, stałą Lebesgue'a. Chociaż jeśli założymy, że zawsze przyjmujemy −1 i 1 jako węzły do ​​interpolacji (co nazywa się kanoniczną konfiguracją węzłów), to taki zbiór jest unikalny i zero-symetryczny. Aby zilustrować tę właściwość, zobaczymy, co się stanie, gdy n = 2 (tj. rozważymy 3 węzły interpolacji, w którym to przypadku właściwość nie jest trywialna). Można sprawdzić, że każdy zbiór (zerowo-symetrycznych) węzłów typu (− a , 0, a ) jest optymalny, gdy 8 / 3 a ≤ 1 (rozważamy tylko węzły w [−1, 1]). Jeśli wymusimy, aby zbiór węzłów był typu (−1, b , 1) , to b musi być równe 0 (spójrz na funkcję Lebesgue'a, której maksimum jest stałą Lebesgue'a). Wszystkie dowolne (tj. zero-symetryczne lub zero-asymetryczne) optymalne zbiory węzłów w [−1,1], gdy n = 2 zostały określone przez F. Schurera, a alternatywnie przez H.-J. Rack i R. Vajda (2014).

Jeśli założymy, że przyjmujemy −1 i 1 jako węzły do ​​interpolacji, to jak pokazuje H.-J. Rack (1984 i 2013), dla przypadku n = 3, znane są jawne wartości optymalnych (unikatowych i zero-symetrycznych) 4 węzłów interpolacji oraz jawna wartość minimalnej stałej Lebesgue'a. Wszystkie arbitralne zestawy optymalne 4 węzłów interpolacji w [1,1], gdy n = 3 zostały wyraźnie określone, na dwa różne, ale równoważne sposoby, przez H.-J. Rack i R. Vajda (2015).

Punkty Padewskie zapewniają kolejny zestaw węzłów o powolnym wzroście (choć nie tak powolnym jak węzły Czebyszewa) iz dodatkową właściwością bycia jednorozpuszczalnym zbiorem punktów .

Czułość wartości wielomianu

Stałe Lebesgue'a pojawiają się również w innym problemie. Niech p ( x ) będzie wielomianem stopnia n wyrażonym w postaci Lagrange'a związanej z punktami wektora t (tzn. wektor u jego współczynników jest wektorem zawierającym wartości ). niech będzie wielomianem uzyskanym przez nieznaczną zmianę współczynników u pierwotnego wielomianu p ( x ) na . Rozważ nierówność:

błąd wartości nie będzie wyższy niż odpowiednia stała Lebesgue'a współczynników. W tym sensie stałą Lebesgue'a można postrzegać jako względny numer warunku operatora odwzorowującego każdy wektor współczynników u na zbiór wartości wielomianu o współczynnikach u w postaci Lagrange'a. Właściwie możemy zdefiniować taki operator dla każdej bazy wielomianu, ale jego numer warunku jest większy niż optymalna stała Lebesgue'a dla najwygodniejszych baz.