果變壞，當(dāng)信號(hào)高度和關(guān)時(shí)估計(jì)性能同樣會(huì)變壞，且對(duì)所設(shè)的參數(shù)有較高的要求，少量的誤差也會(huì)導(dǎo)致算法的失。在ESPRT算法中隨著信噪比的增加，均方誤差越小，DOA估計(jì)效果越好。在陣元數(shù)目不同，其他條件相同的情況下，陣元數(shù)目越大，均方誤差越小，ESPRIT算法的估計(jì)精度越高。在條件相同的情況下，相鄰信號(hào)（以10。為例）的均方差與信噪比關(guān)系隨著角度的垢加而性能降低，，ESPRIT算氵去對(duì)相干信號(hào)的DOA估計(jì)失效。而G算氵去，不僅計(jì)算量比較小，而且保證了算法的精度，但當(dāng)信號(hào)高度相關(guān)時(shí)性能仍然變壞。在G“SE算氵去中隨著信噪比的增加，均方誤差越小，D“估計(jì)效果越好。在陣元數(shù)目不同，其他條件相同的情況下，陣元數(shù)目越大，均方誤差越小，G算法的估計(jì)精度越高。在條曾相同的情況下1相鄰信號(hào)（以釅為例）的均方差與信噪比關(guān)系隨著角度的增師而陸能降低，G算氵去對(duì)相干信號(hào)的D“估計(jì)失效：·3．3實(shí)驗(yàn)33．3·1非平穩(wěn)噪聲下，三種算法對(duì)DOA估計(jì)影響的matiab仿真結(jié)果uu漶問圖34非平穩(wěn)噪聲下s芘算法對(duì)D謐估計(jì)的影響TT非半產(chǎn)下因回厲耳圖35非平穩(wěn)噪聲下ESPRIT算法對(duì)DOA估計(jì)的影響GEESE韓+噪盧下仿真圖36韓平穩(wěn)噪聲下GEESE算法對(duì)DOA估計(jì)的影響3．32色噪聲噪聲下加三種算法對(duì)DOA估計(jì)影響的m仿真結(jié)果圖37色噪聲下“SIC算法對(duì)DOA估計(jì)的影響ESPRT色下懣向仿自圖3．8色噪聲下ESPT算法對(duì)DOA估計(jì)的影響0030卩0舜際角0圖39色噪聲下G“算法對(duì)DOA估計(jì)的影響附錄1第1題仿真程序邑as巳allM爿1245310]；％M為接收機(jī)個(gè)數(shù)N=[15101％為射頻脈沖數(shù)Estill羽tez〔3，(l)，figurefo「s11=060削1err=Cjforcon11100蕘℃1“艸e艸100〕，N(l),sm〗％計(jì)算j個(gè)接收機(jī)和]個(gè)射頻脈沖數(shù)的最小錯(cuò)誤概率allem—allen*sum(En?or)endEsttmate(i十I)=allet22000；endendhodonp屨伍山m0〕漁p乍sum0〗、m〗xlabelCSNR—O、60dB〗2ylabe1Ctitle(ll6llegendC1 110gridonholdoff,endlemor(theslzeofxISnotequalwtthendsu=0,fori=lforFIbendendda=su,%mISthenunlberCfrecewera=raylmdCGatns);%Signasrecewedfori=lmfor .n13endendPeoutll=zeros(m,y,Pcout1),Pcout211),%thestructureofthemacformforFIendend%gettheabsolutevalueofthemacPcoutl1=absQPcout11),end%getthedecestonofthedetectlon2MUSlCforfor14）艸畔@*2，產(chǎn)d*sm〔d。aD／地0M一1以end％％％陣列流形矩陣％$：1儼020嚴(yán)0／t(2)〕*ndn（4，N廾〕竹andi〕蠲〗％仿真信號(hào)〔隨機(jī)信號(hào)）惲10、@til的0以tn（tv0．N一1〗％仿真信號(hào)（正弦信號(hào)〕％s一1儼〔兜0）*(i*w彗0N-l〗％仿真信號(hào)（指數(shù)信號(hào)》％$：1儼020）“t。豐0過．(1)，％仿真信號(hào)（鋸齒汲信號(hào)）％s一2慳即儼（*卩N]》％生成信號(hào)x*葉（11t0尹(randn(h《N廾J*randnO.《〗％師了高斯白噪聲后的陣列接收信號(hào)x+awgn區(qū)〗％加了高斯白噪聲后的陣列接收信號(hào)％Q：drag(v〗％噪聲協(xié)方差矩陣對(duì)角線值相同%v=randn(1,M),％Q到rag@〗％噪聲協(xié)方差矩陣對(duì)角線值不同Qnd囗00％噪聲協(xié)方差矩陣不為對(duì)角陣R斗Q；％數(shù)據(jù)協(xié)方差矩陣％％％％以下是的程序Il-J,V]=eLg(R),I-nJ=U()1M-P);theta=-90050，％％譜峰搜索一1ength(theta)：2艸卩eng山00，forjj=0M1AA(1+jJ)艸KP()j*2劫，1*dn曲婦()D/180，以l)，endms11片s()]上*尹[甫*是D，endnuc：10*適g10@1111巧i血ax(Pmuc〕+，plot(thetaPmusic)xlabeICAngIe\theta/degree)ylabel(P0theta〕劌B》title('MUSIC非相關(guān)源測(cè)向仿真《foniz的13，．fi。nel蜘t',b引d是化n恤《隸書I,color,誕d〗%YITJSIC程序結(jié)束附錄3ESPRIT算法程序一]024，％％快拍數(shù)一10234習(xí)門，，％％信號(hào)到達(dá)角5w爿1的03的02的時(shí)％％信號(hào)頻率M：3％陳元數(shù)gu0％％信號(hào)個(gè)數(shù)，也可以用特征分解的大特征值數(shù)來決定一巧0，％波長(zhǎng)d一10％陣元間距“=15，％％信噪比％％陣列流形矩陣E=2艸$但M〕fopk=lPBOG〕艸xp()j*2，韓din（d“的）／0M-11),end％s：1儼細(xì)r/20氣肱cpt(2)〕氣tandn(),》0十〕〕dn(3bT》％仿真信號(hào)（隨機(jī)信號(hào)〕]儼（兜0〕un@阿0N-l〕，％仿真信號(hào)〔正弦信號(hào)）％$：1儼〔、/20*“的*w豐ON．1〕仿真信號(hào)（指數(shù)信號(hào)）％s：1儼〔/2叻*s“№〔耳0N一]0％仍真信號(hào)（鋸齒波信號(hào)）％$：2*儼（，1N力〗％生成信號(hào)x一+（1／sq2〕尹nd[]0《N）到懺andna《〕〕，x：葉“gn如“驍％加了高斯白噪聲后的陣列接收信號(hào)%v=[111]1111]*tandn％Q：diag回％噪聲協(xié)方差矩陣對(duì)角線值相同Q：di嘔@〕，％噪聲協(xié)方差矩陣對(duì)角線值不同％《ndn0％噪聲協(xié)方差矩陣不為對(duì)角陣R斗0％數(shù)據(jù)協(xié)方差矩陣％％以下是ESPRIT程序Rxx-R(1M一1，1M-I),%%%M1維的自相關(guān)函數(shù)k“0M一1，2M);%%%M1維的互相關(guān)函數(shù)b=[zeros(1tM2〕，eye心]，C囂刪“．mn(elg@二以〕尹eye讎．1〗CxyˉR?'-mm(a1gC?xx〗)*b,a=eig(cxx,cry),％找出最接近]的a值其對(duì)應(yīng)的角度即為al一abab《的一1〗fo「1=1Pc,d]=mm(al),al@〕產(chǎn)1瓿0d片1m0，0衛(wèi)>1disp(Theanglesofsignalsare)%else%disp(Theangeof重alis')O/oend一doa*180／doa=ADO盍：asin@ng！eOb〕01乩DO盍一s00A〗stem(doaDO盍，力：山b引00A實(shí)際角度〕ylabe1(DOA估計(jì)角度．〕刷e(Ek非相關(guān)源測(cè)向仿真《nt丷e，3，。ontweLght'氅01d，》ontnamel沁0巧儼，d〕，Idon0/00/0%%ESPRIT程序結(jié)束附錄4GEESE算法程序1024％％快拍數(shù)doa=[-30030]／180加1，％％信號(hào)利達(dá)角w：0開03開0200吲％％信號(hào)頻率M：8，％％陣元數(shù)g山（〕，％％信號(hào)個(gè)數(shù)，也可以用特征分解的大特征值數(shù)來決定Lambda=50％汲長(zhǎng)d=Lambda/2％陣元間距一]5％％信噪比％％陣列流形矩陣B“e0蠲，for=1PE(l'?〕艸即@*2，1叼L.n@。a00〕能ambda豐0M-IJ),end％％％陣列流形矩陣％1驢0／20》氣1/沖《紗*@nd《3對(duì)片ih?andn（3，N〕〗％仂真信號(hào)〔隨機(jī)信號(hào)〕s：10儼（/2帥*sm00N．一1讓％仿真信號(hào)〔正弦信號(hào)〕％s：]0氣／20XP儼w*[ON一]l)，％仿真信號(hào)（指數(shù)信號(hào)〕％1驢020awt。*0毖J)，％仿真信號(hào)〔鋸齒波信號(hào)）％s：一慳xp(l*@覲1l)〗％生成信號(hào)x=B*s％v：011]1110/0Q=diag(v);%噪聲方差矩陣為對(duì)角陣，且值相同=10000ndl〕0070rag回，％噪聲協(xié)方差年陣為對(duì)角陣，值不同%Q=00傭1ndn，蠲，％噪聲協(xié)方差矩陣不為對(duì)角陣R“+0％數(shù)據(jù)協(xié)方差知陣％％％％以下是GEESE的程序《U，到ig巛〗EI=U(13，5上《EZ：U〔24，5上《),E〕g〔壓1〗DOA=asm(angleE)fpD/pi*180,DOA=s01