压缩感知CS的DOA代码

压缩感知CS的DOA代码压缩感知CS的DOA代码压缩感知CS的DOA代码压缩感知CS的DOA代码压缩感知CS的DOA代码压缩感知CS的DOA代码

%function Farfield_DOA （resignal）
clc;
clear;
close all;
Freq =700;                    %Frequency
snr =40;                      %Source to noise ratio
M =12;                          %The number of microphones
Cspeed = 340;                    % The speed of sound
Lambda = Cspeed/340;             %Wavelength
d = Lambda/2;                     %Array elements spacing
Dirsources = [2035];             %Direction of sound sources
Numsources = length（Dirsources）;      %Number of sound sources
 
for q = 1:Numsources
    a（:q）=exp（-1i * 2 * pi * Freq * d *（0:M-1）‘*sin（Dirsources（q）*pi/180）/Cspeed）;     %
end
 
Energy = [12];                %Sources energy
y = a*Energy‘;                 %Received signal of array
Y = awgn（ysnr‘measured‘）;     %Adding noise
 
grid = 0.1;                       %The value of grid node  
dirmin = -90;                   %The min and max DOA value of the source plane
dirmax = 90;
Dirmicrogrid = dirmin:grid:dirmax;       %Direction between microphone and grid node
Numgrid = length（Dirmicrogrid）;          %Number of grid
 
A = zeros（MNumgrid）;                    %Measurement matrix （1/sqrt（M））*e
for j = 1:Numgrid
    A（:j）=exp（-1i * 2 * pi * Freq * d * （0:M-1）‘ * sin（Dirmicrogrid（j） * pi / 180） / Cspeed）;
end
tic
resignal = CS_OMP（YA） ;       %Recover using OMP
resignal = real（resignal‘）;
toc
 
nonsource = nonzeros（resignal）
 
figure（1）;
 
plot（Dirmicrogridresignal‘k*-‘）                                 %  重建信号
hold on
plot（DirsourcesEnergy‘ks‘）
title（[‘SNR = ‘num2str（snr）]）
xlabel（‘两个信源方向角‘）
ylabel（‘两个信源幅度值‘）
dirction = （find（resignal） -1） *grid -90

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-10-31 06:56  CS_DOA\
     文件        1621  2017-10-31 06:50  CS_DOA\CS_DOA.m
     目录           0  2017-10-31 06:57  __MACOSX\
     目录           0  2017-10-31 06:57  __MACOSX\CS_DOA\
     文件         120  2017-10-31 06:50  __MACOSX\CS_DOA\._CS_DOA.m
     文件         497  2017-10-31 06:48  CS_DOA\CS_OMP.m
     文件         120  2017-10-31 06:48  __MACOSX\CS_DOA\._CS_OMP.m