【正文】 tion [J] .IEEE Trans AP, 1986, 34 (2):276—280.[3] Kundu D .Modified MUSIC algorithm for estimating DOA of signals [J].Signal Processing, 1996, 48(1):8590. [4] HayKin S, Keily J P. Some aspects of arrays signal processing [J].IEEE Proc, 1992, 139(1):100124.[5] Godara L C. Application of antenna arrays to mobile munication, part Ⅱ: Beam forming and directionofarrival consideration [J].Proc, IEEE, 1997, 85(8):11951245.[6] Zhang Q, Wang K M, Statistical analysis of the performance if information theoretic criteria in the detection of the number of signals in array processing [J]. IEEE trans. On ASSP, 1989, 37(10):15571567.[7] H. krim, . Two decades of array signal processing research. [J]. IEEE Magazine, 1996, (7):.[8] . Naidu , . Direction of Arrival estimation in the presence of distributed noise sources: cumulant based approach.[ J].Acoust. Soc , 1995 , Pt .1 .[9] Hasan. DOA and frequency estimation using fast subspace algorithms[J]. Signal Processing,1999 , Vol, 4962 .[10] Rao B D, Hari K V S. Performance analysis of RootMUSIC. IEEE Trans. On ASSP, 1989, 37(12):19391949.[11] 王永良，[M].北京：清華大學(xué)出版社，2004.[12] 劉剛，呂新華，[J].微計(jì)算機(jī)信息，2006年第22卷第43期.[13] .[14] [M].:科學(xué)出版社，1997.[15] [D].廣州:華南師范大學(xué)，2002.[16] [D].華中科技大學(xué)，2007.[17] 陳輝，[J].系統(tǒng)工程與電子技術(shù)，.2004, .[18] 李杰，[J].杭州電子工業(yè)學(xué)院學(xué)報(bào)，,Dec,2001.[19] [D].山東大學(xué)碩士學(xué)位論文，2006.[20] 陳俊仕，江連山，[M].清華大學(xué)出版社，2005.[21] 劉德樹，羅景青，[M].中國(guó)科學(xué)技術(shù)大學(xué)出版社，1997.[22] [M].西北工業(yè)大學(xué)出版社，1989.[23] 王曙，[J].上海交通大學(xué)學(xué)報(bào)，1999.[24] [D].西安電子科技大學(xué)，1994.[25] 郭躍，[J].電子學(xué)報(bào)，2007年09期.[26] 曾慶棟，.[J].科技信息，2008.[27] 張賢達(dá)，.[M].國(guó)防工業(yè)出版社，2000.[28] 張賢達(dá)，現(xiàn)代信號(hào)處理[M].北京：清華大學(xué)出版社，1996.[29] 韓芳明，張守宏，[J].,30(5)：5 85589.[30] 邢超，[J].科技信息2007.[31] [M].清華大學(xué)出版社，2001.[32] 樓順天，李博菡. 基于MATLAB的系統(tǒng)分析與設(shè)計(jì)——信號(hào)處理[M]. 西安: 西安電子科技大學(xué)出版社, 2000.[33] 張志涌，[M].北京航空航天大學(xué)出版社，2010.[34] 王正林，[M].北京：電子工業(yè)出版社，2006.[35] 侯友國(guó)，郭偉，[J]. , ,.[36] 王進(jìn)，趙擁軍，[J].計(jì)算機(jī)工程，.附錄附錄一：MUSIC 算法MATLAB仿真源代碼clc clear all format long %將數(shù)據(jù)顯示為長(zhǎng)整型科學(xué)計(jì)數(shù)N=200。%快拍數(shù)doa=[20 60]/180*pi。 %信號(hào)到達(dá)角w=[pi/4 pi/3]39。%信號(hào)頻率M=10。%陣元數(shù)P=length(w)。 %信號(hào)個(gè)數(shù)lambda=150。%波長(zhǎng)d=lambda/2。%陣元間距 snr=20。%信噪比B=zeros(P,M)。 %創(chuàng)建一個(gè)P行M列的0矩陣for k=1:P B(k,:)=exp(j*2*pi*d*sin(doa(k))/lambda*[0:M1])。 %矩陣賦值end B=B39。 xx=2*exp(j*(w*[1:N]))。 %仿真信號(hào)x=B*xx。 x=x+awgn(x,snr)。%加入高斯白噪聲R=x*x39。 %數(shù)據(jù)協(xié)方差矩陣 [U,V]=eig(R)。 %求R的特征值和特征向量UU=U(:,1:MP)。 %估計(jì)噪聲子空間theta=90::90。 %%譜峰搜索for ii=1:length(theta) AA=zeros(1,length(M))。 for jj=0:M1 AA(1+jj)=exp(j*2*jj*pi*d*sin(theta(ii)/180*pi)/lambda)。end WW=AA*UU*UU39。*AA39。 Pmusic(ii)=abs(1/ WW)。 end Pmusic=10*log10(Pmusic/max(Pmusic))。 %空間譜函數(shù)plot(theta,Pmusic,39。k39。) xlabel(39。角度 \theta/degree39。) ylabel(39。譜函數(shù)P(\theta) /dB39。) title(39。MUSIC算法的DOA估計(jì)譜39。)grid on 附錄二：MUSIC算法DOA估計(jì)與陣元數(shù)的關(guān)系仿真源代碼clc clear all format long %將數(shù)據(jù)顯示為長(zhǎng)整型科學(xué)計(jì)數(shù)N=200。 %快拍數(shù)doa=[20 60]/180*pi。 %信號(hào)到達(dá)角w=[pi/4 pi/3]39。 %信號(hào)頻率M1=10。 %陣元數(shù)M2=50。M3=100。P=length(w)。 %信號(hào)個(gè)數(shù)lambda=150。 %波長(zhǎng)d=lambda/2。 %陣元間距snr=20。 %信噪比B1=zeros(P,M1)。 B2=zeros(P,M2)。B3=zeros(P,M3)。for k=1:P B1(k,:)=exp(j*2*pi*d*sin(doa(k))/lambda*[0:M11])。 %矩陣賦值B2(k,:)=exp(j*2*pi*d*sin(doa(k))/lambda*[0:M21])。 B3(k,:)=exp(j*2*pi*d*sin(doa(k))/lambda*[0:M31])。 end B1=B139。 B2=B239。 B3=B339。 xx=2*exp(j*(w*[1:N]))。 %仿真信號(hào)x1=B1*xx。 x2=B2*xx。 x3=B3*xx。 x1=x1+awgn(x1,snr)。 %加入高斯白噪聲x2=x2+awgn(x2,snr)。x3=x3+awgn(x3,snr)。R1=x1*x139。 %數(shù)據(jù)協(xié)方差矩陣R2=x2*x239。R3=x3*x339。[U1,V1]=eig(R1)。 %求R的特征值和特征向量[U2,V2]=eig(R2)。 [U3,V3]=eig(R3)。 UU1=U1(:,1:M1P)。 。 %估計(jì)噪聲子空間UU2=U2(:,1:M2P)。 UU3=U3(:,1:M3P)。 theta=90::90。 %%譜峰搜索for ii=1:length(theta) AA1=zeros(1,length(M1))。 for jj=0:M11 AA1(1+jj)=exp(j*2*jj*pi*d*sin(theta(ii)/180*pi)/lambda)。end WW1=AA1*UU1*UU139。*AA139。 Pmusic1(ii)=abs(1/ WW1)。 end for ii=1:length(theta) AA2=zeros(1,length(M2))。 for jj=0:M21 AA2(1+jj)=exp(j*2*jj*pi*d*sin(theta(ii)/180*pi)/lambda)。end WW2=AA2*UU2*UU239。*AA239。 Pmusic2(ii)=abs(1/ WW2)。 end for ii=1:length(theta) AA3=zeros(1,length(M3))。 for jj=0:M31 AA3(1+jj)=exp(j*2*jj*pi*d*sin(theta(ii)/180*pi)/lambda)。end WW3=AA3*UU3*UU339。*AA339。 Pmusic3(ii)=abs(1/ WW3)。 end Pmusic1=10*log10(Pmusic1/max(Pmusic1))。 %空間譜函數(shù)Pmusic2=10*log10(Pmusic2/max(Pmusic2))。Pmusic3=10*log10(Pmusic3/max(Pmusic3))。plot(theta,Pmusic1,39。k39。,39。LineWidth39。,)hold onplot(theta,Pmusic2,39。k39。,39。LineWidth39。,) hold onplot(theta,Pmusic3,39。39。,39。LineWidth39。,) hold offxlabel(39。角度 \theta/degree39。) ylabel(39。譜函數(shù)P(\theta) /dB39。) title(39。MUSIC算法的DOA估計(jì)譜39。)grid on 附錄三：MUSIC算法DOA估計(jì)與陣元間距的關(guān)系仿真源代碼clc clear all format long %將數(shù)據(jù)顯示為長(zhǎng)整型科學(xué)計(jì)數(shù)N=200。%快拍數(shù)doa=[20 60]/180*pi。 %信號(hào)到達(dá)角w=[pi/4 pi/3]39。%信號(hào)頻率M=10。%陣元數(shù)P=length(w)。 %信號(hào)個(gè)數(shù)lambda=150。%波長(zhǎng)d=lambda/6。%陣元間距 snr=20。%信噪比B=zeros(P,M)。 %創(chuàng)建一個(gè)P行M列的0矩陣for k=1:P B(k,:)=exp(j*2*pi*d*sin(doa(k))/lambda*[0:M1])。 %矩陣賦值end B=B39。 xx=2*exp(j*(w*[1:N]))。 %仿真信號(hào)x=B*xx。 x=x+awgn(x,snr)。%加入高斯白噪聲R=x*x39。 %數(shù)據(jù)協(xié)方差矩陣 [U,V]=eig(R)。 %求R的特征值和特征向量UU=U(:,1:MP)。 %估計(jì)噪聲子空間theta=90::90。 %%譜峰搜索for ii=1:length(theta) AA=zeros(1,length(M))。 for jj=0:M1 AA(1+jj)=exp(j*2*jj*pi*d*sin(theta(ii)/180*pi)/lambda)。end WW=AA*UU*UU39。*AA39。 Pmusic(ii)=abs(1/ WW)。 end Pmusic=10*log10(Pmusic/max(Pmusic))。 %空間譜函數(shù)plot(theta,Pmusic ,39。k39。,39。linewidth39。,)hold onclc clear all format long %將數(shù)據(jù)顯示為長(zhǎng)整型科學(xué)計(jì)數(shù)N=200。%快拍數(shù)doa=[20 60]/180*pi。 %信號(hào)到達(dá)角w=[pi/4 pi/3]39。%信號(hào)頻率M=10。%陣元數(shù)P=length(w)。 %信號(hào)個(gè)數(shù)lambda=150。%波長(zhǎng)d=lambda/2。%陣元間距 snr=20。%信噪比B=zeros(P,M)。 %創(chuàng)建一個(gè)P行M列的0矩陣for k=1:P B(k,:)=exp(j*2*pi*d*sin(doa(k))/lambda*[0:M1])。 %矩陣賦值end B=B39。 xx=2*exp(j*(w*[1:N]))。 %仿真信號(hào)x=B*xx。 x=x+awgn(x,snr)。%加入高斯白噪聲R=x*x39。 %數(shù)據(jù)協(xié)方差矩陣 [U,V]=eig(R)。 %求R的特征值和特征向量UU=U(:,1:MP)。 %估計(jì)噪聲子空間theta=90::90。 %%譜峰搜索for ii=1:length(theta) AA=zeros(1,length(M))。 for jj=0:M1 AA(1+jj)=exp(j*2*jj*p