Parameterize overlapping - cpsd()
显示 更早的评论
Hello, everyone.
I need to create a script where overlapping is parameterized. I did it this way:
N = length(u); % = length(y)
percentage = 0.5; %range = [0.01 - 0.99]
window = hanning(round(N/4.5));
overlap = round(length(window)*percentage);
[SuyW,fw] = cpsd(u,y,window,overlap,[],fs);
Can anyone tell me if this is correct?
采纳的回答
[edit: correcting my recommendation from overlap=floor(...) to overlap=ceil(...)]
Yes this is reasonable.
You can answer your quesiton by trying some different values for N and observing the resulting values for window length and overlap:
%N = length(u); % = length(y)
N=[500,1000,2000,4000];
percentage = 0.5; %range = [0.01 - 0.99]
for i=1:4
window = hanning(round(N(i)/4.5));
overlap = round(length(window)*percentage);
fprintf('N=%d, winLen=%d, overlap=%d\n',N(i),length(window),overlap);
end
The results look reasonable.
You chose a window length that is 1/4.5 times the signal length, in all cases. This choice will usually allow for 8 half-overlapped windows in the signal. It is possible that, due to rounding, the 8th and final window will not fit. For example, N=98, N=502, N=512, N=1006, N=1024, and others. Therefore I would use
window = hanning(floor(N/4.5));
overlap = ceil(length(window)*percentage);
to assure that the window goes onto the signal 8 times.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Spectral Estimation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
