matlab - 如何使用 Scilab 获得类似 Audacity 的情节

Question

我有一个音频文件，我大胆地分析它，频谱看起来像：

我想用 scilab 得到同样的结果，所以我运行了这段代码：

[y,Fs] = wavread('fileName.wav');
Y_0 = abs(fft(y(1,:),-1));
plot(fftshift(Y_0));

结果是：

我也试过这个：

plot(10*log10(fftshift(Y_0)));

结果不同，但仍然与大胆不一样，知道该怎么做吗？提前感谢您的帮助！

Answer 1

Audacity 的频谱实际上是一个周期图，它是对频谱密度的估计。尝试使用pwelch 函数。这是我的 Audacity 频率分析：

使用 MATLAB 的 pwelch 函数重新创建：

[y, fs] = wavread('inputfile.wav');
pwelch(y, hanning(512))