Aproksymacja

Aproksymacja (łac. approximare – przybliżać) – budowanie rozwiązań przybliżonych, zwłaszcza wtedy, gdy ścisłego rozwiązania nie da się przedstawić dokładnie w postaci analitycznejSzablon:R.

Przykład to zastąpienie pewnej funkcji $f(x)$ inną, zazwyczaj prostszą $\phi (x),$ umożliwiającą efektywne rozwiązanie postawionego problemu. Przykłady takich sytuacji to:

• obliczanie całek oznaczonych z funkcji, które nie dają się scałkować ściśle;
• rozwiązywanie równań różniczkowych – zarówno zwyczajnych, jak i cząstkowych – kiedy poszukuje się niewiadomych funkcji;
• opracowywanie wyników pomiarów znanych tylko na dyskretnym zbiorze punktów, np. w meteorologii.

Aproksymacja może być dokonywana na różne sposoby i dlatego można poszukiwać aproksymacji optymalnej w ściśle określonym sensie.

Sformułowanie uproszczone

Ogólnie rzecz ujmując, aproksymacja polega na przybliżaniu pewnej funkcji $f(x),$ w obszarze $\mathrm{\Omega }$ jej określoności, za pomocą innej, prostszej funkcji przybliżającej $\phi (x),$ określonej w tym samym obszarze, której wartości zależą od pewnej liczby parametrów. Najczęściej jako funkcje $\phi (x),$ stosuje się wielomiany uogólnione w postaci Szablon:Wzór

w której funkcje ${\phi }_i(x)$ tworzą tzw. bazę aproksymacji

$𝐁(x)=[{\phi }_1(x),{\phi }_2(x),\dots {\phi }_n(x)],$

zaś $a_i$ są liczbowymi współrzędnymi funkcji $\phi (x)$ względem przyjętej bazy. Dobór tych współczynników może się odbywać na różne sposoby, przy czym powinien on być taki, aby błąd aproksymacji był jak najmniejszy.

Jednym z praktycznych sposobów budowania aproksymacji w pewnym sensie optymalnej, jest metoda minimalizacji błędu przybliżenia, określonego iloczynem skalarnym różnicy funkcji $\phi (x),f(x)$ Szablon:Wzór

przy czym ten iloczyn może być definiowany na dwa sposobySzablon:R: Szablon:Wzór

Minimalizacja tak określonego błędu wymaga, aby Szablon:Wzór

Opisany powyżej sposób aproksymacji funkcji $f(x)$ za pomocą wielomianu $\phi (x)$ polegał na sformułowaniu i wykorzystaniu konkretnych warunków minimalizacji błędu określonego wzorem (b). Warunki te przybrały postać układu równań (d), w których współczynniki przy niewiadomych $a_i$ określone zostały funkcjonałami $⟨{\phi }_k\cdot {\phi }_i⟩,⟨{\phi }_k\cdot f⟩$ ze względu na funkcje ${\phi }_k.$

Sformułowanie ogólne

Ogólne sformułowanie aproksymacji w przestrzeni liniowej $F$ wymaga określenia warunków, jakie ta aproksymacja ma spełniać. Jeżeli przez $\overline{F}$ oznaczymy podzbiór $(\overline{F}\subset F)$ zbioru $F,$ będący również przestrzenią liniową, to aproksymacja będzie polegać na tym, aby dla każdego elementu $f\in F$ znaleźć taki element $\overline{f}\in \overline{F},$ dla którego zachodzą równości

$l_k(\overline{f})=l_k(f),$ dla $k=1,2,\dots n,$

w których $l_k$ są pewnymi funkcjonałami liniowymi, określającymi warunki dokonywanej aproksymacji.

Zatem zagadnienie aproksymacji wymaga określenia trzech zbiorów:

• funkcji $F,$ funkcji aproksymowanych,
• funkcji $\overline{F},$ funkcji aproksymujących,
• ciągu $l_1,l_2,\dots l_n$ funkcjonałów liniowych.

Najczęściej jako $\overline{F},$ wybiera się zbiór wielomianów uogólnionych o postaci

$\phi (x)=\sum _{i=1}^n{a}_i{\phi }_i(x)$

utworzonej z funkcji bazowych ${\phi }_1,{\phi }_2,\dots {\phi }_n\in F.$ W tym przypadku $\overline{F}$ staje się n-wymiarową podprzestrzenią $F.$

Poszukiwanie elementu $\overline{f}\in \overline{F},$ aproksymującego $f\in F$ polega na zbudowaniu takiego wielomianu $\phi \in \overline{F}$ Szablon:Wzór

który spełnia równości Szablon:Wzór

tworzące układ równań służących do określenia współczynników kombinacji liniowej (e).

Jeżeli za ${\phi }_i(x)\in F$ przyjmiemy dowolne funkcje liniowo niezależne, to macierz układu równań najczęściej będzie bardzo pełna. W celu wygenerowania takiego układu równań, który odznaczałby się macierzą rzadką, budujemy aproksymację w $\overline{f}\in \overline{F}$ zawężoną do interpolacjiSzablon:R, w następującej postaci Szablon:Wzór

w której bazę takiej aproksymacji stanowią funkcje ${\varphi }_i(x)$ o tej własności, że $l_k{\varphi }_i={\delta }_{ki}.$ Stąd wynika Szablon:Wzór

Funkcję ${\varphi }_s(x),s=1,2,\dots n$ otrzymuje się na podstawie kombinacji (e)

${\varphi }_s(x)=\sum _{k=1}^n{a}_i^{(s)}{\phi }_s(x),$

której współczynniki $a_i^{(s)}$ są określone równaniami (f), w których funkcja $f(x)$ zostaje zastąpiona przez kolejne funkcje ${\varphi }_s(x),s=1,2,\dots n.$ Funkcje te nazywane są funkcjami bazowymi Lagrange’a.

Przykład 1

Funkcję $f(x)=x^2$ można aproksymować w przedziale $(0,1)\in R^1,$ funkcją liniową przyjmując, że

$\phi (x)=a_1{\phi }_1(x)+a_2{\phi }_2(x)=a_1+a_{2}x$

i definiując dwa funkcjonały, na przykład w postaci iloczynów skalarnych

$l_k(\phi )=⟨{\phi }_k\cdot \phi ⟩=\underset{0}{\overset{1}{\int }}{\phi }_k(x)\cdot \phi (x)dx,k=1,2.$

Warunki (f) przyjmują postać

$⟨{\phi }_k\cdot {\phi }_1⟩a_1+⟨{\phi }_k\cdot {\phi }_2⟩a_2=⟨{\phi }_k\cdot f⟩,k=1,2.$

Dla obliczenia współczynników $a_1,a_2$ otrzymujemy układ równań

$\left[\begin{array}{cc}1& \frac{1}{2}\\ \frac{1}{2}& \frac{1}{3}\end{array}\right]\left[\begin{array}{c}{a}_1\\ a_2\end{array}\right]=\left[\begin{array}{c}\frac{1}{3}\\ \frac{1}{4}\end{array}\right]\to a_1=-\frac{1}{6},a_2=1.$

Przykład 2

Aproksymację funkcji $f(x)$ określonej w przedziale $[a,b]$ można zastąpić aproksymacją funkcji $F(\xi ),(x=\frac{b+a}{2}+\frac{b-a}{2}\xi )$ w przedziale standardowym $[-1,1].$ Bazę aproksymacji zbudujemy w postaci uogólnionego wielomianu stopnia $n$

$\phi (\xi )=a_0{w}_0(\xi )+a_1{w}_1(\xi )+\dots +a_n{w}_n(\xi ),$

utworzonego z funkcji dowolnych, ale wzajemnie ortogonalnych, spełniających w przedziale standardowym, warunki

$\underset{-1}{\overset{1}{\int }}{w}_k(\xi )w_i(\xi )d\xi =0,$ gdy $i\ne k.$

Funkcjonał $l_k,$ występujący w (f) przyjmiemy w postaci

$l_k(F)=⟨w_k\cdot F⟩=\underset{-1}{\overset{1}{\int }}{w}_k(\xi )F(\xi ),k=0,1,\dots n.$

Mamy również

$\begin{array}{cc}{l}_k(\phi )& =⟨w_k\cdot \phi ⟩\\ & =\underset{-1}{\overset{1}{\int }}{w}_k(\xi )\phi (\xi )\\ & =\underset{-1}{\overset{1}{\int }}{w}_k(\xi )\sum _{i=1}^n{a}_i{w}_i(\xi )d\xi \\ & =\underset{-1}{\overset{1}{\int }}{a}_k{w}_k^2(\xi )d(\xi )\\ & =⟨w_k\cdot w_k⟩a_k,k=0,1,\dots n.\end{array}$

Układ równań (f) redukuje się do najprostszej postaci

$⟨w_k\cdot F⟩=⟨w_k\cdot w_k⟩a_k\to a_k=\frac{⟨w_k\cdot F⟩}{⟨w_k\cdot w_k⟩},k=0,1,\dots n.$

Funkcjonały $⟨w_k\cdot w_k⟩,k=0,1,2,\dots n$ mają przykładowo wartości

• $⟨w_k\cdot w_k⟩=\frac{2}{2k+1},$ gdy funkcje $w_k(x)$ są wielomianami Legendre’aSzablon:R stopnia $k,$
• $⟨w_k\cdot w_k⟩=1,$ gdy $w_0(\xi )=\frac{1}{\sqrt{2}},w_k(\xi )=\cos (k\pi \xi ),k=1,2,\dots ,n,$
• $⟨w_0\cdot w_0⟩=\pi ,⟨w_k\cdot w_k⟩=\frac{\pi }{2},k=1,2,\dots n,$ gdy funkcje $w_k(x)$ są wielomianami CzebyszewaSzablon:R stopnia $k.$

Zdefiniowanie najlepszej aproksymacji

• W przestrzeniach unormowanych

Niech dana będzie przestrzeń liniowa $X$ z normą $\Vert \cdot \Vert$ i niech $V\subset X$ będzie podprzestrzenią liniową $X$ skończonego wymiaru. Zadanie najlepszej aproksymacji polega na znalezieniu takiego $v^{*}\in V$ (elementu najlepszej aproksymacji dla danego $x\in X$), że zachodzi:

$\mathrm{\forall }v\in V\Vert x-v^{*}\Vert ⩽\Vert x-v\Vert$

Należy przez to rozumieć, że element $v^{*}$ jest elementem „najbliższym” do aproksymowanego $x$ spośród wszystkich elementów $v\in V.$

Zadanie najlepszej aproksymacji jest zawsze rozwiązywalne, tzn. dla każdego $x\in X$ istnieje element najlepszej aproksymacji $v^{*},$ ale niekoniecznie jest on jedyny. Należy zauważyć, że element najlepszej aproksymacji zależy od normy, jaka została przyjęta w przestrzeni $X.$

• W przestrzeniach unitarnych

Niech $X$ będzie przestrzenią z iloczynem skalarnym $⟨\cdot ,\cdot ⟩$ i niech norma w $X$ będzie generowana tym iloczynem: $\Vert x\Vert =\sqrt{⟨x,x⟩}.$

Wtedy dla danego $x\in X$ element najlepszej aproksymacji $v^{*}$ jest jedyny i jest określony następująco:

$\mathrm{\forall }v\in V⟨x-v^{*},v⟩=0.$

Zagadnienia aproksymacji funkcji

Aproksymację stosuje się w sytuacjach, gdy nie istnieje analityczna postać funkcji, która pozwalałaby na wyznaczenie wartości dla dowolnego z jej argumentów, a jednocześnie wartości tej nieznanej funkcji są dla pewnego zbioru jej argumentów znane. Z przypadkiem takim mamy do czynienia np. w meteorologii przy sporządzaniu map synoptycznych na podstawie wyników pomiarów terenowych.

Aproksymowanie funkcji może polegać na przybliżaniu jej za pomocą kombinacji liniowej tzw. funkcji bazowychSzablon:Odn. Od funkcji aproksymującej, przybliżającej daną funkcję nie wymaga się, aby przechodziła ona przez jakieś konkretne punkty, tak jak to ma miejsce w interpolacji. Z matematycznego punktu widzenia aproksymacja funkcji $f$ w pewnej przestrzeni Hilberta $H$ jest zagadnieniem polegającym na poszukiwaniu pewnej funkcji $g\in G,$ gdzie $G$ jest podprzestrzenią $H(G\subset H)$ takiej, by odległość (w sensie obowiązującej w $H$ normy) między $f$ a $g$ była jak najmniejsza.

Aproksymacja funkcji powoduje pojawienie się błędów, zwanych błędami aproksymacjiSzablon:Odn. Dużą zaletą aproksymacji w stosunku do interpolacji jest to, że aby dobrze przybliżać, funkcja aproksymująca nie musi być wielomianem wysokiego stopnia (w ogóle nie musi być wielomianem). Przybliżenie w tym wypadku rozumiane jest jako minimalizacja pewnej funkcji błędu. Prawdopodobnie najpopularniejszą miarą tego błędu jest średni błąd kwadratowy, ale możliwe są również inne funkcje błędu, jak choćby błąd średni.

Istnieje wiele metod aproksymacyjnych. Jednymi z najbardziej popularnych są: aproksymacja średniokwadratowa i aproksymacja jednostajna oraz aproksymacja liniowa, gdzie funkcją bazową jest funkcja liniowa.

Funkcja aproksymująca może być przedstawiona w różnej postaci. Najczęściej jest to postać:

• wielomianu (tzw. aproksymacja wielomianowa),
• funkcji sklejanych,
• funkcji matematycznych uzyskanych na drodze statystyki matematycznej (przede wszystkim regresji),
• sztucznych sieci neuronowych.

Aproksymację można formułować również przy rozwiązywaniu zagadnień dwu- i trójwymiarowych.

Zobacz też

• aproksymacja diofantyczna
• aproksymacja punktowa
• krzywa Béziera
• metoda Sheparda
• zaokrąglanie

Bibliografia

• Szablon:Cytuj książkę