《信息论基础实验》指导老师:贺正芸 班级：信息工程 081姓名：卢慈荣学号：湖南工业大学电气与信息工程学院实验一 信道容量的迭代算法程序设计一、实验目的（1）进一步熟悉信道容量的迭代算法；（2）学习如何将复杂的公式转化为程序；（3）掌握 C 语言数值计算程序的设计和调试技术。二、实验要求（1）已知：信源符号个数 r、信宿符号个数 s、信道转移概率矩阵 P。（2）输入：任意的一个信道转移概率矩阵。信源符号个数、信宿符号个数和每个具体的转移概率在运行时从键盘输入。（3）输出：最佳信源分布 P*，信道容量 C。三、信道容量迭代算法 1：procedure CHANNEL CAPACITY(r,s,())2：initialize:信源分布=1/r，相对误差门限，C=—3：repeat4：5：6：C 7：until 8：output P*= ,C9：end procedure -------------------------------------------------------------------------------------------------------四、实验代码/*问题：１ 初始最大容量的设定２ exp 的精确求解 */#include#include#include#include#include#define R 1000#define S 1000#define delta 1e-2#define inf 1e6 using namespace std; double P_i[R],P_ji[S][R],Thi_ij[R][S]; double Pre_C,Now_C; int r,s,Num; double __log2(double a) { return log(a)/log(2); } double __exp(double a) { return pow(2.7045 ,a); } int eps( double a) { if(a>delta||a<-delta) return 1; return 0; } void scan() { int i,j,k; int flag=0; double t; freopen("信道信息.txt","r",stdin); printf("本次信道信息如下：（若要更改信道信息，请终止程序运行后，改变文件信道信息.txt 中的信息）\n"); printf("信源处信息的个数："); scanf("%d",&r); printf("%d\n\n",r); printf("信道接收处信息的种类："); scanf("%d",&s); printf("%d\n\n",s); printf( "信道的分配（每一行表示信源处的一个信息对接收端一处信息的概率分布）:\n"); for(i=0;i
