如何使用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的光谱实际上是一个周期图,它是光谱密度的估计值.尝试使用pwelch功能.这是我的Audacity频率分析:

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

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