%*****單基地雜波模型******
clc
clear
tic
lambda=0.23;d=lambda/2;fr=2434.8;
M=8;N=6;K=8;
V=140;H=12000;V1=10;%V1為風(fēng)速
R0=50e3;????%目標(biāo)所在距離環(huán)斜距
theta0=60*pi/180;%主波束方位角
phi0=asin（H/R0）;%主波束俯仰角
thetap=60*pi/180;??%非正側(cè)面天線與飛行方向夾角
theta_a=1*pi/180;????%水平方位面內(nèi)波束寬度
pt=200e3;
Im=chebwin（M20）;In=chebwin（N20）;
CNR=60;
%***************
%****雷達(dá)系統(tǒng)參數(shù)****
GG=pt*lambda^2/（4*pi）^3;%增益
tau=0.1e-6;%脈沖寬度？
R_cell=150;%距離分辨率
yangben=2*N*K;%
%************
Re=8490e3;c=3e8;
k=1.38e-23;T0=290;Bn=70e6;Fn=3.5;Fn=10^（3.5/10）;
Pn=k*T0*Bn*Fn;
%******
theta=[0:pi/179:pi];
W=length（theta）;
R2=zeros（N*KN*K）;
%***************
earth_leixing=2;%可以用來(lái)控制地表類型
clutter_leixing=2;%控制雜波幅度類型
ops=0;
for?l=1:yangben?%距離門采樣訓(xùn)練樣本
????ops=ops+1;
????cr=zeros（NK）;
????r（l）=R0+（l-yangben/2）*R_cell;
????phi（l）=asin（H/r（l））;
%?????aaa=unifrnd（02*pi11）;
%??????gaosi_noise=Pn*randn（N*K1）.*exp（j*2*pi*rand（N*K1））;
????for?w=1:W???%方位角循環(huán)
????????cos_psi（w）=cos（theta（w））*cos（phi（l））;
????????%*****空間、時(shí)間角頻率
????????ws（lw）=2*pi*d/lambda*cos（theta（w））*cos（phi（l））;
????????wt（lw）=4*pi*V/（lambda*fr）*cos（theta（w）+thetap）*cos（phi（l））;
????????fd（lw）=2*V/lambda*cos（theta（w）+thetap）*cos（phi（l））/fr;
????????%*****天線方向圖
????????F（lw）=sum（In‘.*exp（j*2*pi*d/lambda*[0:N-1]*（cos（theta（w））-cos（theta0））））*...
????????????sum（Im‘.*exp（j*2*pi*d/lambda*[0:M-1]*（cos（theta（w））*cos（phi（l））-cos（theta0）*cos（phi0））））;
????????for?n=1:N
????????????G（n）=sum（Im‘.*exp（j*2*pi*d/lambda*[0:M-1]*（sin（phi（l））-sin（phi0））））;
????????end
????????%****雜波散射單元面積**
????????Ac=r（l）*theta_a*c*tau/2;????%大俯角情況下雜波散射單元面積；theta_a：方位角平面內(nèi)波束寬度。
????????%***雜波散射系數(shù)sigma0*****
????????if?earth_leixing==1?%農(nóng)田
????????????A=0.004;B=pi/2;beta0=0.2;sigma_c0=1;
????????elseif?earth_leixing==2?%丘陵
????????????A=0.0126;B=pi/2;beta0=0.4;sigma_c0=1;
????????elseif?earth_leixing==3?%高山
????????????A=0.04;B=1.24;beta0=0.5;sigma_c0=1;
????????end
????????u=sqrt（c/lambda*1e-9）/4.7;
????????theta_g=asin（H/r（l）+H^2/（2*Re*r（l））-r（l）/（2*Re））;
????????sigma0（l）=A*sigma_c0*sin（theta_g）/lambda+u*cot（beta0）^2*exp（-tan（B-theta_g）^2/tan（beta0）^2）;
????????%***雜波分辨單元雷達(dá)散射截面積
????????sigma_c=sigma0（l）*Ac;
????????%***雜波幅度控制
????????if?clutter_leixing==1???%瑞利分布
????????????al=ra