实验一

代码说明：

说明：  一般信道容量迭代算法 （1）初始化信源分布：P^((0))=(p_(1 ),p_(2 ),…,p_(i ),…,p_r)，置迭代计数器k=0，设信道容量相对误差门限为δ，δ>0，可设C(0)=-∞； （2）φ_ji^((k))=(p_ij p_i^((k)))/(∑p_ij p_i^((k)) ) i=1,…,r;j=1,…,s； （3）p_i^((k+1) )=exp?[∑_j?〖p_ij lnφ_ji^((k)) 〗]/(∑_i?〖{exp?[∑_j?〖p_ij lnφ_ji^((k)) 〗]}〗) i=1,…,r； （4）C^((k+1))=ln∑_i?〖{exp?[∑_j?〖p_ij lnφ_ji^((k)) 〗]}〗； （5）如果， (|C^((k+1) )-C^((k)) |)/C^((k+1)) ≤δ， 转向（7）； （6）否则，置迭代符合 k+1→k，转向（2）； （7）输出p_i^((k+1) )和C^((k+1))； （8）停止。(General channel capacity iterative algorithm)