信息论与编码胡亚洲.docVIP

信息论与编码 实验报告 专业：通信工程 班级：通信三班 学号：120310312 姓名：胡亚洲 实验一 英文文本信息量的计算 一 实验目的 1 通过本实验熟悉Matlab软件编程环境 2 编写M文件实现对英文文本信息量的计算，掌握信源熵的计算方法 二 实验要求 1 了解matlab中M文件的编辑、调试过程 2 编写程序实现英文文本信息量的统计，掌握信源熵及信息量的计算方法 三 实验步骤 ① 查找各个英文字母及空格出现的频率 ② 在Matlab中读取给定的英文文章（该文章为txt格式） ③ 计算英文文章的长度 ④ 统计在该文章中各个字母及空格出现的次数并放入数组N中 ⑤ 计算各个字母和空格的信息量及整篇文章的信息量 ⑥ 计算信源熵 四 给定程序的流程图，以及每条语句的注释 五 自己编写的改进程序，及运行结果记录 wenzhang=textread(english.txt,%c); %读入的英文文章 len=size(wenzhang); length=len(1,1); %文章的长度为length N=zeros(1,52); for i=1:length if double(wenzhang(i,1))96double(wenzhang(i,1))=122 N(1,double(wenzhang(i,1))-96)=N(1,double(wenzhang(i,1)-96))+1; else double(wenzhang(i,1))64double(wenzhang(i,1))=91 N(1,double(wenzhang(i,1))-64)=N(1,double(wenzhang(i,1)-64))+1; end end %统计各字母出现的个数，存入N数组中 p=zeros(1,52); p=N/length; %计算各字母出现的概率 I=zeros(1,52); I=log2(1./p);%计算各字母的信息量 for i=1:52 if p(1,i)==0 I(1,i)=0; end end %出现概率为0的字母的信息量赋值 sumI=sum(I.*N)%计算整篇文章的总的信息量 Ip=I.*p; H=sum(Ip) %计算信源熵 文本内容：huyazhouxinxilunyubianma 运行结果：sumI = 81.7744 H = 3.4073 六实验小结 通过本次实验熟悉了Matlab软件编程环境和一些函数的功能及使用，掌握了信息量、信源熵的计算方法。 实验二 信息熵与图像熵计算 一 实验目的 1 通过本实验复习MATLAB的基本命令，熟悉MATLAB下的基本函数 2 复习信息熵基本定义，能够自学图像熵定义和基本概念 二 实验要求 1 能够写出MATLAB源代码，求信源的信息熵 2 根据图像熵基本知识，设计出MATLAB程序，求出给定图像的图像熵 三 实验步骤 1 求解信息熵过程： ① 输入一个离散信源，并检查该信源是否是完备集。 ② 去除信源中符号分布概率为零的元素。 ③ 根据平均信息量公式，求出离散信源的熵。 2 图像熵计算过程： ① 输入一幅图像，并将其转换成灰度图像。 ② 统计出图像中每个灰度阶象素概率。 ③ 统计出图像中相邻两象素的灰度阶联合分布矩阵。 ④ 根据图像熵和二阶熵公式，计算出一幅图像的熵。 四、实验数据及结果分析： 函数程序 CalEntropy.m % Information Shannon Entropy calculation % jma@hebut.edu.cn, 22/08/2007 % array : Discrete Probabilities Set % H : Output Shannon Entropy function H = CalEntropy(array) % Vector number num = length(array); % Check probabilities