信息论实验最终版.doc

实验一 信道容量的迭代算法程序设计 1、实验目的 （1）进一步熟悉信道容量的迭代算法； （2）学习如何将复杂的公式转化为程序； （3）掌握C语言数值计算程序的设计和调试技术。 2、实验要求 （1）已知：信源符号个数r、信宿符号个数s、信道转移概率矩阵P。 （2）输入：任意的一个信道转移概率矩阵。信源符号个数、信宿符号个数和每个具体的转移概率在运行时从键盘输入。 （3）输出：最佳信源分布P*，信道容量C。 3、信道容量迭代算法 1：procedure CHANNEL CAPACITY(r,s,()) 2：initialize:信源分布=1/r，相对误差门限，C=— 3：repeat 4： 5： 6： C ( 7：until 8：output P*= ,C 9：end procedure ------------------------------------------------------------------------------------------------------- 4、参考代码 自己编写的代码： /*********************************************************************Author: Wan Jiaolong *Date : 2012.04.1 *Copyright: WenjorinLab Purpose: Caculate the capacity of a given channel ********************************************************************* #includestdio.h #includemath.h #includeiostream.h //用宏定义一个DELTA值 #define DELTA 1e-4 float matrixValue( int n, float **a ); float **InverseMatrix(float **A,int n); int main() { int sourseNum,destinationNum,i,j,k,counter; float **matrix=NULL; float **inv_matrix=NULL; float *probilityData=NULL; float *hData=NULL; float *sum=NULL; float *w=NULL; float C,C_old; cout请输入信源数据个数和信道数据个数endl; cinsourseNumdestinationNum; if(sourseNum!=destinationNum) { cout转移矩阵必须是方阵，请再输入一次endl; cinsourseNumdestinationNum; } probilityData=new float [sourseNum]; matrix=new float *[sourseNum]; for(i=0;isourseNum;i++) matrix[i]=new float [sourseNum]; inv_matrix=new float *[sourseNum]; for(i=0;isourseNum;i++) inv_matrix[i]=new float [sourseNum]; cout请输入转移矩阵概率：\n; for(i=0;isourseNum;i++) for(j=0;jsourseNum;j++) cinmatrix[i][j]; while(matrixValue(sourseNum,matrix)==0)//[sourseNum][sourseNum] { cout请再次输入转移矩阵概率：\n; for(i=0;isourseNum;i++) for(j=0;jsourseNum;j++) cinmatrix[i][j]; } //求取行列式的逆矩阵 inv_matrix=InverseMatrix(matrix,sourseNum);//[sourseNum][sourseNum]---- //int matrixSum=matrixValue(sourseNum,matrix[sourseNum][so