麦克风阵列信号处理

在matlab下进行麦克风阵列信号的仿真系统，适用于近场

%%%%%基于均匀线阵的麦克风近场MUSIC算法
%%% 黄祖辉  厦门大学 通信过程系
%%% 2016-12-29
%%% Copyright@2016
clc
clear all
close all
M=16;                          %阵元个数
CenterFrequence=1700;%声波频率
SoundSpeed=340;         %声速
Lamda=SoundSpeed/CenterFrequence;%波长
ArraySpace=Lamda/2;%%%d
R_limit=2*ArraySpace^2/Lamda;%%%r=0.1m
K=2*pi/Lamda;
TestTime=2000;      %%%%快拍数
snr=25;                   %信噪比
%%%%%  signal  perferce

DistanceSignal=[1 2 3];
ThetSignal=[10 20 30];

%%%%%% Dmi  Ar_theta
Dmi=zeros（length（ThetSignal）M）;
for xm=1:M
    for ri=1:length（ThetSignal）
       Dmi（rixm）=sqrt（（DistanceSignal（ri）*sin（ThetSignal（ri）*pi/180）-ArraySpace*（xm-1））^2+（DistanceSignal（ri）*cos（ThetSignal（ri）*pi/180））^2）;
       Ar_theta（rixm）=DistanceSignal（ri）*exp（-j*K*（Dmi（rixm）-DistanceSignal（ri）））/DistanceSignal（ri）;
    end
end
%%%%%
Signal=sin（randn（length（ThetSignal）TestTime）*10）;
X0=Ar_theta‘*Signal;
X=awgn（X0snr‘measured‘）;
R=X*X‘/TestTime;

%%%%%
[emzm]=eig（R）; 
[zm1pos1]=max（zm）; 
 for l=1:length（ThetSignal）
    [zm2pos2]=max（zm1）; 
    zm1（:pos2）=[]; 
    em（:pos2）=[]; 
 end 

 %%%%%
 rangR=0.05:0.05:4.5;
 rangThet=1:90;
 for rr=1:length（rangR）
     for rt=1:length（rangThet）
         for rm=1:M
             D=sqrt（（rangR（rr）*sin（rangThet（rt）*pi/180）-ArraySpace*（rm-1））^2+rangR（rr）*cos（rangThet（rt）*pi/180））;
             B（rm）=rangR（rr）*exp（-j*K*D）/D;
         end
         A=B‘;
         P（rrrt）=（A‘*A）/（A‘*em*em‘*A）;
     end
 end
 
mesh（rangRrangThet20*log10（abs（P）））
xlabel（‘distance（m）‘）;
ylabel（‘thet‘）
zlabel（‘power （DB）‘）
title（‘麦克风阵列：快拍数2000近场MUSIC算法方向和距离参数距离估计‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1701  2016-12-19 10:23  matla源程序\liner_arrymaikf.m

     文件       1153  2016-12-19 10:24  matla源程序\main_2d_doa_music.m

     文件         87  2016-12-19 10:26  matla源程序\readme.txt

     目录          0  2017-01-04 16:22  matla源程序

----------- ---------  ---------- -----  ----

                 2941                    4