gusucode.com > 相关分析法计算相位差 > 相关分析法计算相位差/Untitled15.m

    %相关分析法计算相位差
num=100;
l=0:1:num;
x=sin(2*pi*l/num);
y=cos(2*pi*l/num);
Ix=sum(x.^2)/num;
Iy=sum(y.^2)/num;
Ixy=sum(x.*y)/num;
c=180*acos(2*Ixy/(4*Ix*Iy)^0.5)/pi;
plot(l,x,l,y);
legend('sin(x)','cos(x)');
text(50,0.6,strcat('相位差=',strcat(num2str(c),'度')));