tic;
clc;
t=linspace(0,40*pi,100);  %基于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;

F=fft(A);             %對A和B信號做fft,并且求相位
a=angle(F)*180/pi; 
R=fft(B);
b=angle(R)*180/pi;
Q=fft(Z);              %對A和B信號分別加了躁聲后的實際信號仿真做fft,并且求相位
c=angle(Q)*180/pi;
M=fft(O);  
d=angle(M)*180/pi;

subplot(2,2,1)
plot(c)
xlabel('\omega/\omega_0');
ylabel('angle');
title('A加躁聲后信號的相位 ');
subplot(2,2,2)
plot(d);
xlabel('\omega/\omega_0');
ylabel('angle');
title('B加躁聲后信號的相位 ');
subplot(2,2,3)
plot((b-a));       %求A和B信號相位差
xlabel('\omega/\omega_0');
ylabel('angle');
title('A和B信號相位差 ');
subplot(2,2,4)
plot((d-c));%求A和B信號加入躁聲后的相位差
xlabel('\omega/\omega_0');
ylabel('angle');
title('A和B信號加入躁聲后的相位差');
p1=(b-a)*180/pi
p2=(d-c)*180/pi
toc