stap机载雷达正侧视阵地杂波功率谱capon图

%                杂波建模     正侧阵
%  ======================================================================
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;