一種利用范數(shù)性質(zhì)的無約束盲均衡準則

王大磊， 楊 賓， 吳 瑛， 王秀秀

(1. 信息工程大學 信息系統(tǒng)工程學院，河南 鄭州 450002；2. 西安電子科技大學 綜合業(yè)務網(wǎng)理論及關鍵技術國家重點實驗室，陜西 西安 710071)

在數(shù)字通信中，作為克服碼間干擾(Inter-Symbol Interference，ISI)的重要手段，均衡技術受到廣泛研究[1]．由于盲均衡方法無需訓練序列，可以節(jié)省帶寬資源，且適合多點通信和非協(xié)作接收等場合，因此，涌現(xiàn)出了大量的盲均衡方法．典型的有Bussgang類型，如恒模(CMA)算法[2]、多模算法[3-4]和混合算法[5]等．此類算法基于某種非最小均方誤差(MSE)的代價函數(shù)，采用隨機梯度法求極值，計算量小，實現(xiàn)簡單．然而，由于這類代價函數(shù)大都為非凸的，往往會局部收斂，導致不能充分消除ISI．

Shtrom和Fan依據(jù)向量范數(shù)的性質(zhì)，提出了一類迫零盲均衡算法Shtrom-Fan算法(SFA)[6]，此類代價函數(shù)在信道與均衡器聯(lián)合響應域不存在假峰，能夠保證全局收斂．然而，SFA算法需要固定均衡器的某一個系數(shù)或?qū)ζ錃w一化約束，致使均衡器輸出與發(fā)送信號的增益比不為1，需要自動增益控制，且這種約束方法對均衡器收斂性能的影響并沒有嚴格的數(shù)學證明．之后，Abrar等[7]同樣利用向量范數(shù)性質(zhì)，結(jié)合星座圖特點，提出了一種自適應的盲均衡算法βCMA．然而，該算法的參數(shù)是依據(jù)Bussgang的特性得到的，只反映了信號的部分統(tǒng)計特征，致使算法收斂后的穩(wěn)態(tài)誤差較大．基于以上考慮，筆者依據(jù)向量范數(shù)性質(zhì)，結(jié)合文獻[8]中將約束方程轉(zhuǎn)換為無約束方程的方法，給出了一種新的無約束盲均衡準則，證明了在無噪聲情況下，在代價函數(shù)的局部最優(yōu)解處可實現(xiàn)理想均衡，且能夠?qū)崿F(xiàn)發(fā)送信號幅值的恢復．可采用批數(shù)據(jù)處理和在線自適應迭代兩種方法更新均衡器系數(shù)，無需自動增益控制，實現(xiàn)簡單，同時穩(wěn)態(tài)誤差更?。?/p>

1 系統(tǒng)模型

2 新的無約束盲均衡準則

2.1 SFA算法及βCMA算法

其中，γ為發(fā)送符號的最大幅值．然后，通過拉格朗日乘子法將約束條件引入到目標函數(shù)，得出一種自適應的迭代算法βCMA，即

wn+1=wn+μφ(yn)*xn，

(1)

其中，φ(yn)為誤差函數(shù)，可依據(jù)Bussgang特性求得．當|yn|<γ時，φ(yn)=yn；當 |yn|=γ時，φ(yn)=0；當 |yn|>γ時，φ(yn)=-βyn，β由發(fā)送信號的平均能量以及星座圖決定．盡管該算法能夠恢復發(fā)送信號的幅值，但由于其參數(shù)是依據(jù)Bussgang特性得到的，只反映了信號的部分統(tǒng)計特征，致使算法收斂后的穩(wěn)態(tài)誤差較大．

2.2 新的無約束盲均衡準則

其中，p，q和ζ為實數(shù)，且0

(2)

其中，f: [0，∞]→R1，為一個分段連續(xù)實函數(shù)；g(x)=x2+f(x)，是在0≤x<1區(qū)間內(nèi)單調(diào)遞增，在x>1 區(qū)間單調(diào)遞減的函數(shù)，并在x=1 處取得最大值．代價函數(shù)建立后，需要考察它的收斂性和極值點的分布問題．下面給出關于該代價函數(shù)收斂性的兩個定理．

定理1 當且僅當{si}序列中只有1個非零元素，且幅度為1，即s=δ(i-k) exp(jθ)時，J(s)取得最大值．

證明 由向量范數(shù)性質(zhì)可得

定理2J(s)不存在假峰，即任何局部極大值解均對應理想均衡條件．

證明J(s)的平穩(wěn)點滿足

(3)

式(3)的解為具有M個相同幅度元素的向量，且幅度aM滿足

2aM+f′(MaM)=0 ．

(4)

式(5)中，當α≠0時，不等式嚴格成立．因此，sM在M≥2時，不是局部極大值點．

因此，sM在M≥2時，也不是局部極小值點．從而這些點都是鞍點．綜上，J(s)平穩(wěn)點中只有解向量集s1是局部極值點，同時也是全局最優(yōu)解．證畢．

定理1保證了代價函數(shù)能夠取得最大值，定理2保證了代價函數(shù)的收斂性，采用隨機梯度法時能夠?qū)ふ业阶畲笾担?/p>

為得到算法的具體形式，考慮一個二次函數(shù)，g(x)=-a(x-1)2+a，a>0，該函數(shù)滿足 0≤x<1 區(qū)間內(nèi)單調(diào)遞增，x>1 區(qū)間單調(diào)遞減的性質(zhì)．此時，代價函數(shù)為

(7)

由于該代價函數(shù)是定義在s域的，而s不能直接得到，要想將代價函數(shù)應用于盲均衡，必須將代價函數(shù)轉(zhuǎn)換為均衡器輸出{yn}的函數(shù)．在無噪聲情況下，有

(8)

(9)

其中，γ=max{|an|}．將式(8)和式(9)代入式(7)中，得到

(10)

該代價函數(shù)同大多數(shù)盲均衡算法[6，8，10]一樣，是在理想的無噪聲情況下得到的．然而實際中信道噪聲是存在的，對該代價函數(shù)的優(yōu)化不能保證實現(xiàn)完全均衡．下面對有噪聲情況下算法的性能進行討論．由于s=wHH，則wH=sH-1，此時，式(8)修正為

(11)

(12)

3 算法實現(xiàn)

算法的兩種具體實現(xiàn)方式: 基于一段數(shù)據(jù)的批處理方法和在線自適應迭代方法．

3.1 批數(shù)據(jù)處理方法

采用隨機梯度下降法更新均衡器權值:

(14)

(15)

將式(16)和式(17)代入方程式(15)，得到

(18)

歸結(jié)起來，批數(shù)據(jù)處理算法步驟如下:

步驟1 收集一段長度為N的信道輸出數(shù)據(jù)，計算均衡器輸入向量的協(xié)方差矩陣的期望E[Xn]，并初始化均衡器．

步驟2 計算均衡器輸出序列{yn}．

步驟3 根據(jù)式(14)和式(18)計算新的均衡器權值，重復迭代直至收斂．

將式(14)和式(18)組成的批數(shù)據(jù)處理方法命名為批數(shù)據(jù)處理的無約束盲均衡算法(Batch Processing Unconstrained Equalization Algorithm，BP-UEA)．

3.2 在線自適應處理方法

(19)

(20)

其中，誤差函數(shù)Φ(yn)為

(21)

(22)

(23)

將式(14)、式(20)和式(23)組成的自適應迭代方法稱為簡化的自適應無約束盲均衡算法(Simplified Adaptive Processing Unconstrained Equalization Algorithm，SAP-UEA)．

4 仿真與分析

為了驗證文中算法的性能，可選取同樣基于向量范數(shù)性質(zhì)的βCMA和SFA算法進行性能對比．在仿真中，如未特殊說明，條件統(tǒng)一設置如下: 均衡器輸入信號的疊加噪聲為高斯白噪聲，信噪比為 30 dB；均衡器階數(shù)為15階，初始化時中心抽頭為1，其他系數(shù)為0；AP-UEA和SAP-UEA算法中的參數(shù)p和a分別設為10和1，實驗次數(shù)為100次．

4.1 算法有效性的驗證

對文中算法的3種實現(xiàn)方法BP-UEA、AP-UEA和SAP-UEA的有效性進行了仿真驗證．信道采用文獻[12]中圖2所示的語音通信信道，信號采用16QAM調(diào)制的獨立同分布信號源．圖1(a)為BP-UEA算法的剩余ISI隨迭代次數(shù)的變化曲線，在步長相同時，N越大，收斂速度越快，穩(wěn)態(tài)誤差越?。畧D1(b)為AP-UEA算法和SAP-UEA算法剩余ISI隨迭代次數(shù)的變化曲線．可以看出，兩種自適應算法的性能幾乎無差別．相比自適應方法，批數(shù)據(jù)處理算法BP-UEA需要的樣點個數(shù)少，適合短時、數(shù)據(jù)量少的場合．

圖1 BP-UEA，AP-UEA和SAP-UEA算法的剩余ISI曲線

為驗證不同的初始化條件對均衡器性能的影響，對3種處理算法BP-UEA、AP-UEA及SAP-UEA進行如下仿真: 仿真條件同上，均衡器系數(shù)初始化為1的位置依次改變，得到剩余ISI隨不同初始化位置改變的曲線如圖1(c)所示．可以看出，盡管均衡器系數(shù)初始化為1的位置不同，使得穩(wěn)態(tài)ISI不一致，但均衡器均能收斂．

為驗證算法在中等信噪比下的性能，以SAP-UEA算法為例進行如下仿真: 仿真條件同上，信噪比分別設為 15 dB，20 dB，25 dB，經(jīng)過 6 000 次迭代，糾正相偏后的均衡輸出結(jié)果如圖2所示．

圖2 SAP-UEA均衡輸出星座圖

從圖2可以看出，隨著信噪比的降低，算法性能逐漸下降．然而，在中等至較高信噪比條件下，能夠得到較理想的均衡輸出結(jié)果．

4.2 與βCMA算法的性能比較

對自適應算法SAP-UEA與βCMA算法進行比較與分析．信號為8APSK調(diào)制的獨立同分布信號源，信道仍采用文獻[12]中的復信道．圖3為剩余ISI隨迭代次數(shù)的變化曲線．可以看出，SAP-UEA算法相比βCMA算法具有更低的穩(wěn)態(tài)誤差，其原因可以從兩者的誤差函數(shù)Φ(yn)及φ(yn)看出．圖4為兩者的誤差函數(shù)曲線，其中，γ=2．當 |yn|?γ時，|Φ(yn)|? |φ(yn)|．因此，SAP-UEA算法具有更快的收斂速度；當 |yn|=γ，|Φ(yn)|=0 時，φ(yn)的符號會隨機變化，導致βCMA算法的穩(wěn)態(tài)誤差較大．

4.3 與SFA算法的性能比較

對SAP-UEA算法與SFA算法的性能進行對比．發(fā)射信號為獨立同分布的16QAM調(diào)制信號，信道采用在文獻[13-14]所用的復信道，其脈沖響應如圖5所示．圖6為剩余ISI隨迭代次數(shù)變化的曲線．可以看出，SAP-UEA算法相比SFA算法的穩(wěn)態(tài)誤差更?。@是由于SFA算法采用了均衡器中心抽頭系數(shù)固定為1的約束所造成的，這種約束在避免代價函數(shù)收斂至0的同時，也限制了均衡器系數(shù)的調(diào)整．而SAP-UEA算法的代價函數(shù)對均衡器沒有約束，在均衡器系數(shù)的更新過程中，都可以使其向最優(yōu)方向調(diào)整．

圖3 SAP-UEA與βCMA算法的剩余ISI曲線圖4 誤差函數(shù)曲線

圖5 信道的脈沖響應圖6 SAP-UEA與SFA算法的剩余ISI曲線

5 結(jié) 束 語

基于向量范數(shù)的性質(zhì)，給出了一種新的無約束盲均衡代價函數(shù)．并對代價函數(shù)的收斂性，極值點的分布進行了分析和討論；證明了在無噪聲情況下，在代價函數(shù)的每個局部最優(yōu)解處可實現(xiàn)理想均衡，保證了采用隨機梯度法迭代時能夠使代價函數(shù)收斂；并給出了適應不同場合的批數(shù)據(jù)處理和在線自適應處理兩種方法．通過分析和仿真實驗說明，該方法在收斂速度和穩(wěn)態(tài)誤差等方面優(yōu)于文獻[6-7]中的方法．

[1] 趙良, 汪春梅, 孫洪波. 逐幸存序列算法引導的快速盲均衡算法[J]. 西安電子科技大學學報, 2013, 40(2): 201-206.

Zhao Liang, Wang Chunmei, Sun Hongbo. Fast Blind Equalization Algorithm Based on PSP[J]. Journal of Xidian University, 2013, 40(2): 201-206.

[2] Godard D N. Self-Recovering Equalization and Carrier Tracking in Two-Dimensional Data Communication Systems[J]. IEEE Transactions on Communications, 1980, 28(11): 1867-1875.

[3] 李進, 馮大正, 劉文娟. 快速Q(mào)AM信號多模盲均衡算法[J]. 電子與信息學報, 2013, 35(2): 273-279.

Li Jin, Feng Dazheng, Liu Wenjuan. A Fast Multimodulus Blind Equalization Algorithm for QAM Signal[J]. Journal of Electronics & Information Technology, 2013, 35(2): 273-279.

[4] Yuan J T, Chao J H, Lin T C. Effect of Channel Noise on Blind Equalization and Carrier Phase Recovery of CMA and MMA[J]. IEEE Transactions on Communications, 2012, 60(11): 3274-3285.

[5] Silva M T M, Arenas-Garcia J. A Soft-Switching Blind Equalization Scheme via Convex Combination of Adaptive Filters[J]. IEEE Transactions on Signal Processing, 2013, 61(5): 1171-1182.

[6] Shtrom V, Fan H. New Class of Zero-Forcing Cost Functions in Blind Equalization[J]. IEEE Transactions on Signal Processing, 1998, 46(10): 2674-2683.

[7] Abrar S, Nandi A K. Adaptive Minimum Entropy Equalization Algorithm[J]. IEEE Communications Letters, 2010, 14(10): 966-968.

[8] Shalvi O, Weinstein E. New Criteria for Blind Deconvolution of Nonminimum Phase Systems (Channels)[J]. IEEE Transactions on Information Theory, 1990, 36(2): 312-321.

[9] Verbu S, Verdu S, Kenndy R, et al. Convex Cost Functions in Blind Equalization[J]. IEEE Transactions on Signal Processing, 1994, 42(8): 1952-1960.

[10] Chen M, Tuqan J, Ding Z. A Quadratic Programming Approach to Blind Equalization and Signal Separation[J]. IEEE Transactions on Signal Processing, 2009, 57(6): 2232-2244.

[11] Han H D, Ding Z. Steepest Descent Algorithm Implementation for Multichannel Blind Signal Recovery[J]. IET Communications, 2012, 6(18): 3196-3203.

[12] Picchi G, Prati G. Blind Equalization and Carrier Recovery Using a ‘Stop-and-Go’ Decision-Directed Algorithm[J]. IEEE Transactions on Communications, 1987, 35(9): 877-887.

[13] Kennedy R A, Ding Z. Blind Adaptive Equalizers for Quadrature Amplitude Modulated Communication Systems Based on Convex Cost Functions[J]. Optical Engineering, 1992, 31(6): 1189-1199.

[14] Ding Z, Luo Zhiquan. A Fast Linear Programming Algorithm for Blind Equalization[J]. IEEE Transactions on Communications, 2000, 48(9): 1432-1436.