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

xxxx学院综合性实验报告 专业：xxxx 年级/班级：xxxx级xxxx—xxxx学年第-一学期 课程名称 信息论与编码 指导教师 本组成员 学号姓名 实验地点 实验时间 项目名称 二元哈夫曼编码 实验类型 综合性 一、实验口的 用Matlab语言编程实现霍夫曼（Huffman）编码。 二、 实验仪器或设备 一台计算机。 MATLAB 软件。 三、 总体设计 将信源消息符号按其出现的概率大小依次排列，pl$p22…2pq 取两个概率最小的字母分别配以0和1两个码元，并将这两个概率相加 作为一个新字母的概率，从而得到只包含q-1个符号的新信源SU 对重排后的缩减信源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-l个概率；第（n, 2）元素放置合并元素在 木列中位置； ??? 如果输入n个概率，则需要缩减n-l次后，剩下t=2两个概率； 所以编码表可以定义为n行n-l列。 编码表生成后，从最后一列开始编码。 首先最后一列两个元素分别编码为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); %b=sort(p, descend，); b=f1iplr(sort(p)); for i=l:m-1 %p=sort (p, descend); p=fliplr (sort (p)); n=size(p, 2); a(:, i)二p; p(n-i)=p(n-i + l)+p(n-i); for j=n-i+l:m p(j)=o； end end %p中元索个数 %用于般后的输岀 %对输入概率矩阵P进行编码前的准备，用矩阵a记录哈弗曼编码的过程 n=m; for i=l:m~2 x=a(n-i+l, i)+a(n-i, i); r=find(a(:, i+l)==x); a(n, i+l)=r (end); end %矩阵a从第二列开始，每一列的最后一个元索用于保存前次求和后元索所处位置 a %- %输出编码衣 %根据a対所输入概率矩阵p进行哈州曼编码 codcl二sym( [0 1]); %给最后一列的元素编码 h二codel; t 二3; d二 1; for j二n-2:-1:1 for i=l:t-2 %控制编码个数，初始为3 %从倒数第二列开始依次对各列元素编码 if il a(i, j)=a(i-l, j) d=d+l; else d=l; end a (a (n, j+1), j+1)二T ; %对已经编码的赋值为T temp=a(:,j+1); x二find(temp==a(i, j)); h(i)=codel (x(d)); end y=a (n, j+1); h(t-l) = [char (codel(y)), O]; h(t) = [char(codel (y)), 1]; t=t+l; codel二h; end %排序后的原概率序列 h %编码结果 结果： ? P二[0. 2 0. 19 0. 18 0. 17 0. 15 0. 1 0. 01]; huffman(p) a =0. 2000 0. 2000 0. 2600 0. 3500 0. 3900 0.6100 0. 1900 0. 1900 0. 2000 0. 2600 0. 3500 0. 3900 0. 1800 0. 1800 0. 1900 0. 2000 0. 2600 0 0. 1700 0. 1700 0. 1800 0. 1900 0 0 0. 1500 0. 1500 0. 1700 0 0 0 0. 1000 0. 1100 0 0 0 0 0.0100 6. 0000 1.0000 1.0000 1.