Sigma aproximace

Animované srovnání aproximace a (zeleně - špatně viditelné). Je vidět, že při aproximaci Fourierovou řadou se v okolí diskontinuit projevuje Gibbsův jev, zatímco σ-aproximace konverguje hladce. σ-aproximace je v matematice úprava Fourierovy sumace, která značně omezuje Gibbsův jev projevující se oscilacemi v místě diskontinuit.

Součet řady s periodou T lze při použití σ-aproximace zapsat takto:

:s(\theta) = \frac{1}{2} a_0 + \sum_{k=1}^{m-1} \operatorname{sinc} \frac{k}{m} \cdot \left[a_{k} \cos \left( \frac{2 \pi k}{T} \theta \right) + b_k \sin \left( \frac{2 \pi k}{T} \theta \right) \right],

při vyjádření normalizovanou funkcí sinc

: \operatorname{sinc} x = \frac{\sin \pi x}{\pi x}.

Člen

:\operatorname{sinc} \frac{k}{m}

je Lanczosův σ-faktor, díky kterému je eliminována velká část Gibbsova jevu. Gibbsův jev není odstraněn úplně, ale použitím druhé nebo třetí mocniny výrazu jej lze ve většině extrémních případů výrazně utlumit.

Vysvětlení

Lanczosova myšlenka je utlumit Fourierovy koeficienty vysokého řádu, které způsobují lokální divergenci řady. Studuje tedy případy, kdy se derivace Fourierovy řady může lokálně výrazně měnit. +more Pro částečný součet funkce rozvinuté na Fourierovu řadu tvaru.

:f_m(x) = \sum_{k=-(m-1)}^{m-1} c_{k} {\rm e}^{\mathrm{i}kx}, definuje :\rho_m(x) = \sum_{k=0}^{+\infty} c_{m+k} {\rm e}^{\mathrm{i}kx}. Pak lze zbytek Fourierovy řady zapsat ve tvaru :R_m(x) = f(x)-f_m(x) = {\rm e}^{\mathrm{i}mx} \rho_m(x) + {\rm e}^{-\mathrm{i}mx}\rho_{-m}(x).

Lanczos si všiml, že v obecném případě má \rho_m(x) tvar hladké nosné modulované vysokou frekvencí, takže derivace zbytku je :R_m'(x) = \mathrm{i}m ({\rm e}^{\mathrm{i}mx} \rho_m(x) - {\rm e}^{-\mathrm{i}mx}\rho_{-m}(x)) + {\rm e}^{\mathrm{i}mx} \rho_m'(x) + {\rm e}^{-\mathrm{i}mx}\rho_{-m}'(x), což způsobuje, že při velkých hodnotách zbytek řady nekonverguje „dost rychle“. Proto definuje upravený diferenciální operátor: :\mathrm{D}_m f(x) = \frac{f\left(x+\frac{\pi}{m}\right)-f\left(x-\frac{\pi}{m}\right)}{\frac{2\pi}{m}}, který dobře konverguje k operátoru derivace pro velká m, což dává :\mathrm{D}_m R_m(x) = - {\rm e}^{\mathrm{i}mx}\mathrm{D}_m \rho_{m}(x) - {\rm e}^{-\mathrm{i}mx}\mathrm{D}_m \rho_{-m}(x), a funkce \rho_m(x), \rho_{-m}(x) jsou dostatečně hladké, takže hodnoty jejich derivací nemají velký vliv na chybu aproximace. +more Všimneme-li si, že :\mathrm{D}_m ({\rm e}^{\mathrm{i}kx}) = \frac{\sin (\frac{\pi k}{m})}{\frac{\pi}{m}} \times \mathrm{i}k {\rm e}^{\mathrm{i}kx}, vidíme, že použití tohoto diferenciálního operátoru odpovídá vynásobení Fourierových koeficientů faktorem σ.

Odkazy

Reference

Související články

Lanczosovo převzorkování

Literatura

Kategorie:Fourierovy řady Kategorie:Numerická matematika