信息论实验二报告.docVIP

实 验 报 告 课程名称 信息论与编码 实验项目 霍夫曼编码、率失真函数的计算 实验仪器 计算机 系 别 信息与通信工程学院 专 业 电子信息工程 班级/学号 电信1103/2011010550 学生姓名 张 琪 实验日期 成 绩 指导教师 实验项目：霍夫曼编码、率失真函数的计算 实验项目性质：验证性实验 课程：《信息论与编码》 计划学时：4学时 一、实验目的 掌握使用计算机进行霍夫曼编码、计算率失真函数的方法。 二、实验内容 利用matlab语言编程实现霍夫曼编码、计算率失真函数。 三、实验原理 1．参照教材（焦瑞莉等编著《信息论基础教程》）第111页表4-8理解霍夫曼编码的过程 2．参照教材（焦瑞莉等编著《信息论基础教程》）第148页式（6-7）理解率失真函数的定义 四、实验方法与实验步骤 （一）参照教材第223页信源熵程序建立文件huffman.m 文件huffman.m的代码： % huffman编码生成器 % % 函数说明： % % [W，L．q]=huffman(P)为huffman编码函数 % % P为信源的概率矢量，W为编码返回的码字 % % L为编码返回的平均码字长度，q为编码效率 % % ***************************************% function [W,L,q]=huffman(P) if(length(find(P=0))~=0) error(Not a prob.vector.negative component); % 判断是否符合概率分布条件 end if(abs(sum(P)-1)10e-10) error(Not a ponent do not add up to 1) % 判断是否符合概率和为1 end n=length(P); % 计算输入元素个数 p=P; mark=zeros(n-1,n); % mark为n-1行、n列矩阵，用来记录每行最小两概率叠加后概率排列次序 % 1) 确定概率大小值的排列，得到mark矩阵。 for i=1:n-1 [p,num]=sort(p); % 对输入元素排序并纪录 mark(i,:)=[num(1:n-i+1),zeros(1,i-1)]; p=[p(1)+p(2),p(3:n),1]; end % 2)生成一个n-1行、n1(n×n)列矩阵table，每行可看做n个段， % 每段长为n,记录一个码字(每个码字的长度不会超过n)。 for i=1:n-1 table(i,:)=blanks(n*n); end % 3)计算各个元素码字，循环n 2次，决定矩阵table % 从倒数第二行开始到第行的每段的码字值，到编码表格table table(n-1,n)=l; % 小值赋1 table(n-1,2*n)=0; % 大值赋0 for i=2:n-1 table(n-i,1:n-1)=table(n-i+1,n*(find(mark(n-i+1,:)==1))-(n-2)... :n*(find(mark(n-i+1,:)==1))); % 按mark的记录依次赋值 table(n-i,n)=1; table(n-i,n+1:2*n-1)=table(n-i,1: n-1); table(n-i,2*n)=0; for j=1:i-1 table(n-i,(j+1)*n+1:(j+2)*n)=table(n-i+1,... n*(find(mark(n-i+1,:)==j+1)-1)+1:n*find(mark(n-i+1,:)==j+1)); % 按mark的记录依次赋值 end end % 4)得到编码后的码字 for i=1:n W(i,1:n)=tabl