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