日日行,不怕千万里;时时学,不怕千万卷
分类:
2008-03-28 16:59:17
%相关分析法计算相位差 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),'度'))); |
这个程序大家一看应该会懂,具体的公示推导可以查找相关的信号处理书籍,这里省略