实验二霍夫曼编码.docxVIP

信息论与编码实验报告 院系: 专业 学号 姓名 哈尔滨理工大学荣成校区 电子信息工程 霍夫曼编码 信息论与编码第二次实验报告 —、实验目的 根据霍夫曼编码的原理，川MATLAB设计进行霍夫曼编码的程序，并得岀正确的结果。 二、实验仪器或设备 1、 计算机； 2、 MATLAB 2013b。 三、二元霍夫曼编码原理 1、 将信源消息符号按其出现的概率大小依次排列，p} p2…… pn ; 2、 取两个概率最小的字母分别配以0和1两个码元，并将这两个概率相加作为一个新 字母的概率，从而得到只包含q-1个符号的新信源S1； 3、 对垂诺后的缩减信源SIS?以递减次序排序，两个概率最小符号重复步骤2的过程。 4、 不断继续上述过程，直到最后两个符号配以0和1为1卜， 5、 从最后一级开始，向前返回得到各个信源符号所对应的码元序列， 即相应的码字。 function [h,1]=huffmancode(P) %输入概率数组P%表示p的数组长度%对卩进行冒泡排序（从大到小） %输入概率数组P %表示p的数组长度 %对卩进行冒泡排序（从大到小） for i=l:n-1 for j=i+l:n if P(i)=P(j) p=P (i); P(i)=P(j)； P(j)=p; end disp （ 1概率分布T Q=P； m=zeros(n-1,n); for i=l:n-1 咎初始化全零数组 [Q,b]=sort (Q) ; %sort换数是对Q进行升序排列，返回值1显示排序后位置的变动信息 m (iz :) = [b(l:n-i + l),zeros(l,i-l)]; Q=[Q(l)+Q(2),Q(3:n),1]; end for i=l:n-1 %最小两数相加与第三个数比较 c(i, :)=blanks (n*n); end %blanks是空格函数，创建空格串 c(n-1z n) = 111; 概率排在下面 c(n-l,2*n)=,0,; for i=2:n-1 %以下计算各个元素码字，概率相同时下一次输出 c(n-izl:n-l)=c(n-i + l,n*(find(m(n-i + 1, :)==]))-(n-2) :n*(find(m(n-i + 1,:) ==1))); c(n-i,n)=111; 咎对每次相加的两数进行0、1编码 c (n-if n+1:2*n-l)=c(n-iz1:n-1); c(n-i,2*n)=101; for j=l:i-1 %对每次和加的两数进行0、1编码 c(n-i, (j + 1)*n + l: (j + 2)*n)=c(n-i + l,n*(find(m(n-i + lz :)==j + l)-l)+l:n*find (m(n-i+l,:)==j+1)); end end for i=l:n h(i,l:n)=c(l,n* (find (m (lz :) ==i) -1)+1: find (m (lz :) ==i) *n) ; % 从右向左读码 11 (i) =length (find (abs (h (iz :) ) ~=32) ) ; %输出码 K end l=sum(P.*11); hl=-log2(P); H=P*(hl 1); disp ( *信息爛1 ) $平均码长 %自信息量 %求爛 H disp (，二兀崔夫曼编码平均码长?) 1 $显示平均码1 丄 disp( *编码效率T G=H/1 disp「二元霍夫曼编码1 ) %显示编码效率G 咎显示码字 五、实验结果 在 command window 中输入的概率数组为 p二[0.4 0.18 0.1 0」0.07 0.06 0.05 0.04]后输出 结果如下 huffmancode 输入概率:[0.4 0.18 0.1 0.1 0.07 0.06 0.05 0.04] 概率分布 信息爛 H = 2.5524 二元霍夫曼编码平均码长 1 = 2.6100 编码效率 G = 0.9779 二元崔夫曼编码 ans = 1 001 0000 Oil 0100 0101 00010 00011