- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
离散信源的熵信道容量实验报告
计算离散信源熵、离散信道容量
实验任务和目的
实验任务:
(1)简要总结信源的熵、信道容量的物理意义,概念;
(2)写出离散信源熵、离散信道容量计算的基本步骤,画出实现离散信源熵、离散信道容量计算的程序流程图;
(3)讨论信源的熵的大小与前后符号之间相关性的关系,讨论信道容量与信源先验概率及信道转移概率的关系。
实验目的:
掌握信源的熵、信道容量的物理意义,概念;熟练掌握离散信源熵、离散信道容量的计算方法步骤;利用Matlab编写离散信源熵、离散信道容量的计算程序;验证程序的正确性。
实验过程和结果
实验过程
1、简要总结信源的熵、信道容量的物理意义,概念。
信源熵的物理意义是指信源中的各个符号的平均不确定性;熵是信源符号的平均信息量,是信源符号的平均不确定度。
信道容量 概念:在信道可以传输的基本前提下,对信源的一切可能的概率分布而言,信道能够传输的最大(接收)熵速率称为信道容量。
意义:求出了某个信道的信道容量,也就找到了信源的最佳概率分布。从而指导人们改造信源,使之最大可能地利用信道的传输能力。
2、写出离散信源熵、离散信道容量计算的基本步骤,画出实现离散信源熵、离散信道容量计算的程序流程图;
离散信源熵的计算步骤:
信道容量的计算步骤:
3、(1)讨论信源的熵的大小与前后符号之间相关性的关系,讨论信道容量与信源先验概率及信道转移概率的关系。
信源的相关性是信源符号间的依赖程度的度量。由于信源输出符号间的依赖关系也就是 信源的相关性使信源的实际熵减小。信源输出符号间统计约束关系越长,信源的实际熵越小。当信源输出符号间彼此不存在依赖关系且为等概率分布时,信源的实际熵等于最大熵。
(2)信道容量与信源先验概率及信道转移概率的关系。
信道容量是信道的一个参数,反映了信道所能传输的最大信息量,其大小与信源无关。对不同的输入概率分布,互信息一定存在最大值。我们将这个最大值定义为信道的容量。一但转移概率矩阵确定以后,信道容量也完全确定了。尽管信道容量的定义涉及到输入概率分布,但信道容量的数值与输入概率分布无关。我们将不同的输入概率分布称为试验信源,对不同的试验信源,互信息也不同。其中必有一个试验信源使互信息达到最大。这个最大值就是信道容量。
实验结果
计算离散信源熵:
计算离散信道容量:
实验总结和心得
通过本次的试验,让我加深了对信息论与编码这门课知识点的认识,学会了怎么样去计算信源熵以及信道容量,怎样去把它具体实现出来。在此次实验过程中,让我认识到了做实验之前预习是很重要的,通过查找资料可以解决一些我们不能解决的问题的方法。还有就是不懂的问题及时问老师,能解决的问题随堂解决。
附录(代码)
计算信道容量:
clc;clear?all;?
N?=?input(输入信源符号X的个数N=);?
?M?=?input(输出信源符号Y的个数M=);??
p_yx=zeros(N,M);???%程序设计需要信道矩阵初始化为零?
fprintf(输入信道矩阵概率\n)?
for?i=1:N
?for?j=1:M?
?p_yx(i,j)=input(p_yx=);%输入信道矩阵概率??????
if?p_yx(i)0?
?error(不符合概率分布)??????
end??????
end?
end?
for?i=1:N???%各行概率累加求和?????
?s(i)=0;?????
?for?j=1:M?
s(i)=s(i)+p_yx(i,j);?????
?end?
end?
for?i=1:N?????%判断是否符合概率分布
?if?(s(i)=0.999999||s(i)=1.000001)???
?error(不符合概率分布)?
end
?end?
b=input(输入迭代精度:);%输入迭代精度?for?i=1:N?
?p(i)=1.0/N;???%取初始概率为均匀分布
?end?
for?j=1:M???%计算q(j)??????
q(j)=0;????
for?i=1:N?
q(j)=q(j)+p(i)*p_yx(i,j);???
?end?
end?????????
for?i=1:N?????%计算a(i)??????d(i)=0;?
?for?j=1:M?
if(p_yx(i,j)==0)????????????d(i)=d(i)+0;??????????else?
?d(i)=d(i)+p_yx(i,j)*log(p_yx(i,j)/q(j));???????
?end????
end?
?a(i)=exp(d(i));?end?u=0;?
for?i=1:N???%计算u??????
u=u+p(i)*a(i);?end?
IL=log2(u);?%计算IL?
IU=log2(max(a));%计算IU?n
原创力文档


文档评论(0)