资源简介

压缩感知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

评论

共有 条评论