信息论实验报告2.docVIP

实 验 报 告 课程名称 信息论与编码 实验项目霍夫曼编码、率失真函数的计算 指导教师 学 院 光电信息与通信工程 _ 专 业 电子信息工程 班级/学号 学生姓名 实验日期 ______ ____ 成 绩 _____________________ 实验目的 掌握使用计算机进行霍夫曼编码、计算率失真函数的方法。 实验内容 利用matlab语言编程实现霍夫曼编码、计算率失真函数。 实验原理 用huffman编码方法得到的码并非是唯一的，因为首先每次对辅助信源最后两个概率的最小的符号，用0和1码符号可以是任意的，所以可以得到不同的码。但他们只是码字具体形式不同，而其码长不变，平均码长也不变，所以没有本质差别。率失真函数表示当失真不超过D时传输所需的最小平均互信息量。它的单位就是平均互信息量的单位。 源程序代码 function H=entropy(P,r) if ~isempty(find(P=0,1)) error(Not aprob.Vector,negative component); end if(abs(sum(P)-1)10e-10) error(Not aprob cector,component do not add up to 1); end H=(sum(-P.*log2(P)))/(log2(r)+eps); Exercise4_8_1.m close all;clear all;clc; % Information Theory exp2 % author：Wang Haixi % time:29-Nov-2011 23:05:51 % since:1.0 % See also huffman, RateDF. disp(==========P117练习4.8(1)==========); huffman([0.49,0.14,0.14,0.07,0.07,0.04,0.02,0.02,0.01]); disp(==========P161练习6.2=============); RateDF([0.5,0.5],[0,1,0.25;1,0,0.25],-100); % huffman编码生成器 % % 函数说明： % % [W,L,q]=huffman(P) 为huffman编码函数 % % P为信源的概率矢量，W为编码返回的码字 % % L为编码返回的平均码字长度，q为编码效率 % %***************************************** % function [W,L,q]=huffman(P) if (~isempty(find(P=0, 1))) error(Not a prob.vector,negative component); % 判断是否符合概率分布条件 end if (abs(sum(P)-1)10e-10) error(Not a prob.vector,component 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)=1; % 小值赋1 table(n-1,2*n)=0; % 大值赋0 for i=2:n-1 table(n-i,1:n