MATLAB从零开始实现短时傅里叶变换STFT

一、基础目标

在MATLAB中从零开始实现短时傅里叶变换（STFT）是一项非常有益的工作，它能让你深入理解这种经典时频分析技术的核心思想与实现细节。下面将梳理其基本原理、具体的实现步骤、演示代码以及一些重要的注意事项。

二、短时傅里叶变换的核心原理

STFT的基本思想非常直观：假定一个非平稳信号在很短的时间间隔内是近似平稳的。通过一个沿时间轴滑动的窗函数，将长的非平稳信号分割为一系列短的信号片段（称为“帧”），然后对每一帧信号分别进行傅里叶变换，从而得到信号频率成分随时间变化的规律。

其数学定义清晰地表达了这一过程。对于连续信号x ( t ) x(t)x(t)和窗函数w ( t ) w(t)w(t)，在时间t tt处的STFT为：

X ( t , f ) = ∫ − ∞ ∞ x ( τ ) w ( τ − t ) e − j 2 π f τ d τ \begin{aligned} X(t,f)=\int_{-\infty}^{\infty}x(\tau)w(\tau-t)e^{-j2\pi f\tau}d\tau \end{aligned}X(t,f)=∫−∞∞​x(τ)w(τ−t)e−j2πfτdτ​

对于离散信号x [ n ] x[n]x[n]，相应的离散STFT为：

X [ m , k ] = ∑ n = 0 N − 1 x [ n ] w [ n − m H ] e − j 2 π k n / N X[m,k]=\sum_{n=0}^{N-1}x[n]w[n-mH]e^{-j2\pi kn/N}X[m,k]=∑n=0N−1​x[n]w[n−mH]e−j2πkn/N

其中m mm是时间帧索引，k kk是频率索引，H HH是帧移（hop size），N NN是窗长度同时也是FFT点数。

STFT的核心在于时频分辨率权衡。根据Heisenberg不确定性原理，时间分辨率和频率分辨率无法同时无限提高。窗函数的选择和参数设置直接决定了分析的特性：

• 窗长越长，频率分辨率越高，但时间分辨率越差。
• 窗长越短，时间分辨率越高，但频率分辨率越差。

三、从零实现STFT的步骤与代码

第一步：实现基础STFT函数

以下是一个从零实现的STFT函数，它不直接依赖MATLAB的stft或spectrogram函数，而是通过循环完成分帧、加窗和FFT的过程。

function[S,f,