Interpolacja trygonometryczna – metoda przybliżania funkcji za pomocą wielomianu trygonometrycznego (szeregu Fouriera). Taka interpolacja daje szczególnie dobre rezultaty przy przybliżaniu funkcji okresowych, gdyż metody używające klasycznych wielomianów, pozbawionych okresowości, powodują duże błędy interpolacji.
Przypadkiem szczególnym jest sytuacja, gdy punkty węzłowe są równoodległe. W takim przypadku najlepszym rozwiązaniem jest dyskretna transformata Fouriera.
Metoda ogólna
Opracowano na podstawie materiału źródłowego.
Założeniem każdej interpolacji jest spełnienie warunków: gdzie:
Wtedy:
Dla nieparzystej ilości punktów węzłowych:
Dla parzystej ilości punktów węzłowych:
Dla obu powyższych przypadków:
Przykład
Punkty węzłowe z przykładu i funkcja interpolująca przez nie przechodząca
Dokonać interpolacji punktów za pomocą wielomianu trygonometrycznego:
Rozwiązanie
Ilość punktów interpolowanych: (parzyste)
Stopień:
Odpowiedź
Wielomian zespolony
Problem staje się bardziej naturalny jeśli sformujemy go w dziedzinie zespolonej. Możemy wtedy zapisać zależność na wielomian trygonometryczny w postaci:
Redukuje to problem interpolacji trygonometrycznej do interpolacji wielomianowej na okręgu jednostkowym. Dowód i jednoznaczność interpolacji trygonometrycznej staje się więc wtedy równoważnym odpowiednim założeniom dla interpolacji wielomianowej.
![{\displaystyle \Theta (x)={\frac {A_{0}}{2}}+\sum _{k=1}^{m}[A_{k}\cdot \cos(k\cdot x)+B_{k}\cdot \sin(k\cdot x)].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11d00aafeedbd13dd900f40fa71eb334962d88d5)
![{\displaystyle \Theta (x)={\frac {A_{0}}{2}}+\sum _{k=1}^{m-1}[A_{k}\cdot \cos(k\cdot x)+B_{k}\cdot \sin(k\cdot x)]+{\frac {A_{m}}{2}}\cdot \cos(m\cdot x).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ae58e574b62eebe1242bc519efd093317cf7e5a3)
![{\displaystyle A_{j}={\frac {2}{n}}\sum _{k=0}^{n-1}[f(x_{k})\cdot \cos(j\cdot x_{k})],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ef390399bc7f9d69bb0127cb5e48a8064222f379)
![{\displaystyle B_{j}={\frac {2}{n}}\sum _{k=0}^{n-1}[f(x_{k})\cdot \sin(j\cdot x_{k})].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47dc5c92ea673e075e6412c127ebe196cb96793a)
![{\displaystyle A_{1}={\frac {2}{4}}\sum \limits _{k=0}^{3}f_{k}\cdot \cos(1\cdot x_{k})={\frac {1}{2}}\left[1\cdot \cos(0)+3\cdot \cos \left({\frac {\pi }{2}}\right)-2\cdot \cos(\pi )-1\cdot \cos \left({\frac {3}{2}}\pi \right)\right]={\frac {3}{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c5bb6408523b921114c2f2137e338827def902c1)
![{\displaystyle A_{2}={\frac {2}{4}}\sum \limits _{k=0}^{3}f_{k}\cdot \cos(2\cdot x_{k})={\frac {1}{2}}[1\cdot \cos(0)+3\cdot \cos(\pi )-2\cdot \cos(2\pi )-1\cdot \cos(3\pi )]=-{\frac {3}{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56a35fe44f83c03388cd4b3856ee1031001b2cb4)
![{\displaystyle B_{1}={\frac {2}{4}}\sum \limits _{k=0}^{n-1}f_{k}\cdot \sin(1\cdot x_{k})={\frac {1}{2}}[1\cdot 0+3\cdot 1-2\cdot 0-1\cdot (-1)]=2}](https://wikimedia.org/api/rest_v1/media/math/render/svg/592f508d46390f9f4672686c09e011b0f04de2e6)
