35、傅里叶分析与多项式插值：原理、应用与实现

傅里叶分析与多项式插值：原理、应用与实现

1. 傅里叶分析相关内容

1.1 离散傅里叶变换（DFT）结果解读

在对信号进行离散傅里叶变换（DFT）后，我们可以得到一系列系数。以某信号为例，其DFT结果的部分系数如下表所示：
| Index | k | Frequency | Period | Real | Imaginary |
| — | — | — | — | — | — |
| 1 | 0 | constant | | 5 | 0 |
| 2 | 1 | 6.25 | 0.16 | 0 | 0 |
| 3 | 2 | 12.5 | 0.08 | 0.5 | 0 |
| 4 | 3 | 18.75 | 0.053333 | 0 | 0.5 |
| 5 | 4 | 25 | 0.04 | 0 | 0 |
| 6 | -3 | 31.25 | 0.032 | 0 | -0.5 |
| 7 | -2 | 37.5 | 0.026667 | 0.5 | 0 |
| 8 | -1 | 43.75 | 0.022857 | 0 | 0 |

从表中可以看出，第一个系数对应信号的均值。同时，快速傅里叶变换（FFT）检测到了12.5 Hz和18.75 Hz的信号。需要注意的是，我们标记了奈奎斯特频率，该频率以下的部分值是冗余的，它们只是奈奎斯特频率以下结果的镜像。

若去除常数项，我们可以绘制DFT的实部和虚部与频率的关系图，具体代码如下：

nyquist = fs/2;
fmin = 1/T