clear?all;clc;
M=5;%陣列個數(shù)
P=1;%信號個數(shù)
N=1024;%快拍數(shù)
C=3.0E+8;
f0=1.0E+9;%載波頻率
d=（C/f0）/2;%兩個相鄰陣元距離，不能大于載波的半波長
Radius=2*d;%均勻圓陣半徑
Bearings=[40];%入射角方位角
Pitchings=[60];%入射信號俯仰角


AmpCoeff=[100];%陣元放大系數(shù)
f=[500];%基帶信號頻率
fs=4*max（f）+1;%采樣頻率


AmpCoeff=repmat（AmpCoeff[1?N]）;
S=AmpCoeff.*sin（2*pi*f*[1:N]/fs）;
Noise=randn（MN）;
Bearings=Bearings*pi/180;
Pitchings=Pitchings*pi/180;


TimeDelay=d*[0:M-1].‘*sin（Bearings）/C;%線陣的延遲

A=exp（-1*sqrt（-1）*2*pi*f0*TimeDelay）;

X=A*S+50*Noise;

Rxx=（X*X‘）/N;??%correlation?matrix

[EigenVectorsEigenValues]=eig（Rxx）;

Lemda=diag（EigenValues）;??%計算矩陣特征值
[SortedLemdaIndex]=sort（Lemda）;
Index=flipud（Index）;%將特征值降序排列


NoiseSubspace（1:M1:M-P）=EigenVectors（1:MIndex（P+1:M））;


Delta=[-90:1:90];
for??k=1:length（Delta）
????sum_spec（k）=0;
end

for?i=1:1:180
????a（1:）=exp（-1*sqrt（-1）*2*pi*f0*Radius*cos（2*pi*0/8-Delta*（pi/180））*cos（i*（pi/180））/C）;
????a（2:）=exp（-1*sqrt（-1）*2*pi*f0*Radius*cos（2*pi*1/8-Delta*（pi/180））*cos（i*（pi/180））/C）;
????a（3:）=exp（-1*sqrt（-1）*2*pi*f0*Radius*cos（2*pi*2/8-Delta*（pi/180））*cos（i*（pi/180））/C）;
????a（4:）=exp（-1*sqrt（-1）*2*pi*f0*Radius*c