1、 沒實(shí)驗(yàn)一 繪制二進(jìn)熵函數(shù)曲線（2個(gè)學(xué)時(shí)）一、實(shí)驗(yàn)?zāi)康模?. 掌握Excel的數(shù)據(jù)填充、公式運(yùn)算和圖表制作2. 掌握Matlab繪圖函數(shù)3. 掌握、理解熵函數(shù)表達(dá)式及其性質(zhì)二、實(shí)驗(yàn)要求：1. 提前預(yù)習(xí)實(shí)驗(yàn)，認(rèn)真閱讀實(shí)驗(yàn)原理以及相應(yīng)的參考書。2. 在實(shí)驗(yàn)報(bào)告中給出二進(jìn)制熵函數(shù)曲線圖三、實(shí)驗(yàn)原理：1. Excel的圖表功能2. 信源熵的概念及性質(zhì)四、實(shí)驗(yàn)內(nèi)容：用Excel或Matlab軟件制作二進(jìn)熵函數(shù)曲線。具體步驟如下：1、啟動(dòng)Excel應(yīng)用程序。2、準(zhǔn)備一組數(shù)據(jù)p。在Excel的一個(gè)工作表的A列（或其它列）輸入一組p，取步長為0.01，從0至100產(chǎn)生101個(gè)p（利用Excel填充功能）。3、

2、取定對(duì)數(shù)底c，在B列計(jì)算H(x) ,注意對(duì)p=0與p=1兩處，在B列對(duì)應(yīng)位置直接輸入0。Excel中提供了三種對(duì)數(shù)函數(shù)LN(x),LOG10(x)和LOG(x,c)，其中LN(x)是求自然對(duì)數(shù)，LOG10(x)是求以10為底的對(duì)數(shù)，LOG(x,c)表示求對(duì)數(shù)。選用c=2,則應(yīng)用函數(shù)LOG(x,2)。在單元格B2中輸入公式：=-A2*LOG(A2,2)-(1-A2)*LOG(1-A2,2)雙擊B2的填充柄，即可完成H(p)的計(jì)算。4、使用Excel的圖表向?qū)В瑘D表類型選“XY散點(diǎn)圖”，子圖表類型選“無數(shù)據(jù)點(diǎn)平滑散點(diǎn)圖”，數(shù)據(jù)區(qū)域用計(jì)算出的H(p)數(shù)據(jù)所在列范圍，即$B$1:$B$101。在“系列

3、”中輸入X值(即p值)范圍，即$A$1:$A$101。在X軸輸入標(biāo)題概率，在Y軸輸入標(biāo)題信源熵。實(shí)驗(yàn)二：香農(nóng)編碼軟件實(shí)現(xiàn)（2個(gè)學(xué)時(shí)）1、實(shí)驗(yàn)?zāi)康模?）了解香農(nóng)編碼的基本原理及其特點(diǎn)；（2）熟悉掌握香農(nóng)編碼的方法和步驟；（3）掌握C語言或者M(jìn)atlab編寫香農(nóng)編碼的程序。2、實(shí)驗(yàn)報(bào)告要求（1）簡(jiǎn)要總結(jié)香農(nóng)編碼的基本原理與特點(diǎn)（2）寫出香農(nóng)編碼的基本步驟，畫出實(shí)現(xiàn)香農(nóng)編碼的程序流程圖（3）實(shí)現(xiàn)香農(nóng)編碼的Matlab或者C源程序3、實(shí)驗(yàn)內(nèi)容（1）根據(jù)香農(nóng)編碼的方法和步驟，用香農(nóng)編碼編寫程序（2）用編寫的源程序驗(yàn)證書中例題的正確性。實(shí)驗(yàn)三：Huffman編碼軟件實(shí)現(xiàn)（2個(gè)學(xué)時(shí)）1、實(shí)驗(yàn)?zāi)康模?）進(jìn)一步

4、熟悉Huffman編碼過程；（2）掌握C語言遞歸程序的設(shè)計(jì)和調(diào)試技術(shù)（或者使用Matlab）。2、實(shí)驗(yàn)要求（1）輸入：信源符號(hào)個(gè)數(shù)r、信源的概率分布P；（2）輸出：每個(gè)信源符號(hào)對(duì)應(yīng)的Huffman編碼的碼字。3、實(shí)驗(yàn)內(nèi)容（1）算法 1、從鍵盤輸入組成信源C的字符個(gè)數(shù)N； 2、從鍵盤輸入信源C和組成信源的字符所對(duì)應(yīng)的概率數(shù)組P； 3、用函數(shù)來對(duì)信源進(jìn)行二進(jìn)制編碼；先對(duì)P按從大到小進(jìn)行排序，與此同時(shí)要把C中相應(yīng)的字符的位置做相應(yīng)的調(diào)換；用數(shù)組來記錄編碼：在進(jìn)行記錄編碼時(shí)是從數(shù)組的最后一個(gè)開始存儲(chǔ)的，而且，每進(jìn)行一次編碼所記錄下來的兩個(gè)編碼是按從數(shù)組的最后一個(gè)元素開始服從countm-k-j、cou

5、ntm-k-j-1，其中k表示編碼所進(jìn)行的次數(shù)，j表示每次編碼都只有;最后用函數(shù)來輸出編碼。（2）部分偽代碼：（一）節(jié)點(diǎn)信息結(jié)構(gòu)體struct HuffNode int weight;/信源符號(hào)的概率 int parent;int lchild;int rchild;（二）算法void Huffman(int weight, int n, HuffNode hn, HuffCode hc) for(i = 0; i != 2*n - 1; +i) /create Huffman Node,step 1for(i = 0; i != n-1; +i) /create Huffman Node,

6、step 2for(j = 0; j != n+i; j+) if(hnj.weight < min1 && hnj.parent = 0)else if(hnj.weight < min2 && hnj.parent = 0)else ; 在此逆序存儲(chǔ)Huffman編碼 int tempmaxlen; for(i = 0; i != n; +i)int parent = hni.parent;while(hnchild.parent != 0) 4、實(shí)驗(yàn)報(bào)告（1）簡(jiǎn)要總結(jié)Huffman編碼的原理與特點(diǎn)（2）寫出Huffman編碼的基本步驟，畫出實(shí)現(xiàn)

7、Huffman編碼的程序流程圖（3）給出Huffman編碼的源程序，并給出實(shí)驗(yàn)過程中的測(cè)試結(jié)果（4）總結(jié)實(shí)驗(yàn)過程遇到的問題及解決方法實(shí)驗(yàn)四 循環(huán)碼的編碼和譯碼程序設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康模?1通過實(shí)驗(yàn)了解循環(huán)碼的工作原理。 2了解生成多項(xiàng)式g(x)與編碼、譯碼的關(guān)系。 3了解碼距d與糾、檢錯(cuò)能力之間的關(guān)系。4分析(7.3)循環(huán)碼的糾錯(cuò)能力。二、實(shí)驗(yàn)要求：1、編、譯碼用上述的計(jì)算法程序框圖編寫。2、計(jì)算出所有的碼字集合，可糾錯(cuò)誤圖樣E（x）表和對(duì)應(yīng)的錯(cuò)誤伴隨式表。3、考查和分析該碼檢、糾一、二位錯(cuò)誤的能力情況。4、整理好所有的程序清單，變量名盡量用程序框圖所給名稱，并作注釋。5、出示軟件報(bào)告.三、實(shí)驗(yàn)設(shè)

8、計(jì)原理 1、循環(huán)碼編碼原理設(shè)有一（，）循環(huán)碼，碼字Cn-1CrCr-1C0，其中r=n-k。碼字多項(xiàng)式為： C(x)= Cn-1xn-1+Cn-2xn-2+C1x+C0。 碼字的生成多項(xiàng)式為： g（x）=gr-1xr-1gr-2xr-2+g1x+g0 待編碼的信息多項(xiàng)式為： m（x）=mK-1xK-1+m0 xn-k.m（x）=Cn-1xn-1+Cn-Kxn-K 對(duì)于系統(tǒng)碼有： Cn-1=mK-1，Cn-2=mK-2，Cn-K=Cr=m0 設(shè)監(jiān)督多項(xiàng)式為： r（x）=Cr-1Xr-1+C1x+C0 根據(jù)循環(huán)碼的定義，則有： C（x）=xn-Km（x）+r（x）=q（x）.g(x) Xn-Km(

9、x)=q(x).g(x)+r(x) r(x)=Rg(x)xn-Km(x)即監(jiān)督多項(xiàng)式是將多項(xiàng)式xn-Km(x)除以g(x)所得的余式。編碼過程就是如何根據(jù)生成多項(xiàng)式完成除法運(yùn)算求取監(jiān)督多項(xiàng)式的過程。 設(shè)循環(huán)碼(7.3)碼的字多項(xiàng)式為: C(x)=C6x6+C5x5+C4x4+C3x3+C2x2C1x+C0 (n=7) 生成多項(xiàng)式為: g(x)=x4+x2+x+1 信息多項(xiàng)式為: m(x)=m2x2+m1x+m0 (k=3), 設(shè)m(x)=x2+x 監(jiān)督多項(xiàng)式為: r(x)= Cr-1Xr-1+C1x+C0根據(jù)循環(huán)碼的定義：生成多項(xiàng)式的倍式均是碼字，編碼實(shí)際上是做xn-km(x)除以g(x)的

10、運(yùn)算求得r(x)。編碼程序框圖見圖4.1(a)左,二進(jìn)制多項(xiàng)式除法示意圖見圖4.1(b)。 2、譯碼原理設(shè)R(x)為接收碼字多項(xiàng)式，E(x)為錯(cuò)誤圖樣多項(xiàng)式，S(x)為伴隨式，則根據(jù)循環(huán)碼的性質(zhì)有： S(x)=g(x)R(x)=g(x)E(x) 當(dāng)R(x)=C(x)時(shí)，有E(x)=0，S(x)=0 當(dāng)R(x)不等于C(x)時(shí)，有E(x)為非0，S(x)為非0 111 .商數(shù) g(x): 10111| 1100000 .xrm(x) + 10111 .第一步 11110 + 10111 .第二步 10010 + 10111 .第三步 101 .余式:x2+1 編碼步驟： 、n-k=r=7-3=4

11、,用x4乘m(x)的運(yùn)算實(shí) 際上相當(dāng)于在信息碼110后附上4個(gè) 0,變?yōu)?100000 、用xrm(x)=x4(x2+x)=x6+x5除以g(x), 如圖(a)所示,得到監(jiān)督余式r(x)=x2+1。 、編出相應(yīng)的發(fā)送碼字為： C(x)=xrm(x)+r(x)C=1100000+101=11001014、按上述步驟，將得到下述碼表：信息位 監(jiān)督位000 0000001 0111010 1110011 1001100 1011101 1100110 0101111 0010除法子程序NNYYC xnkm(x),DC,r=n-kGg(x)系數(shù)設(shè)循環(huán)變量B=KC的第B+r位=0?C+GCG右移一位BB

12、+1B=0?DC+D碼字D輸出 (a)編碼計(jì)算程序框圖 (b)二進(jìn)制多項(xiàng)式除法示意圖 圖4.1編碼計(jì)算程序框圖及多項(xiàng)式除法示意圖譯碼過程如下:計(jì)算每一種可能被糾的錯(cuò)誤圖樣E(x)的伴隨式， Si(x)=g(x)E(x) 本地作數(shù)據(jù)表存儲(chǔ)好。根據(jù)已接收碼字多項(xiàng)式R(x)，計(jì)算相應(yīng)的伴隨式： S(x)=g(x)R(x)將實(shí)際接收碼字求出的S(x)與本地存儲(chǔ)的各Si(x)相比較，查出兩者相等的那個(gè)唯一匹配的Si(x)，找出并得到相應(yīng)的錯(cuò)誤圖樣E(x)。 (4) 糾錯(cuò): C(x)=R(x)+E(x)否則由S(x)找不出唯一匹配的Si(x)，則報(bào)出錯(cuò)信息，表示出現(xiàn)不可糾錯(cuò)的錯(cuò)誤圖樣，即碼元出錯(cuò)的個(gè)數(shù)超出該循環(huán)碼的糾錯(cuò)能力。 譯碼流程圖4.2所示:YN初始化可糾錯(cuò)誤圖樣種類總數(shù)=N,i=0Gg(x)系數(shù)CEi(x)系數(shù),D0調(diào)除法子程序求Rg(x)Ei(x)SiCi=i+1i=N?將所有的EI(x)系數(shù)=Ei及對(duì)應(yīng)SI造表存儲(chǔ)待用YNNNNYY初始化Gg(x)系數(shù),i=0,N賦值輸入接收矢量R(x)的系數(shù)R,CR(x)系數(shù)調(diào)除法子程序求Rg(