高速率低密度奇偶檢驗碼的線性規劃譯碼方法技術

本發明專利技術公開了一種高速率低密度奇偶檢驗碼的線性規劃譯碼方法，解決了現有技術的LDPC碼譯碼方法譯碼速率慢的問題。本發明專利技術實現方法的步驟：譯碼參數初始化；計算初始變量節點消息；計算校驗節點消息殘余；更新校驗節點消息；更新變量節點消息；重新計算校驗節點消息殘余；譯碼終止判決；譯碼結果輸出。本發明專利技術通過在每次譯碼迭代中僅更新一個對譯碼收斂作用最大的校驗節點的校驗節點消息以及與該校驗節點相鄰的變量節點的變量節點消息，使得譯碼速率大大提高。

【技術實現步驟摘要】
高速率低密度奇偶檢驗碼的線性規劃譯碼方法
本專利技術屬于通信
，更進一步涉及信道編碼
中一種高速率低密度奇偶檢驗(Low-DensityParity-Check,LDPC)碼的線性規劃譯碼方法。本專利技術可實現基于ADMM算法的低密度奇偶檢驗LDPC碼高速率線性規劃譯碼。
技術介紹
低密度奇偶檢驗LDPC碼譯碼復雜度低，有逼近香農限的良好性能，被廣泛應用于現代通信領域中，如深空通信、無線通信等，并被各種現代通信標準采納，如802.11n、802.16e、10GBASE-T等，是近年信道編碼領域令人矚目的研究熱點。低密度奇偶檢驗LDPC碼譯碼方法分為兩類：消息傳遞類方法和線性規劃譯碼方法。消息傳遞類方法，如置信度傳播(BeliefPropagation,BP)等，具有實現簡單、譯碼復雜度較低等優勢，但是存在易受短環影響、難以數學分析等缺點。線性規劃譯碼方法基于凸優化理論，具有最大似然保證的特點，便于數學分析，但是由于沒有充分利用低密度奇偶檢驗LDPC碼校驗矩陣稀疏的特點，譯碼復雜度較高。山東大學擁有的專利技術“一種應用于LDPC碼的自適應線性規劃譯碼算法”(授權公告日：2013年6月5日，授權公告號：CN102148619B)中公開了一種自適應低密度奇偶檢驗LDPC碼譯碼方法。該專利申請中的譯碼方法通過自適應地添加對當前錯誤解有效的奇偶校驗，有目的的收緊線性規劃的可行域范圍，因此相對于傳統的低密度奇偶檢驗LDPC碼線性規劃譯碼，整個譯碼過程中所用到的奇偶校驗大大減少。該方法雖然解決低密度奇偶檢驗LDPC碼的線性規劃譯碼方法譯碼速率低的問題，但是仍然存在的不足之處是，譯碼過程中需要利用標準線性規劃問題求解方法多次求解線性規劃問題，因此該譯碼方法雖然能提高譯碼速率，但提高有限。Barman等人在“DecompositionmethodsforlargescaleLPdecoding”(IEEEInternationalConferenceonDataEngineering,Hannover,2011:253-260.)提出了一種基于交替方向乘子法(AlternatingDirectionMethodofMultiplier,ADMM)的低密度奇偶檢驗LDPC碼線性規劃譯碼方法。該方法通過重新設計適用于低密度奇偶檢驗LDPC碼譯碼的線性規劃問題求解方法，使譯碼方法的譯碼效率相對于利用標準線性規劃問題求解方法的譯碼方法大大提高。該方法存在的不足是：基于ADMM的低密度奇偶檢驗LDPC碼線性規劃譯碼方法每次迭代中所有校驗節點的校驗節點消息和變量節點的變量節點消息均需要更新，且更新校驗節點的校驗節點信息需要執行時間復雜度較高的歐幾里得投影運算，譯碼速度慢且計算量大，譯碼效率不高。
技術實現思路
本專利技術的目的在于針對上述已有技術的不足，提出一種高速率低密度奇偶檢驗LDPC碼線性規劃譯碼方法，可實現低密度奇偶檢驗LDPC碼線性規劃譯碼方法的快速收斂，提高譯碼效率。為了實現上述目的，本專利技術方法的思路是：改變原有譯碼方法每次迭代中所有校驗節點信息和變量節點信息均需要更新的更新策略，通過對校驗節點信息殘余進行計算，每次僅更新最大校驗消息殘余的校驗節點的校驗節點消息以及與最大校驗消息殘余的校驗節點相關的變量節點的變量節點消息，從而減少不必要的消息更新。本專利技術方法的實現步驟如下：(1)譯碼參數初始化：(1a)按照下式，計算低密度奇偶檢驗碼的碼字比特對應的代價參數：γi＝log(Pr(ri|ci＝0)/Pr(ri|ci＝1))其中，γi表示低密度奇偶檢驗碼的碼字比特i對應的代價參數，i的取值范圍為[1，n]，n表示低密度奇偶校驗碼的長度，ri表示加性高斯白噪聲信道輸出的第i位信息，ci表示發送的第i位信息，Pr(·)表示括號內表示的事件發生概率；(1b)設置當前迭代次數k＝0，設置最大迭代次數N＝1000、容差值ε＝10-5；(2)按照初始變量節點消息的公式，計算每個變量節點的初始變量節點消息；(3)按照下式，計算每個校驗節點的校驗節點消息殘余：r(Lj→i)＝Tjx1其中，r(Lj→i)表示校驗節點j對相鄰變量節點i的校驗節點信息殘余，j表示校驗節點，i表示變量節點，變量節點i的索引范圍是與校驗節點j相鄰的變量節點的索引集，Tj表示校驗節點j對應的轉換矩陣，x1表示初始變量節點消息向量；(4)更新校驗節點消息：(4a)從校驗節點中找出最大校驗節點信息殘余的校驗節點；(4b)按照下式，計算第k+1次迭代時與最大校驗節點消息殘余的校驗節點對應的輔助向量：其中，表示第k+1次迭代時與最大校驗節點消息殘余的校驗節點m對應的輔助向量，Pm表示最大校驗節點消息殘余的校驗節點m對應的校驗多胞體，表示在最大校驗節點信息殘余的校驗節點m對應的校驗多胞體Pm上做歐幾里得投影操作，Tm表示最大校驗節點信息殘余的校驗節點m對應的轉換矩陣，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量；(4c)按照下式，計算第k+1次迭代時與最大校驗節點消息殘余的校驗節點對應的拉格朗日乘子向量：其中，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量，Tm表示最大校驗節點信息殘余的校驗節點m對應的轉換矩陣，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的輔助向量；(4d)按照下式，計算第k+1次迭代時最大校驗節點信息殘余的校驗節點消息：其中，表示第k+1次迭代時最大校驗節點信息殘余的校驗節點m向相鄰變量節點i發送的校驗節點信息，m表示最大校驗節點信息殘余的校驗節點，i表示變量節點，變量節點i的索引范圍是與最大校驗節點信息殘余的校驗節點m相鄰的變量節點的索引集，表示第k+1次迭代時與最大校驗節點消息殘余的校驗節點m對應的輔助向量zm中與第i個變量節點對應的向量位，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量ym中與第i個變量節點對應的向量位；(4e)將最大校驗節點信息殘余的校驗節點的校驗節點消息殘余歸零；(5)按照變量節點消息的公式，計算與最大校驗節點消息殘余的校驗節點相鄰的變量節點的變量節點消息；(6)按照下式，計算與變量節點消息發生改變的變量節點相鄰的校驗節點的校驗節點消息殘余：其中，表示第k+1次迭代時校驗節點j對相鄰變量節點i的校驗節點消息殘余，j表示校驗節點，i表示變量節點，變量節點i的索引范圍是與校驗節點j相鄰的變量節點的索引集，Tj表示校驗節點j對應的轉換矩陣，xk+2表示第k+2次迭代時的變量節點消息向量，表示第k+1次迭代時與第個j校驗節點對應的拉格朗日乘子向量，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與第個j校驗節點對應的拉格朗日乘子向量；(7)判斷是否滿足譯碼終止條件，若是，則執行步驟(8)，否則，執行步驟(4)；(8)將變量節點消息向量作為譯碼結果輸出。本專利技術與現有技術相比具有以下優點：第一，由于本專利技術在每次迭代中僅更新一個校驗節點的校驗節點消息以及與該校本文檔來自技高網...

【技術保護點】
高速率低密度奇偶檢驗碼的線性規劃譯碼方法，包括如下步驟：(1)譯碼參數初始化：(1a)按照下式，計算低密度奇偶檢驗碼的碼字比特對應的代價參數：γi＝log(Pr(ri|ci＝0)/Pr(ri|ci＝1))其中，γi表示低密度奇偶檢驗碼的碼字比特i對應的代價參數，i的取值范圍為[1，n]，n表示低密度奇偶校驗碼的長度，ri表示加性高斯白噪聲信道輸出的第i位信息，ci表示發送的第i位信息，Pr(·)表示括號內表示的事件發生概率；(1b)設置當前迭代次數k＝0，設置最大迭代次數N＝1000、容差值ε＝10?5；(2)按照初始變量節點消息的公式，計算每個變量節點的初始變量節點消息；(3)按照下式，計算每個校驗節點的校驗節點消息殘余：r(Lj→i)＝Tjx1其中，r(Lj→i)表示校驗節點j對相鄰變量節點i的校驗節點信息殘余，j表示校驗節點，i表示變量節點，變量節點i的索引范圍是與校驗節點j相鄰的變量節點的索引集，Tj表示校驗節點j對應的轉換矩陣，x1表示初始變量節點消息向量；(4)更新校驗節點消息：(4a)從校驗節點中找出最大校驗節點信息殘余的校驗節點；(4b)按照下式，計算第k+1次迭代時與最大校驗節點消息殘余的校驗節點對應的輔助向量：zmk+1=Πpm(Tmxk+1+ymk)]]>其中，表示第k+1次迭代時與最大校驗節點消息殘余的校驗節點m對應的輔助向量，Pm表示最大校驗節點消息殘余的校驗節點m對應的校驗多胞體，表示在最大校驗節點信息殘余的校驗節點m對應的校驗多胞體Pm上做歐幾里得投影操作，Tm表示最大校驗節點信息殘余的校驗節點m對應的轉換矩陣，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量；(4c)按照下式，計算第k+1次迭代時與最大校驗節點消息殘余的校驗節點對應的拉格朗日乘子向量：ymk+1=Tmxk+1+ymk-zmk+1]]>其中，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量，Tm表示最大校驗節點信息殘余的校驗節點m對應的轉換矩陣，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的輔助向量；(4d)按照下式，計算第k+1次迭代時最大校驗節點信息殘余的校驗節點消息：Lm→ik+1=(zmk+1)i-(ymk+1)i]]>其中，表示第k+1次迭代時最大校驗節點信息殘余的校驗節點m向相鄰變量節點i發送的校驗節點信息，m表示最大校驗節點信息殘余的校驗節點，i表示變量節點，變量節點i的索引范圍是與最大校驗節點信息殘余的校驗節點m相鄰的變量節點的索引集，表示第k+1次迭代時與最大校驗節點消息殘余的校驗節點m對應的輔助向量zm中與第i個變量節點對應的向量位，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量ym中與第i個變量節點對應的向量位；(4e)將最大校驗節點信息殘余的校驗節點的校驗節點消息殘余歸零；(5)按照變量節點消息的公式，計算與最大校驗節點消息殘余的校驗節點相鄰的變量節點的變量節點消息；(6)按照下式，計算與變量節點消息發生改變的變量節點相鄰的校驗節點的校驗節點消息殘余：r(Lj→ik+1)=(Tjxk+2+yjk+1)i-(Tjxk+1+yjk)i]]>其中，表示第k+1次迭代時校驗節點j對相鄰變量節點i的校驗節點消息殘余，j表示校驗節點，i表示變量節點，變量節點i的索引范圍是與校驗節點j相鄰的變量節點的索引集，Tj表示校驗節點j對應的轉換矩陣，xk+2表示第k+2次迭代時的變量節點消息向量，表示第k+1次迭代時與第個j校驗節點對應的拉格朗日乘子向量，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與第個j校驗節點對應的拉格朗日乘子向量；(7)判斷是否滿足譯碼終止條件，若是，則執行步驟(8)，否則，執行步驟(4)；(8)將變量節點消息向量作為譯碼結果輸出。...

【技術特征摘要】
1.高速率低密度奇偶檢驗碼的線性規劃譯碼方法，包括如下步驟：(1)譯碼參數初始化：(1a)按照下式，計算低密度奇偶檢驗碼的碼字比特對應的代價參數：γi＝log(Pr(ri|ci＝0)/Pr(ri|ci＝1))其中，γi表示低密度奇偶檢驗碼的碼字比特i對應的代價參數，i的取值范圍為[1，n]，n表示低密度奇偶校驗碼的長度，ri表示加性高斯白噪聲信道輸出的第i位信息，ci表示發送的第i位信息，Pr(·)表示括號內表示的事件發生概率；(1b)設置當前迭代次數k＝0，設置最大迭代次數N＝1000、容差值ε＝10-5；(2)按照初始變量節點消息的公式，計算每個變量節點的初始變量節點消息；所述初始變量節點消息的公式如下：其中，表示第i個變量節點的初始變量節點消息，Π[0,1]表示在區間[0,1]內做歐幾里得投影操作，γi表示與第i個變量節點對應的初始化后的譯碼參數，ρ表示針對具體信道、具體碼字，根據誤碼率大小優化得到的拉格朗日參數，di表示與變量節點i相鄰的校驗節點的個數；(3)按照下式，計算每個校驗節點的校驗節點消息殘余：r(Lj→i)＝Tjx1其中，r(Lj→i)表示校驗節點j對相鄰變量節點i的校驗節點信息殘余，j表示校驗節點，i表示變量節點，變量節點i的索引范圍是與校驗節點j相鄰的變量節點的索引集，Tj表示校驗節點j對應的轉換矩陣，x1表示初始變量節點消息向量；(4)更新校驗節點消息：(4a)從校驗節點中找出最大校驗節點信息殘余的校驗節點；(4b)按照下式，計算第k+1次迭代時與最大校驗節點消息殘余的校驗節點對應的輔助向量：其中，表示第k+1次迭代時與最大校驗節點消息殘余的校驗節點m對應的輔助向量，pm表示最大校驗節點消息殘余的校驗節點m對應的校驗多胞體，表示在最大校驗節點信息殘余的校節點m對應的校驗多胞體pm上做歐幾里得投影操作，Tm表示最大校驗節點信息殘余的校驗節點m對應的轉換矩陣，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量；(4c)按照下式，計算第k+1次迭代時與最大校驗節點消息殘余的校驗節點對應的拉格朗日乘子向量：其中，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量，Tm表示最大校驗節點信息殘余的校驗節點m對應的轉換矩陣，xk+1表示第k+1次迭代時變量節點消息向量，表示第k次迭代時與最大校驗節點信息殘余的校驗節點m對應的拉格朗日乘子向量，表示第k+1次迭代時與最大校驗節點信息殘余的校驗節點m對應的輔助向量；(4d)按照下式，...

【專利技術屬性】
技術研發人員：慕建君，王彪，魏浩源，范慶輝，焦曉鵬，
申請(專利權)人：西安電子科技大學，
類型：發明
國別省市：陜西;61