線性分組碼編碼與譯碼

線性分組碼編碼與譯碼第1頁，共31頁，2023年，2月20日，星期二0000010100111001011101110000001001100101101100100011011011011111消息碼字碼許用碼字禁用碼字編碼效率漢明重量漢明距離最小漢明距離糾檢錯能力群子群分元陪集0001101100000010011001011011消息碼字基本概念第2頁，共31頁，2023年，2月20日，星期二0000001001100101101100100011011011011111000010100010011110100010101100101111111000010010111000011001001100111110100111010001101010100011100000111011101010111100分元陪集劃分第3頁，共31頁，2023年，2月20日，星期二0000010100111001011101110000001001100101101100100011011011011111消息碼字碼許用碼字禁用碼字編碼效率漢明重量漢明距離最小漢明距離糾檢錯能力群子群分元陪集域GF(2)上的矢量空間子空間矢量張成的子空間基底維數(shù)零化空間矩陣行空間0001101100000010011001011011消息碼字GF(2)基本概念第4頁，共31頁，2023年，2月20日，星期二線性分組碼長度為n，有2k個碼字的碼組，當且僅當這2k個碼字是GF(2)上n維矢量空間的一個k維子空間時，稱為(n,k)線性分組碼，簡稱(n,k)碼。由于k維子空間是在模2加法下運算的，構(gòu)成了一個加法交換群，所以線性分組碼也稱為群碼。碼率R=k/n，就是傳輸效率。最小漢明距離dmin等于非零碼字的最小重量。系統(tǒng)碼n-k

kk

n-k碼字信息位與輸入信息序列相同，并且與校驗位分開第5頁，共31頁，2023年，2月20日，星期二生成矩陣線性分組碼是GF(2)上n維空間中的一個k維子空間，因此它可以由k個線性無關(guān)n維矢量完全確定。由于這組矢量的所有線性組合給出了碼C中的任一個碼字，稱生成碼C。C中任何一組基底所構(gòu)成的矩陣G稱作碼C的生成矩陣第6頁，共31頁，2023年，2月20日，星期二生成矩陣對于任何一個給定的信息序列，可指定作為相應(yīng)的碼字。G矩陣的每一行都是一個碼字矢量，分別對應(yīng)信息位為(10…0),(010…0)…(00…01)時的碼字。第7頁，共31頁，2023年，2月20日，星期二生成矩陣(n,k)分組碼實際上就是這k個線性無關(guān)的碼字矢量張成的k維子空間，這k個矢量組成的矩陣就是生成矩陣。確定(n,k)碼的生成矩陣的問題實際上就是確定n重矢量空間中k維子空間的k個線性無關(guān)的碼字矢量的問題，也就是尋找基底的問題。(n,k)碼的n重矢量空間中可以有多個k維子空間，產(chǎn)生不同的碼組，即有不同的基底。(n,k)碼的n-重矢量空間中的一個k維子空間的基底可以有多個，因此可以有不同的生成矩陣G，但都產(chǎn)生相同的碼組。第8頁，共31頁，2023年，2月20日，星期二基底的線性組合等效于G的行初等變換，可以產(chǎn)生一組新的基底。利用這一點，可使生成矩陣具有如下“系統(tǒng)形式”，稱之為典型生成矩陣。典型生成矩陣即：G=[IkQ]，Q為k×r矩陣，Ik為k×k單位陣。非系統(tǒng)碼與系統(tǒng)碼并無本質(zhì)區(qū)別，它的生成矩陣可以通過行初等變換轉(zhuǎn)變?yōu)橄到y(tǒng)形式，這個過程叫做系統(tǒng)化。系統(tǒng)化并不會改變碼集，其糾錯能力完全等價。第9頁，共31頁，2023年，2月20日，星期二

例題設(shè)二元（5,3）碼，其生成矩陣為將其化為標準形式？第10頁，共31頁，2023年，2月20日，星期二一致校驗矩陣與任何一個(n,k)碼的碼空間C相對應(yīng)，一定存在一個對偶空間D，它的每個矢量都與C中的每個矢量正交，其維數(shù)為n-k。事實上，若找出生成空間D的基底（n-k個）用這n-k個矢量同樣可以生成包含個碼字的(n,n-k)線性分組碼，我們稱其(n,k)碼的對偶碼，生成矩陣為第11頁，共31頁，2023年，2月20日，星期二一致校驗矩陣由對偶空間的定義知，有對任意的即H可以檢驗一個n重是否為碼字，稱H為碼C的一致校驗矩陣。第12頁，共31頁，2023年，2月20日，星期二典型一致校驗矩陣系統(tǒng)碼的一致校驗矩陣為即H=[PIr],

其中，Ir為r×r單位陣，P為r×k矩陣。第13頁，共31頁，2023年，2月20日，星期二一致校驗矩陣與生成矩陣之間的關(guān)系

由于生成矩陣每一行都是一個碼字，因此應(yīng)當滿足一致校驗矩陣所規(guī)定的校驗關(guān)系，即應(yīng)當有：GHT=0

或者HGT=0

因此H與G互為正交矩陣，說明G和H的行空間互為零化空間。第14頁，共31頁，2023年，2月20日，星期二一致校驗矩陣與生成矩陣之間的關(guān)系對于系統(tǒng)碼，上式可以寫成[IkQ][PIr]T=0得[PT+Q]=0所以PT=Q或者P=QT即P矩陣與Q矩陣互為轉(zhuǎn)置矩陣。對于系統(tǒng)碼，已知校驗矩陣H就可以確定典型生成矩陣G，反之，已知生成矩陣也就可以確定校驗矩陣。第15頁，共31頁，2023年，2月20日，星期二例題【例】設(shè)二元(7,4)碼的生成矩陣為求其一致校驗矩陣？【例】設(shè)二元(5,3)碼，其生成矩陣為求其一致校驗矩陣？第16頁，共31頁，2023年，2月20日，星期二線性分組碼編碼線性分組碼的編碼過程分為兩步：把信息序列按一定長度分成若干信息碼組，每組由k位組成；編碼器按照預定的線性規(guī)則，把信息碼組變換成n重(n>k)碼字。信息碼組長k位，有2k個不同的信息碼組，則有2k個碼字與它們一一對應(yīng)。第17頁，共31頁，2023年，2月20日，星期二一致校驗矩陣編碼設(shè)c=[c0,c1,c2,c3,c4,c5,c6]，其中，[c0,c1,c2,c3]為信息位，[c4,c5,c6]為校驗位。由HCT=0可知校驗方程為：c4=c0+c2+c3c5=c0+c1+c2c6=c1+c2+c3信息碼元m=[1101]則編得的碼字c=[1101000]第18頁，共31頁，2023年，2月20日，星期二生成矩陣編碼若信息碼元m=[1101]，則有c=mG=[1101000]。第19頁，共31頁，2023年，2月20日，星期二譯碼準則設(shè)發(fā)送碼字為c=(c0,c1,……,cn-1)，由于信道干擾產(chǎn)生差錯，反映到接收碼字上可以用一個二元矢量e表示，e=(e0,e1,……,en-1)，稱為錯誤圖樣，其中，ei=1表明相應(yīng)位有錯，ei=0表明相應(yīng)位無錯。這時接收碼字可以表示為r=c+e=(c0+e0，c1+e1,……cn-1+en-1)譯碼器就是從接收碼字r得到發(fā)送碼字的估計值，或者說從接收碼字中確定錯誤圖樣e，然后由c=r-e得到發(fā)送碼字的估計值。如果估計正確則譯碼正確，否則譯碼錯誤。如何得到發(fā)送碼字的估計值，根據(jù)什么準則？第20頁，共31頁，2023年，2月20日，星期二譯碼準則最大后驗概率譯碼最大似然譯碼最小距離譯碼對于給定的接收矢量，計算它與M個可能的發(fā)送碼字之間的距離，從中選擇能使距離達到最小的碼字作為判決結(jié)果。第21頁，共31頁，2023年，2月20日，星期二若信道是對稱DMC信道，其轉(zhuǎn)移概率為1-p和p/(q-1)，則

則對數(shù)似然函數(shù)為最大似然譯碼準則可簡化為：若對所有的，有則判定最小距離譯碼最小距離譯碼第22頁，共31頁，2023年，2月20日，星期二伴隨式設(shè)碼C的一致校驗矩陣為H，接收矢量為r，定義稱s為接收矢量r的伴隨式。由伴隨式的定義可知s=rHT=(c+e)HT=cHT+eHT=eHT可以看出伴隨式完全由e決定，它充分反映了信道的干擾情況。如果伴隨式s≠0，接收碼字一定有錯誤；如果伴隨式s=0，譯碼器認為接收碼字無錯誤。第23頁，共31頁，2023年，2月20日，星期二譯碼步驟由接收碼字r計算伴隨式sT=HrT若s=0，則譯碼器認為接收碼字沒錯，否則有錯，并由s計算錯誤圖樣e由錯誤圖樣進行譯碼，估計發(fā)送的碼字c=r-e=r+e

其中最困難的是確定錯誤圖樣，即錯誤定位。如何進行錯誤定位？第24頁，共31頁，2023年，2月20日，星期二譯碼思路1根據(jù)sT=HrT=HeT列出線性方程組（含有n-k個相互獨立的方程），通過求解線性方程得到e。上式有n個未知數(shù)e0,e1,…,en-1，有n-k個方程，因此有多解。（伴隨式s是一個r=n-k維矢量，共有個，而錯誤圖樣e是n維矢量，共有個，因此，s與e不存在一一對應(yīng)關(guān)系。）最終根據(jù)譯碼準則選取其中一個，常常選取重量最輕的為錯誤圖樣e的估計值，從而得到發(fā)送碼字的估計值，體現(xiàn)最小距離譯碼的思想。第25頁，共31頁，2023年，2月20日，星期二譯碼思路2(n,k)分組碼的2k個碼字，是n維矢量空間Vn中的一個k維子空間，它在GF(2)上是一個子群。利用分元陪集的方法，可以利用該子空間的2k元素，生成Vn中的所有2n個元素。陪集首c0c1c2…c2k-1e1c1+e1c2+e1c2k-1+e1e2c1+e2c2+e2c2k-1+e2……………e2r-1c1+e2r-1c2+e2r-1c2k-1+e2r-1第26頁，共31頁，2023年，2月20日，星期二陪集首c0c1c2…c2k-1e1c1+e1c2+e1c2k-1+e1e2c1+e2c2+e2c2k-1+e2……………e2r-1c1+e2r-1c2+e2r-1c2k-1+e2r-1標準陣列譯碼表譯碼按以下規(guī)則進行：收到碼字R必然在這個表中，如果落在表中某一列，譯碼器就譯成第一行的相應(yīng)碼字。第27頁，共31頁，2023年，2月20日，星期二非常直觀的譯碼方法，選擇重量最輕的禁用碼字作為陪集首，體現(xiàn)最小距離譯碼思想。同一陪集中元素的伴隨式都相同，并且，陪集首與伴隨式矢量有一一對應(yīng)的關(guān)系。根據(jù)這種關(guān)系，可以將譯碼表簡化。標準陣列譯碼第28頁，共31頁，2023年，2月20日，星期二陪集首c0c1c2…c2k-1e1c1+e1c2+e1c2k-1+e1e2c1+e2c2+e2c2k-1+e2……………e2r-1c1+e2r-1c2+e2r-1c2k-1+e2r-1標準陣列譯碼表伴隨式s0s1s2…s2r-1第29頁，共31頁，2023年，2月20日，星期二非常直觀的譯碼方法，選擇重量最輕的禁用碼字作為陪集首，體現(xiàn)最小距離譯碼思想。同一陪集中元素的伴隨式都相同，并且，陪集首與伴隨式矢量有一一對應(yīng)的關(guān)系。根據(jù)這種關(guān)系，可以將譯碼表簡化。當n-k=r較小時(小于