r.m

计算太阳黑子的活动周期

%求信号自相关函数
function R=r(x)
N=length(x);
for m=1:N                     
    S=0; 
    for n=1:N-m
       S=S+x(n)*x(n+m-1);
    end
    R(m)=S/N;
end