wav格式的音频信号的时域和时频谱图分析

首先，录入文件

[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]);

最终显示图：