用MATLAB写一段代码,绘制音频文件的时域波形图和频谱波形图
时间: 2024-02-06 15:12:52 浏览: 100
好的,以下是使用MATLAB绘制音频文件的时域波形图和频谱波形图的代码:
```matlab
% 读取音频文件
[y, Fs] = audioread('audio_file.wav');
% 定义时间向量和频率向量
t = 0:1/Fs:(length(y)-1)/Fs;
N = length(y);
f = (-Fs/2:Fs/N:Fs/2-Fs/N);
% 计算信号的傅里叶变换
Y = fftshift(fft(y));
% 绘制时域波形图
subplot(2,1,1);
plot(t, y);
xlabel('Time (s)');
ylabel('Amplitude');
title('Time Domain Waveform');
% 绘制频谱波形图
subplot(2,1,2);
plot(f, angle(Y));
xlabel('Frequency (Hz)');
ylabel('Phase');
title('Frequency Domain Waveform');
```
在该代码中,我们首先使用`audioread`函数读取音频文件,并将其存储在`y`中,采样率存储在`Fs`中。
然后,我们定义了时间向量`t`和频率向量`f`,并使用`fftshift`函数将信号的傅里叶变换移动到中心。
接着,我们使用`subplot`函数创建一个2x1的图形区域,并在第一个子图中绘制音频文件的时域波形图。我们使用时间向量`t`和信号向量`y`来绘制时域波形图。
然后,我们在第二个子图中绘制音频文件的频谱波形图。我们使用`fft`函数计算信号的傅里叶变换,并使用`angle`函数来计算相位谱。我们使用频率向量`f`来绘制频谱波形图。
希望这个代码可以帮到你!
阅读全文