1

pwelch用来绘制功率谱密度。我想使用格式

pwelch=(x,window,noverlap,nfft,fs,'onesided')

但在 x 轴上有一个对数刻度。

我也试过

[P,F]=(x,window,noverlap,nfft,fs);
plot(F,P)

但它没有给出与上面相同的结果图。所以,

semilogx(F,P) 

不是一个好的解决方案。

4

1 回答 1

2

好的,所以首先,我从未听说过此功能或此方法。但是,我能够生成与使用输出参数生成的函数相同的图。我从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;生成线性版本(相同的精确图,减去标签):

图1

或者

semilogx(bar,10*log10(foo)); grid on;对于 x 轴上的对数刻度。 图2

我不喜欢 x 比例是线性采样但对数显示(这是一个词对吗?),但它看起来还不错。

够好了?

于 2013-10-15T22:52:00.583 回答