relation.m

用matlab 实现强噪声中弱信号的提取

function y = relation(x1,x2,x1Equ,x2Equ,x1Var,x2Var,taoMax,N)

for tao = 1:(taoMax)
    t = 1;
    sum = 0;
    while t + tao - 1 < N
        sum = x1(t) * x2(t + tao - 1) + sum;
        t = t + 1;
    end

    y(tao) = ((sum / (N - tao)) - x1Equ * x2Equ) / (sqrt(x1Var) * sqrt(x2Var));
end