%-----------------------------------------%
%???????????xulingji?2011.7.7?????????????%
%?雙邊相關變換（TCT）估計寬帶相干信號源方位?%
%-----------------------------------------%
clear?all
close?all
clc
tic

jay=sqrt（-1）;
M=16;???????????????%陣元數
N=15000;????????????%仿真產生的快拍數
N_use=12800;????????%實際用的快拍數
n1=1000;????????????%數據起始點??????????
f0=2000;????????????%以f0為聚焦頻率
f_l=3*f0/5;?????????%寬帶信號上邊界頻率
f_u=f0;?????????????%寬帶信號下邊界頻率
delta_F=40;
F=f_l:delta_F:f_u;
B=（f_u-f_l）;????????%絕對帶寬
B_delta=（f_u-f_l）/（（f_l+f_u）/2）;??%相對帶寬（一般取50%）%%%這樣子難道頻率上下界都要有規律？
Fs=5*f0;
c=1500;
lambda=c/f0;????????%f0對應的波長
d=1/2*lambda;???????%陣元間距
rsj=[-1.5?1];???????
Num=length（rsj）;
Q1=-15;
Q2=15;
deta_theta=0.02;
Angles=Q1:deta_theta:Q2;
z=0:M-1;
SNR=15;?

L=100;???????????????
k1=-L:L;
index=1;
nfft=256;
delta_f=Fs/nfft;
f1=round（f_l/delta_f+1）;????%頻率對應于FFT中的點數
f2=round（f_u/delta_f+1）;
freq=f1:f2;
freq_hz=（freq-1）*Fs/nfft;???%實際分析的頻率

Q=512;????????????????????????
wn1=800*2/Fs;wn2=2200*2/Fs;???
b=fir1（Q[wn1wn2]）;
no1=randn（1N）;
no2=randn（1N）;
no3=randn（MN）;
s1=filter（b1no1）;
s2=filter（b1no1）;
noise=filter（b1no3[]2）;??%噪聲
t0=10*1/Fs;
for?i=1:M
????tao1（i）=d*sin（rsj（1）*pi/180）/c*z（i）*Fs;????????%入射信號1延遲點數?
????tao2（i）=d*sin（rsj（2）*pi/180）/c*z（i）*Fs+t0*Fs;??%入射信號2延遲點數?
end
for?j=1:M
????Coef1=sinc（k1+tao1（j））;
????x1（j:）=filter（Coef11s1）;
????Coef2=sinc（k1+tao2（j））;
????x2（j:）=filter（Coef21s2）;
end
S=x1+x2+（10^（-SNR/20））*noise;
S=S（:n1:n1+N_use-1）;
%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?
Section=1;??????????????????????%分段數的起始值
Rk=zeros（MMlength（freq））;
while?index+nfft-1?????X=fft（S（:index:index+nfft-1）nfft