王蘭勛，佟婧麗，孟祥雅

(河北大學電子信息工程學院，河北保定 071002)

在實際數(shù)字通信信道中，為了抗擊噪聲的影響，提高數(shù)據(jù)傳輸可靠性，需要進行信道編碼技術(shù)，即在信息序列中增加冗余碼元，使接收序列具有糾檢錯能力。對信息截獲方而言，如何根據(jù)接收的序列盲識別出信道編碼體制和參數(shù)是當前編碼研究領(lǐng)域中一個嶄新的研究課題，具有重要的現(xiàn)實意義和應(yīng)用價值［1－2］。

目前，根據(jù)公開發(fā)表文獻知，大部分信道編碼識別的研究主要集中在卷積碼的盲識別上，針對線性分組碼盲識別研究的文獻相對較少。文獻［3］根據(jù)碼重分布距離公式估計二進制線性分組碼的碼長，通過矩陣化簡得到生成矩陣，在低誤碼率下識別效果較好。文獻［4］提出了碼重分布信息熵算法識別線性分組碼碼長，適用于較高誤碼環(huán)境，但沒有考慮截獲序列非同步的情況。文獻［5］提出了利用比特頻率檢測法識別碼長和起始點，更適合于較低誤碼率的環(huán)境且所用數(shù)據(jù)量較大。文獻［6］提出了利用漢明距離識別碼長和起始點，所需數(shù)據(jù)量較少，但適應(yīng)于誤碼率較低的線性分組碼。文獻［7］提出了一種基于統(tǒng)計顯著水平的快速識別法，首次解決了BCH縮短碼的識別問題。文獻［8］提出一種利用歐幾里得算法識別最大公因式的方法，該方法無繁瑣計算。文獻［7－8］只針對BCH碼而言，針對性較強。文獻［9］提出一種通過解調(diào)輸出的軟判決序列來求解含錯方程的方法，所需數(shù)據(jù)量較大。

上述算法或者是數(shù)據(jù)量較大、針對性較強，不能識別同步點，或者是在低誤碼條件下才能達到較好的識別效果。針對以上方法的不足，本文提出碼重相似度算法識別碼長和同步點(本文碼字起始點簡稱同步點)。在此基礎(chǔ)上，通過計算碼字深度值，選取非零不同深度值對應(yīng)的碼字擺成矩陣形式，即為生成矩陣。最后對碼重相似度識別方法進行仿真分析，并討論該方法的容錯性，與其他方法進行了比較，仿真結(jié)果表明文中方法在較高誤碼率為0.01條件下能夠有效識別中短碼，容錯性較好，且所需數(shù)據(jù)量較少。

1 線性分組碼的相關(guān)定義及其描述

定義1:(n，k)線性分組碼是GF(q)上的n維線性空間中的k維子空間，由qk個碼字集合構(gòu)成，而k個獨立向量構(gòu)成的一組基底完全可以線性表示這qk個碼字，將這組基底寫成矩陣形式，即為生成矩陣G，G表示形式不唯一，但都是k行n列矩陣［10］。若G有如下形式:G=［Ik，P］，稱G為典型生成矩陣。對于循環(huán)碼，生成多項式的向量形式可由典型生成矩陣的第k行的第k列到第n列表示。本文以二進制線性分組碼為例，即q=2。

定義2:一個碼字的重量等于該碼字中非零元素的個數(shù)。對于(n，k)二進制線性分組碼，其碼字重量為碼字中含“1”的個數(shù)，也稱碼重［10］。設(shè)qi是 (n，k)分組碼中碼重為i的碼字個數(shù)，則重量分布表示為{q0，q1，…，qn}。碼重分布概率Pi是重量為i的碼字個數(shù)在碼字總數(shù)中出現(xiàn)的概率。

定義3:令c=(b1，b2，…，bn)是二元域上碼長為n的線性分組碼C的一個碼字，定義D(c)=(b2－b1，…，bnbn－1)，稱D為c的微分運算［11］。碼字c的深度是滿足Dz(c)=(0n－z)的最小非負整數(shù)z，若這樣的z不存在，則c的深度為n。其中Dz(c)表示對c進行z次微分運算。記Dj以j為深度的碼字個數(shù)，集合 {D0，D1，…，Dn}稱為碼c的深度分布。下標集{j|1＜j＜n，Dj≠0}稱為碼c的深度譜。

定理:有限域GF(q)上，任意一個(n，k)線性碼c的深度譜中都有且僅有k個非零值，深度互不相同的任意多個碼字向量是線性無關(guān)的［11］。

2 碼重相似度盲識別方法

線性分組碼碼字之間存在較強的線性約束關(guān)系，相比于隨機序列，并不是所有任意組合的碼字都會出現(xiàn)，導致碼重分布不平衡，其碼重分布概率與隨機序列的碼重分布概率之間差異很大，低碼率分組碼的這種差異程度更明顯，因此本文主要針對低碼率線性分組碼進行識別。將兩序列的碼重分布概率看成兩個向量，采用向量間的相似度衡量兩個碼重分布概率之間的差異程度。本文根據(jù)隨機序列與實際序列碼重分布概率間差異最大這一特性，提出一種基于碼重相似度識別線性分組碼碼長和同步點的方法。

在統(tǒng)計學中，積矩相關(guān)ρW1W2的平方稱為判定系數(shù)R，用判定系數(shù)來衡量實際與隨機序列分布的相似度，則W1與W2的相似度表達式為

式中:COV(W1，W2)表示W(wǎng)1與W2的協(xié)方差;D(W1)與D(W2)分別為W1與W2的方差，即

將式(3)～(5)帶入式(2)，將式(2)結(jié)果帶入式(1)得判定系數(shù)即相似度

同步點已知，識別碼長步驟如下:

1)初始化實際接收序列的碼長n。

2)接收序列按碼長n分為M個碼字。n的變化范圍設(shè)為1～S，S為最大可能的值。S的選取要根據(jù)實際應(yīng)用需求綜合考慮。

4)利用式(6)求R。判斷R值是否最小，若最小，則識別碼長或碼長整數(shù)倍，當首次出現(xiàn)谷值時對應(yīng)碼長為真實碼長;否則，n=n+1，轉(zhuǎn)到步驟2)，以此循環(huán)，直到R最小為止。

碼長已知，識別同步點步驟如下:

1)初始化接收序列的同步點m，從任意m位開始，以確定的碼長n劃分為M個碼字。

2)同識別碼長步驟3)。

3)利用式(6)求R。判斷R值是否最小，若最小，則識別同步點，否則，m=m+1，轉(zhuǎn)到步驟1)，以此循環(huán)，直到R最小為止，m的變化范圍是1～n+1。

3 深度分布性質(zhì)識別生成矩陣

用上述方法識別出碼長或同步點后，利用線性分組碼的深度分布特性識別G。由定義1知，只需找到k個不相關(guān)的向量組成矩陣形式，即為G。由定理知，深度值互不相同的非零碼字向量作為行向量生成的矩陣是G，根據(jù)此結(jié)論識別G。將G進行初等行變換，得到典型生成矩陣。對于循環(huán)碼，進而得到生成多項式。

該方法關(guān)鍵點是識別k值。(n，k)線性碼深度譜中有且僅有k個非零值，當已知碼長和同步點后，求出M個碼字中每個碼字的深度值，記為j，j∈{0，1，…，n}，將深度為j的碼字個數(shù)相加，結(jié)果為Dj，則深度分布為{D0，D1，…，Dn}，深度譜為{j|1＜j＜n，Dj≠0}，由定理知選出深度譜中所有不同深度值分別對應(yīng)的任意一個碼字，擺成矩陣形式，即為G。此方法簡單，過程易于實現(xiàn)。

4 仿真驗證及分析

4.1 已知同步點識別碼長

當 BSC 信道無誤碼時，采用(6，3)，(15，5)，(31，11)，(63，18)，(127，50)5 種線性分組碼進行實驗，用 MATLAB產(chǎn)生一段隨機序列進行編碼。取10 000個碼元進行仿真，結(jié)果如圖1所示。

圖1 5種低碼率線性分組碼的碼長識別

經(jīng)圖1仿真分析，當遍歷的碼長是真實碼長或碼長整數(shù)倍時，碼組內(nèi)有完整的線性約束關(guān)系，不同碼重的碼字分布是非等概的，與隨機序列碼重分布相差最大，相似度最小，即判定系數(shù)最小，此時，相似度相對于臨近點其值變化較大，因此當首次出現(xiàn)低谷時對應(yīng)的碼長為真實碼長。當遍歷的碼長不是真實碼長或碼長整數(shù)倍時，碼組內(nèi)不具有完整的線性約束關(guān)系，不同碼重的碼字分布趨于等概率，接近隨機序列碼重分布，此時實際序列與隨機序列相似度較高。經(jīng)仿真分析該方法能正確識別碼長n。

當BSC信道有誤碼時，以(15，5)線性分組碼為實驗對象，碼元個數(shù)為10 000個，誤碼率分別為Pe=0.001，0.01，0.05，0.08，仿真結(jié)果如圖2 所示。

圖2 不同誤碼率對應(yīng)的仿真圖

從圖2可知，碼重相似度識別算法受誤碼率的影響。常規(guī)誤碼率為0.001～0.01，在常規(guī)誤碼率下可以正確識別碼長，之后隨著誤碼率的增加，在正確碼長或碼長整數(shù)倍處判定系數(shù)變化相對緩慢，在碼長整數(shù)倍處可能不會出現(xiàn)谷值，當誤碼率增大到一定程度，出現(xiàn)失真，不能識別碼長。如在Pe=0.08時正確識別率很低，由給出的仿真圖勉強識別碼長。

4.2 已知碼長識別同步點

設(shè)BSC信道有誤碼，以誤碼率為Pe=0.03的(6，3)，(15，5)，(31，11)和Pe=0.01 的(63，18)4 種線性分組碼為實驗對象，碼元個數(shù)為10 000個，在碼長n已知條件下，從任意起點開始，順序遍歷n個點，在各點處均以碼長n劃分序列，并統(tǒng)計判定系數(shù)，尋找判定系數(shù)為最小值時對應(yīng)的點，即為同步點。用MATLAB進行仿真，結(jié)果如圖3所示。

圖3 碼字同步點識別仿真圖

計算同步點由m=1遍歷到m=n+1對應(yīng)的判定系數(shù)，并作分析，其中m=1表示起點，沒有移位，m=n+1表示從起點向右移n位，由于碼長為n，所以m=1與m=n+2對應(yīng)的判定系數(shù)是一樣的，因此僅比較起點和移位n個起點即可。由圖3可知，m分別為5，11，17，22處判定系數(shù)最小，即碼重相似度最低，碼字內(nèi)線性約束關(guān)系最強，此時對應(yīng)的m值為同步點;而在其他位置，隨機性比較強，與隨機序列相似度較高。因此在誤碼率較高的條件下，運用碼重相似度方法也能準確識別碼字同步點。

4.3 深度分布性質(zhì)識別生成矩陣

假設(shè)截獲序列已知碼長和同步點，以(15，5)二進制線性分組碼為例，取200個碼字。計算每個碼字的深度值，統(tǒng)計深度值相同的碼字個數(shù)。通過MATLAB編程，產(chǎn)生深度分布直方圖如圖4所示，圖4中橫坐標為深度值j，縱坐標為深度值相同的碼字個數(shù)之和Dj。

圖4對應(yīng)的深度分布為{4 10 0 0 0 0 0 0 0 0 0 0 10 22 57 97}，深度譜為{1，12，13，14，15}，根據(jù)定理知，此線性碼k=5即生成矩陣有5行，從這k個非零不同深度值對應(yīng)的碼字中，各取一個碼字，組成5×15的矩陣形式，即為生成矩陣，對其初等行變換得到典型生成矩陣。任取k=1對應(yīng)碼字為(111111111111111)，k=12對應(yīng)碼字為(101011001000111)，k=13對 應(yīng) 碼 字 為(100110111000010)，k=14對 應(yīng) 碼 字 為(111010110010001)，k=15對 應(yīng) 碼 字 為(111000010100110)，將此組成的矩陣形式如圖5所示。

圖4 (15，5)線性分組碼深度分布仿真圖

圖5 矩陣形式

從圖5可知，G3為典型生成矩陣。本節(jié)實驗對象(15，5)也屬于循環(huán)碼，G3最后一行為生成多項式的系數(shù)，即循環(huán)碼的生成多項式為g(x)=x10+x8+x5+x4+x2+x+1。利用深度分布識別生成矩陣的方法簡單，過程易于實現(xiàn)。

5 容錯性能比較分析

根據(jù)上述分析，碼重相似度識別碼長和同步點方法是有效的。下面以(7，4)，(15，5)，(31，16)，(63，16)這 4種線性分組碼為實驗對象，分別對識別碼長和同步點兩種情況進行容錯性分析。

針對本文和文獻［5］兩種識別同步點的方法在相同環(huán)境下比較分析，對1 000個碼字在不同誤碼率下進行1 000次蒙特卡洛仿真實驗，各種線性碼的識別率與誤碼率曲線如圖6、圖7所示:圖6為碼重相似度識別同步點的仿真圖，明顯看出，當誤碼率為0.18 時，(7，4)，(15，5)分組碼的同步點識別率超過了90%;誤碼率為0.06時，(31，16)分組碼的同步點識別率高達90%;誤碼率為0.016時，(63，16)識別率超過90%;因此本文算法能在高誤碼率下識別同步點。圖7為比特頻率檢測法識別同步點，其容錯性較好，但本文算法相對于比特頻率檢測法容錯性更有優(yōu)勢。

圖6 碼重相似度算法識別碼字同步點

圖7 比特頻率檢測方法識別碼字同步點

根據(jù)文獻［4］知碼重信息熵算法容錯性強于碼重分布距離［3］算法，在此僅對文獻［4－5］和本文提出的識別碼長算法進行比較。分別取1 000個4種長度的線性分組碼碼字為實驗對象，進行仿真，圖8為各種算法的誤碼適應(yīng)能力［13］曲線圖，圖中的誤碼適應(yīng)能力是以正確識別率90%以上為門限，不難看出，本文算法與以往算法相比，誤碼適應(yīng)能力有一定的提升，在適應(yīng)誤碼率1%以上識別碼長的正確率高達90%以上，容錯性較好。

圖8 不同算法誤碼適應(yīng)能力

6 結(jié)論

根據(jù)線性分組碼碼重分布不平衡與隨機序列差異較大的特性，提出了一種碼重相似度算法識別碼長和同步點。在此基礎(chǔ)上，利用深度分布特性，識別生成矩陣，利用高斯消去法產(chǎn)生典型生成矩陣，針對循環(huán)碼，識別生成多項式，實現(xiàn)了線性分組碼的盲識別。最后在不同誤碼率下對低碼率線性分組碼進行大量仿真實驗，并與其他算法進行對比，碼重相似度識別算法能在誤碼率為0.01情況下識別低碼率線性分組碼碼長和同步點，具有較好的容錯性，此方法的誤碼適應(yīng)能力也有一定的提高。利用深度分布特性識別出生成矩陣，此方法簡單，隨著碼長增加數(shù)據(jù)量不斷擴大，此識別生成矩陣的方法對無誤碼情況下的中短碼及低碼率的線性分組碼識別效果較好。運用深度分布特性識別高誤碼率下、高碼率的線性分組碼參數(shù)成為下一步研究的重點。

:

［1］解輝，黃知濤，王豐華.信道編碼盲識別技術(shù)研究進展［J］.電子學報，2013，41(6):1166－1176.

［2］王偉祥，劉玉君.基于信道編碼的信息隱藏技術(shù)研究［J］.電視技術(shù)，2006，30(3):284－286.

［3］昝俊軍，李艷斌.低碼率二進制線性分組碼的盲識別［J］.無線電工程，2009，39(1):19－24.

［4］陳金杰，楊俊安.基于碼重信息熵低碼率線性分組碼的盲識別［J］.電路與系統(tǒng)學報，2012，17(1):41－47.

［5］陳金杰，楊俊安.一種對線性分組碼編碼參數(shù)的盲識別方法［J］.電路與系統(tǒng)學報，2013，18(2):248－254.

［6］閆郁翰，湯建龍.基于漢明碼距離的二進制線性分組碼盲識別方法［J］.通信對抗，2011(4):20－23.

［7］張永光，鄭仕鏈.BCH碼的參數(shù)識別研究［J］.西安電子科技大學學報，2013，40(5):92－98.

［8］王蘭勛，李丹芳，汪洋.二進制本原BCH碼的參數(shù)盲識別［J］.河北大學學報，2012，32(4):416－420.

［9］于沛東，李靜，彭華.一種利用軟判決的信道編碼識別新算法［J］.電子學報，2013(2):301－306.

［10］王新梅，肖國鎮(zhèn).糾錯碼——原理與方法［M］.修訂版.西安:西安電子科技大學出版社，2001.

［11］李超，耿晉.有限域上線性碼的深度分布［D］.長沙:國防科學技術(shù)大學，2006.

［12］王磊，胡以華，王勇，等.基于碼重分布的系統(tǒng)循環(huán)碼識別方法［J］.計算機工程與應(yīng)用，2012，48(7):150－153.

［13］解輝，王豐華，黃知濤，等.基于頻譜預處理的RS碼盲檢測識別方法［J］.宇航學報，2013，34(1):128－132.