clear?all;clf;close?all;clc;
Fs=2048;
t=0:1/Fs:1;
x1=2*cos（2*pi*40*t）;
nfft=2048;
R11=xcorr（x1‘unbiased‘）;
g11=fft（R11nfft）;
%g11=fftshift（R11）;
G11=abs（g11）;
index=0:round（nfft/2-1）;
f=index*Fs/nfft;
x2=2*cos（2*pi*40*（t+0.2））;
R22=xcorr（x2‘unbiased‘）;
g22=fft（R22nfft）;
%g22=fftshift（R22）;
G22=abs（g22）;
G11=1./G11;
plot_G11=10*log10（G11（index+1））;
plot_G22=10*log10（G22（index+1））;
figure（1）plot（f1）;
xlabel（‘頻率‘）
ylabel（‘幅度dB‘）
%title（‘正交相關‘）
figure（2）plot（fabs（plot_G11）‘-g‘）;
xlabel（‘頻率‘）
ylabel（‘幅度（dB）‘）
%title（‘Rogh‘）
G11=abs（g11）;
plot_G11=10*log10（G11（index+1））;
plot_G11=10*log10（G11（index+1））;
scot=1./（sqrt（G11.*G22））;
SCOT=10*log10（scot（index+1））;
%scot=sqrt（plot_G11+plot_G22）;
%SCOT=scot.^-1;
figure（3）plot（fabs（SCOT））;
xlabel（‘頻率‘）
ylabel（‘幅度（dB）‘）
%title（‘SCOT‘）
R12=xcorr（x1x2‘unbiased‘）;
g12=fft（R12nfft）;
%g12=fftshift（R12）;
G12=1./（abs（g12））;
plot_G12=10*log10（G12（index+1））;
phat=（plot_G12）;
PHAT=phat;
figure（4）plot（f