戚 林 ，郝士琦 ，王 磊

(1.電子工程學(xué)院 脈沖功率激光技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，安徽 合肥 230037；2.電子工程學(xué)院 安徽省電子制約技術(shù)重點(diǎn)實(shí)驗(yàn)室，安徽 合肥 230037)

在空間數(shù)據(jù)通信中，由于噪聲干擾及信號(hào)衰落等因素的影響，通信質(zhì)量得不到保障。為此，通常采用信道編碼技術(shù)來(lái)降低傳輸中的誤碼率，提高通信質(zhì)量。國(guó)際空間數(shù)據(jù)系統(tǒng)咨詢(xún)委員會(huì)CCSDS(Consultative Committee for Space Data Systems)根據(jù)空間數(shù)據(jù)通信的特點(diǎn)，制定了適合空間數(shù)據(jù)傳輸?shù)男诺谰幋a標(biāo)準(zhǔn)，其中規(guī)定的刪除卷積碼是通過(guò)刪除(2，1，6)卷積碼的碼字中某些特指位置的碼元而得到的，具有編碼方式靈活、帶寬利用率較高、編碼冗余較低等特點(diǎn)，在空間數(shù)據(jù)通信中得到了廣泛的應(yīng)用。美國(guó)國(guó)防部的天基紅外系統(tǒng)SBDRS-High、法國(guó)國(guó)防部的衛(wèi)星星座ESSAIM、英國(guó)國(guó)防部的戰(zhàn)術(shù)光學(xué)衛(wèi)星TopSae等[1]都采用了CCSDS的編碼標(biāo)準(zhǔn)。

在空間數(shù)據(jù)信號(hào)截獲分析領(lǐng)域中，實(shí)現(xiàn)CCSDS標(biāo)準(zhǔn)的刪除卷積碼的盲解碼，主要是要知道刪除卷積碼的刪除圖案和刪除位置。參考文獻(xiàn)[2]中介紹了一種對(duì)(2，1，6)卷積碼的識(shí)別和碼字同步方法，可以將這種方法推廣到刪除卷積碼中。參考文獻(xiàn)[3]介紹了一種對(duì)刪除卷積碼的刪除圖案的識(shí)別方法，但不能識(shí)別碼頭。本文對(duì)此方法加以改進(jìn)，通過(guò)理論推導(dǎo)得出刪除卷積碼的等效校驗(yàn)矩陣，再以等效校驗(yàn)矩陣為先驗(yàn)知識(shí)，來(lái)識(shí)別刪除卷積碼的參數(shù)，并進(jìn)一步提出了刪除卷積碼的盲解碼方法。

1 CCSDS刪除卷積碼編譯碼原理

1.1 編碼原理

CCSDS選擇的刪除卷積碼是以碼率為1/2，約束長(zhǎng)度為 7，生成多項(xiàng)式矩陣為 G(D)=[g(1，1)，g(1，2)]。 (其中，g(1，1)=1+D+D2+D3+D6，g(1，2)=1+D2+D3+D5+D6)的(2，1，6)卷積碼作為源碼，再按照刪除圖案分別對(duì)C1、C2兩組碼序列刪除后得到的[4]。CCSDS標(biāo)準(zhǔn)的刪除卷積碼的編碼器框圖如圖1所示，刪除圖案如表1所示。

圖1 CCSDS標(biāo)準(zhǔn)刪除卷積碼的編碼器框圖

表1 不同碼率的刪除圖案

刪除圖案中“1”表示這一位不刪除，“0”表示這一位刪除。如當(dāng)碼率r=2/3時(shí)，如果輸入序列M=1101100111，源碼的輸出序列為 C1=1001101100，C2=1111000100，經(jīng)過(guò)刪除圖案刪除后，源碼的輸出序列是C1=10110，C2=1111000100，刪除卷積碼的輸出序列 C=111011100101000。

1.2 譯碼原理

刪除卷積碼的譯碼過(guò)程是：將接收的碼序列安裝發(fā)送端的刪除圖案，在刪除圖案為“0”的位置插入“0”，然后輸入(2，1，6)卷積碼的 Viterbi譯碼器進(jìn)行譯碼[5]，刪除卷積碼的譯碼框圖如圖2所示。在分析截獲的刪除卷積碼時(shí)，并不知道刪除圖案，同時(shí)也不能確定截獲的碼段的起始位置即為發(fā)送碼字的起始位置，進(jìn)而即使在知道刪除圖案的情況下也不能明確插入比特的位置。因此，實(shí)現(xiàn)刪除卷積碼的盲解碼的關(guān)鍵是識(shí)別刪除圖案和定位刪除位置。

圖2 刪除卷積碼譯碼器框圖

2 CCSDS刪除卷積碼的盲解碼方法

2.1 刪除卷積碼的等效校驗(yàn)矩陣

[6]中，通過(guò)將碼率為(n0-1)/n0的刪除卷積碼等效成碼率為(n0-1)/n0的卷積碼，按照求解卷積碼校驗(yàn)矩陣的方法，求得刪除卷積碼的等效校驗(yàn)矩陣，并且證明每個(gè)碼率為(n0-1)/n0的刪除卷積碼只有唯一的等效校驗(yàn)矩陣。采用這種方法，可以求得CCSDS標(biāo)準(zhǔn)的刪除卷積碼的等效校驗(yàn)矩陣，如表2所示。由于等效校驗(yàn)矩陣與各種刪除圖案是一一對(duì)應(yīng)的，因此，可以將對(duì)刪除圖案的識(shí)別轉(zhuǎn)化為對(duì)等效校驗(yàn)矩陣的識(shí)別。

表2 各碼率的等效校驗(yàn)矩陣

2.2 刪除圖案和碼頭識(shí)別

建立碼率為 (n0-1)/n0的刪除卷積碼識(shí)別的模型如圖3所示。模型中在時(shí)間t=i時(shí)輸出n0個(gè)碼字，(n0-1)/n0卷 積 碼 的 校 驗(yàn) 多 項(xiàng) 式 矩 陣 為 ：H(D)=[H(1，1)(D)，H(1，2)(D)，…，H(1，n0)(D)]，其中，H(1，i)(D)=h0i+h1iD+h2iD2+…h(huán)LiDL，（其中L為編碼存儲(chǔ)）。由卷積碼的性質(zhì)[7]可知：C(D)HT(D)=0，則：

將其轉(zhuǎn)化為矩陣的形式：

圖3 刪除卷積碼識(shí)別模型

可以將上式看成C·H=0的形式，矩陣可以通過(guò)接收到的碼序列進(jìn)行構(gòu)造，因此通過(guò)高斯消元法[8]求解式（1）即可得到矩陣H，進(jìn)而得到刪除卷積碼的等效校驗(yàn)矩陣，然后再按表2中等效校驗(yàn)矩陣與刪除圖案的關(guān)系，就能識(shí)別出刪除圖案。

如果接收到的碼段不能同步，假設(shè)接收到的碼段列是：

在分析時(shí)，按照已同步的碼序列構(gòu)造矩陣，得：

將矩陣代入式（1）中，解得：

由此可以看出，以c0m開(kāi)頭的接收碼段，所求解出的H′正好是以c00開(kāi)頭的碼段所解出的H向前周期循環(huán)m×(L+1)位。反過(guò)來(lái)看，在分析截獲的碼段時(shí)，解出的H′與表2中的 H對(duì)比，發(fā)生了m×(L+1)位周期循環(huán)，即可判斷接收的碼字是以c0m開(kāi)頭。

2.3 盲解碼算法

綜合上述，CCSDS標(biāo)準(zhǔn)卷積碼的盲解碼算法步驟為：

(1)理論計(jì)算CCSDS標(biāo)準(zhǔn)卷積碼各碼率的校驗(yàn)矩陣H，將其作為先驗(yàn)知識(shí)。

(2)初始化，取 n0=8，L≥6，構(gòu)造矩陣 C′。

(3)計(jì)算式(1)，求得接收的碼段的等效校驗(yàn)矩陣H′。

(4)將H′與(1)求得的碼率為7/8的 H進(jìn)行比較。如果不匹配，則返回(2)，再取 n0=6、4、3，直到求得的 H′與(1)求得的對(duì)應(yīng)碼率的H匹配為止。此時(shí)，接收碼段的刪除圖案即為H對(duì)應(yīng)的刪除圖案。

(5)將 H′與 H進(jìn)行比對(duì)，可以得到 H′的循環(huán)量為m×(L+1)，則接收碼段的第(n0-m)位即為碼頭，再按照(4)識(shí)別出來(lái)的刪除圖案插入虛擬比特，然后輸入(2，1，6)卷積碼的Viterbi譯碼器，最終實(shí)現(xiàn)CCSDS標(biāo)準(zhǔn)的刪除卷積碼的盲解碼。

3 仿真結(jié)果分析

選取碼率為r=2/3的CCSDS標(biāo)準(zhǔn)的刪除卷積碼來(lái)驗(yàn)證上述算法。假設(shè)信息序列M=1011100011…經(jīng)過(guò)編碼后得到r=2/3刪除卷積碼C=111000101011111…。如果截獲到的數(shù)據(jù)為C′=00101011111…，應(yīng)用MATLAB進(jìn)行仿真。由于并不知道編碼存儲(chǔ)，取較大的L=8，n0從8開(kāi)始取值，取到n0=3時(shí)，計(jì)算結(jié)果如圖4所示。

圖4 仿真結(jié)果

圖4中倒數(shù)第二列即為計(jì)算出的校驗(yàn)矩陣H′100010100011100100010110110。 將 H′平均分成 3組，每組有 9 位，H′=[100010100 011100100 010110110]， 很顯然與r=2/3的刪除卷加碼的校驗(yàn)矩陣H=[1011011 1000101 0111001]相匹配，且剛好發(fā)生了9位循環(huán)，這樣接收到的碼段的第2位即為碼頭。此時(shí)的刪除矩陣P=[1 1 0 1]，插入虛擬碼字后為 C′=000100101011101…，將其輸入 (2，1，6)卷積碼的 Viterbi譯碼器即可實(shí)現(xiàn)解碼。由式(1)可知，構(gòu)造矩陣C需要截獲連續(xù)(L+1)×(n0+1)×n0-n0位沒(méi)有誤碼的碼段，對(duì)于CCSDS標(biāo)準(zhǔn)的刪除卷積碼，需要的最長(zhǎng)無(wú)誤碼位數(shù)是496位，而數(shù)字通信系統(tǒng)的誤碼率一般為10-3～10-4。因此，可以通過(guò)依次截取獲得的碼序列來(lái)構(gòu)造矩陣C，再運(yùn)用上述方法，必然能找到無(wú)誤碼碼段的正確解。本方法在高斯信道下對(duì)該盲解碼算法進(jìn)行了仿真，得到不同碼率下的盲解碼誤碼率如圖5所示。當(dāng)不能識(shí)別傳輸?shù)拇a率時(shí)，設(shè)盲解碼的誤碼率為0.5，從仿真結(jié)果可以看出，每種碼率的刪除卷積碼在能識(shí)別的臨界信噪比處，誤碼率值有個(gè)階躍，這是由該方法需要搜索一定長(zhǎng)的連續(xù)無(wú)誤碼碼段導(dǎo)致的，當(dāng)信噪比較低時(shí)不能找到滿(mǎn)足條件的碼段。當(dāng)SNR＞5 dB時(shí)，能對(duì)所有CCSDS標(biāo)準(zhǔn)的刪除卷積碼進(jìn)行盲解碼，而一般的數(shù)字通信系統(tǒng)要求信噪比要達(dá)到10 dB，因此該方法在實(shí)際中具有較高的實(shí)用性。

本文分析了CCSDS標(biāo)準(zhǔn)的刪除卷積碼的編譯碼原理，給出了該類(lèi)刪除卷積碼的等效校驗(yàn)矩陣，并利用等效校驗(yàn)矩陣與各種碼率間唯一對(duì)應(yīng)的特點(diǎn)，提出了一種盲解碼的方法。仿真結(jié)果表明，該方法在信噪比較低的信道中也能實(shí)現(xiàn)盲解碼，稍加修改還可以推廣到其他碼率為(n0-1)/n0的刪除卷積碼的盲解碼上，具有較高的實(shí)用價(jià)值。

參考文獻(xiàn)

[1]譚維熾，顧瑩琦.空間數(shù)據(jù)系統(tǒng)[M].北京：中國(guó)科學(xué)技術(shù)出版社，2004.

[2]眭惠巧.基于校驗(yàn)矩陣的卷積碼識(shí)別和碼字同步[J].無(wú)線電通信技術(shù)，2008，34(1)：26-28.

[3]韓國(guó)賓.刪除卷積碼的識(shí)別技術(shù)[D].成都：電子科技大學(xué)，2009.

[4]CCSDS 131.0-B-1 Blue Book.TM synchronization and channel coding[S].CCSDS Press，2003.

[5]CAIN J B.Punctured convolutional codes of rate (n-1)/n and simplified maximum likelihood decoding[J].IEEE Trans.on IT， 1979(1)：97-100.

[6]陳發(fā)新.刪除卷積碼生成矩陣及最簡(jiǎn)信息恢復(fù)式的求法[J].無(wú) 線 通 信 技 術(shù) ，2009(2)：5-10.

[7]王新梅，肖國(guó)鎮(zhèn).糾錯(cuò)碼—原理與方法[M].西安：西安電子科學(xué)技術(shù)大學(xué)出版社，2001.

[8]陳大新.矩陣?yán)碚揫M].上海：上海交通大學(xué)出版社，1997.