%使用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/pia=fft(B);[F,I]=max(abs(a));            c=angle(i*real(a(I))-imag(a(I)))*180/pia=fft(Z);[F,I]=max(abs(a));            %對A和B信號分別加了躁聲后的實際信號仿真做fft,并且求相位 d=angle(i*real(a(I))-imag(a(I)))*180/pia=fft(O);[F,I]=max(abs(a));            e=angle(i*real(a(I))-imag(a(I)))*180/pitoc