b2.m

相关信号的相位差检测

%使用fft检测相位

tic;
clear all;
clc;
sample=4096;
t=linspace(0,4*pi,sample);  %基于fft求相位的算法
A=sin(t+pi/3);
B=sin(t+3*pi/4);
N=length(t);
tp=5;                                      
Ws=N/tp;
R=0.0001*rand(1,N);
P=0.0001*rand(1,N);
Z=A+R;
O=B+P;

a=fft(A);
[F,I]=max(abs(a));             %对A和B信号做fft,并且求相位
b=angle(i*real(a(I))-imag(a(I)))*180/pi
a=fft(B);
[F,I]=max(abs(a));            
c=angle(i*real(a(I))-imag(a(I)))*180/pi
a=fft(Z);
[F,I]=max(abs(a));            %对A和B信号分别加了躁声后的实际信号仿真做fft,并且求相位 
d=angle(i*real(a(I))-imag(a(I)))*180/pi
a=fft(O);
[F,I]=max(abs(a));            
e=angle(i*real(a(I))-imag(a(I)))*180/pi

toc