clc;clear;
%%%%%%%%%%%%%%用Compressive?Sensing方法利用估計MIMO雷達(dá)的DOA
M=30;????????????????????????????????????%發(fā)射陣元個數(shù)
N=10;?????????????????????????????????????%接收陣元個數(shù)
L=256;???????????????????????????????????%快拍數(shù)目
K=3;?????????????????????????????????????%稀疏度
SNR=-20;??????????????????????????????????%信噪比
%f=8.62e12;???????????????????????????????%載波頻率
%c=3e8;???????????????????????????????????%傳播速度
%lamda=c/f;???????????????????????????????%波長
dt=0.5;????????????????????????????%發(fā)射陣元間距取半個波長
dr=0.5;????????????????????????????%接收陣元間距取半個波長
%%%%%有3個目標(biāo)方位2°3°4°
theat0=[-5?2?7];??????????????????????????%目標(biāo)與陣列的夾角
%d=[1e9?2e8?3e10];????????????????????????%目標(biāo)距離坐標(biāo)原點即第一個發(fā)射天線的距離
beata=[5?3?4];???????????????????????????%目標(biāo)的反射系數(shù)
Q=length（theat0）;
s0=hadamard（L）;??????????????????????????%產(chǎn)生一組正交信號
for?m=1:M
????s（:m）=s0（:m）;??????????????????????%M個發(fā)射信號源
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]‘）;%發(fā)射陣列引導(dǎo)矢量
????b（:q）=exp（-j*2*pi*dr*sin（theat0（q）*pi/180）*[0:N-1]‘）;%接收陣列引導(dǎo)矢量
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]‘）;??????%發(fā)射陣列引導(dǎo)矢量
????bb（:p）=exp（-j*2*pi*dr*sin（theat（p）*pi/180）*[0:N-1]‘）;??????%接收陣列引導(dǎo)矢量
????psi（:p）=reshape（（bb（:p）*aa（:p）.‘*s.‘）.‘N*L1）;??%稀疏域的基矩陣???????????????????
end

KK