tct宽带非相干信号目标方位估计matlab程序

%-----------------------------------------%
%           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 < size（S2）  %构造出频域相关矩阵
    X=fft（S（:index:index+nfft-1）nfft