periodogram for 200 Hz sample rate returns 129 array length

Question

sam 2023-1-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1899915-periodogram-for-200-hz-sample-rate-returns-129-array-length

评论： sam 2023-1-25

why periodogram for 200 Hz sample rate, for 1 second, returns an array lenght of 129, if called without parameters. why it's not 100

signal = linspace(0,3000,200)'; %this is 1 second of my signal (200 because that's the sample rate)
[x,y] = periodogram(signal,[],[],200);
x
x = 129×1
0e+06 *

2500
1035
2091
0750
0438
0316
0222
0152
0112
0093

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Paul 2023-1-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1899915-periodogram-for-200-hz-sample-rate-returns-129-array-length#answer_1155475

在 MATLAB Online 中打开

Hi sam,

signal = linspace(0,3000,200)'; %this is 1 second of my signal (200 because that's the sample rate)
[x,y] = periodogram(signal,[],[],200);
numel(signal)
ans = 200

Because signal is only 200 elements and nfft is set for default, nfft will be (according to periodogram)

nfft = max(256,2^nextpow2(length(signal)))
nfft = 256

The default for "freqrange" is 'onesided', in which case with nfft = 256, x will have length nfft/2 + 1 = 129 (also on periodogram)

3 个评论
显示 1更早的评论隐藏 1更早的评论

Paul 2023-1-24

编辑：Paul 2023-1-24

在 MATLAB Online 中打开

I don't know what you mean by "100 Hz and not 129 bins," as in I don't know how you're relating Hz and bins. In the original code, the periodogram was computed as

signal = linspace(0,3000,200)'; %this is 1 second of my signal (200 because that's the sample rate)
[x,y] = periodogram(signal,[],[],200);

The resulting frequency vector y has

numel(y)
ans = 129

129 bins and the final frequency in y is

y(end)
ans = 100

100 Hz, which is what happens when nfft is even and freqrange is onesided. With nfft even, numel(y) will always be odd when freqrange is onesided.

In the new code

[x,y] = periodogram(signal,[],199,200);

the frequency vector has 100 bins

numel(y)
ans = 100

but the final frequency point is just short of 100 Hz

y(end)
ans = 99.4975

which is what happens when nfft is odd and freqrange is onesided.

sam 2023-1-25

I beleive I need to do more reading on nfft, I wanted to end up with a frquency vector of 100 bins only not 129, that's why I thought using 199.

the part that confuses me is if changing the nfft might affect my periodogram results. I guess this is beyond the scope of this post and MATLAB itself.

thank you.

Sam

请先登录，再进行评论。