shiyan.m

%*****单基地杂波模型******
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为风速
R0=50e3;    %目标所在距离环斜距
theta0=60*pi/180;%主波束方位角
phi0=asin（H/R0）;%主波束俯仰角
thetap=60*pi/180;  %非正侧面天线与飞行方向夹角
theta_a=1*pi/180;    %水平方位面内波束宽度
pt=200e3;
Im=chebwin（M20）;In=chebwin（N20）;
CNR=60;
%***************
%****雷达系统参数****
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;%可以用来控制地表类型
clutter_leixing=2;%控制杂波幅度类型
ops=0;
for l=1:yangben %距离门采样训练样本
    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   %方位角循环
        cos_psi（w）=cos（theta（w））*cos（phi（l））;
        %*****空间、时间角频率
        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：方位角平面内波束宽度。
        %***杂波散射系数sigma0*****
        if earth_leixing==1 %农田
            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）;
        %***杂波分辨单元雷达散射截面积
        sigma_c=sigma0（l）*Ac;
        %***杂波幅度控制
        if clutter_leixing==1   %瑞利分布
            al=ra