%壓縮感知應用于DOA估計的SMV模型下的算法
%BY?胡德順
%2015.4.17

close?all;
clear?all;
clc;
tic

%創建信號的數據采集矩陣
k=5;
m=100;
n=1;
theta=[-40?-20?11?50?80]*pi/180;
d=0.5;
snr=10;
A=exp（-j*2*pi*d*（0:m-1）‘*sin（theta））;
S=sqrt（10.^（snr/10））*exp（j*2*pi*randn（k1））;
X=A*S+（1/sqrt（2））*（randn（m1）+j*randn（m1））;

%等角度劃分，產生間距為1度的測量矩陣
Theta=（-90:0.1:90）*（pi/180）;
A1=exp（-j*2*pi*d*（0:m-1）‘*sin（Theta））;
%size（A1）

%采用OMP算法進行重構，進而得知角度信息
%測量矩陣為A1，觀測信號為X，信號稀疏度為k。
S0=ones（size（Theta））;
r0=X;
%error=;
ompin