信息论与编码综合性实验报告.docxVIP

xxxx学院 xxxx学院综合性实验报告 专业：xxxx 年级/班级： xxxx级 xxxx—xxxx学年第一学期 课程名称 信息论与编码 指导教师 本组成员 学号姓名 实验地点 实验时间 项目名称 二元哈夫曼编码 实验类型 综合性 一、实验目的 用Matlab语言编程实现霍夫曼（Huffman）编码。 二、实验仪器或设备 一台计算机。 MATLAB软件。 三、总体设计 将信源消息符号按其出现的概率大小依次排列，p1≥p2≥…≥pq 取两个概率最小的字母分别配以0和1两个码元，并将这两个概率相加作为一个新字母的概率，从而得到只包含q-1个符号的新信源S1。 对重排后的缩减信源S1重新以递减次序排序，两个概率最小符号重复步骤(2)的过程。 不断继续上述过程，直到最后两个符号配以0和1为止。 从最后一级开始，向前返回得到各个信源符号所对应的码元序列，即相应的码字。 先生成的编码表 0.4 0.4 0.4 0.4 0.2 0.2 0.4 0.6 0.2 0.2 0.2 0 0.1 0.2 0 0 0.1 2 2 2 编码表每列为缩减信源概率 第一列为输入n个信源概率； 第二列为第一次缩减后t=n-1个概率；第（n，2）元素放置合并元素在本列中位置； … 如果输入n个概率，则需要缩减n-1次后，剩下t=2两个概率； 所以编码表可以定义为n行n-1列。 编码表生成后，从最后一列开始编码。 首先最后一列两个元素分别编码为0和1； 对其中某j列，其前t-2个元素编码为j+1列非合并概率移植过来。（问题,如何知道j=1列中非合并元素的位置？） 第t-1元素和第t个元素为j+1列合并概率编码分别加0和加1。（合并概率的位置已经在本列最后一个元素中标出） 四、实验步骤 %输入一个概率矩阵p，判断其合理性，如果合理则通过函数调用对其进行哈弗曼编码 function huffman(p) if sum(p)~=1 %判断所输入概率矩阵p的元素的和是否为1 display(please input the right numbers!) return end %如果所输入概率矩阵合理，则进行哈弗曼编码 m=size(p,2); %p中元素个数 %b=sort(p,descend); b=fliplr(sort(p)); %用于最后的输出 for i=1:m-1 %对输入概率矩阵p进行编码前的准备，用矩阵a记录哈弗曼编码的过程 %p=sort(p,descend); p=fliplr(sort(p)); n=size(p,2); a(:,i)=p; p(n-i)=p(n-i+1)+p(n-i); for j=n-i+1:m p(j)=0; end end n=m; for i=1:m-2 %矩阵a从第二列开始，每一列的最后一个元素用于保存前次求和后元素所处位置 x=a(n-i+1,i)+a(n-i,i); r=find(a(:,i+1)==x); a(n,i+1)=r(end); end a %输出编码表 %--------------------------------------------------------- %根据a对所输入概率矩阵p进行哈弗曼编码 code1=sym([0 1]); %给最后一列的元素编码 h=code1; t=3; %控制编码个数，初始为3 d=1; for j=n-2:-1:1 %从倒数第二列开始依次对各列元素编码 for i=1:t-2 if i1 a(i,j)==a(i-1,j) d=d+1; else d=1; end a(a(n,j+1),j+1)=-1; %对已经编码的赋值为-1 temp=a(:,j+1); x=find(temp==a(i,j)); h(i)=code1(x(d)); end y=a(n,j+1); h(t-1)=[char(code1(