本發(fā)明提供一種基于EM算法的LDPC碼譯碼噪聲方差的估計方法，它是一種基于期望最大(即EM)算法的低密度奇偶校驗(即LDPC)碼譯碼噪聲方差的估計方法，屬于通信技術領域。

背景技術：

信道編碼和譯碼是通信系統(tǒng)中保證信息正確傳輸?shù)闹匾h(huán)節(jié)。對于編碼領域常用的碼字LDPC碼或渦輪碼(即Turbo)碼，其理論性能優(yōu)異，但是其優(yōu)異的性能建立在一些前提條件下，包括理想的同步以及一些準確的參數(shù)信息，任何非理想條件都會導致譯碼性能惡化。以LDPC碼為例，置信傳播譯碼算法為其主要譯碼方法。在譯碼的開始，通常需要求解接收信號相應比特后驗概率的對數(shù)似然比來初始化迭代譯碼器，而計算后驗概率的對數(shù)似然比需要已知接收信號的噪聲方差。然而，噪聲方差的估計精度直接影響了譯碼性能。因此，在信道譯碼時需要合理的估計噪聲方差。

傳統(tǒng)的方差估計需要在利用線性變換去除調(diào)制信息后進行。但當噪聲較大時，非線性變化導致信噪比進一步降低，需要通過擴大觀察區(qū)間來提升噪聲方差估計準確性，同時也增大了計算量。若已知譯碼環(huán)節(jié)接收數(shù)據(jù)，可用極大似然法估計模型的參數(shù)。但是，當模型含有隱變量時，就不能簡單地使用這些估計方法。信道編碼后的0、1序列通過二進制相移鍵控(BPSK)調(diào)制后變?yōu)?、-1的發(fā)送序列，發(fā)送序列通過高斯白噪聲信道即為接收序列，接收序列的概率分布是一種高斯混合模型。每一個接收信息所對應的發(fā)送碼字可能是1或-1，因此，發(fā)送序列是高斯混合模型的隱變量。EM算法是一種經(jīng)典的機器學習方法，用于含有隱變量的概率模型的參數(shù)估計，1977年由鄧普斯特(即Dempster)等人提出。因此，可通過EM算法對譯碼環(huán)節(jié)的噪聲方差進行估計，保證正確譯碼。

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

(一)發(fā)明的目的

本發(fā)明的目的是提出一種基于EM算法的LDPC碼譯碼噪聲方差的估計方法，它能夠在已知譯碼環(huán)節(jié)接收序列的情況下，較準確地計算信道噪聲方差，輔助迭代譯碼，減小由方差不準確導致的譯碼性能損失。

(二)技術方案

為實現(xiàn)上述目的，本發(fā)明一種基于EM算法的LDPC碼譯碼噪聲方差的估計方法，其步驟如下：

步驟一：獲得信道譯碼環(huán)節(jié)的接收序列；

步驟二：確定方差估計范圍和初值；

步驟三：基于EM算法估計噪聲方差；

步驟四：將信道噪聲方差估計結果應用于譯碼初始化和迭代環(huán)節(jié)；

其中，在步驟一中所述的“獲得信道譯碼環(huán)節(jié)的接收序列”，其作法如下：

在發(fā)送端，將信源數(shù)據(jù)通過LDPC碼信道編碼獲得編碼序列u，并進行BPSK調(diào)制獲得發(fā)送序列c，在發(fā)送序列上加入方差為σ²的高斯白噪聲n，獲得信道譯碼環(huán)節(jié)的接收序列r。

其中，在步驟二中所述的“確定方差估計范圍和初值”，其作法如下：

根據(jù)LDPC碼碼率、接收信息信噪比與方差的對應關系計算方差的合理估計范圍，在該范圍內(nèi)選取隨機數(shù)作為方差估計初值。

其中，在步驟三中所述的“基于EM算法估計噪聲方差”，其作法如下：

首先設置EM算法迭代的初值，該初值為高斯混合模型的主要參數(shù)，主要包括噪聲方差、噪聲均值和接收序列來自均值為1或-1的高斯模型概率；其次執(zhí)行EM算法的E步：根據(jù)模型參數(shù)計算分模型對接收序列的響應度；執(zhí)行EM算法的M步：根據(jù)響應度重計算模型參數(shù)。將E步和M步反復迭代直至模型參數(shù)中的方差收斂。

其中，在步驟四中所述的“將信道噪聲方差估計結果應用于譯碼初始化和迭代環(huán)節(jié)”，其作法如下：

根據(jù)接收序列和步驟三中的噪聲方差估計值計算接收序列相應比特后驗概率的對數(shù)似然比，從而初始化迭代譯碼器；執(zhí)行LDPC碼置信傳播譯碼算法的迭代過程，完成信道譯碼。

通過以上步驟，本發(fā)明實現(xiàn)了對LDPC碼信道譯碼環(huán)節(jié)所需的信道噪聲方差的估計，能有效地避免錯誤方差帶來的譯碼性能衰減。相比傳統(tǒng)方差估計方法，本發(fā)明方法估計精度高，復雜度較低，對碼字長度和信道噪聲方差變化的魯棒性好，方差估計結果較準確。

(三)優(yōu)點和功效

基于上述步驟，本發(fā)明的一種基于EM算法估計LDPC碼譯碼噪聲方差的方法可以達到以下目的：

一：本發(fā)明解決了實際工程中噪聲方差難以估計的問題，可以較準確地估計信道噪聲方差。

二：本發(fā)明將噪聲方差估計很好地應用于LDPC碼信道譯碼，防止由于噪聲方差不準確而帶來的譯碼性能衰減。

三：本發(fā)明在LDPC碼長度變化和噪聲方差變化時均能正確估計噪聲方差，魯棒性較好。

四：本發(fā)明的方差估計結果不僅適用于LDPC碼譯碼，還適用于其他同樣需要噪聲方差的譯碼過程，如Turbo碼的譯碼，應用范圍較廣泛。

總之，本發(fā)明能夠解決通信系統(tǒng)中噪聲方差的估計問題，進而輔助信道譯碼過程，是正確譯碼的必要條件。該發(fā)明魯棒性較好，在碼字長度不同、噪聲方差大小不同的情況下，均能準確地估計信道噪聲方差，可廣泛應用于LDPC碼和Turbo碼譯碼。

附圖說明

圖1是本發(fā)明所述方法的流程圖。

圖2是信道編譯碼及其噪聲方差估計的系統(tǒng)框圖。

圖3是EM算法的流程圖。

圖4是接收序列的概率密度分布函數(shù)。

圖5是不同信噪比下的方差估計值。

圖6是初始方差為正確方差的1/2時，長528的LDPC碼采用本發(fā)明的譯碼性能。

圖7是初始方差為正確方差的1/2時，長1200的LDPC碼采用本發(fā)明的譯碼性能。

圖中序號、符號、代號說明如下：

EM表示期望最大算法；LDPC碼表示低密度奇偶校驗碼；BPSK表示二進制相移鍵控；u表示LDPC碼編碼后的序列；c表示通過BPSK調(diào)制的發(fā)送序列；n表示均值為0、方差為常數(shù)的高斯白噪聲序列；r表示接收序列；Optimum Decision Threshold表示接收序列的硬判決門限；E_b/N₀表示接收序列的信噪比，單位為分貝(即dB)。

具體實施方式

下面結合附圖和具體實施例對本發(fā)明的一種基于EM算法的LDPC碼譯碼噪聲方差的估計方法進行更進一步的介紹。

如圖1所示，本發(fā)明提供一種基于EM算法的估計LDPC碼譯碼噪聲方差的方法，其具體實施步驟如下：

第一步：獲得信道譯碼環(huán)節(jié)的接收序列；

信道編譯碼及其噪聲方差估計的系統(tǒng)框圖如圖2所示，將信源輸入信息通過LDPC碼信道編碼獲得編碼序列u＝(u₁,u₂,…,u_K)，并進行BPSK調(diào)制獲得發(fā)送序列c＝(c₁,c₂,…,c_N)，發(fā)送序列c將編碼序列u中的0和1轉換為1和-1，具體計算方法如下：

c_k＝1-2·u_k

接著，在發(fā)送序列上按位疊加均值為0、方差為σ²的高斯白噪聲n＝(n₁,n₂,…,n_N)，從而獲得信道譯碼環(huán)節(jié)的接收序列r＝(r₁,r₂,…,r_N)。其中，K表示信源輸入的未編碼信息的長度，N表示編碼后碼字長度。

第二步：確定方差估計范圍和初值；

為合理地進行方差估計，需先確定方差估計范圍，再選擇方差估計的初值。噪聲方差與接收序列的信噪比E_b/N₀的關系如下：

其中信噪比E_b/N₀的單位為分貝(即dB)，R表示碼率，且滿足R＝K/N。由上式可知，當碼率確定時，接收序列信噪比越大，噪聲方差越小。假設信噪比最小值為0dB，方差估計范圍的上限為：

σ_max²＝0.5R

因此，方差估計范圍為(0,0.5R]。在該范圍下去任意隨機數(shù)可作為方差估計的初始值，即σ₀²∈(0,0.5R]。

第三步：基于EM算法估計噪聲方差；

基于EM算法估計信道噪聲方差的流程圖如圖3所示。

首先，設置EM算法迭代的初值，該初值為高斯混合模型的主要參數(shù)。接收序列是一組隨機序列，其概率分布符合高斯混合模型，概率密度函數(shù)如圖4所示，可表示為：

其中，r為接收序列，α₁和α₂代表單個接收碼元來自均值為1和-1的高斯模型的概率，α_k≥0且α₁+α₂＝1；φ(r|θ_k)是高斯概率密度函數(shù)，可表示為：

因此，當前高斯混合模型的參數(shù)主要包括噪聲方差噪聲均值μ_k和模型概率α_k?？蓪?shù)設置如下：μ₁＝1，μ₂＝-1，α₁＝α₂＝0.5，

設置完EM算法初值后，開始EM算法的迭代環(huán)節(jié)。首先，執(zhí)行E步：根據(jù)當前模型參數(shù)，計算分模型k對接收序列每個元素r_j的響應度，即當前模型參數(shù)下第j個接收數(shù)據(jù)r_j來自第k個模型的概率

接著，執(zhí)行M步：計算新一輪迭代的模型參數(shù)：

其中k＝1,2。上述的E步和M步的推導過程見《統(tǒng)計學習方法》書中的163-165頁。

每執(zhí)行完一次M步，計算方差估計值：

循環(huán)執(zhí)行E步和M步，直到方差估計值收斂時停止循環(huán)迭代，輸出噪聲方差估計值。

第四步：將信道噪聲方差估計結果應用于譯碼初始化和迭代環(huán)節(jié)；

根據(jù)接收序列和噪聲方差估計值進行LDPC碼譯碼初始化，即計算接收序列相應比特后驗概率的對數(shù)似然比：

執(zhí)行LDPC碼置信傳播譯碼算法的迭代過程，更新校驗節(jié)點和變量節(jié)點，直至硬判決結果滿足校驗方程或迭代次數(shù)達到上限時，結束譯碼，輸出譯碼結果。

為驗證本專利提出算法的有效性、合理性、可行性及科學性，利用本方法對實際的LDPC碼接收序列進行方差估計。圖5是不同信噪比下的方差估計值；圖6是初始方差為正確方差的1/2時長528的LDPC碼采用本發(fā)明的譯碼性能；圖7初始方差為正確方差的1/2時長1200的LDPC碼采用本發(fā)明的譯碼性能。試驗中各參數(shù)設置如下：LDPC碼(528,264)和(1200,600)的碼率為1/2；LDPC碼譯碼采用經(jīng)典的和積算法，最大迭代次數(shù)為50。

由試驗結果可以看出，本文算法可以有效地估計信道噪聲方差，保證LDPC碼譯碼的正確執(zhí)行；從圖5可以看出，基于EM算法的方差估計可以在已知接收序列的前提下較準確地估計高斯白噪聲方差，隨信噪比的改變，方差估計的準確性較高。從圖6可以看出，相比采用錯誤的初始方差，基于EM算法的方差估計可有效地降低譯碼性能衰減。結合圖6和圖7可以看出，在誤碼率為10^-5時，長度為1200的LDPC碼采用本方法相比直接使用初始方差譯碼所提升的編碼增益為0.3dB，其對譯碼性能的提升效果優(yōu)于長度為528的碼字。因此，碼字越長，方差估計越準確，能更有效地在方差未知時實現(xiàn)最優(yōu)譯碼性能。