multiple fisher criteria for periodic components
Show older comments
i have following code
clear all;
clc;
close all;
fs=100;
ts=1/fs;
t=0:ts:2.93;
x=24*sin(2*pi*20*t)+randn(size(t));
[Pxx,f]=periodogram(x,[],[],fs);
[maxval,index] = max(Pxx);
fisher_g = Pxx(index)/sum(Pxx);
N = length(Pxx);
upper = floor(1/fisher_g);
for nn = 1:3
I(nn) = ...
(-1)^(nn-1)*nchoosek(N,nn)*(1-nn*fisher_g)^(N-1);
end
pval = sum(I)
%end program;
pval is
pval =
1.2667e-70
what about when we more then one periodic component? please help me
Answers (0)
Categories
Find more on Pulse and Transition Metrics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!