信息论绘制熵函数matlab,实验1绘制二进熵函数曲线串联信道容量曲线.doc

实验1绘制二进熵函数曲线串联信道容量曲线.doc

信息论与编码 实验 1 绘制二进熵函数曲线、串联信道容量曲线 学院:物理与电子学院 班级:电信1105班 姓名: 学号:140411072一、实验内容 用 Excel 或 Matlab 软件制作二进熵函数曲线、串联信道容量曲线。 二、实验目的 1. 掌握 Excel 的数据填充、公式运算和图表制作 2. 掌握 Matlab 绘图函数 3. 掌握、理解熵函数表达式及其性质 三、实验要求 1. 提前预习实验，认真阅读实验原理以及相应的参考书。 2. 认真高效的完成实验，实验中服从实验室管理人员以及实验指导老师的管理。 3. 认真填写实验报告。 四、实验原理 1. Excel 的图表功能(略) 2. 信源熵的概念及性质             ) ( 1 ) ( 1 . log ) ( . ) ( 1 log 1 log ) ( log ) ( ) ( 1 0 , 1 1 0 ) ( 2 1 Q H P H Q P H b n X H a p H p p p p x p x p X H p p p x x X P X i i i                                      3. 串联信道的信道容量(图 1) 。 p  1 p  1 p p 0 1 1 1 0 p  1 p  1 p p 0 1 1 1 0 p  1 p  1 p p 0 1 1 1 0 I II III 图 1 三个二元对称信道(BSC)的串联 串联 信道转移矩阵 信道容量 I          p p p p x y p 1 1 )] / ( [ ) ( 1 p H C I  I,II                                        2 2 2 2 1 1 2 1 2 11 1 1 1)] / ( )][ / ( [ )] / ( [ p p p p p p p p p p p p p p p p y z p x y p x z p     p p H C II I    1 2 1 , I,II,III                                                            p p p p p p p p p p p p p p p p p p p p p p p p z w p x z p x w p 1 3 1 1 3 1 3 1 3 11 1 1 1 2 1 2 1)] / ( )][ / ( [ )] / ( [ 2 3 2 3 2 3 2 3 2 2 2 2     2 3 , , 1 3 1 p p p H C III II I     五、实验步骤： 1. 信息熵函数 H.m 文件 function y = H(p) y=-(p.*log(p)+(1-p).*log(1-p))./log(2); end 2 matlab 命令 p=0.01:0.01:0.99; y1=H(p); y2=H(2.*p.*(1-p)); y3=H(p.^3+3.*p.*(1-p).*(1-p)); c1=1-y1;c2=1-y2;c3=1-y3; subplot(3,1,1);plot(p,c1);xlabel( p );ylabel( C1 ); subplot(3,1,2);plot(p,c2);xlabel( p );ylabel( C2 ); subplot(3,1,3);plot(p,c3);xlabel( p );ylabel( C3 );六、实验心得通过本次信息与编码实验，了解了二进制上函数的曲线以及串联信道容量曲线，使自己对书上知识 点的把握更加直观，收获很大，更深刻的了解了曲线产生的原理，比书上了解的更加具体，使自己受益 匪浅，增加了学习兴趣。