%%%%%基于均勻線陣的麥克風(fēng)近場MUSIC算法
%%%?黃祖輝??廈門大學(xué)?通信過程系
%%%?2016-12-29
%%%?Copyright@2016
clc
clear?all
close?all
M=16;??????????????????????????%陣元個數(shù)
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;??????%%%%快拍數(shù)
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（‘麥克風(fēng)陣列：快拍數(shù)2000近場MUSIC算法方向和距離參數(shù)距離估計(jì)‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????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