MSK信號(hào)檢測(cè)識(shí)別的FPGA實(shí)現(xiàn)
采用MSK 調(diào)制的跳頻通信具有主瓣能量集中、旁瓣衰落滾降快、頻譜利用率高和抗干擾能力強(qiáng)等優(yōu)點(diǎn),在軍事通信中應(yīng)用廣泛。如美軍現(xiàn)役的聯(lián)合戰(zhàn)術(shù)信息分發(fā)系統(tǒng)采用的通信信號(hào),工作帶寬969~1 206 MHz,跳頻速率為70000 多跳/ s, 單個(gè)頻點(diǎn)駐留時(shí)間約為13 s,信號(hào)持續(xù)時(shí)間* s, 總共有51個(gè)間隔為3 MHz 的信道,碼速率為5 MHz。已知在該工作頻段內(nèi)主要還存在單頻、窄帶調(diào)幅和線性調(diào)頻等信號(hào)。為了準(zhǔn)確截獲并識(shí)別目標(biāo)信號(hào),針對(duì)此信號(hào)環(huán)境設(shè)計(jì)了一種MSK 信號(hào)檢測(cè)識(shí)別方法,并使用FPGA進(jìn)行了設(shè)計(jì)實(shí)現(xiàn)。
本文引用地址:http://m.butianyuan.cn/article/150603.htm 1 算法設(shè)計(jì)
1.1 寬帶跳頻信號(hào)實(shí)時(shí)檢測(cè)算法
用現(xiàn)代技術(shù)來實(shí)現(xiàn)寬帶數(shù)字化接收的一個(gè)實(shí)用的方法是通過信道化技術(shù),實(shí)現(xiàn)信道化通常的方法是采用快速傅里葉變換(FFT)。 在近代電信裝備和各類控制系統(tǒng)中,濾波器應(yīng)用極為廣泛;在所有的電子部件中,使用最多,技術(shù)最復(fù)雜要算濾波器了。濾波器的優(yōu)劣直接決定產(chǎn)品的優(yōu)劣,所以,對(duì)濾波器的研究和生產(chǎn)歷來為各國(guó)所重視。
設(shè)計(jì)法更容易實(shí)現(xiàn),因?yàn)镕FT所需要的運(yùn)算量更少。
某個(gè)由FFT運(yùn)算輸出的頻率分量,可以看成輸入信號(hào)與某個(gè)脈沖函數(shù)的卷積。為了處理一個(gè)連續(xù)的輸入信號(hào)。必須在不同時(shí)刻對(duì)各段數(shù)據(jù)進(jìn)行FFT處理。通常,起始點(diǎn)記為n = 0,數(shù)據(jù)段可以滑動(dòng)M點(diǎn),相應(yīng)的FFT可以寫成:
M 的值必須隨著輸入信號(hào)連續(xù)變化,這種運(yùn)算也叫作短時(shí)傅里葉變換(STFT)。
FFT的長(zhǎng)度和重疊點(diǎn)數(shù)是非常重要的參數(shù),這些參數(shù)與最小脈寬和頻率分辨率
在電視工業(yè)中,分辨率是用清晰度來度量,單位是電視線(TVLine)。
有關(guān),它們決定了接收機(jī)的靈敏度。若FFT的長(zhǎng)度為N,信號(hào)的采樣頻率為f s,那么經(jīng)FFT計(jì)算后,信號(hào)的頻率分辨率為:
數(shù)據(jù)重疊點(diǎn)數(shù)決定了時(shí)間分辨率和處理的最短脈寬,數(shù)據(jù)重疊率越高,則時(shí)間分辨率越高。
本設(shè)計(jì)中使用STFT 的方法實(shí)現(xiàn)一個(gè)粗測(cè)頻引導(dǎo)數(shù)字接收機(jī),為覆蓋整個(gè)跳頻帶寬,采用700 MHz采樣率對(duì)目標(biāo)信號(hào)進(jìn)行采樣,粗測(cè)頻引導(dǎo)精度在1 MHz以內(nèi),因此FFT長(zhǎng)度選擇為1 024 點(diǎn),數(shù)據(jù)重疊率50%,保證時(shí)間分辨率在1 s 之內(nèi)。
利用粗測(cè)頻接收單元的檢測(cè)結(jié)果去引導(dǎo)一個(gè)精測(cè)單元,可以對(duì)目標(biāo)信號(hào)進(jìn)行更為精確的測(cè)量和識(shí)別。精測(cè)單元采用數(shù)字正交下變頻的實(shí)現(xiàn)方法,數(shù)字本振頻率隨粗測(cè)引導(dǎo)結(jié)果而設(shè)置,變頻后的帶寬依據(jù)目標(biāo)信號(hào)而確定。1.2 MSK 信號(hào)識(shí)別算法
MSK 信號(hào)可以寫成:
式中,Tb 為碼元周期,θk 是第k 個(gè)碼元的相位常數(shù),取值為nπ, Pk 為二進(jìn)制雙極性碼元,取值為 1,所以MSK 信號(hào)相位分段線性變化,每個(gè)碼元周期內(nèi)相對(duì)前*元載波相位上升或下降π/2
根據(jù)三角函數(shù)展開,設(shè)θk 起始參考值為0,得:
又根據(jù)Ik= 1, Qk= 1, 令f L= f c-14Tb, f H=f c+14Tb,MSK 信號(hào)經(jīng)過平方環(huán)可得:
可知,MSK 信號(hào)經(jīng)過平方運(yùn)算后,含有2f L 和2fH兩個(gè)離散頻率分量,反映在FFT幅度譜上,在這2 個(gè)頻點(diǎn)上存在2 個(gè)明顯的譜峰,2 個(gè)譜峰的距離為1 個(gè)碼元速率,并且距離2 倍載頻處均為碼元速率的一半。這些特征是此頻段內(nèi)其他信號(hào)不具備的,可以根據(jù)這些特征對(duì)MSK 進(jìn)行有效的識(shí)別。因此,將經(jīng)過正交下變頻的信號(hào)進(jìn)行平方運(yùn)算,確定信號(hào)出現(xiàn)后對(duì)其平方的結(jié)果進(jìn)行FFT處理計(jì)算幅度譜,再對(duì)譜峰點(diǎn)進(jìn)行分析即可完成MSK 信號(hào)的識(shí)別。
2 算法的FPGA 實(shí)現(xiàn)
2.1 算法實(shí)現(xiàn)
跳頻MSK 信號(hào)檢測(cè)識(shí)別的實(shí)現(xiàn)框圖如圖1 所示,采用FPGA 實(shí)現(xiàn),包括短時(shí)傅里葉變換(STFT) 粗測(cè)頻引導(dǎo)、數(shù)字正交下變頻、平方運(yùn)算和幅度譜分析等主要模塊。STFT 粗測(cè)頻引導(dǎo)在寬帶條件下進(jìn)行實(shí)時(shí)的信號(hào)檢測(cè)和頻率粗測(cè),測(cè)量的結(jié)果引導(dǎo)數(shù)字正交下變頻模塊,對(duì)信號(hào)進(jìn)行變頻、濾波和抽取,得到低采樣速率的零中頻數(shù)據(jù),平方運(yùn)算模塊對(duì)零中頻數(shù)據(jù)進(jìn)行平方處理,在確定存在信號(hào)后,對(duì)平方運(yùn)算模塊的輸出進(jìn)行FFT運(yùn)算得到信號(hào)的幅度譜,通過幅度譜分析模塊得到最終的識(shí)別結(jié)果。
評(píng)論