Przejdź do zawartości

Metody numeryczne fizyki/Interpolacja

Z Wikibooks, biblioteki wolnych podręczników.
Metody numeryczne fizyki
Metody numeryczne fizyki
Interpolacja

Licencja
Autor: Mirosław Makowiecki
Absolwent UMCS Fizyki Komputerowej Uniwersytetu Marii Curie-Skłodowskiej w Lublinie
Email: miroslaw(kropka)makowiecki(małpa)gmail(kropka)pl
Dotyczy: książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami.
Użytkownika książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami nie zwalnia z odpowiedzialności prawnoautorskiej nieprzeczytanie warunków licencjonowania.
Umowa prawna: Creative Commons: uznanie autorstwa, na tych samych warunkach, z możliwością obowiązywania dodatkowych ograniczeń.
Autor tej książki dołożył wszelką staranność, aby informacje zawarte w książce były poprawne i najwyższej jakości, jednakże nie udzielana jest żadna gwarancja, czy też rękojma. Autor nie jest odpowiedzialny za wykorzystanie informacji zawarte w książce, nawet jeśli wywołaby jakąś szkodę, straty w zyskach, zastoju w prowadzeniu firmy, przedsiębiorstwa lub spółki bądź utraty informacji, niezależnie czy autor (a nawet Wikibooks) został powiadomiony o możliwości wystąpienie szkód. Informacje zawarte w książce mogą być wykorzystane tylko na własną odpowiedzialność.


Następny rozdział: Aproksymacja.

Podręcznik: Metody numeryczne fizyki.

Zagadnienie interpolacji

[edytuj]

Załóżmy, że mamy n+1 punktów x0,x1,...,xn, co nazywamy węzłami interpolacji, w którym dla każdego xi jest przyporządkowana pewna wartość, która jest pewną wartością funkcji, czyli yi, co ogólnie możemy powiedzieć:

(1.1)

Przeprowadźmy pewną krzywą przechodzącej przez punkty (1.1), co nazywamy interpolacją, czyli znając te nasze punkty możemy policzyć wartość funkcji pomiędzy tymi punktami.

Zagadnienie interpolacji przy pomocy wielomianów

[edytuj]

Będziemy szukali zagadnienia interpolacyjnego w postaci wielomianu przechodzącego poprzez te punkty, co ten wielomian jest stopnia n-tego, który to piszemy:

(1.2)

Zakładamy, że węzły znajdują się w dowolnym miejscu w przedziale <a,b>. Jeśli wykorzystamy wzór (1.1), wtedy możemy napisać n+1 równań, które to zapisane są w postaci układów równań, z którego będziemy wyznaczać parametry a0,a1,...,an.

(1.3)

Możemy wyznaczyć wyznacznik główny układu równań (1.3), który jest wyznacznikiem Vandermonde'a pisanych w postaci:

(1.4)

Układ równań (1.3) nazywamy układem Cramera, których to współczynniki ai wyliczamy ze wzoru Cramera w następującej postaci:

(1.5)

Z twierdzenia Cramera istnieje taki wielomian (1.2) spełniający warunki (1.1).

Interpolacja według Lagrange'a

[edytuj]

Jeśli wzór na ai (1.5) podstawimy do wzoru na wielomian n-tego stopnia, w których poszczególne wielomiany bazowe są co najwyżej n-tego stopnia (1.2), otrzymujemy wtedy wzór na wielomian Wn(x):

(1.6)
  • gdzie wielomiany Φ0, Φ1,...,Φn są wielomianami co najwyżej n-tego stopnia. Jeśli założymy, że węzły są rozłożone w dowolny sposób w danym ściśle określonym przedziale, to dla argumentu x=xi funkcję (1.6) możemy napisać też w postaci zwartej i przy pomocy oznaczenia yi=F(xi), gdzie ta ostatnia jest funkcją interpolowaną w punktach xi:
(1.7)

Z równania (1.7) wynika następująca własność, którą piszemy:

(1.8)

Teraz skonstruujmy wielomian Φj(x), który w punktach x0,x1,...,xj-1,xj+1,..,xn jest określony według (1.8) dla wspomnianych węzłów jako:

(1.9)

Równość (1.9), która przestawia funkcję Φj(x), w której możemy policzyć parametr λ, która zależy od węzłów przestawionych powyżej, z własności funkcji Φ(x) przestawionych w punkcie (1.8) możemy powiedzieć:

(1.10)

Z wartości λ otrzymanego z tożsamości (1.10) jako odwrotności jedynki i wyrażenia zależnego od różnic pewnych węzłów i podstawieniu jego do równania (1.9), wtedy otrzymujemy równość na funkcję Φj:

(1.11)

Zatem całkowity wielomian interpolacyjny Wn(x) przestawionego w punkcie (1.6), na podstawie wyliczonej funkcji Φj(x) napisaną w punkcie (1.11), piszemy wedle schematu:



(1.12)

Oznaczmy przez oznaczenie ω(x), który zależy od iloczynów różnic argumentu x i j-tego węzła w problemie interpolacyjnym:

(1.13)

Zatem wielomian Wn(x) (1.12) według przestawienia ωn (1.13) jako wielomianu n+1-tego stopnia przepisujemy w postaci:

(1.14)
  • gdzie ω'n jest wielomian pochodzący od ωn(x) w punkcie x=xj, który jest zapisany bez użycia czynnika (x-xj), która jest miejscem zerowym funkcji ωn.

Otrzymany wzór (1.12) spełnia wszystkie warunki interpolacji i nazywamy go wielomianem interpolacyjnym Lagrange'a.

Metoda interpolacyjna Aitkena

[edytuj]

Zbudujmy teraz wielomian stopnia pierwszego w postaci wyznacznika, które w punktach xi, xj przyjmuje wartości yi i yj, którego definicja jest przestawiona dla dwóch węzłów interpolacyjnych:

(1.15)

W oparciu o wielomian pierwszego stopnia (1.15), który zbudujemy jako wielomian drugiego stopnia, który w punktach xi, xj,xk przyjmuje wartości yi, yj, yk, który jest w postaci:

(1.16)

Ogólnie rzecz można powiedzieć, że wielomian n-tego stopnia przyjmujący dla argumentów x1, x2,...,xn, określony koleino dla odpowiednich węzłów, której odpowiadają określone wartości y1, y2,...,,yn, że jest to wielomian interpolacyjny przestawionych następującym wzorem interpolacyjnym dla wspomnianych węzłów:

(1.17)

Wielomian otrzymany wychodząc od (1.15) i (1.16) piszemy jako wielomian (1.17) o stopniu n, który spełnia wszystkie warunki interpolacyjne dla węzłów interpolacyjnych oznaczonej ogólnie xi.

Oszacowanie błędu dla problemu interpretacyjnego Lagrange'a

[edytuj]

Zdefiniujmy sobie wielomian na przedziale <a,b>, który jest określony wzorem:

(1.18)

Musimy dalej określić funkcję ε(x) mając na uwadze, że pochodna f(x) ma pochodne do rzędu n+1 włącznie, zatem teraz wprowadzamy funkcje pomocniczą z pewną stałą K:

(1.19)

Z powodów oczywistych przy definicji φ(u) funkcja dla punktów x0, x1,..,xn, które są węzłami w omawianym problemie interpolacyjnym, przyjmuje wartość zerową w tychże punktach. Współczynnik K bywa tak zwykle dobierany, by funkcja φ(u) miała również pierwiastek w punkcie , zatem tożsamość (1.19), z której wyliczymy K, jest wyrażana:

(1.20)

Mianownik w równości (1.19) jest zawsze nierówny zero, dla , tzn. . Funkcja φ(u) ma n+2 miejsc zerowych w przedziale <a,b> dla punktów . Jeśli skorzystamy z twierdzenia Rolle'a możemy powiedzieć, że pierwsza pochodna funkcji (1.20) ma n+1 miejsc zerowych funkcji, każda pochodna pierwszego rzędu równa zero należy do odpowiednich przedziałów, w których każdy taki przedział posiada jedną pochodną równą zero. W każdym tak opisanej w przedziale, jeśli nań go na małe odcinki podzielimy między punktami interpolacyjnymi, w ten sposób otrzymujemy, że druga pochodna według naszego wcześniej wspomnianego wzoru przyjmuje wartość zerową w pewnym punkcie wewnątrz tego przedziału, a tych miejsc zerowych w całym omawianym przedziale jest n, wtedy dochodzimy do wniosku:

(1.21)

Jeśli będziemy kontynuować nasze rozważania dochodzimy do wniosku, ze mamy conaj mniej jeden taki punkt, który mieści się w przedziale (min(x,x0,),max(x,xn)), dla którego zachodzi , natomiast jeśli mamy , ale i też n+1 pochodna funkcji ω, którą to piszemy wzorem , zatem n+1-ta pochodna funkcji (1.16), która na tej podstawie i wcześniejszych rozważań, zapisujemy wedle schematu:

(1.22)

I jeśli weźmiemy miejsce zerowe funkcji φ(u), które jest w punkcie ξ dla (n+1) według wcześniejszych określeń jest równa zero, to wtedy z tożsamości (1.22) wynika wzór na parametr K dla ściśle określonych węzłów:

(1.23)

Jako ostatni krok należy wykorzystać równość (1.20) i wyliczoną stałą K przestawioną w punkcie (1.23), wtedy mamy następujący wzór na różnicę funkcji, z którego chcemy napisać interpolację Wn(x) i funkcji, którą chcemy interpolować, tzn.: f(x), wtedy możemy napisać tożsamość używając tutaj zamiast symbol x:

(1.24)

Przyjmijmy teraz jako maksimum funkcji f(n+1) w przedziale <a,b>, jeśli w ogóle w tym przedziale istnieje jego maksimum, które oznaczmy je przez Mn+1, tą wartość funkcji oznaczamy jako:

(1.25)

Bezwzględny błąd opisywania funkcji f(x) w postaci wielomianów Wn(x) jest to wzór napisanych na podstawie (1.25) i (1.24), który piszemy:

(1.26)

Problem najlepszego wyboru węzłów interpolacji

[edytuj]

Problem najlepszego wyboru wezłów interpolacji po raz pierwszy przestawił Czybyszew, który to wielomiany Czybyszewa przestawiamy wedle wzoru na przedziale <-1,1>:

(1.27)

Jak sie okazuje, że wielomiany Czybyszewa dla pierwszych jego wyrazów przestawiamy:

....

A jego kolejne wyrazy są przestawione wzorem na wielomian Czybyszewa, który to związek iteracyjny przestawiony jest w zależności od liczby naturalnej n jako wskaźnika i od argumentu x:

(1.28)

Wielomian Czybyszewa (1.27) jak można wykazać ma pierwiastki w punktach xn, które to są zależne od liczby naturalnej m i n, które to pierwiastki po podstawieniu wzoru poniższego do (1.27) daje nam wartość zerową funkcji Tn(x), zatem:

(1.29)

A dowodem tego, że liczby (1.29) są pierwiastkami równania (1.27) jest dokonanie do niego podstawienia wspomnianych liczb, tzn. xn, co po dokonaniu obliczeń otrzymujemy zero:

(1.30)

Parząc na kolejne wyrazy przedstawione w punkcie (1.28), to one mają przy najwyższej potędze wartość 2n-1, zatem szukamy wielomianu, który przy najwyższej potędze będzie miał wartość jeden, zatem na podstawie tego możemy napisać wielomian, który to (1.27) ma pierwiastki w punktach (1.29), zatem na podstawie tego możemy napisać równość po podzieleniu wielomianu Tn(x) przez liczbę 2n, wtedy:

(1.31)

zatem wielomian na podstawie obliczeń przeprowadzonych w punkcie (1.31), który to w punktach (1.29) ma pierwiastki, możemy przestawić według następującego sposobu:

(1.32)

co wówczas na podstawie, że wielomian Czybyszewa (1.27), który ma ograniczoną wartość leżącą między wartościami <-1,1>, co możemy napisać:

(1.33)

Ograniczenia przedziału błędów między szukanym wielomianem Wn(x), a funkcją interpolowaną f(x) (1.26), możemy przedstawić wykorzystując fakt ograniczenia wielomianu ωn(x), zatem ten nasz błąd bezwzględny przy liczeniu wartości funkcji przy pomocy funkcji interpolującej względem samej funkcji interpolowanej możemy przestawić według:

(1.34)

Wzór interpolacyjny Newtona dla dowolnych odstępu argumentu

[edytuj]

Określmy teraz funkcję f(x), określone dla tablicy x0, x1,..,xn, które są węzłami interpolacji, a natomiast punkty odpowiadająca tym węzłom są to wartości f(x0), f(x1), f(x2),.. f(xn), zatem zdefiniujmy teraz ilorazy różnicowe:

(1.35)

Iloraz różnicowy drugiego rzędu, który to piszemy przy pomocy ilorazów różnicowych pierwszego rzędu, wygląda:

(1.36)

Dla dowolnego ilorazu różnicowego n-tego rzędu, który to piszemy przy pomocy ilorazów różnicowych n-1 rzędu, mamy:

(1.37)

Napiszmy teraz wielomian interpolacyjny Wn, który to spełnia następujący warunek interpolacyjny poniżej, czyli wartość funkcji f w punkcie xi jest równa funkcji, której jest funkcją interpolacyjną Wn, też w tym samym punkcie, którą to piszemy wedle schematu ogólnie:

(1.38)

Zapiszmy teraz wielomian interpolacyjny Wn(x) rozpisując ją jako tożsamość, która jest sumą różnic wielomianu Wi(x) i Wi-1(x), co do całego wielomianu dodamy funkcję W0(x), wtedy taki przebieg rozważań sugeruje, że ona jest funkcją Wn(x), którą to możemy zapisać do postaci:

(1.39)

Oczywiste jest, że wielomian Wk-Wk-1(x) jest różnicą rozważanych wielomianów, który jest iloczynem stałej A przez iloczyn różnić argumentu x i argumentu xi dla całej gamy węzłów od i=0 do i=k-1, wtedy:

(1.40)
  • gdzie A jest pewna stałą, aby je obliczyć należy wykorzystując wzór (1.38) i podstawić xk do (1.40) wykorzystując fakt (1.14) na wielomian interpolacyjny Lagrange'a, wtedy na podstawie tego możemy powiedzieć:
(1.41)

Wzór (1.41) z oczywistych powodów na podstawie tożsamości na wielomian ωk, według wzoru (1.13), przyjmuje postać:

(1.42)

Na podstawie obliczeń przeprowadzonych w punkcie (1.42) różnica wielomianów dla k i k-1, czyli według wzoru (1.40), jest zapisana dla dowolnego argumentu x tejże funkcji:

(1.43)

Jeśli chcemy zapisać (1.37) w postaci dogodnej używając ilorazów różnicowych, to iloraz różnicowy rzędu n przedstawiamy:



(1.44)

Wzór (1.34) możemy udowodnić metodą indukcji zupełnej, którego dla n=1 ten etap dowodu piszemy w postaci poniżej, jak się przekonamy, że wspomniany wzór pokrywa się ze wzorem powyżej dla przypadku n=1, zatem na podstawionych tych wypowiedzeń mamy:

(1.45)

Teraz dokonajmy następnej części dowodu i sprawdźmy, czy z części dowodu dla n=k-1 wynika dowód n=k, wtedy:








(1.46)

W przeprowadzonych obliczeń w punkcie (1.46) możemy uzyskać taką prawidłowość, że wyrazy stojące przy f(xi) oraz f(xi+k) występują tylko jeden raz, Natomiast składniki zawierające czynnik f(xi) dla j=i+1,i+2,..,i+k-1 występują parami, zatem rozpisując te pary, możemy powiedzieć:





(1.47)

Na podstawie obliczeń przeprowadzonych w punkcie (1.47) i (1.46) udowodniliśmy za pomocą indukcji matematycznej, że wyrazy pod sumą w (1.43) możemy przestawić je jako ilorazy różnicowe do następującej postaci:

(1.48)

Przy dalszych rozważaniach skorzystamy ze wzoru (1.43), i (1.48) dla przypadku i=0 , a także z W0(x0)=f(x0), co według przedstawienia (1.14) można zapisać wielomian Wn(x):

(1.49)

Wzór interpolacyjny (1.49) możemy zapisać w postaci interpolacyjnej, który ma nazwę wzoru interpolacyjnego Netwona dla dowolnych odstępów dla węzłów lub wzoru interpolacyjnego Newtona z ilorazem różnicowym.

Różnice progresywne i wsteczne

[edytuj]

Wezły w sieci w poprzednich rozważaniach mogły być całkowicie dowolne, tzn. węzły: x0, x1, ...,xn. Tym razem wybierzemy węzły, które są napisane:

(1.50)

Różnice progresywne

[edytuj]

Niech będą określone wartości funkcji f(x) w punktach x0,x1,...,xn określone wzorami według punktu (1.50), zatem różnica progresywną pierwszego rodzaju piszemy według schematu:

(1.51)

Aby zdefiniować różnice progresywne wyższych rzędów należy określić tożsamość, która jest zdefiniowana przy pomocy kilku operatorów różnic:

(1.52)

Przykładem różnicy drugiego rzędu nazywamy różnicę przy wykorzystaniu wzoru (1.52) określoną wzorem:


(1.53)

Ogólnie rzecz biorąc możemy obrać wartości y0, y1,..,yn, wtedy n-ta różnica elementu ciągu yi określamy:

(1.54)

Co wyrażenie (1.54) możemy napisać w formie bardzo uproszczonej, używając definicji sumy, w postaci:

(1.55)

Teraz wzór (1.54) udowodnimy przy pomocy indukcji matematycznej, co dla n=1 powyższy wzór przechodzi w tożsamość:

(1.56)

Co wzór (1.54) na podstawie jej szczególnego przypadku (1.55), czyli dla n=1, czyli on się pokrywa ze wzorem (1.51). Teraz udowodnijmy wzór (1.55), ze jeśli ten nasz wzór jeśli jest prawdziwy dla n, to jest prawdziwy dla n+1, wtedy możemy powiedzieć:







(1.57)

Co kończy dowód twierdzenia (1.54) na podstawie obliczeń przeprowadzonych w punkcie (1.57) i (1.56). Następny bardzo ważnym twierdzeniem jest twierdzenie podane przy pomocy wzoru na element yk przy pomocy k-tych rzędów różnic progresywnych:

(1.58)

Aby udowodnić wzór (1.58) należy wykorzystać twierdzenie o indukcji matematycznej. Wedle twierdzenia od indukcji zupełnej sprawdźmy dla k=0, zatem dla tego nasz wzór wspomniany przechodzi w wzór, który jest zwykłą tożsamością. Sprawdźmy, czy jeśli nasz wzór (1.58) jest prawdziwy dla k, to czy z niego wynika prawdziwość dla k+1, w tym celu napiszmy do niego wyrażenie:




(1.59)

Pewne dowody ważnych twierdzeń na różnicach progresywnych

[edytuj]

Najpierw podamy pewne trzy twierdzenia bez dowodu, których to dowody są trywialne, oto twierdzenia:

(1.60)
(1.61)
(1.62)

Teraz podamy następne twierdzenia, którego dowód w cale nie jest trywialny, zatem to pierwsze twierdzenie piszemy wedle schematu:

(1.63)

A jego dowód przeprowadzamy w sposób następujący:


(1.64)

Teraz podamy następne twierdzenia, którego dowód w cale nie jest trywialny, zatem to pierwsze twierdzenie piszemy wedle następującego schematu:

(1.65)

A jego dowód, który przeprowadzamy w sposób następujący:


(1.66)

Różnice wsteczne

[edytuj]

Różnice wsteczne przyjęło się oznaczać jako pierwszą różnice progresywną w punkcie xi-1, którego zapis tejże rozważanej różnicy piszemy:

(1.67)

Katą różnicę wsteczną przyjęło się rozważać jaką różnicę według (1.67), którą piszemy:

(1.68)

Twierdzenie, które będziemy podawać i udowodniać poniżej jest to twierdzenie:

(1.69)

Twierdzenie (1.69) udowadniamy przy pomocy indukcji zupełnej, zatem dla k=0 wspomniane twierdzenie można udowodnić, co jest trywialne, zatem nasze twierdzenie dla tego k jest prawdziwe. Załóżmy teraz, że nasze twierdzenie (1.69) jest prawdziwe dla k i udowodnimy, że jest ono słuszne dla k+1, zatem wtedy możemy powiedzieć:

(1.70)

Równanie interpolacyjne Newtona dla jednakowych różnic argumentów

[edytuj]

W punkcie (1.49) napisaliśmy wzór interpolacyjnego Newtona, w której występują pewne wielkości (1.39), dla węzłów ogólnie o różnych odstępach, które tutaj będziemy pisali ten wzór dla jednakowych ostępów dla punktów xi=x0+ih, gdzie i=0,1,2,..,n. Na podstawie tego napiszemy teraz ilorazy różnicowe, które to są dla naszego przypadku:





(1.71)

Na podstawie obliczeń przeprowadzonych w punkcie (1.71) ilorazy różnicowe o numerze "k" definiujemy jako iloraz różnicy k-tej zmiennej y0 przez iloczyn k silni i k-tej potęgi stałej:

(1.72)

Wzór (1.71) udowadniamy przy pomocy twierdzenia o indukcji zupełnej, zatem dla k=0 wzór (1.72) przyjmuje postać, którą jest wartość funkcji y0, zatem wtedy . Następnym krokiem jest stwierdzenie, że wzór (1.72) jest prawdziwy, i udowodniamy tą równość dla k+1, zatem do dzieła.

(1.73)

Jeśli podstawimy wzór (1.72), który jest ilorazem różnicowym, do równości (1.49). który jest ilorazem różnicowym, w ten sposób otrzymujemy następującą tożsamość na wzór interpolacyjny Newtona wykorzystując definicję różnicy dowolnego rzędu i wiedząc, że różnica między dowolnymi a sąsiednimi wyrazami jest równa h:

(1.74)

Wprowadźmy teraz nową zmienną q, którą to definiujemy jako iloraz różnicy x i x0 przez liczbę h, który jest odległością między punktami sąsiednimi, to zmienną q piszemy:

(1.75)

Ale też zachodzą następujące tożsamości (1.75), które są ilorazami różnicy zmiennej x i węzła o numerze i-tym przez stałą odległość pomiędzy punktami, które to piszemy:




(1.76)

Wedle obliczeń przeprowadzonych w punkcie (1.76), który to wynika z definicji zmiennej q zdefiniowanej w punkcie (1.75), wtedy możemy napisać tożsamość (1.74):


(1.77)

Jeśli znów skorzystamy ze wzoru (1.74) i wzoru na różnice wsteczna w zależności od różnic progresywnych (1.69) i zmieniając numeracje z xi na xi-n bo interesuje nas zachowanie funkcji w końcowej części przedziału, tzn. że na pewno numerujemy argumenty w taki sposób by ich numeracja była od tyłu, możemy powiedzieć:


(1.78)

Jeśli przyjmować będziemy xi=x0+ih, dla i=0,-1,-2,..-n, a także przyjmować będziemy, że zmienna q jest ilorazem różnicy węzła zerowego i zmiennej x przez stałą odległość pomiędzy punktami:

(1.79)

To wtedy wykonując podobne obliczenia jak w punkcie (1.76), to wtedy wzór (1.78), na podstawie definicji zmiennej q (1.79), piszemy w następującej postaci:

(1.80)

Wzór (1.80) nazywamy drugim wzorem interpolacyjnym Newtona na interpolację wstecz.

Interpolacja przy pomocy funkcji sklejanych

[edytuj]

Wstęp do funkcji sklejanych

[edytuj]

Określmy sobie teraz n+1 punktów x0, x1, x2,..xn dla której zachodzi:

(1.81)

Funkcję s(x) będziemy określać jako funkcję sklejaną stopnia (), jeżeli w każdym z tych przedziałów funkcja jest stopnia m, a także wartości tychże funkcji należą do wartości liczb zespolonych. Na każdym z tych przedziałów (xi,xi+1) i=0,1,2,..,n-1 określmy teraz funkcję sklejaną w postaci:

(1.82)

Ponieważ współczynników cim dla ściśle określonego jest m+1, a ponieważ tych naszych przedziałów jest n określonych według (1.81), zatem liczba wszystkich współczynników cik jest n(m+1). Jeśli teraz określimy pochodną funkcji s(x), to liczba współczynników określonej przez "i" jest ich m, bo stała ci0 w tej pierwszej pochodnej znika, a w każdym węźle xi dla i=1,2,..,n określamy według m(n-1) warunków brzegowych, zatem na podstawie tego możemy powiedzieć, że funkcja s(x) zależy od n(m+1)-m(n-1)=nm+n-mn+m=n+m zmiennych . Określmy teraz funkcję s(x) jako funkcję potęgową, dla której z oczywistych powodów zachodzi xi=x0+ih dla i=-1,0,1,..,n,n+1, określmy dodatkowo i=-3,-2,-1,n,n+1,n+2,n+3. Zatem na podstawie tego możemy zbudować pewną funkcję sklejaną, który to przepisujemy wzorami:

(1.83)
(Rys. 1.1) Splot funkcji Φ3i(x)

Wykresem ogólnym funkcji (1.83) jest funkcja podana obok na prawej jego stronie tego modułu. Teraz określmy funkcję w punktach xj-2, xj-1,xi, xi+1,xi+2 i wartości tej funkcji (1.83) w tymże punktach określamy przy pomocy tabelki poniżej.

0 1 4 1 0
0 0 0
0 0

Każdą funkcję s(x) należącej do pewnego zbioru możemy przestawić w bazie funkcji , którą to wspomniana funkcja jest w postaci kombinacji liniowej funkcji bazy, co ten zapis piszemy w postaci wzoru:

(1.84)

Interpolacja funkcji interpolowanej przy pomocy funkcji sklejanej stopnia trzeciego

[edytuj]

Obierzmy sobie teraz funkcję interpolacyjną należącej do zbioru funkcji sklejanych stopnia trzeciego, które dla węzłów sieci xi wartość funkcji i funkcji interpolującej dla tego argumentu wynosi f(xi) piszemy:

(1.85)

Jak wiemy, że funkcja trzeciego stopnia, której to przestawiamy wzorem (1.82) dla m=3 ma n+3 niezależnych parametrów, w tym wliczając (1.85), w ten sposób otrzymujemy dwa parametry swobodne, dodatkowe warunki jakie funkcja powinna spełniać są dwa warunki nałożone na pierwszą pochodną, którą to możemy określić jako:

(1.86)
(1.87)

Lub dwa warunki określone przez wzory na drugich pochodnych:

(1.88)
(1.89)

Na funkcję s(x) możemy określić też przez inne warunki, tzn. warunki określoności funkcji f(x), o okresie b-a, którą to tak składamy by za pomocą funkcji sklejanych tak by było można je sprowadzić na przedział jednowymiarowy nieskończony , tzn. warunki określoności możemy napisać dla i-tej pochodnej funkcji s(x) w postaci wzoru:

(1.90)

Określmy teraz funkcję s''(x), która w punkcie xj przyjmuje wartość stałą funkcji równą Mj:

(1.91)

Zgodnie z określeniem drugiej pochodnej funkcji s(x) w punkcie xj (1.91), to określmy drugą pochodną funkcji x dla dowolnego argumentu należącego do <xj-1,xj> w postaci:

(1.92)

Możemy z całkować funkcję (1.92) względem argumentu x, w ten sposób otrzymujemy tożsamość, którą to piszemy wzorem poniżej:

(1.93)

Zatem na ostatnim wniosku możemy otrzymamy funkcje s(x), która wynika z całkowania funkcji (1.93) w naszym wyżej określonym przedziale:

(1.94)

Na funkcję s(x) (1.94) możemy nałożyć warunki interpolacji (1.85), w ten sposób wzór (1.94) dla tych węzłów możemy napisać:

(1.95)
(1.96)

Z warunku (1.95) wyznaczamy stałą Bj, którą schematycznie piszemy:

(1.97)

Z warunku (1.96) możemy wyznaczyć wzór na stałą Aj, którą to piszemy po podstawieniu do niego wzoru na stałą Bj (1.97), otrzymujemy:

(1.98)

Piszemy teraz warunek, która mówi, że s'(x) jest funkcja ciągłą na przedziale obustronnie zamkniętym <a,b>, zatem w tym celu należy obliczyć granicę jednostronne w punktach xj+0, xj-0. Możemy policzyć te pochodne przy wykorzystaniu tożsamości na pewne stałe, czyli stałe zdefiniowane według wzoru Aj (1.98), ale najpierw przedstawmy pochodną s(x) w punkcie xj-0:


(1.99)

i pochodną funkcji s(x) w punkcie xj+0:


(1.100)

I zażądamy by funkcja sklejana s(x) była funkcją ciągłą w punkcie xj, tzn, pochodna lewostronna (1.99) tej funkcji jest równa pochodnej prawostronnej tej samej funkcji (1.100), tzn. by spełniała warunek:

(1.101)

Dalej wykorzystamy warunek (1.101), który przedstawia ciągłość funkcji s'(x) w punkcie xj, następnie do tego wzoru podstawiamy formułę na pochodną lewostronną (1.99) i prawostronną (1.100), zatem takim razie możemy otrzymać warunek, który będziemy dalej przekształcać w tym samym punkcie:


(1.102)

Wzór (1.102) mnożymy obustronnie przez , w ten sposób możemy dojść do następnego celu:

(1.103)

Wprowadźmy teraz oznaczenia, które będą nam potrzebne do dalszego kroku obliczeń:

(1.104)
(1.105)
(1.106)

Na podstawie oznaczeń (1.104), (1.105), (1.106) równanie (1.103) piszemy wedle przepisu:

(1.107)

Warunek (1.86) możemy przedstawić podstawiając w funcji s(x) (1.100) taki argument by było "x" dążące do x0 prawostronnie:


(1.108)

Ze wzoru (1.108) otrzymujemy dwie bardzo ważne tożsamości na d0, które zapisujemy w zależności od współczynników M0 i M1 dla pierwszej tożsamości i w zalezności od y0, y1 i h1 w drugiej tożsamości:

(1.109)
(1.110)

Warunek (1.87) możemy przedstawić podstawiając za funkcję s(x) (1.100) dla i=n-1 i dla "x" dążącego lewostronnie do xn:,


(1.111)

Z tożsamości (1.111) otrzymujemy dwie bardzo ważne tozsamości na dn, które zapisujemy w zależności od współczynników Mn i Mn-1 dla pierwszej tożsamości i w zależności od yn, yn-1 i hn w drugiej tożsamości:

(1.112)
(1.113)

Na podstawie tożsamości (1.107) , (1.109) i (1.112) możemy napisać równanie macierzowe jako:

(1.114)

Jeśli będziemy przyjmować funkcje okresową s(x) o okresie b-a, to wtedy na podstawie (1.90) możemy powiedzieć, że wielkość o wskaźniku n+1 jest to wielkość o wskaźniku 1, a wielkość o wskaźniku n jest to wielkość o wskaźniku 0. Zatem na podstawie tego i wcześniejszych rozważań tożsamość (1.107) możemy napisać:

(1.115)

Wzory na λn według (1.104), μn według (1.105) i według (1.106) piszemy tutaj na podstawie określoności funkcji s(x) wedle schematów:

(1.116)
(1.117)
(1.118)

Na podstawie równania (1.107) i (1.115) możemy napisać tożsamość macierzową wynikającej z przyjęcia okresowości funkcji s(x) rugując stałą M0 (M0=Mn):

(1.119)

Można obliczyć stałe Mi z warunku macierzowego (1.114) dla funkcji nieokresowej s(x) lub dla funkcji okresowej s(x) według warunku macierzowego (1.119), zatem przy pomocy obliczonych współczynników Mi możemy policzyć wartość s(x) w punkcie "x" ze wzoru (1.94).