首先,录入文件
[wave,f]=audioread('文件路径及名称.wav');
wave存储为音频信号的时域波形幅度值,f为采样频率
n=length(wave);
n为文件总采样点数
t=(1:n)/(f*60);
将时间轴变为以分钟为单位,这步可自行改变
画时域波形图
subplot(211);plot(t,wave);
%axis([1 n/(f*60) -0.01 0.01]);
上步为调节横纵轴的显示范围
xlabel('Time(min)');ylabel('Amplitude');
画时频谱图
subplot(212);
[S, F, T] = spectrogram(wave, hanning(1024), 512, 1024, f);
tt=T/60;
将时间轴变为以分钟为单位
Ff=F/1000;
将频率轴变为以kHz为单位
imagesc(tt, Ff, log10(abs(S)));set(gca, 'YDir', 'normal');
xlabel('Time(min)');
ylabel('Frequency(kHz)');
%axis([0 10 0 16]);
最终显示图: