www.pudn.com > fourierCyclic.rar > calcyclic.m

function scf=calcyclic(tfr); 
% calculate the scf of the signal by time-frequence expression 
% 
% tfr is a matrix of [f,t] 
 
[xrow,xcol]=size(tfr); 
 
scf=zeros(xrow,xcol); 
X=zeros(1,xrow); 
X1=zeros(1,xrow); 
X2=zeros(1,xrow); 
 
%for cyclic frequence loop 
for a=1:xrow 
    %for time loop 
    %let X = 0 
    X=zeros(1,xrow); 
 
    %set X1 and X2 move length 
    a1=floor(a/2); 
    if a/2==.5, 
        a2=a1+1; 
    else 
        a2=a1; 
    end 
 
    for i=1:xcol 
        X1=[tfr((xrow-a1+1):xrow,i)' tfr(1:(xrow-a1),i)']; 
        X2=[tfr((a2+1):xrow,i)' tfr(1:a2,i)']; 
        X=X+(X1).*conj((X2)); 
    end 
 
    if a==64 
        %plot(1:512,X1); 
        %figure; 
        %plot(1:512,X2); 
        %figure; 
        %plot(1:512,X); 
        b=a; 
    end     
    %sum up the short time frequence 
    scf(:,a+1)=X'./xcol; 
end 
 
end