MUSIC算法在DSP實(shí)現(xiàn)

　　隨著信號處理技術(shù)的發(fā)展，陣列信號處理逐漸表現(xiàn)出優(yōu)越性，并被廣泛應(yīng)用，以提高通信系統(tǒng)對信號處理的效率。毫米波數(shù)字陣列近年來一直是人們研究的熱點(diǎn)。波達(dá)方向（Direction of　Arrival， DOA）估計是陣列信號處理的一個基本問題，它以譜估計技術(shù)為基礎(chǔ)，通過求解空間譜函數(shù)得到目標(biāo)準(zhǔn)確方位，對空間目標(biāo)精確定位。多重信號分類（Multiple Signal Classification，MUSIC）法是經(jīng)典的超分辨率空間譜估計方法。本文主要基于毫米波段進(jìn)行MUSIC算法的DSP實(shí)現(xiàn)。

　　近年來，在MUSIC算法優(yōu)化領(lǐng)域有學(xué)者分別提出了root-MUSIC算法、降維MUSIC算法等優(yōu)化算法?，F(xiàn)有文獻(xiàn)中優(yōu)化方法主要集中在協(xié)方差矩陣簡化或者譜峰搜索優(yōu)化等MUSIC算法本身的簡化上，但是具體實(shí)現(xiàn)過程中，MUISC算法的DSP實(shí)現(xiàn)仍有優(yōu)化空間。

　　本文DSP實(shí)現(xiàn)主要基于TMS320VC5502（C-5502）。C5502的函數(shù)庫中不包含特征值的求解也不能直接進(jìn)行復(fù)數(shù)運(yùn)算。本文采用一種QR變換法求解特征值。該方法把協(xié)方差矩陣轉(zhuǎn)換成上H矩陣并進(jìn)行QR分解。最后所得的對角陣的對角元素即是扔方差矩陣的特征值。同時本文把復(fù)矩陣轉(zhuǎn)換為實(shí)矩陣計算特征值，對算法進(jìn)行優(yōu)化，減少計算量。最后本文通過優(yōu)化前后對比和MALTAB結(jié)果驗證了本文優(yōu)化算法的可行性和實(shí)用性。