matlab傅立叶变换频谱分析

以下是使用Matlab进行傅立叶变换频谱分析的步骤：

1.读取音频文件并绘制时域波形图

[y, Fs] = audioread('audio.wav'); % 读取音频文件
t = (0:length(y)-1)/Fs; % 时间向量
subplot(2,1,1); % 绘制时域波形图
plot(t, y);
xlabel('Time (s)');
ylabel('Amplitude');
title('Time Domain Waveform');

2.进行傅立叶变换并绘制频域波形图

Y = fft(y); % 进行傅立叶变换
L = length(y); % 信号长度
P2 = abs(Y/L); % 双侧频谱
P1 = P2(1:L/2+1); % 单侧频谱
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L; % 频率向量
subplot(2,1,2); % 绘制频域波形图
plot(f, P1);
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Frequency Domain Waveform');

3.绘制频谱图

spectrogram(y, 256, 250, [], Fs, 'yaxis'); % 绘制频谱图

注意事项：
- 在进行傅立叶变换之前，需要对信号进行预处理，例如去除直流分量、加窗等操作。
- 在绘制频谱图时，可以使用spectrogram函数，该函数可以绘制出信号在时间和频率上的变化情况。