我pwelch
用来绘制功率谱密度。我想使用格式
pwelch=(x,window,noverlap,nfft,fs,'onesided')
但在 x 轴上有一个对数刻度。
我也试过
[P,F]=(x,window,noverlap,nfft,fs);
plot(F,P)
但它没有给出与上面相同的结果图。所以,
semilogx(F,P)
不是一个好的解决方案。
好的,所以首先,我从未听说过此功能或此方法。但是,我能够生成与使用输出参数生成的函数相同的图。我从help
文本中运行了示例。
EXAMPLE:
Fs = 1000; t = 0:1/Fs:.296;
x = cos(2*pi*t*200)+randn(size(t)); % A cosine of 200Hz plus noise
pwelch(x,[],[],[],Fs,'twosided'); % Uses default window, overlap & NFFT.
这产生了这个情节:
然后我做了:plot(bar,10*log10(foo)); grid on;
生成线性版本(相同的精确图,减去标签):
或者
semilogx(bar,10*log10(foo)); grid on;
对于 x 轴上的对数刻度。
我不喜欢 x 比例是线性采样但对数显示(这是一个词对吗?),但它看起来还不错。
够好了?