如何使用Scilab获得类似大胆的类似情节

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的光谱实际上是periodogram，它是光谱密度的估计值。尝试使用pwelch function。这是我的Audacity频率分析：

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

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