%????????????????雜波建模?????正側陣
%??======================================================================
clc;????close?all;??clear?all;
j=sqrt（-1）;
M=8;
N=8;
K=10;???????????????????%?采樣快拍數
lambda=0.32;

%?lambda=0.12;
d=lambda/4;
Va=120;?????????????????%?飛機速度120m/s
H=6000;?????????????????%?飛機高度
Vc=0;???????????????????%?目標方位向速度
Vr=100;?????????????????%?目標距離向速度
prf=1500;???????????????%?PRF：脈沖重復頻率
c=3*10^8;
Ru=c/（2*prf）;???????????%?最大不模糊距離100km
Rmax=100000;????????????%?雷達遠地斜距?100km
Rmin=9000;??????????????%?近地距離??飛機正下方開始


theta0=90/180*pi;???????%?天線指向的方位角
phi0=30/180*pi;?????????%?天線指向的高低角
theta=[0:180]/180*pi;???%?方位角方位（0~180°）
distance=15;????????????%?距離環寬度
L=200;??????????????????%?距離環個數：L>2MN
R0?=?H./sin（phi0）;
R?=?R0?+?distance*（-round（L/2）+1:round（L/2））;??%?各個距離環對應的斜距?[-k/2:k/2]?認為平行????????????1*200?????????????
phi=asin（H./R）;?????????%?各個距離環對應的高低角

SCR=-30;????????????????%信雜比
CNR=60;?????????????????%雜噪比
%%%%%%-----------天線方向圖-------------
dB=16;
Im=chebwin（MdB）;???%天線陣列電流強度，加窗????????????????????????????
In=chebwin（NdB）;???%天線陣行電流強度，加窗
F_row=sum（diag（Im）*exp（j*2*pi*d/lambda*（0:M-1）‘*（sin（phi）-sin（phi0）））1）;??????%天線陣列合成電流幅度加窗???1*91
F=zeros（length（phi）length（theta））;?????????????????????????????????????????????????????????%天線方向圖二維矩陣
for?i=1:length（phi）
????F_col=sum（diag（In）*exp（j*2*pi*d/lambda*（0:N-1）‘*（cos（theta）*cos（phi（i））-cos（theta0）*cos（phi0）））1）;???????%線陣天線在方位角維上的分布圖，電流幅度加窗????-cos&s為指向目標的加權項
????F（i:）=F_col;%*F_row（i）;?????%i對應高低角的位置同一個i中高低角不變遍歷方位角???F為原始發送到各點的信號
end

%%%%%%%%%%------------雜波反射------------%%%目標所在的方位角在天線指向的位置????????并假定目標位于中間距離環上?
clutter1=zeros（N*K1）;??????????%用于存儲每一次快拍（距離環）的列矢量（雜波＋信號）　???臨時存放
clutter2=zeros（N*K1）;??????????%用于存儲每一次快拍（距離環）的列矢量（雜波）
signal=zeros（N*K1）;????????????%目標信號矩陣
clutter1_matrix=zeros（N*KL）;???%用于存儲所有距離環的雜波＋信號?矩陣
clutter2_matrix=zeros（N*KL）;???%用于存儲除去動目標的所有距離環雜波矩陣?（包含動目標所在距離環的雜波）
Rcs_clutter_signal=zeros（N*KN*KL）;????%三維矩陣?用于存儲每一個距離環一次快拍的相關矩陣（未平均）（包含信號）
Rc_clutter=zeros（N*KN*KL）;????????????%三維矩陣?用于存儲每一個距離環一次快拍的相關矩陣（未平均）（不包含信號）
Rcs=zeros（N*K）;?????????????????????????%包含動目標距離環估計的協方差矩陣
Rc=zeros（N*K）;??????????????????????????%包含動目標距離環但不包含動目標信號估計的協方差矩陣?（雜波＋噪聲?協方差矩陣）
for?i=1:L
????i;
????B_coe=randn（1length（theta））.*exp（j*2*pi*randn（1length（theta）））;????????%第i個距離環上的雜波反射系數幅度服從正態分布，相位服從正態分布
????B_coe_target=randn（11）*exp（j*2*pi*randn（11））;
????for?k=1:K????????
????????for?n=1:N
%????????????Wt=4*pi*Va/lambda/prf*?cos（theta+theta1）*cos（phi（i））;
????????????a=4*pi*Va/lambda/prf;
????????????Wt=4*pi*Va/lambda/prf*?cos（theta）*cos（phi（i））;???????????????%第k次快拍時的時域角頻率???i距離環的多普勒角頻率????????
????????????Ws=4*pi*d/lambda*cos（theta）*cos（phi（i））;????????????????????%第n個天線的空域角頻率
????????????s=exp（j*（k-1）*Wt+j*（n-1）*Ws）/R（i）^2;????????????????????????