💻MATLAB快速傅里叶变换绘制频谱⚡️

在信号处理领域，快速傅里叶变换（FFT）是分析信号频率特性的强大工具！今天用 MATLAB 实现一个简单的 FFT 示例，帮助大家轻松绘制信号频谱✨。

首先，我们需要准备一段时域信号。假设是一个正弦波叠加了噪声的声音信号，代码如下：

```matlab

fs = 1000;% 采样频率

t = 0:1/fs:1-1/fs;% 时间向量

f = 50; % 信号频率

x = sin(2pift) + 0.5 randn(size(t));% 正弦波+噪声

```

接下来，调用 `fft` 函数进行频域转换，并绘制频谱图：

```matlab

X = fft(x); % 计算 FFT

P2 = abs(X)/length(x);% 双边频谱

P1 = P2(1:length(x)/2+1);% 单边频谱

P1(2:end-1) = 2P1(2:end-1);

f = fs(0:(length(x)/2))/length(x);% 频率轴

% 绘制频谱

plot(f, P1)

title('单边幅度谱')

xlabel('频率 (Hz)')

ylabel('幅值')

```

运行后，你会看到清晰的频谱图，50Hz 的信号跃然纸上🎉。FFT 不仅能用于音频分析，还能广泛应用于图像处理、通信等领域。快来试试吧！💬

Matlab FFT 频谱分析