实验二-计算信道容量.doc

《实验二-计算信道容量.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验二计算信道容量一、实验题目1、已知：信源符号个数r、信宿符号个数s、信道转移概率矩阵P。2、输入：任意的一个信道转移概率矩阵。信源符号个数、信宿符号个数和每个具体的转移概率在运行时从键盘输入。3、输出：信道容量C。二、实验目的1、理解和掌握信道容量的概念和物理意义；2、理解计算离散信道容量的迭代算法。三、算法设计四、程序分析1、信道容量：2、当正向传输的信道容量和反向传输的信道容量在误差范围内时表示此时信道稳定，该信道容量即为所求。3、计算正、反向传输的信道容量的迭代算法公式如下：五、程序代码/***********************

2、***********************************************************************【名称】：信道容量迭代法【函数】：计算a[i]函数、计算cap_result函数、计算cap_max函数、主函数【思路】：利用迭代法，使cap_result初步逼近cap_max，当误差小于e时，cap_result即为信道容量******************************************************************************************

3、****/#include#include#includeintX_num,Y_num;//X_num为信源个数，Y_num为信宿个数intn=1;doublee;//迭代法精度误差doublePXi[50];//输入符号的概率P(xi)数组doubleP[50][50];//信道转移概率矩阵doublea[50];doublecap_result;doublecap_max;//迭代法参数a[i]、C、CCdoubleresult;//信道容量/**********************

4、****************************************************函数名：doubleCalculate_a(intk,doublePXi[])功能：计算并输出迭代法所需的参数a[k]参数：返回值：a[k]，也就是P(yi

5、xj)*log(P(yi

6、xj)/P(yi))/log2对i从0到Y_num-1求和**************************************************************************/doubleCalculate_a(intk,doub

7、lePXi[]){inti,j;doublesum1=0;for(j=0;j

8、xi);}sum1+=P[k][j]*log(P[k][j]/sum);//求和得到a[k]}printf("a[%d]=%ft",k,sum1);//输出a[k]returnexp(sum1);}/*****************

9、*********************************************************函数名：doubleCalculate_cap_result(doublePXi[],doublea[])功能：计算并输出迭代法所需的参数cap_result参数：返回值：cap_result，也即是用来逼近最值的参数**************************************************************************/doubleCalculate_cap_result(doubleP

10、Xi[],doublea[]){inti;doublesum=0;for(i=0;i

11、x参数：返回值：cap_max，也即a[k]的最大值*//**********************************************