实验一信道容量的迭代算法

1、陕西科技大学实验报告班级： 学号： 姓名： 实验组别： 实验日期： 2012-12-18 报告日期： 2012-12-18 成绩： 报告内容：（目的和要求、原理、步骤、数据、计算、小结等）实验名称 ：信道容量的迭代算法一、实验目的1） 进一步熟悉信道容量的迭代算法。2） 学习如何将复杂的公式转化为程序。3） 掌握高级语言数值计算程序的设计和调试技术。二、实验要求（1）已知：信源符号个数 r 、信宿符号个数 s 、信道转移概率矩阵 P。 （2）输入：任意的一个信道转移概率矩阵。信源符号个数、信宿符号个数和每 个具体的转移概率在运行时从键盘输入。（ 3）输出：信道容量 C。三、实验原理（流程图）

2、：1: procedure ChannelCapacity( r, s, ( pji )12: initialize: 信源分布 pi 1, 相对误差门限 , C r3: repeat4:pi pjirpi pjii15:exp pji log ijj1pirrj1s6:exppji log iji 1 j1rC log expi1j1pjilog ij7: until8: output C9: end procedure流程图 :四、实验数据(源代码) ：#include #include #define MAX 100 double Calculate_a(int k,double pa)

3、; double Calculate_C1(double pa,double a); double Calculate_C2(double a);int r,s; double pba33; void main() int i,j;double C1,C2,E; double aMAX,pa3; E=0.000001;r=3;word 文档 可自由复制编辑 s=3;pa0=0.2;pa1=0.3; pa2=0.5; pba00=0.5; pba01=0.33333333; pba02=0.16666666; pba10=0.16666666; pba11=0.5; pba12=0.333333

4、333; pba20=0.33333333; pba21=0.16666666; pba22=0.5;dofor(i=0;i=E)double sum=0; for(i=0;ir;i+) sum+=pai*ai;for(i=0;ir;i+) pai=pai*ai/sum;elseprintf( 最佳信源概率 :n); for(i=0;i=E);printf( 信道容量为 :%lfn,C1/log(2); double Calculate_a(int k,double pa)int i,j;double temp,sum2=0; for(j=0;js;j+) double sum1=0; for

5、(i=0;ir;i+)word 文档 可自由复制编辑 sum1+=pai*pbaij; temp=pbakj/sum1; temp=log(temp);sum2+=pbakj*temp; return exp(sum2);double Calculate_C1(double pa,double a) int i;double sum=0; for(i=0;ir;i+) sum+=pai*ai; return log(sum);double Calculate_C2(double a)int i;double max=a0; for(i=0;ir;i+) if(maxai) max=ai; return log(max);五、程序运行结