《信息论基础实验》指导老师:贺正芸 班级:信息工程 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 -------------------------------------------------------------------------------------------------------四、实验代码/*问题:1 初始最大容量的设定2 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) {