压缩感知DOA估计

clc;clear;
%%%%%%%%%%%%%%用Compressive Sensing方法利用估计MIMO雷达的DOA
M=30;                                    %发射阵元个数
N=10;                                     %接收阵元个数
L=256;                                   %快拍数目
K=3;                                     %稀疏度
SNR=-20;                                  %信噪比
%f=8.62e12;                               %载波频率
%c=3e8;                                   %传播速度
%lamda=c/f;                               %波长
dt=0.5;                            %发射阵元间距取半个波长
dr=0.5;                            %接收阵元间距取半个波长
%%%%%有3个目标方位2°3°4°
theat0=[-5 2 7];                          %目标与阵列的夹角
%d=[1e9 2e8 3e10];                        %目标距离坐标原点即第一个发射天线的距离
beata=[5 3 4];                           %目标的反射系数
Q=length（theat0）;
s0=hadamard（L）;                          %产生一组正交信号
for m=1:M
    s（:m）=s0（:m）;                      %M个发射信号源
end
n=sqrt（10^（-SNR/10）/2）;
for q=1:Q
    a（:q）=exp（-j*2*pi*dt*sin（theat0（q）*pi/180）*[0:M-1]‘）;%发射阵列引导矢量
    b（:q）=exp（-j*2*pi*dr*sin（theat0（q）*pi/180）*[0:N-1]‘）;%接收阵列引导矢量
end

noise=n*randn（NL）+j*n*randn（NL）;        %N个接收通道L次快拍的噪声
r=zeros（NL）;
for q=1:Q
    r=r+beata（q）*b（:q）*a（:q）.‘*s.‘;
end
r=r+noise;                                %接收的信号
rr=reshape（r.‘N*L1）;
%%%%%CS定理
theat=-10:0.1:10;P=length（theat）;           %将角度空间稀疏化
for p=1:P
    aa（:p）=exp（-j*2*pi*dt*sin（theat（p）*pi/180）*[0:M-1]‘）;      %发射阵列引导矢量
    bb（:p）=exp（-j*2*pi*dr*sin（theat（p）*pi/180）*[0:N-1]‘）;      %接收阵列引导矢量
    psi（:p）=reshape（（bb（:p）*aa（:p）.‘*s.‘）.‘N*L1）;  %稀疏域的基矩阵                   
end

KK