music--doa.txt

演示传统波束形成及ＭＵＳＩＣ方法的ＤＯＡ估计．信号为宽带信号源．

%WIDEEX1 用传统的方法及MUSIC方法估计宽带信号的DOA.
%
%Start:981212 Svante Bjoklund(svabj).
%
********************************************
*
%--------------------------------------------%
%参数
%--------------------------------------------%
lambda=0.1;           %波长.
D=0.45*lambda;        %阵元间距
T=16;                 %快拍数
K=12;                 %数字通道数
theta=d2r([25 40]).'; %目标角度,目标数由角度个数给出
%phi=zeros(size(theta));%目标角度
%SNR=[10 5]';          %每个阵元的信噪比dB
%alpha=d2r([0 16]).';
%dalpha=d2r([34-13]).';
%dist=Inf*ones(size(theat));
%tgtModel='const';
%noiseModel='rndnw';
%noiseModel='rndn';
MMu=2;
SNR=0;
%-------------------------------------------%
%运行命令.
%-------------------------------------------%
ant=defant('isotropULA',[K,D]);
%[sigOut,sourceSigTime,noiseSigTime]=simwidebandkernel(srcType,antenna,carrierFreq,bandwidth,no Times,noTrials,srcSignals,srcDoas,srcPowers,srcCorr,orthSig,propSpeed)
antenna=ant;
noSrc=2;
propSpeed=speedoflight;%[m/s].
carrierFreq=propSpeed/lambda%[Hz].(3 GHz)
bandwidth=2e9%[Hz].
noTimes=T;
noTrials=1;
srcSignals=[];
srcDoas=theta.';
srcPowers=ones(noSrc,1);
noiseType='randn';
noisePower=srcPowers(1)/(10^(SNR/10))
srcCorr=[];
orthSig=[];
spect1=sdoaspc('music',sig,[],MMu);
R=ecorrm(sig);
spect2=sdoaspc('cdf',R);
splot2(spect1,",spect2)
legend('Music','Conventional beamforming')
plotxline(r2d(srcDoas),'m','--')
%printfm wideex1a1.eps[]12