如何用matlab比较两个音频信号

Question

我有两个音频信号要使用 Matlab 进行比较，我的问题是我可以清晰地听到它们之间的区别，但是当使用函数 pwelch 比较它们的 PSD 时，我看不到太多它们之间的区别，知道如何将它们与不同的方法进行比较！提前致谢 ！ PS：1.我已经在DSP.stackexchange中问过问题，但没有答案！ 2.我不是要代码，这就是为什么我没有放我的！

更新

在@Bas Swinckels 的回答之后，我正在使用结果的代码标准，但我仍然看不到描述这个的方式？

  clear;
clc;
[x0,Fs] = audioread('Nonoise.wav');
x0 = x0(:,1);
[Y0,G] = pwelch(x0,hamming(512));
plot(G,10*log10(Y0));
grid on 

[x50,Fs] = audioread('50% noise.wav');
x50 = x50(:,1);
[Y50,G] = pwelch(x50,hamming(512));
hold on ;
plot(G,10*log10(Y50),'r');


[x100,Fs] = audioread('100% noisy.wav');
x100 = x100(:,1);
[Y100,G] = pwelch(x100,hamming(512));
hold on ;
plot(G,10*log10(Y100),'g');

%% spectrogram  
spectrogram(x0,hann(64),32,64,Fs)
figure();
spectrogram(x50,hann(64),32,64,Fs)
figure();
spectrogram(x100,hann(64),32,64,Fs)

结果如下：

谁能告诉我我添加到源中的噪声或效果如何影响我的频谱？

Answer 1

您可以尝试绘制两个信号的spectrogram，如下所示：

spectrogram(x, hann(nfft), nfft/2, nfft, fsample)

当您使用pwelch 计算 PSD 时，您只会获得平均频谱，并且会丢失所有时间信息：它不会显示信号在开始时或结束时是否较大。计算频谱图也与您的耳朵所感知的非常相似，它们测量一种频谱作为时间的函数。因此，如果您能听到差异，它应该在具有足够 SNR 的频谱图中显示为差异。

当您制作频谱图时，您应该尝试使用频率和时间分辨率，因为它们遵循某种海森堡定律。缩短 FFT 窗口（通过选择较小的 nfft）会给您提供良好的时间分辨率，但会降低频率分辨率，反之亦然。