本發(fā)明涉及MPSK信號解調(diào)的頻率同步技術(shù)，特別涉及一種適用于MPSK解調(diào)的快速精確頻率同步方法。

背景技術(shù)：

在無線數(shù)字通信系統(tǒng)中，載波頻率偏移是影響系統(tǒng)解調(diào)性能指標(biāo)的重要因素，為了降低誤碼率，提高通信系統(tǒng)的解調(diào)性能，必須在接收端進(jìn)行估計(jì)和校正。針對MPSK的頻率估計(jì)，一般采用FFT的方法進(jìn)行，但由于資源限制，一般FFT點(diǎn)數(shù)不宜超過8192點(diǎn)，過大的FFT不僅資源占用多，處理延時(shí)也大，影響了頻率補(bǔ)償?shù)男Ч?#xff1b;而FFT點(diǎn)數(shù)過小，又無法達(dá)到所需的頻率精度，使后續(xù)的鎖相環(huán)等帶寬過大，增加了復(fù)雜度，降低了通信系統(tǒng)的穩(wěn)定性。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的在于提供一種適用于MPSK解調(diào)的快速精確頻率同步方法，以解決現(xiàn)有的對MPSK的頻率估計(jì)在采用FFT方法時(shí)所存在的受資源限制，特別是在FFT點(diǎn)數(shù)過大時(shí)處理延時(shí)較大導(dǎo)致多普勒變化率累積，影響頻率補(bǔ)償；在FFT點(diǎn)數(shù)小時(shí)，無法滿足頻率精度需求的問題。

為實(shí)現(xiàn)上述目的，本發(fā)明提供了一種適用于MPSK解調(diào)的快速精確頻率同步方法，包括以下步驟：

S1：獲取MPSK接收信號，對所述MPSK接收信號采用M次方復(fù)乘去調(diào)制的方式進(jìn)行去調(diào)制處理，恢復(fù)得到調(diào)制載波信號，對所述調(diào)制載波信號進(jìn)行K點(diǎn)的FFT處理，并對FFT處理的結(jié)果進(jìn)行求模后選取其中的最大值得到粗估計(jì)頻率值；其中，K為2的整數(shù)次冪值；

S2：根據(jù)所述粗估計(jì)頻率值對所述MPSK接收信號的數(shù)據(jù)序列進(jìn)行頻率粗補(bǔ)償，并對頻率粗補(bǔ)償后的數(shù)據(jù)序列再次進(jìn)行M次方復(fù)乘去調(diào)制，然后進(jìn)行窄帶濾波處理，再進(jìn)行N倍抽取，最后通過求模后選取其中的最大值得到精細(xì)估計(jì)頻率值；

S3：將所述精細(xì)估計(jì)頻率值與所述粗估計(jì)頻率值進(jìn)行相加，得到最終的DDS頻率字，并采用該DDS頻率字對輸入的MPSK接收信號進(jìn)行精確頻率補(bǔ)償。

較佳地，采用第一級狀態(tài)機(jī)、第二級狀態(tài)機(jī)及第三級狀態(tài)機(jī)分別執(zhí)行所述步驟S1、S2及S3。

較佳地，采用第一級狀態(tài)機(jī)將獲取的所述MPSK接收信號中的復(fù)數(shù)數(shù)據(jù)存儲在RAM中，存儲的長度為FFT的長度，并由所述第一級狀態(tài)機(jī)進(jìn)行調(diào)制處理及求模后選取其中的最大值得到所述粗估計(jì)頻率值。

較佳地，所述第一級狀態(tài)機(jī)將所述粗估計(jì)頻率值發(fā)送至所述第二級狀態(tài)機(jī)，所述第二級狀態(tài)機(jī)讀取RAM中的MPSK接收信號的復(fù)數(shù)數(shù)據(jù)并進(jìn)行頻率粗補(bǔ)償、調(diào)制、窄帶濾波處理、抽取及求模后選取其中的最大值，得到精細(xì)估計(jì)頻率值。

較佳地，所述第二級狀態(tài)機(jī)將所述精細(xì)估計(jì)頻率值發(fā)送至所述第三級狀態(tài)機(jī)，所述第三級狀態(tài)機(jī)采用計(jì)算的DDS頻率字對外部輸入的MPSK接收信號進(jìn)行精確頻率補(bǔ)償。

較佳地，窄帶濾波處理時(shí)，窄帶濾波的帶寬小于采樣帶寬的1/N。

較佳地，所述FFT處理的K點(diǎn)的值為1024，且所述FFT處理選用基-4FFT運(yùn)算器執(zhí)行。

本發(fā)明采用新型的分時(shí)復(fù)用兩級FFT技術(shù)，可以用資源占用較小的1024點(diǎn)FFT實(shí)現(xiàn)65536點(diǎn)FFT的頻率分辨率和效果，大幅減少了程序資源占用，降低了處理時(shí)延，減少了多普勒頻率變化率對頻率估計(jì)的影響，提升了精確頻率估計(jì)的效果，對解調(diào)中的頻率精確同步的實(shí)現(xiàn)有著重要的實(shí)際意義。

本發(fā)明方法適用于快速精確同步系統(tǒng)，也適用于連續(xù)或者突發(fā)信號的MPSK解調(diào)系統(tǒng)。

附圖說明

圖1為本發(fā)明方法整體流程圖；

圖2為本發(fā)明優(yōu)選實(shí)施例使用的基-4FFT運(yùn)算器結(jié)構(gòu)示意圖；

圖3為本發(fā)明優(yōu)選實(shí)施例的QPSK解調(diào)的快速精確頻率同步方法流程圖。

具體實(shí)施方式

為更好地說明本發(fā)明，茲以一優(yōu)選實(shí)施例，并配合附圖對本發(fā)明作詳細(xì)說明，具體如下：

如圖1所示，本實(shí)施例提供的一種適用于MPSK解調(diào)的快速精確頻率同步方法，包括以下步驟：

S1：獲取MPSK接收信號，對所述MPSK接收信號采用M次方復(fù)乘去調(diào)制的方式進(jìn)行去調(diào)制處理，恢復(fù)得到調(diào)制載波信號，對所述調(diào)制載波信號進(jìn)行K點(diǎn)的FFT處理，并對FFT處理的結(jié)果進(jìn)行求模后選取其中的最大值得到粗估計(jì)頻率值；這里的K為2的整數(shù)次冪值。

具體地，本實(shí)施例中假設(shè)在接收端進(jìn)行載波同步處理時(shí)，用于載波頻偏估計(jì)的MPSK接收信號基帶表達(dá)式如下：

其中，ak為獨(dú)立同分布的等概率數(shù)據(jù)，fe為待估計(jì)得的載波頻率偏移誤差，T是符號周期，θ0是未知的載波相位，ωk為加性高斯白噪聲。

首先，對輸入星座圖為2π/M旋轉(zhuǎn)對稱形式的MPSK復(fù)數(shù)基帶信號rk采用M次方復(fù)乘去調(diào)制的方式進(jìn)行去調(diào)制處理，恢復(fù)出調(diào)制載波。對去調(diào)制后的信號進(jìn)行K點(diǎn)的FFT處理(其中K取2的整數(shù)次冪值，方便FFT的蝶形運(yùn)算)，并求模選大找出載波頻率中心位置，此時(shí)，FFT頻率估計(jì)范圍為采樣頻率的1/M，頻率分辨率為采樣頻率的1/(M×K)倍。FFT的運(yùn)算過程為：

其中，rk為接收基帶信號，X(k)為FFT處理后的數(shù)據(jù)序列。

其中為了減少FFT運(yùn)算延時(shí)，FFT長度不宜選的過長，本實(shí)施例中的FFT處理的K點(diǎn)的值為1024，且FFT處理選用基-4FFT運(yùn)算器執(zhí)行，基-4FFT運(yùn)算器的結(jié)構(gòu)如圖2所示。

S2：根據(jù)所述粗估計(jì)頻率值對所述MPSK接收信號的數(shù)據(jù)序列進(jìn)行頻率粗補(bǔ)償，并對頻率粗補(bǔ)償后的數(shù)據(jù)序列再次進(jìn)行M次方復(fù)乘去調(diào)制，然后進(jìn)行窄帶濾波處理，再進(jìn)行N倍抽取，最后通過求模后選取其中的最大值得到精細(xì)估計(jì)頻率值；

具體地，根據(jù)步驟S1中FFT結(jié)果進(jìn)行求模選大(求模后選取其中的最大值)得到粗估計(jì)頻率值，對原始數(shù)據(jù)信號序列進(jìn)行DDS頻率粗補(bǔ)償。

然后，對粗頻率補(bǔ)償后的數(shù)據(jù)序列再次進(jìn)行M次方復(fù)乘去調(diào)制，復(fù)乘后的序列進(jìn)行窄帶濾波處理，再進(jìn)行N倍抽取，以縮小頻率范圍，進(jìn)行精細(xì)頻率估計(jì)和補(bǔ)償。此時(shí)窄帶濾波在粗同步后進(jìn)行，由于前端頻率粗同步已經(jīng)把頻率誤差控制在采樣頻率的1/(M×K)倍以內(nèi)，只需使用一個(gè)帶寬小于后端抽取倍數(shù)N的窄帶低通濾波器即可，階數(shù)可選32階，以濾除M次方復(fù)乘后的噪聲和雜散，并降低N倍抽取后信號的混疊影響。經(jīng)過N倍抽取，精細(xì)頻率同步FFT的鑒頻分辨率縮小為采樣頻率的1/(M×K×N)倍，抽取后的K/N點(diǎn)數(shù)據(jù)后面補(bǔ)齊(K-K/N)個(gè)零，以湊齊K點(diǎn)數(shù)據(jù)序列進(jìn)行精確同步的K點(diǎn)FFT計(jì)算，FFT輸出求模選大后得到精細(xì)同步頻率估計(jì)結(jié)果。

其中，窄帶濾波處理時(shí)，窄帶濾波的帶寬小于采樣帶寬的1/N。

S3：將所述精細(xì)估計(jì)頻率值與所述粗估計(jì)頻率值進(jìn)行相加，得到最終的DDS頻率字，并采用該DDS頻率字對輸入的MPSK接收信號進(jìn)行精確頻率補(bǔ)償。

具體地，將步驟S2中FFT輸出求模選大后的精細(xì)同步頻率估計(jì)結(jié)果加上粗同步頻率估計(jì)結(jié)果，合并產(chǎn)生出最終頻率估計(jì)的DDS頻率字，對外部輸入的MPSK接收信號進(jìn)行一次性的精確頻率補(bǔ)償。

其中，本實(shí)施例中分別采用第一級狀態(tài)機(jī)、第二級狀態(tài)機(jī)及第三級狀態(tài)機(jī)來執(zhí)行上述的步驟S1、S2及S3的流程，使用3級狀態(tài)機(jī)實(shí)現(xiàn)頻率估計(jì)和補(bǔ)償，可以使RAM、FFT、DDS、復(fù)數(shù)乘法器等資源在不同狀態(tài)機(jī)狀態(tài)下可以分時(shí)復(fù)用，節(jié)省了程序的資源占用。

具體地，采用第一級狀態(tài)機(jī)將獲取的所述MPSK接收信號中的復(fù)數(shù)數(shù)據(jù)存儲在RAM中，存儲的長度為FFT的長度，并由所述第一級狀態(tài)機(jī)進(jìn)行調(diào)制處理及求模后選取其中的最大值得到所述粗估計(jì)頻率值。

第一級狀態(tài)機(jī)將所述粗估計(jì)頻率值發(fā)送至所述第二級狀態(tài)機(jī)，所述第二級狀態(tài)機(jī)讀取RAM中的MPSK接收信號的復(fù)數(shù)數(shù)據(jù)并進(jìn)行頻率粗補(bǔ)償、調(diào)制、窄帶濾波處理、抽取及求模后選取其中的最大值，得到精細(xì)估計(jì)頻率值。

所述第二級狀態(tài)機(jī)將所述精細(xì)估計(jì)頻率值發(fā)送至所述第三級狀態(tài)機(jī)，所述第三級狀態(tài)機(jī)采用計(jì)算的DDS頻率字對外部輸入的MPSK接收信號進(jìn)行精確頻率補(bǔ)償。

具體地，本實(shí)施例以MPSK中的QPSK解調(diào)為實(shí)施例(M為調(diào)制相位點(diǎn)數(shù)，QPSK中M＝4)，設(shè)采樣速率為100MHz，調(diào)制信號符號速率為25Msps，采樣頻率為4倍符號速率。如圖3所示，為本發(fā)明方法實(shí)現(xiàn)過程示意圖，用3級狀態(tài)機(jī)完成程序資源的復(fù)用和信號處理，3級狀態(tài)機(jī)為流水線處理結(jié)構(gòu)。

步驟(1)：第1級狀態(tài)機(jī)(State 1)

首先，對輸入的復(fù)數(shù)(I、Q路)數(shù)據(jù)用RAM_1進(jìn)行存儲，存儲長度為FFT的長度，如這里選用1024點(diǎn)長度。然后用M次方法對RAM讀出的1024點(diǎn)數(shù)據(jù)進(jìn)行M次復(fù)數(shù)乘法，以完成去調(diào)制，恢復(fù)調(diào)制載波，這里針對QPSK解調(diào)，M取4。進(jìn)行4次方復(fù)乘后，頻率范圍縮小為采樣頻率的1/4，即頻率檢測范圍為±12.5MHz，經(jīng)過1024點(diǎn)FFT處理，并求模選大求出載波頻率位置，頻率分辨率為25MHz/1024＝24.41KHz。求模選大的結(jié)果A1送給狀態(tài)機(jī)2進(jìn)行處理。

步驟(2)：第2級狀態(tài)機(jī)(State 2)

接收第1級狀態(tài)機(jī)求模選大結(jié)果A1進(jìn)行DDS粗補(bǔ)償頻率字產(chǎn)生，然后再次利用RAM_1讀出原存儲的1024點(diǎn)數(shù)據(jù)結(jié)合DDS進(jìn)行I、Q路頻率粗同步結(jié)果補(bǔ)償，頻率誤差在補(bǔ)償后控制在±24.41KHz以內(nèi)。此時(shí)對粗頻率補(bǔ)償后的數(shù)據(jù)序列繼續(xù)進(jìn)行M次方復(fù)乘，復(fù)乘后的序列需要進(jìn)行窄帶濾波，再進(jìn)行64倍抽取，以縮小頻率范圍，進(jìn)行精細(xì)頻率估計(jì)和補(bǔ)償。此時(shí)的窄帶濾波因?yàn)槭窃诖滞胶筮M(jìn)行，頻率誤差已經(jīng)控制在0.1％以內(nèi)，只需使用一個(gè)要求不高的帶寬為1％(帶寬需小于后端抽取倍數(shù))采樣率的窄帶低通濾波器即可，階數(shù)可選32階，以濾除M次方復(fù)乘后的噪聲和雜散，并降低64倍抽取后信號的混疊影響。經(jīng)過64倍抽取，FFT的鑒頻范圍縮小為±12.5M/64＝±195.3KHz，大于和覆蓋粗同步補(bǔ)償后的±24.41KHz頻偏范圍。1024點(diǎn)窄帶濾波數(shù)據(jù)經(jīng)64倍抽取后，抽出點(diǎn)數(shù)為16個(gè)點(diǎn)，后面補(bǔ)齊1024-16＝1008個(gè)0，以湊齊1024點(diǎn)進(jìn)行精確同步的1024點(diǎn)FFT計(jì)算，此時(shí)FFT的頻率分辨率為195.3KHz×2/1024＝381Hz。求模選大后的結(jié)果B1送給狀態(tài)機(jī)3進(jìn)行處理。

步驟(3)：第3級狀態(tài)機(jī)(State 3)

接收狀態(tài)機(jī)2的FFT求模選大精確頻率結(jié)果B1，并加上粗同步結(jié)果，產(chǎn)生出最終精確頻率估計(jì)的DDS頻率字，對外部輸入的信號進(jìn)行一次性的精確頻率補(bǔ)償，得到同步結(jié)果C1。其中“最大指數(shù)乘法器”完成FFT求模選大后的最大值位置乘以DDS相位步進(jìn)的乘法功能，產(chǎn)生出送給DDS補(bǔ)償頻率的相位控制字，例如DDS相位累積器位寬為32位，前端采用1024點(diǎn)的FFT，則相位步進(jìn)為232/1024，FFT最大值在第860個(gè)點(diǎn)的位置，則經(jīng)最大指數(shù)乘法器后的輸出結(jié)果為860×232/1024，此結(jié)果作為后端DDS的輸入相位增量控制字。“上一狀態(tài)最大相位增量值”即為粗頻率估計(jì)的最大指數(shù)乘法器輸出值，方法同上。“IQ補(bǔ)償乘法器”即是對原輸入I、Q路數(shù)據(jù)(RAM1輸出)進(jìn)行正交頻率補(bǔ)償，以抵消頻率差。

其中，此三級狀態(tài)機(jī)所用的資源，如RAM、1024點(diǎn)FFT、DDS、復(fù)數(shù)乘法器等在不同狀態(tài)機(jī)狀態(tài)下可以分時(shí)復(fù)用，以降低程序資源的占用率。估計(jì)出的頻率精度一次可以達(dá)到(采樣頻率÷(M×1024×64))的精確效果，且頻率估計(jì)范圍仍保持(采樣頻率÷M)的大帶寬。

該方法采用兩級頻率估計(jì)方法，完成頻率的粗估計(jì)和精細(xì)估計(jì)，其中兩級頻率估計(jì)的資源占用模塊可以分時(shí)復(fù)用，第一級使用常規(guī)的FFT頻率估計(jì)和補(bǔ)償方法，第二級在第一級粗補(bǔ)償?shù)幕A(chǔ)上進(jìn)行M次方復(fù)乘后的窄帶濾波和N倍抽取，以縮小頻率估計(jì)的帶寬，并降低雜散和噪聲，抽取后的結(jié)果經(jīng)過補(bǔ)零處理后再次進(jìn)行FFT精細(xì)頻率估計(jì)和補(bǔ)償，最終完成MPSK解調(diào)的快速精確頻率同步。本發(fā)明的頻率同步方法資源占用率小，頻率估計(jì)精度高，可以適用于連續(xù)和突發(fā)MPSK調(diào)制信號的頻率快速精確估計(jì)和補(bǔ)償。

以上所述，僅為本發(fā)明的具體實(shí)施方式，但本發(fā)明的保護(hù)范圍并不局限于此，任何本領(lǐng)域的技術(shù)人員在本發(fā)明揭露的技術(shù)范圍內(nèi)，對本發(fā)明所做的變形或替換，都應(yīng)涵蓋在本發(fā)明的保護(hù)范圍之內(nèi)。因此，本發(fā)明的保護(hù)范圍應(yīng)以所述的權(quán)利要求的保護(hù)范圍為準(zhǔn)。

總結(jié)

以上是生活随笔為你收集整理的fft 估计载波频率程序_一种适用于MPSK解调的快速精确频率同步方法与流程的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。