博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
三角插值的 Fourier 系数推导
阅读量:5994 次
发布时间:2019-06-20

本文共 1769 字,大约阅读时间需要 5 分钟。

给定 $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

你可能感兴趣的文章
Codeforces Round #196 (Div. 2) A. Puzzles 水题
查看>>
Can only modify an image if it contains a bitmap
查看>>
【软件分析与挖掘】Multiple kernel ensemble learning for software defect prediction
查看>>
[.net 面向对象程序设计进阶] (1) 开篇
查看>>
JavaScript权威设计--JavaScript类型,值,变量(简要学习笔记三)
查看>>
Oracle一个中文汉字占用几个字节
查看>>
汇编开发环境
查看>>
Android开发环境搭建
查看>>
git reset --hard 回滚以后 以后怎么再回去?
查看>>
【转】测试趋势之我的观点
查看>>
使用poi读写Excel
查看>>
静态和动态链接
查看>>
500 OOPS: vsftpd: cannot locate user specified in 'chown_username':whoever
查看>>
解锁redis锁的正确姿势
查看>>
Kylin如何进行JDBC方式访问或者调用
查看>>
学习jQuery
查看>>
oracle spm使用1
查看>>
修改用户所在组,修改文件的所有者,明明是自己的文件什么不能解压?
查看>>
JXL组件生成报表报错(一)
查看>>
atom中vue高亮支持emmet语法
查看>>