multiple fisher criteria for periodic components

2 次查看（过去 30 天）

dato datuashvili 2015-2-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/174599-multiple-fisher-criteria-for-periodic-components

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;