給定 $k$ 個互不相同的復數 $x_0,\cdots,x_{k-1}$,以及 $k$ 個復數$y_0,\cdots,y_{k-1}$.我們知道存在唯一的復系數 $k-1$ 次多項式
$$
\mathcal{P}_{k-1}(x)=\xi_0+\xi_1x+\cdots+\xi_{k-1}x^{k-1}
$$
使得
$$
\mathcal{P}_{k-1}(x_0)=y_0,\cdots,\mathcal{P}_{k-1}(x_{k-1})=y_{k-1}.
$$
其中 $\xi_0,\cdots,\xi_{k-1}\in\mathbf{C}$.這個結論是范德蒙行列式不為0的一個簡單推論.特別的,我們令 $x_i=\omega^i$,其中 $\omega=e^{\frac{2\pi i}{k}}$,我們就得到了三角插值多項式.為了確定三角插值多項式的系數,我們使用 Cramer 法則.我們知道
$$
\begin{cases}
\xi_0+\xi_1\omega^0+\cdots+\xi_{k-1}\omega^0=y_0\\
? \xi_0+\xi_1\omega^1+\cdots+\xi_{k-1}\omega^{k-1}=y_1\\
\xi_0+\xi_1\omega^2+\cdots+\xi_{k-1}\omega^{2(k-1)}=y_2\\
\vdots\\
\xi_0+\xi_1\omega^{k-1}+\cdots+\xi_{k-1}\omega^{(k-1)(k-1)}=y_{k-1}.
\end{cases}
$$
因此
$$
\xi_i=\frac{\begin{vmatrix}
??? \omega^{0}&\omega^0&\cdots&y_{0}&\cdots&\omega^{0}\\
??? \omega^0&\omega^1&\cdots&y_{1}&\cdots&\omega^{k-1}\\
??? \vdots&\vdots&\cdots&\vdots&\cdots&\vdots\\
\omega^0&\omega^{k-1}&\cdots&y_{k-1}&\cdots&\omega^{(k-1)(k-1)}\\
? \end{vmatrix}}{\begin{vmatrix}
??? \omega^{0}&\omega^0&\cdots&\omega^{0}\\
\omega^0&\omega^1&\cdots&\omega^{k-1}\\
\vdots&\vdots&\cdots&\vdots\\
\omega^0&\omega^{k-1}&\cdots&\omega^{(k-1)(k-1)}\\
? \end{vmatrix}}.
$$
設
$$
\begin{pmatrix}
? y_0\\
y_1\\
\vdots\\
y_{k-1}\\
\end{pmatrix}=\alpha_0 \begin{pmatrix}
? \omega^{0}\\
\omega^{0}\\
\vdots\\
\omega^0\\
\end{pmatrix}+\cdots+\alpha_i \begin{pmatrix}
? \omega^{0}\\
\omega^{i}\\
\vdots\\
\omega^{i(k-1)}\\ ?
\end{pmatrix}+\cdots+\alpha_{k-1}\begin{pmatrix}
? \omega^{0}\\
\omega^{k-1}\\
\vdots\\
\omega^{(k-1)(k-1)}
\end{pmatrix},
$$
則 $\xi_i=\alpha_i$.于是我們只用求 $\alpha_i$ 即可.易得
$$
\alpha_i=\frac{1}{k}\begin{pmatrix}
? y_0\\
y_1\\
\vdots\\
y_{k-1}\\
\end{pmatrix}\cdot \begin{pmatrix}
? \omega^0\\
\omega^{-i}\\
\vdots\\
\omega^{-i(k-1)}\\
\end{pmatrix}.
$$

轉載于:https://www.cnblogs.com/yeluqing/p/3827412.html

總結

以上是生活随笔為你收集整理的三角插值的 Fourier 系数推导的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。