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

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


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

SCR=-30;????????????????%信雜比
CNR=60;?????????????????%雜噪比
%%%%%%-----------天線方向圖-------------
dB=16;
Im=chebwin（MdB）;???%天線陣列電流強(qiáng)度，加窗????????????????????????????
In=chebwin（NdB）;???%天線陣行電流強(qiáng)度，加窗
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為指向目標(biāo)的加權(quán)項(xiàng)
????F（i:）=F_col;%*F_row（i）;?????%i對(duì)應(yīng)高低角的位置同一個(gè)i中高低角不變遍歷方位角???F為原始發(fā)送到各點(diǎn)的信號(hào)
end

%%%%%%%%%%------------雜波反射------------%%%目標(biāo)所在的方位角在天線指向的位置????????并假定目標(biāo)位于中間距離環(huán)上?
clutter1=zeros（N*K1）;??????????%用于存儲(chǔ)每一次快拍（距離環(huán)）的列矢量（雜波＋信號(hào)）　???臨時(shí)存放
clutter2=zeros（N*K1）;??????????%用于存儲(chǔ)每一次快拍（距離環(huán)）的列矢量（雜波）
signal=zeros（N*K1）;????????????%目標(biāo)信號(hào)矩陣
clutter1_matrix=zeros（N*KL）;???%用于存儲(chǔ)所有距離環(huán)的雜波＋信號(hào)?矩陣
clutter2_matrix=zeros（N*KL）;???%用于存儲(chǔ)除去動(dòng)目標(biāo)的所有距離環(huán)雜波矩陣?（包含動(dòng)目標(biāo)所在距離環(huán)的雜波）
Rcs_clutter_signal=zeros（N*KN*KL）;????%三維矩陣?用于存儲(chǔ)每一個(gè)距離環(huán)一次快拍的相關(guān)矩陣（未平均）（包含信號(hào)）
Rc_clutter=zeros（N*KN*KL）;????????????%三維矩陣?用于存儲(chǔ)每一個(gè)距離環(huán)一次快拍的相關(guān)矩陣（未平均）（不包含信號(hào)）
Rcs=zeros（N*K）;?????????????????????????%包含動(dòng)目標(biāo)距離環(huán)估計(jì)的協(xié)方差矩陣
Rc=zeros（N*K）;??????????????????????????%包含動(dòng)目標(biāo)距離環(huán)但不包含動(dòng)目標(biāo)信號(hào)估計(jì)的協(xié)方差矩陣?（雜波＋噪聲?協(xié)方差矩陣）
for?i=1:L
????i;
????B_coe=randn（1length（theta））.*exp（j*2*pi*randn（1length（theta）））;????????%第i個(gè)距離環(huán)上的雜波反射系數(shù)幅度服從正態(tài)分布，相位服從正態(tà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次快拍時(shí)的時(shí)域角頻率???i距離環(huán)的多普勒角頻率????????
????????????Ws=4*pi*d/lambda*cos（theta）*cos（phi（i））;????????????????????%第n個(gè)天線的空域角頻率
????????????s=exp（j*（k-1）*Wt+j*（n-1）*Ws）/R（i）^2;????????????????????????