漢明碼編譯碼

1、漢明碼編譯碼一設計思想漢明碼是一種常用的糾錯碼，具有糾一位錯誤的能力。本實驗使用M at lab平臺，分別 用程序語言和simulink來實現(xiàn)漢明碼的編譯碼。用程序語言實現(xiàn)就是從原理層面，通過產(chǎn) 生生成矩陣，錯誤圖樣，伴隨式等一步步進行編譯碼。用simul ink實現(xiàn)是用封裝好的漢明 碼編譯碼模塊進行實例仿真，從而驗證程序語言中的編譯碼和誤碼性能分析結果。此外，在 結合之祈信源編碼的基礎上，還可實現(xiàn)完整通信系統(tǒng)的搭建。二實現(xiàn)流程1. 漢明碼編譯碼圖1漢明碼編譯碼框圖1）根據(jù)生成多項式，產(chǎn)生指定的生成矩陣G2）產(chǎn)生隨機的信息序列3）由C = MG得到碼字4）進入信道傳輸5）計算S=RHr得到伴隨

2、式6）得到解碼碼流7）得到解碼信息序列2. 漢明碼誤碼性能分析誤碼率（SER）是指傳輸前后錯誤比特數(shù)占全部比特數(shù)的比值。 誤幀率（FER）是指傳輸前后錯誤碼字數(shù)占全部碼字數(shù)的比值。 通過按位比較、按幀比較可以實現(xiàn)誤碼率和誤幀率的統(tǒng)計。3. 構建完整通信系統(tǒng)圖2完整通信系統(tǒng)框圖三結論分析1.漢明碼編譯碼編寫了 GUI界面方便呈現(xiàn)過程和結果。圖3漢明碼編譯碼演示GUI界面 以產(chǎn)生（7, 4）漢明碼為例說明過程的具體實現(xiàn)。1）根據(jù)生成多項式.產(chǎn)生指定的生成矩陣G用H, G, n, k = hammgen（3, ,D"3+D+1,）函數(shù)得到系統(tǒng)碼形式的校驗矩陣H、G以及碼字長度n和信息位數(shù)k

3、1H= 00110 100 110 1G =1110 010 10 00 00 01 00 12）產(chǎn)生隨機的信息序列M0 0 10M= 0 10 00 1113）由C = MG得到碼字01 00011c =01 1010000 101114）進入信道傳輸假設是BSC信道.錯誤轉移概率設定為 傳輸后接收端得到的碼流為_0 0 0 0 0 1 rR= 11101000 0 10 111.紅邑表示錯誤比特。5）計算S=RHr得到伴隨式'o 1 r5= 1000 0 1錯誤圖樣0000001000001000001000001000001000001000001000000伴隨式1011110

4、11110001010100查表可知第一行碼字錯誤圖樣為0100000,第二行碼字錯誤圖樣為1000000, 第三行碼字錯誤圖樣為0000001。進行即可得到糾錯解碼的碼字C2o6）得到解碼碼流0 1 10 10 0'C2 =00000001110 0 10刀得到解碼信息序列_0 1 0 0_M2= 00000 0 10可以看出解碼信息序列與原信息序列一樣，體現(xiàn)了漢明碼的糾錯能力。2 性能分析1) BSC信道仿真設置BSC錯誤轉移槪率Pe從0到1變化，步進為.在每個Pe值進行1000次蒙 特卡洛仿真，得到圖4所示誤碼率隨Pe變化曲線圖和圖5所示誤幀率隨Pe變化曲 線圖。圖中綠線為BSC

5、信道誤碼率，紅線為設定Pe值，藍線為Hamming碼解碼誤碼率。 由圖線可以看出仿真的BSC信道誤碼率與Pe 一致。在PeC時，Hanging碼的解碼誤碼率隨著BSC信道錯誤傳輸概率Pe的減小而減 小o Hamming碼的解碼誤碼率顯著下降，約為Pe的1/2。Hamming碼的糾1位錯是 到了很好的效果。Pe吋，Hamming碼的解碼誤碼率大于Pe。這是因為在Pe時，傳一個碼字錯 誤比特數(shù)近似為2,而Hamming碼只能糾一位錯，兩位同時出錯時會糾成另一個碼 字，這樣就可能增加誤比特數(shù)，使得“越糾越錯”。Pe時，請況恰好相反。0.8FER with BSCFER with HammingTb

6、WcikjpocoTb kVcrtjpacnl7oWcr hpKQ2To WcOcspacoJBSC信道仿真框圖0.6QTULI0.40.20IIIIIIIII00.10.20.30.40.50.60.70.80.91Pe圖7誤幀率隨Pe變化曲線圖可以看出隨著Pe增加.BSC傳輸誤幀率和Hamming譯碼誤幀率成S曲線上升達 到1o Hamming譯碼誤幀率要低于BSC傳輸誤幀率，體現(xiàn)了其糾錯能力使得碼字錯誤 減少這一效果。與誤碼率的圖對比可以發(fā)現(xiàn)，誤幀率要比誤比特率爲。為了進一步驗證結果的正確性，進行了 simulink仿真。圖8用伯勞利二進制發(fā)生器產(chǎn)生隨機序列，進行漢明碼編碼，進入BSC信

7、道傳輸, 之后進行漢明碼譯碼，用Error Rate Calculation模塊統(tǒng)計誤碼率，結果如下：oroMotz山ErrorRate with SNRPe圖9 simul ink仿真BER隨Pe變化曲線圖與程序實現(xiàn)仿真的結果幾乎一樣。2) AWGN信道仿真AWGN信道仿真直接用simul ink實現(xiàn)。HjTMMvqCiYOdafRaI： 0rr*<wwrT i IToWoij«<x1圖10 AWGN信道仿真框圖設置系統(tǒng)的數(shù)字調制方式為2FSK,設定AWGN信道的SNR從0到8dB以1dB步 進變化，得到誤碼率統(tǒng)計圖。ErrorRate with SNR10°

8、a)e(roL山Hamming code error rate FSK error rate1°“。 12345678SNRindB圖11 simul ink仿真BER隨SNR變化曲線圖圖中綠線為2FSK調制誤碼率，是由于AWGN帶來的。藍線為漢明碼解碼后誤碼 率?？梢钥闯?，漢明碼能夠很好的降低誤碼率。在SNR達到5dB時錯誤槪率降低為.3完整通信系統(tǒng)的構建以傳輸圖片為例，信道設置為BSC信道。在不加入漢明碼和加入漢明碼兩種情 況下觀察傳輸后圖像的情況。結果如下表所示。BSC錯誤轉 移概率Pe0BSC傳輸圖像'3沁二尸- 信信加漢明碼誤 碼率00加漢明碼傳 輸圖像4 上

9、9;N *MT宿信信由結果可以看出，加入信道編碼后，當BSC錯誤轉移槪率Pe后，圖像恢復性 能有明顯的改罷。這體現(xiàn)了漢明碼雖然只有糾一位錯的能力，但由于一般信道的Pe 不會很大，其糾錯的實用性和效果還是很好的。四思考題解答1 采用循環(huán)Hamm in g碼在硬件實現(xiàn)中的優(yōu)點與普通的線性分組碼譯碼電路相比，循環(huán)漢明碼不需要存儲伴隨式及錯誤圖樣，顯 著的節(jié)省了寄存器的使用，起到簡化電路的作用。碼如何改進可提高糾檢錯性能可以在H校驗矩陣基越上進行擴展，最后一行為全1行，最后一列矢量為00-1 這樣任何3列是線性無關的，出用4,進行奇偶校驗，糾錯能力為仁 檢錯能力為2。即001 1首先在（7,4）漢明碼

10、的基礎上進行擴展,'1110100o'H=011101001101001011111111下面通過實例的方式說明擴展H校驗矩陣的檢錯性能。簡化譯碼表如下:得到（8, 4）擴展?jié)h明碼的生成矩陣H。錯誤圖樣0000000100000010000001000000100000010000001000000100000010000000伴隨式00010011010110010111110111111011錯碼個數(shù)原序列錯碼序列伴隨式查表結果10000 00000000 01010100有21110 10001110 10100011無30001 01110010 01011001有40

11、100 11101100 00100111有說明對于錯碼個數(shù)為1的，既可以檢錯也可以糾錯； 錯碼個數(shù)為2的，可以檢錯，但不能糾錯；錯碼個數(shù)大于2的，被認為是錯碼個數(shù)為1,糾成其他碼字。附錄clear al IH,G,n,k = harm)gen(3, D+D+D;%H, G, n, k = harnngen (3, 1 D"4+D"2+D+1');%H, G, n, k = hammgen (4, 1 D'4+D+11); %H, G, n, k = hawngen(5J5+2+)；桃產(chǎn)生校驗矩陣E= zeros (1, n) ;f I ipl r (eye

12、 (n, n);%產(chǎn)生錯誤圖樣一共是n+1個S=mod (H*E', 2) ； %生成錯誤圖樣的伴隨式%暫產(chǎn)生消息序列%二進制隨機矩陣%M= randi(0,1,1,k);%產(chǎn)生 4 位消息列nm=3;M= randi (0,1,nm, k) ;%產(chǎn)生 4 位消息列%消息序列code=mod (M*G, 2) ;%對消息序列編碼EBSC信道進入Pe=;for j=1:nmfor i=1:ncode_bsc(j, i)=mod(code(j, i) + (unidrnd (round (1/Pe) =1 ),2);%模2加得到傳輸后的 綸碼delta(j, i)=code_bsc(j,

13、i)-code(j, i);%作差來計算錯誤位置end endep=find(delta"=0);%error position display (Iength (ep), ' BSC 錯誤位數(shù)') d i sp I ay (I ength (ep) / (nm*n), 1 BSC 誤比特率');Scode=mod (code_bsc*H1,2 )' %Scode=1 10' errow2=0;for i=1:nmif sum(code_bsc (i, :)-code (i, :)=0 errow2=errow2+1;endenddisplay

14、 (er row2, 'BSC 錯誤碼字數(shù)');% display (er r ow2/nm, * BSC 誤碼率J;for m=1:nmfor i=1:n+1if S(:, i)=Scode (:, m)j=i；endend %找到對應的伴隨式的位置dcode (m, :) =mod (code_bsc(m, :)+E(j, :), 2);ender= I ength (f i nd (dcode-code"=0) ;%計算誤比特的個數(shù) enta=er/(nm*n);display (code,'信息序列碼字)display(code.bsc, 'BSC傳輸后的信息序列碼字')di splay (dcode,'解碼后的信息序列')% display (errow2/nm,'解碼后誤碼率');%