Nejlepší polynomiální aproximace

Abstract

U některých matematických úloh nedokážeme nalézt řešení exaktně, ale pouze přibližně metodami numerické matematiky. Klíčovým konceptem je pak nahrazení neboli aproximace. Mezi nejčastěji používané metody patří například polynomiální interpolace, nicméně chyba této metody může být v některých bodech intervalu velmi vysoká. Aproximace, která v každém bodě daného intervalu minimalizuje maximum chyby, se nazývá nejlepší stejnoměrná aproximace (best approximation). Obecně existuje právě jedno řešení úlohy nalezení nejlepší polynomiální aproximace, ale nelze ji vypočítat přesně. K přibližnému řešení se lze dopracovat například použitím Remezova algoritmu. Překvapivě dobré (near-best approximation) řešení nám také nabízí interpolace v Čebyševových uzlech. Cílem naší práce je zejména nastudovat a naimplementovat Remezův algoritmus a otestovat jeho funkčnost. Remezův algoritmus nakonec porovnáme s Lagrangeovou interpolací v Čebyševových uzlech.

For some mathematical problems we are not able to find solutions exactly, but only approximately by using methods of numerical mathematics. The key concept is then replacement or in other words approximation. One of the most used methods is, for example, polynomial interpolation, however at some points of the interval the error of this method might be significant. The approximation which minimizes the maximum error at every point of the interval is called the best approximation. Generally, there is a unique solution to the problem of finding the best polynomial approximation, but it is not possible to calculate it exactly. The approximate solution can be obtained, for example, by using the Remez algorithm. The interpolation at Chebyshev nodes also yields surprisingly good results, generally being near-best approximations. The goal of our thesis is namely to study and implement the Remez algorithm and test its functionality. Finally, we compare the Remez algorithm with Lagrange interpolation at Chebyshev nodes.