信息论与编码实验报告-信源熵值的计算.doc

《信息论与编码实验报告-信源熵值的计算.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、个人收集整理勿做商业用途实验报告课程名称:信息论与编码姓名：系：专业：年级：学号：指导教师：职称：年月日个人收集整理勿做商业用途实验一信源熵值的计算一、实验目的1进一步熟悉信源熵值的计算2熟悉Matlab编程二、实验原理熵（平均自信息）的计算公式MATLAB实现：；或者流程：第一步：打开一个名为“nan311"的TXT文档,读入一篇英文文章存入一个数组temp,为了程序准确性将所读内容转存到另一个数组S，计算该数组中每个字母与空格的出现次数(遇到小写字母都将其转化为大写字母进行计数），每出现一次该字符的计数器+1;第二步：计算信源总大小计算出每个字母和空格出现

2、的概率；最后,通过统计数据和信息熵公式计算出所求信源熵值（本程序中单位为奈特nat）。程序流程图：个人收集整理勿做商业用途三、实验内容1、写出计算自信息量的Matlab程序2、已知:信源符号为英文字母（不区分大小写）和空格。输入：一篇英文的信源文档。输出：给出该信源文档的中各个字母与空格的概率分布，以及该信源的熵。四、实验环境MicrosoftWindows7Matlab6.5五、编码程序#include”stdio。h”＃include〈math。h〉#include〈string。h〉#defineN1000intmain（void){chars[N]；i

3、nti，n=0；floatnum［27］={0}；doubleresult=0,p［27]=｛0｝；FILE＊f;char＊temp=newchar［485];f=fopen（"nan311.txt"，"r”）;while(!feof(f）){fread(temp,1，486，f)；｝fclose（f);s［0］=*temp;for（i=0；i

4、m［s［i]-97］++；elseif（s［i]>=’A'＆＆s[i]<='Z’)个人收集整理勿做商业用途num[s[i]—65]++；｝printf（”文档中各个字母出现的频率：");for（i=0;i〈26；i++）{p［i］=num［i]/strlen(s）;printf（”%3c:%ft"，i+65,p[i］）;n++;if（n==3)｛printf("”);n=0；}｝p[26］=num[26]/strlen（s)；printf（”空格:％ft”,p[26]）；printf（”"）；for(i=0；i〈27；i++）｛if（p［i]！

5、=0）result=result+p[i］*log（p[i］）;}result=-result；printf(”信息熵为:％f",result);printf(""）;return0；｝六、求解结果其中nan311.txt中的文档如下：Thereisnohatewithoutfear.Hateiscrystallizedfear,fear’sdividend，fearobjectivized。Wehatewhatwefearandsowherehateis，fearislurking.Thuswe个人收集整理勿做商业用途hatewhatthreatenso

6、urperson,ourvanityandourdreamsandplansforourselves.Ifwecanisolatethiselementinwhatwehatewemaybeabletoceasefromhating。七、实验总结通过这次实验，我们懂得了不必运行程序时重新输入文档就可以对文档进行统计，既节省了时间而且也规避了一些输入错误.在实验中，我们进一步了解到信源熵的计算，理论和实践的结合让我们对这个知识点了解的更加深刻了。