function?[?doa?]?=?TctSpectrum（SensorLd_Xanglesnr?）
bw1=1e7;????????????????????????????????????%?信號帶寬
f0=1e8;?????????????????????????????????????%固定頻率
T1=1e-5;????????????????????????????????????%采樣長度
dt1=T1/L;?????????????????????????????????%采樣間隔
gc?=?3e+8;?????????????????????????????????%光速
snr=5;??????????????????????????????????%信噪比
degrad=pi/180;
p=1;???????????????????????????????????????%?信號源個數(shù)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%構(gòu)造采樣序列
t1=（0:L-1）*dt1;???????????????????????????%時域采樣點
sig1=exp（j*pi*bw1*t1.*t1/T1+j*2*pi*f0*t1）;???????????????????%產(chǎn)生線性調(diào)頻信號
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%線性調(diào)頻信號頻率
F0=f0+bw1/（T1）*t1;
le=L/8;?????????????????????????%寬帶信號分為8個窄頻段
for??i=1:8
????f（i）=F0（i*le）;
????F（1i）=f（i）;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%構(gòu)造采樣數(shù)據(jù)矩陣
Lambda=gc/f（8）;???????????????????????????????%?波長
d=d_X*Lambda;???????????????????????????????%陣元間距離
tmp2=d*[0:Sensor-1?]‘;
???????for??i=1:8
???????????v1i=sig1（le*（i-1）+1:le*i）;
???????????s（i:）=v1i;
???????end
???????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
???????
???????%構(gòu)造方向矩陣
???????A=[exp（-j*2*pi*tmp2*F/gc*sin（angle*degrad））];???%參考頻率點的陣列流型信號源決定列數(shù)，陣元數(shù)決定行數(shù)8行8列（信號分為8段）
???????
???????%?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
???????%信號接收
???????S=zeros（8Sensor）;
???????P=zeros（SensorSensor*8）;
?
???????for?i=1:8
???????????xi=A（:i）*s（i:）;
???????????xi=awgn（xisnr）;??????????????%接收信號
???????????Xi=fft（xi.‘）;
???????????Xi=Xi.‘;??????????????????????%FFT
???????????Ri=Xi*Xi‘;??????????????????%計算協(xié)方差矩陣
???
???????????temp=eig（Ri）;
???????????temp=sort（temp）;
???????????Pi=Ri-[sum（temp）-sum（temp（Sensor-p+1:Sensor））]/（Sensor-p）*eye（SensorSensor）;???%計算協(xié)方差矩陣
???????????[UiSiVi]=svd（Pi）;??????????????%奇異值分解
???????????S（i:）=sum（Si）;
???????????P（:Sensor*（i-1）+1:Sensor*i）=Pi;
???????????X（Sensor*（i-1）+1:Sensor*i:）=Xi;
???????end
???????
???????%?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%