信息论与编码课程设计哈夫曼编码的分析与实现

1、教师评语:吉林建筑大学电气与电子信息工程学院信息理论与编码课程设计报告设计题目：哈夫曼编码的分析与实现专业班级：电子信息工程101学生姓名：学号：指导教师：吕卅王超设计时间：2013.11.182013.11.29成绩评阅教师日期一、设计的作用、目的信息论与编码是一门理论与实践密切结合的课程，课程设计是其实践性教学环节之一，同时也是对课堂所学理论知识的巩固和补充。其主要目的是加深对理论知识的理解，掌握查阅有关资料的技能，提高实践技能，培养独立分析问题、解决问题及实际应用的能力。通过完成具体编码算法的程序设计和调试工作，提高编程能力，深刻理解信源编码、信道编译码的基本思想和目的，掌握编码的基本原

2、理与编码过程，增强逻辑思维能力，培养和提高自学能力以及综合运用所学理论知识去分析解决实际问题的能力，逐步熟悉开展科学实践的程序和方法二、设计任务及要求通过课程设计各环节的实践，应使学生达到如下要求：1,理解无失真信源编码的理论基础，掌握无失真信源编码的基本方法；2 .掌握哈夫曼编码/费诺编码方法的基本步骤及优缺点；3 .深刻理解信道编码的基本思想与目的，理解线性分组码的基本原理与编码过程；4,能够使用MATLAB或其他语言进行编程，编写的函数要有通用性。三、设计内容一个有8个符号的信源X,各个符号出现的概率为：,X1;产，X2,x3,x4,x5,xx7x8.1P(X)_、0

3、.10.070.060.050.04J编码方法：先将信源符号按其出现的概率大小依次排列，并取概率最小的字母分别配以0和1两个码元(先0后1或者先1后0,以后赋值固定)，再将这两个概率相加作为一个新字母的概率，与未分配的二进制符号的字母重新排队。并不断重复这一过程，直到最后两个符号配以0和1为止。最后从最后一级开始，向前返回得到各个信源符号所对应的码元序列，即为对应的码字。哈夫曼编码方式得到的码并非唯一的。在对信源缩减时，两个概率最小的符号合并后的概率与其他信源符号的概率相同时，这两者在缩减中的排序将会导致不同码字，但不同的排序将会影响码字的长度，一般讲合并的概率放在上面，这样可获得较小的码方差

4、。四、设计原理4.1 哈夫曼编码步骤(1)将信源消息符号按照其出现的概率大小依次排列为pl_p2_pn(2)取两个概率最小的字母分别配以0和1两个码元，并将这两个概率相加作为一个新的概率，与未分配的二进制符号的字母重新排队。(3)对重新排列后的两个最小符号重复步骤(2)的过程。(4)不断重复上述过程，知道最后两个符号配以0和1为止。(5)从最后一级开始，向前返回得到的各个信源符号所对应的码元序列，即为相应的码字。4.2 哈夫曼编码特点哈夫曼编码是用概率匹配的方法进行信源匹配方法进行信源。它的特点是：(1)哈夫曼的编码方法保证了概率大的符号对应于短码，概率小的符号对应于长码，充分应用了短码。(2

5、)缩减信源的最后两个码字总是最后一位不同，从而保证了哈夫曼编码是即时码。(3)哈夫曼编码所形成的码字不是唯一的，但编码效率是唯一的，在对最小的两个速率符号赋值时可以规定大的为“1”，小得为“0”，如果两个符号的出现概率相等时，排列时无论哪个在前都可以，所以哈夫曼所构造的码字不是唯一的，对于同一个信息源，无论上述的顺序如何排列，他的平均码长是不会改变的，所以编码效率是唯一的。(4)只有当信息源各符号出现的概率很不平均的时候，哈夫曼编码的效果才明显。(5)哈夫曼编码必须精确的统计出原始文件中每个符号出现频率，如果没有这些精确的统计将达不到预期效果。哈夫曼编码通常要经过两遍操作，第一遍进行统计，第二

6、遍产生编码，所以编码速度相对慢。另外实现电路复杂，各种长度的编码的编译过程也是比较复杂的，因此解压缩的过程也比较慢。(6)哈夫曼编码只能用整数来表示单个符号，而不能用小数，这很大程度上限制了压缩效果。哈夫曼所有位都是合在一起的，如果改动其中一位就可以使其数据变得面目全非。五、设计步骤5.1 以框图形式画出哈夫曼编码过程(哈夫曼编码要求构建哈夫曼树)表1哈夫曼编码信源编码概率编码过程码字码长X0.601.090.230.37赤.4/11X3(0.1870.1900.2310013X30.10.1卜10.13/0.1

7、810113X40.10.0700.091,0.06100004X50.0701004X60.0601014X70.05J000105X80.041000115哈夫曼树:给定n个实数w1,w2,wn(n2),求一个具有n个结点的二叉数，使其带权路径长度最小。所谓树的带权路径长度，就是树中所有的叶结点的权值乘上其到根结点的路径长度(若根结点为0层，叶结点到根结点的路径长度为叶结点的层数)。树的带权路径长度为WPL=(W1*L1+W2*L2+W3*L3+.+Wn*Ln),N个权值Wi(i=1,2,.n)构成一棵有N个叶结点的二叉树，相应的叶结点的路径长度为Li(i=1,2,

8、.n)。可以证明哈夫曼树的WPL是最小的。(1)根据与n个权值w1,w2川口对应的n个结点构成具有n棵二叉树的森林F=T1,T2北,其中第i棵二叉树Ti(1in)都只有一个权值为wi的根结点，其左、右子树均为空。(2)在森林F中选出两棵根结点的权值最小的树作为一棵新树的左、右子树，且置新树的根结点的权值为其左、右子树上根结点权值之和。(3)从F中删除构成新树的那两棵，同时把新树加入F中。(4)重复第(2)和第(3)步，直到F中只含有一棵为止，此树便为Huffman图1哈夫曼树5.2 计算平均码长、编码效率、冗余度。平均码长为：_8K=Zp(xi)Ki=0.4X1+0.18X3+0.1X3+0.

9、1X4+0.07X4+0.06X4+i10.05X5+0.04X5=2.61(码元/符号)信源嫡为：nH(X)=、p(xi)logp(xi)=-(0.4log0.4+0.1810g0.18+0.1log0.1+i10.1log0.1+0.07+log0.07+0.0610g0.06+0.0510g0.05+0.0410g0.04)=2.55bit/符号信息传输速率为：r=HXl=2.55=0.977Et/码元K2.61编码效率为：H(X)2.55-T=-=0.977K2.61冗余度为：=1-=1-0.977=0.023六、哈夫曼编码的实现6.1 软件介绍VisualC+6.0,简称VC或者VC

10、6.0,是微软推出的一款C+编译器，将高级语言”翻译为机器语言(低级语言)”的程序。VisualC+是一个功能强大的可视化软件开发工具。自1993年Microsoft公司推出VisualC+1.0后，随着其新版本的不断问世，VisualC+已成为专业程序员进行软件开发的首选工具。VisualC+6.0由Microsoft开发，它不仅是一个C+编译器，而且是一个基于Windows操作系统的可视化集成开发环境(integrateddevelopmentenvironment,IDE)。VisualC+6.0由许多组件组成，包括编辑器、调试器以及程序向导AppWizard、类向导ClassWizar

11、d等开发工具。这些组件通过一个名为DeveloperStudio的组件集成为和谐的开发环境。Microsoft的主力软件产品。VisualC+是一个功能强大的可视化软件开发工具。VisualC+6.0以拥有语法高亮”，自动编译功能以及高级除错功能而著称。比如，它允许用户进行远程调试，单步执行等。还有允许用户在调试期间重新编译被修改的代码，而不必重新启动正在调试的程序。其编译及创建预编译头文件(stdafx.h)、最小重建功能及累加连结(link)著称。这些特征明显缩短程序编辑、编译及连结的时间花费，在大型软件计划上尤其显著。(1) DeveloperStudio这是一个集成开发环境，我们日常工

12、作的99%都是在它上面完成的，再加上它的标题赫然写着“MicrosoftVisualC+,所以很多人理所当然的认为，那就是VisualC+了。其实不然，虽然DeveloperStudio提供了一个很好的编辑器和很多Wizard,但实际上它没有任何编译和链接程序的功能，真正完成这些工作的幕后英雄后面会介绍。我们也知道，DeveloperStudio并不是专门用于VC的，它也同样用于VB,VJ,VID等VisualStudio家族的其他同胞兄弟。所以不要把DeveloperStudio当成VisualC+,它充其量只是VisualC+的一个壳子而已。这一点请切记！(2) MFC从理论上来讲，MFC

13、也不是专用于VisualC+,BorlandC+,C+Builder和SymantecC+同样可以处理MFC。同时，用VisualC+编写代码也并不意味着一定要用MFC,只要愿意，用VisualC+来编写SDK程序，或者使用STL,ATL,一样没有限制。不过，VisualC+本来就是为MFC打造的，VisualC+中的许多特征和语言扩展也是为MFC而设计的，所以用VisualC+而不用MFC就等于抛弃了VisualC+中很大的一部分功能。但是，VisualC+也不等于MFC。(3) PlatformSDK这才是VisualC+和整个VisualStudio的精华和灵魂，虽然我们很少能直接接触到

14、它。大致说来，PlatformSDK是以MicrosoftC/C+编译器为核心(不是VisualC+,看清楚了)，配合MASM,辅以其他一些工具和文档资料。上面说到DeveloperStudio没有编译程序的功能，那么这项工作是由谁来完成的呢？是CL,是NMAKE,和其他许许多多命令行程序，这些我们看不到的程序才是构成VisualStudio的基石。6.2 编程/*哈夫曼编码*#include#include#include#include#include#includeusingnamespacestd;structHuffman_InformationSourcecharInformati

15、onSign10；doubleProbability;charCode10;intCodeLength;structHuffNodecharInformationSign10；doubleProbabilityHuffNode*LeftSubtree,*middleSubtree,*RightSubtree,*Next;charCode10intCodeLengthclassCHuffman_2public:CHuffman_2()ISNumber=0AvageCodeLength=0.0InformationRate=0.0;CodeEfficiency=0.0Redundancy=0.0C

16、Huffman_2()DestroyBTree(HuffTree);voidHuffman_Input();voidHuffman_Sort();voidHuffman_Tree();voidHuffman_Coding()；voidHuffman_CodeAnalyzing()；voidHuffman_Display();voidDestroyBTree(HuffNode*TreePointer)；private:vectorISarray;intISNumber；doubleAvageCodeLengthdoubleInformationRate；doubleCodeEfficiencyH

17、uffNode*HuffTree;private：voidHuffman_Code(HuffNode*TreePointer);；/输入信源信息voidCHuffman_2:Huffman_Input()(Huffman_InformationSourcetemp1=x1,0.40,”,0;ISarray.push_back(temp1);Huffman_InformationSourcetemp2=x2,0.18,”,0;ISarray.push_back(temp2);Huffman_InformationSourcetemp3=x3,0.10,”,0;ISarray.push_back(

18、temp3);Huffman_InformationSourcetemp4=x4,0.10,”,0;ISarray.push_back(temp4);Huffman_InformationSourcetemp5=x5,0.07,”,0;ISarray.push_back(temp5);Huffman_InformationSourcetemp6=x6,0.06,0;ISarray.push_back(temp6);Huffman_InformationSourcetemp7=x7,0.05,0;ISarray.push_back(temp7);Huffman_InformationSource

19、temp8=x8,0.04,”,0;ISarray.push_back(temp8);ISNumber=ISarray.size();/按概率从大到小”排序voidCHuffman_2二Huffman_Sort()Huffman_InformationSourcetemp;intI,j;for(i=0;iISNumber-1;i+)for(j=i+1;jISNumber;j+)if(ISarrayi.ProbabilityInformationSign,ISarray0.InformationSign);ptr1-Probability=ISarray0.Probability;strcpy(

20、ptr1-Code,ISarray0.Code);ptr1-LeftSubtree=NULL;ptr1-middleSubtree=NULL;ptr1-RightSubtree=NULL;ptr1-Next=NULL;HuffTree=ptr1;for(i=1;iInformationSign,ISarrayi.InformationSign);ptr2-Probability=ISarrayi.Probability;strcpy(ptr2-Code,ISarrayi.Code);ptr2-LeftSubtree=NULL;ptr2-middleSubtree=NULL;ptr2-Right

21、Subtree=NULL;ptr2-Next=ptr1;ptr1=ptr2;HuffTree=ptr1;intk;ints;k=ceil(double)(ISNumber-3)/(3-1);s=3+k*(3-1)-ISNumber;if(s=1)ptr2=ptr1-Next;ptr4=newHuffNode;strcpy(ptr4-InformationSign,*);ptr4-Probability=ptr1-Probability+ptr2-Probability;strcpy(ptr4-Code,);ptr4-LeftSubtree=NULL;ptr4-middleSubtree=ptr

22、1;ptr4-RightSubtree=ptr2;HuffTree=ptr2-Next;temp1=HuffTree;while(temp1&(ptr4-Probabilitytemp1-Probability)temp2=temp1;temp1=temp1-Next;ptr4-Next=temp1;if(temp1=HuffTree)HuffTree=ptr4;elsetemp2-Next=ptr4;ptr1=HuffTree;while(ptr1-Next)/合并概率最小的结点ptr2=ptr1-Next;ptr3=ptr2-Next;ptr4=newHuffNode;strcpy(ptr

23、4-InformationSign,*);ptr4-Probability=ptr1-Probability+ptr2-Probability+ptr3-Probability;strcpy(ptr4-Code,);ptr4-LeftSubtree=ptr1;ptr4-middleSubtree=ptr2;ptr4-RightSubtree=ptr3;HuffTree=ptr3-Next;temp1=HuffTree;while(temp1&(ptr4-Probabilitytemp1-Probability)(temp2=temp1;temp1=temp1-Next;ptr4-Next=te

24、mp1;if(temp1=HuffTree)HuffTree=ptr4;elsetemp2-Next=ptr4;ptr1=HuffTree;/释放：ptr1=NULL;ptr2=NULL;ptr3=NULL;ptr4=NULL;temp1=NULL;temp2=NULL;strcpy(HuffTree-Code,);HuffTree-CodeLength=0;/生成哈夫曼码voidCHuffman_2二Huffman_Code(HuffNode*TreePointer)(if(TreePointer=NULL)return;chartempstr10=;if(!TreePointer-Left

25、Subtree&!TreePointer-middleSubtree&!TreePointer-RightSubtree)for(inti=0;iInformationSign)=0)(strcpy(ISarrayi.Code,TreePointer-Code);ISarrayi.CodeLength=TreePointer-CodeLength;return;return;if(TreePointer-LeftSubtree)(strcpy(tempstr,TreePointer-Code);strcat(tempstr,2);strcpy(TreePointer-LeftSubtree-C

26、ode,tempstr);TreePointer-LeftSubtree-CodeLength=TreePointer-CodeLength+1;Huffman_Code(TreePointer-LeftSubtree);if(TreePointer-middleSubtree)(strcpy(tempstr,TreePointer-Code);strcat(tempstr,1);strcpy(TreePointer-middleSubtree-Code,tempstr);TreePointer-middleSubtree-CodeLength=TreePointer-CodeLength+1

27、;Huffman_Code(TreePointer-middleSubtree);if(TreePointer-RightSubtree)(strcpy(tempstr,TreePointer-Code);strcat(tempstr,0);strcpy(TreePointer-RightSubtree-Code,tempstr);TreePointer-RightSubtree-CodeLength=TreePointer-CodeLength+1;Huffman_Code(TreePointer-RightSubtree);voidCHuffman_2:Huffman_Coding()(H

28、uffman_Code(HuffTree);/编码结果voidCHuffman_2:Huffman_CodeAnalyzing()(for(inti=0;iISNumber;i+)AvageCodeLength+=ISarrayi.Probability*ISarrayi.CodeLength;intL=1;intm=2;/InformationRate=AvageCodeLength/L*(log(m)/10g(2);doubleHx=0;for(intj=0;jISNumber;j+)Hx+=-ISarrayj.Probabi1ity*1og(ISarrayj.Probabi1ity)/1

29、og(2);CodeEfficiency=Hx/InformationRate;Redundancy=1-CodeEfficiency;voidCHuffman_2:Huffman_Display()(cout码字：endl;for(inti=0;iISNumber;i+)(coutISarrayi.InformationSign:ISarrayi.Codeendl;coutendl;cout“平均码长：AvageCodeLengthendlendl;cout编码效率：CodeEfficiencyendlendl;cout冗余度：RedundancyendlLeftSubtree);Destr

30、oyBTree(TreePointer-middleSubtree);DestroyBTree(TreePointer-RightSubtree);deleteTreePointer;TreePointer=NULL;voidmain()(CHuffman_3YYY;YYY.Huffman_Input();YYY.Huffman_Sort();YYY.Huffman_Tree();YYY.Huffman_Coding();YYY.Huffman_CodeAnalyzing();YYY.Huffman_Display();6.3 运行结果及分析m,C:DOCUMENTSANDSETTINGSAD

31、MINISTRATOR桌面.码字；：iX2r001* X3iOil* X4*：OOOOX5:0100X6,:0101:00010* XW.:00011平均码长，2.61编码效率：0,977011冗余度，0.0229S9Pressanykeytocontinue-图2运行结果运行结果分析：从运行结果上可以看出，该结果与理论计算一致，并且可以看出哈夫曼编码的特点：(1)哈夫曼的编码方法保证了概率大的符号对应于短码，概率小的符号对应于长码，充分应用了短码。(2)缩减信源的最后两个码字总是最后一位不同，从而保证了哈夫曼编码是即时码。(3)每次缩减信源的最长两个码字有相同的码长。这三个特点保证了所得的哈

32、夫曼码一定是最佳码。因此哈夫曼是一种应用广泛而有效的数据压缩技术。利用哈夫曼编码进行通信可以大大提高信道利用率，加快信息传输速度，降低传输成本。数据压缩的过程称为编码，解压的过程称为译码。进行信息传递时，发送端通过一个编码系统对待传数据（明文）预先编码，而接受端将传来的数据（密文）进行译码。要求数据这样一个简单的哈夫曼编码译码器。在进行哈夫曼编码时，为了得到码方差最小的码，应使合并的信源符号位于缩减信源序列尽可能高的。七、体会及建议在这次课程设计中，通过对程序的编写，调试和运行，使我更好的掌握了哈夫曼树等数据结构方面的基本知识和各类基本程序问题的解决方法，熟悉了各种调用的数据类型，在调试和运行

33、过程中，加深我对程序运行的环境了解和熟悉的程度，同时也提高了我对程序调试分析的能力和对错误纠正的能力。这次信息论与编码的程序设计，对于我来说是一个挑战。我对数据结构的学习在程序的设计中也有所体现。课程设计是培养学生综合运用所学知识，发现问题、提出问题、分析问题和解决问题的过程，锻炼学生的逻辑思维能力和实践能力，是对学生实际工作能力的具体训练和考察过程。在整个课程程序中，我们充分应用和调用各个程序模块，从而部分实现了此次程序设计的所应该有的功能。就是我在课程设计是比较成功的方面，而在这个过程中，让我感觉收获最大的就是我们都能利用这次课程设计学到很多我们在课本上没有的知识，充分的发挥了我们的主动性，使我们学会了自主学习，和独立解决问题的能力。很多程序在结构上是独立的，但是本此设计的程序功能不是零散的，它有一个连接是的程序是一个整体，达到这种统一体十分重要，因为这个输出连接是贯穿始终的。通过翻阅相关书籍，在网上查询资料学习有关哈夫曼编码的实现过程，我实在是获益匪浅，更加深刻的感觉到哈夫曼编码能够大大提高通信的效率，对于我们通信工程专业来说