信道容量的迭代算法.docVIP

班级： 学号：姓名： 实验组别： 实验日期： 报告日期： 成绩： 报告内容：（目的和要求、原理、步骤、数据、计算、小结等） 实验名称：信道容量的迭代算法 一、实验目的 进一步熟悉信道容量的迭代算法。 学习如何将复杂的公式转化为程序。 （3） 掌握高级语言数值计算程序的设计和调试技术。 二、实验要求 （1）已知：信源符号个数r、信宿符号个数s 、信道转移概率矩阵P。 （2）输入：任意的一个信道转移概率矩阵。信源符号个数、信宿符号个数和每个具体的转移概率在运行时从键盘输入。 （3）输出：信道容量C。 三、实验原理（流程图）： 输入 是 否 终止 实验数据（源代码）： #includeiostream using namespace std; #includemath.h #define MAX 50 double Calculate_a(int k,double pa[]); double Calculate_C1(double pa[],double a[]); double Calculate_C2(double a[]); int r,s; double pba[MAX][MAX]; void main() {int i,j; double C1,C2,E; double a[MAX],pa[MAX]; E=0.000001; cout请输入信源符号个数r:\n; cinr; cout请输入信宿符号个数s:\n; cins; cout请输入精确度E:\n; cinE; cout请输入信源P[ai]:\n; for(i=0;ir;i++) cinpa[i]; cout请输入信道转移概率矩阵P[bj][ai]:\n; for(i=0;ir;i++) for(j=0;js;j++) cinpba[i][j]; do{for(i=0;ir;i++) {int i,j;double temp,sum2=0; for(j=0;js;j++) {double sum1=0; for(i=0;ir;i++) {sum1+=pa[i]*pba[i][j]; } temp=pba[k][j]/sum1; temp=log(temp); sum2+=pba[k][j]*temp; } A[i]=exp(sum2);} double sum=0; for(i=0;ir;i++) sum+=pa[i]*a[i]; return log(sum);} double Calculate_C2(double a[]) {int i;double max=a[0]; for(i=0;ir;i++) if(maxa[i]) max=a[i]; C1=log(max); double max=a[0]; for(i=0;ir;i++) if(maxa[i]) max=a[i]; C2=log(max); if(C2-C1=E) {double sum=0; for(i=0;ir;i++) sum+=pa[i]*a[i]; for(i=0;ir;i++) pa[i]=pa[i]*a[i]/sum;} else{ cout最佳信源概率:\n; for(i=0;ir;i++) coutpa[i];} } while(C2-C1=E); cout信道容量为:C1/log(2);} 程序运行结果： 实验小结： 在做本次实验过程中，我深刻体会到细节决定成败，同时也发觉自己知识学的太浅薄，没有从本质上把握住要点，运用知识也有些生硬。自己给自己上了一堂课，在今后学习当中，要注意去深度思考，把所学知识有机联系起来，从而把书由薄变厚，再总结进而把书由厚变薄，掌握本质，理解要领。此次实验进一步加强了我的编程能力和严谨的逻辑思维能力，受益颇多。 附页 4 陕西科技大学实验报告 C(n+1,n)—(n+1,n) C=C(n+1,n)