中如果初始頻偏估計為正值，在計算候 選修正值時選擇加法，否則選擇減法：
[0272]
[0273]
[0275] 根據(jù)三組候選相位值分別計算如下三組相位差值：
[0274]
計算如下絕對差值Di: ，確認(rèn) f 對應(yīng)絕對差值最小的修訂相位值
作為最終修訂的相位值，即令：
[0276] 其中，所述判決模塊，當(dāng)判斷滿足下式時，則判決修正后的相位為正常值：
其中，S為判斷相位是否異常的可調(diào)判決閾值。
[0277] 其中，所述第二確定模塊，可以分別根據(jù)所述兩組相位滑窗緩存中 相位值計算平均相位值ave_anglel(n)、ave_angle2(n)是通過下式實(shí)現(xiàn)的：
;根據(jù)所述平均相位 值確定最大似然頻偏估值e(n)是通過下式實(shí)現(xiàn)的：
[0278]
[0279] 其中，所述處理模塊，對緩存的最大似然頻偏估值進(jìn)行平滑濾波處理得到當(dāng)前時
刻的最大似然頻偏估值是可以通過下式實(shí)現(xiàn)的 其中，P為所述頻偏 t.ML 估計結(jié)果滑窗緩存的指定長度。
[0280] 本領(lǐng)域普通技術(shù)人員可以理解上述方法中的全部或部分步驟可通過程序來指令 相關(guān)硬件完成，所述程序可以存儲于計算機(jī)可讀存儲介質(zhì)中，如只讀存儲器、磁盤或光盤 等?？蛇x地，上述實(shí)施例的全部或部分步驟也可以使用一個或多個集成電路來實(shí)現(xiàn)。相應(yīng) 地，上述實(shí)施例中的各模塊/單元可以采用硬件的形式實(shí)現(xiàn)，也可以采用軟件功能模塊的 形式實(shí)現(xiàn)。本發(fā)明不限制于任何特定形式的硬件和軟件的結(jié)合。
[0281] 以上僅為本發(fā)明的優(yōu)選實(shí)施例，當(dāng)然，本發(fā)明還可有其他多種實(shí)施例，在不背離本 發(fā)明精神及其實(shí)質(zhì)的情況下，熟悉本領(lǐng)域的技術(shù)人員當(dāng)可根據(jù)本發(fā)明作出各種相應(yīng)的改變 和變形，但這些相應(yīng)的改變和變形都應(yīng)屬于本發(fā)明所附的權(quán)利要求的保護(hù)范圍。
【主權(quán)項】
1. 一種基于差分相位的頻偏估計方法，包括： 接收到導(dǎo)頻符號時，更新接收導(dǎo)頻數(shù)據(jù)緩存； 從所述接收導(dǎo)頻數(shù)據(jù)緩存中提取當(dāng)前接收的導(dǎo)頻符號，確定初始頻偏估計，根據(jù)所述 初始頻偏估計確定相位差值； 對所述接收導(dǎo)頻數(shù)據(jù)緩存中所存儲的相鄰兩個導(dǎo)頻符號進(jìn)行預(yù)處理，分別計算所述兩 個導(dǎo)頻符號對應(yīng)的相位； 根據(jù)所述初始頻偏估計、所述相位差值對所述兩個導(dǎo)頻符號對應(yīng)的相位進(jìn)行修正； 如判決修正后的相位為正常值，則分別對兩組相位滑窗緩存進(jìn)行更新； 分別根據(jù)所述兩組相位滑窗緩存中相位值計算平均相位值，根據(jù)所述平均相位值確定 最大似然頻偏估值； 根據(jù)所述最大似然頻偏估值更新頻偏估計結(jié)果滑窗緩存，對緩存的最大似然頻偏估值 進(jìn)行平滑濾波處理得到當(dāng)前時刻的最大似然頻偏估值。2. 如權(quán)利要求1所述的方法，其特征在于：所述更新接收導(dǎo)頻數(shù)據(jù)緩存是通過以下方 式實(shí)現(xiàn)的：獄f用于緩存當(dāng)前時刻接收到的導(dǎo)頻符號R(n)，Rl用于緩存上一時刻R；；1中緩存的導(dǎo) 頻符號；表示第q根接收天線 接收到的第η個正交頻分復(fù)用（OFDM)時域符號塊，Nt為多輸入多輸出（ΜΙΜΟ)系統(tǒng)的發(fā)射 天線數(shù)，Nk為MMO系統(tǒng)的接收天線數(shù)，代表第η個導(dǎo)頻符號對應(yīng)的頻偏矩陣， Eixi =e/2~/wi -代表載波頻偏矩陣，N為OFDM子載波數(shù)，為歸一化頻偏值，Ln代表第η個導(dǎo)頻符號第一 個抽樣時刻所對應(yīng)的序號，為對應(yīng)的發(fā)送導(dǎo)頻符號，I代表第P根發(fā)送天線發(fā)送的第η個頻域?qū)ьl 符號矩陣； IeiJjvw為離散傅里葉變換矩陣，其中第（l，m)個元素為#代表DFT矩陣的前L列，即且有3. 如權(quán)利要求2所述的方法，其特征在于：所述確定初始頻偏估計是通過以下方式 實(shí)現(xiàn)的：所述根據(jù)所述初始頻偏估計確定相位差值是通過下面方式實(shí)現(xiàn)的：其中，Lp代表相鄰導(dǎo)頻符號間隔長度。4. 如權(quán)利要求3所述的方法，其特征在于：對所述接收導(dǎo)頻數(shù)據(jù)緩存中所存儲的相鄰 兩個導(dǎo)頻符號Rf、進(jìn)行預(yù)處理，分別計算所述兩個導(dǎo)頻符號對應(yīng)的相位anglel (n)、 angle2(n)，包括： 對所述接收導(dǎo)頻數(shù)據(jù)緩存中所存儲的相鄰兩個導(dǎo)頻符號Rf、R，進(jìn)行預(yù)處理，得到 兩個修訂測量信號若發(fā)送端發(fā)送的兩個相鄰頻域?qū)ьl符號相同，則若發(fā)送端發(fā)送的頻域?qū)ьl符號為長期演進(jìn)標(biāo)準(zhǔn)所規(guī)定的導(dǎo)頻圖樣，則其中，B=vec (A)，BiWec (Ah)，vec (Anxn)代表取對角陣Anxn元素組成的行向量， A=diag[f，...，eXN-1)]，α = π。5. 如權(quán)利要求4所述的方法，其特征在于：所述根據(jù)所述初始頻偏估計、所述相位差值 對所述兩個導(dǎo)頻符號對應(yīng)的相位進(jìn)行修正，包括： 你捉下才;+管和仿折碁的願撒timpQ .其中，round(·)代表四舍五入運(yùn)算； 按以下方法計算anglel(n), angle2(n)的三組相位候選修正值，其中如果初始頻偏估計 為正值，在計算候選修正值時選擇加法，否則選擇減法：根據(jù)三組候選相位值分別計算如下三組相位差值：確認(rèn)對應(yīng)絕對差值最小的修訂相位值作為最終修訂的相位值，即 令：6. 如權(quán)利要求5所述的方法，其特征在于：當(dāng)判斷滿足下式時，則判決修正后的相位為 正常值：,其中，S為判斷相位是否異常的可調(diào)判決閾值。7. 如權(quán)利要求6所述的方法，其特征在于：所述分別根據(jù)所述兩組相位滑窗緩存中相 位值計算平均相位值ave_anglel(n)、ave_angle2 (n)是通過下式實(shí)現(xiàn)的：其中M為所述 相位滑窗緩存的指定長度； 所述根據(jù)所述平均相位值確定最大似然頻偏估值ε (η)是通過下式實(shí)現(xiàn)的： 8. 如權(quán)利要求7所述的方法，其特征在于：所述對緩存的最大似然頻偏估值進(jìn)行平滑 濾波處理得到當(dāng)前時刻的最大似然頻偏估值是通過下式實(shí)現(xiàn)的：其中，P為所述頻偏估計結(jié)果滑窗緩存的指定長度。9. 一種基于差分相位的頻偏估計的裝置，其特征在于，包括： 更新模塊，用于接收到導(dǎo)頻符號時，更新接收導(dǎo)頻數(shù)據(jù)緩存； 第一確定模塊，用于從所述接收導(dǎo)頻數(shù)據(jù)緩存中提取當(dāng)前接收的導(dǎo)頻符號，確定初始 頻偏估計，根據(jù)所述初始頻偏估計確定相位差值； 預(yù)處理模塊，用于對所述接收導(dǎo)頻數(shù)據(jù)緩存中所存儲的相鄰兩個導(dǎo)頻符號進(jìn)行預(yù)處 理，分別計算所述兩個導(dǎo)頻符號對應(yīng)的相位； 修正模塊，用于根據(jù)所述初始頻偏估計、所述相位差值對所述兩個導(dǎo)頻符號對應(yīng)的相 位進(jìn)行修正； 判決模塊，用于如判決修正后的相位為正常值，則分別對兩組相位滑窗緩存進(jìn)行更 新； 第二確定模塊，用于分別根據(jù)所述兩組相位滑窗緩存中相位值計算平均相位值，根據(jù) 所述平均相位值確定最大似然頻偏估值； 處理模塊，用于根據(jù)所述最大似然頻偏估值更新頻偏估計結(jié)果滑窗緩存，對緩存的最 大似然頻偏估值進(jìn)行平滑濾波處理得到當(dāng)前時刻的最大似然頻偏估值。10. 如權(quán)利要求9所述的裝置，其特征在于： 所述更新模塊，更新接收導(dǎo)頻數(shù)據(jù)緩存是通過以下方式實(shí)現(xiàn)的：1用于緩存所述當(dāng)前時刻接收到的導(dǎo)頻符號R (n)， 用于緩存上一時刻中緩存的導(dǎo)頻符號，表示第q根接收天線接收到的第η個正交頻 分復(fù)用（OFDM)時域符號塊，Nt為多輸入多輸出（ΜΙΜΟ)系統(tǒng)的發(fā)射天線數(shù)，Nk 為MMO系統(tǒng)的接收天線數(shù)，代表第η個導(dǎo)頻 符號對應(yīng)的頻偏矩陣，代表載波頻偏矩陣，N為OFDM子載波數(shù)，ε為歸一化頻偏值，Ln代表第η個導(dǎo)頻符號第 一個抽樣時刻所對應(yīng)的序號，對應(yīng)的發(fā)送導(dǎo)頻符號，代表第P根發(fā)送天線 發(fā)送的第η個頻域?qū)ьl符號矩陣；為離散傅里葉變換矩陣，其中第（l，m)個元 素為代表DFT矩陣f的前L列，即11. 如權(quán)利要求10所述的裝置，其特征在于： 所述第一確定模塊，確定初始頻偏估計是通過以下方式實(shí)現(xiàn)的：初始頻偏估計確定相位差值是通過下面方式實(shí)現(xiàn)的：_ ,其中，Lp代表相鄰導(dǎo)頻 符號間隔長度。12. 如權(quán)利要求11所述的裝置，其特征在于： 所述預(yù)處理模塊，對所述接收導(dǎo)頻數(shù)據(jù)緩存中所存儲的相鄰兩個導(dǎo)頻符號Rif、Rl1 進(jìn)行預(yù)處理，分別計算所述兩個導(dǎo)頻符號對應(yīng)的相位anglel(n)、angle2(n)，包括：對 所述接收導(dǎo)頻數(shù)據(jù)緩存中所存儲的相鄰兩個導(dǎo)頻符號Rf、RiJ進(jìn)行預(yù)處理，得到兩 個修訂測量信號若發(fā)送端發(fā)送的兩個相鄰頻域?qū)ьl符號相同，則:若發(fā)送端發(fā)送的頻域?qū)ьl符號為長期 演進(jìn)標(biāo)準(zhǔn)所規(guī)定的導(dǎo)頻圖樣，則，其 中，B=vec⑷，BiWec (Ah)，vec (Anxn)代表取對角陣Anxn元素組成的行向量，A=diag[ej°°，… ，#(，，α = π。13. 如權(quán)利要求12所述的裝置，其特征在于： 所沭修iF樽塊，具體用于依據(jù)下式計算相位折疊的圈數(shù)times :,其中，round (.)代表四舍五入運(yùn)算；按以下方法計算 anglel(n), angle2(n)的三組相位候選修正值，其中如果初始頻偏估計為正值，在計算候 選修正值時選擇加法，否則選擇減法：根據(jù)三組候選相位值分別計算如下三組相位差值：計算如下絕對差值Di :?確認(rèn) 對應(yīng)絕對差值最小的修訂相位值乍為最終修訂的相位值，即令：t14. 如權(quán)利要求13所述的裝置，其特征在于： 所述判決模塊，當(dāng)判斷滿足下式時，則判決修正后的相位為正常值：，其中，S為判斷相位是否異常的可調(diào)判決閾值。15. 如權(quán)利要求14所述的裝置，其特征在于： 所述第二確定模塊，分別根據(jù)所述兩組相位滑窗緩存中相位值 計算平均相位值ave_anglel(n)、ave_angle2(n)是通過下式實(shí)現(xiàn)的：，其中M為所述相位 滑窗緩存的指定長度；根據(jù)所述平均相位值確定最大似然頻偏估值ε (η)是通過下式實(shí)現(xiàn)16.如權(quán)利要求15所述的裝置，其特征在于：所述處理模塊，對緩存的最大似然頻偏估值進(jìn)行平滑濾波處理得到當(dāng)前時刻的最大似 然頻偏估值是通過下式實(shí)現(xiàn)的： ，其中，P為所述頻偏估計結(jié)果滑窗緩 cML 存的指定長度。
【專利摘要】本發(fā)明公開了一種基于差分相位的頻偏估計方法及裝置，該方法包括：接收到導(dǎo)頻符號時更新接收導(dǎo)頻數(shù)據(jù)緩存；從接收導(dǎo)頻數(shù)據(jù)緩存中提取當(dāng)前接收的導(dǎo)頻符號確定初始頻偏估計，進(jìn)而確定相位差值；對接收導(dǎo)頻數(shù)據(jù)緩存中所存儲的相鄰兩個導(dǎo)頻符號進(jìn)行預(yù)處理，分別計算出對應(yīng)的相位；對所述兩個導(dǎo)頻符號對應(yīng)的相位進(jìn)行修正；如判決修正后的相位為正常值，則分別對兩組相位滑窗緩存進(jìn)行更新；分別根據(jù)兩組相位滑窗緩存中相位值計算平均相位值，進(jìn)而確定最大似然頻偏估值；根據(jù)最大似然頻偏估值更新頻偏估計結(jié)果滑窗緩存，對緩存的最大似然頻偏估值進(jìn)行平滑濾波處理得到當(dāng)前時刻的最大似然頻偏估值。本發(fā)明具有頻偏估計精確，估值波動范圍小的優(yōu)點(diǎn)。
【IPC分類】H04L27/26, H04L25/02
【公開號】CN104980375
【申請?zhí)枴緾N201410134439
【發(fā)明人】陳慶春, 江海, 喬靜, 李斌, 丁遠(yuǎn)晴, 何志謙, 呂穎麗
【申請人】中興通訊股份有限公司
【公開日】2015年10月14日
【申請日】2014年4月3日
【公告號】WO2015149199A1