一般信道容量的计算matlab,DMC信道容量迭代计算的matlab实现

DMC 信道容量迭代计算的matlab 实现--通信与信息系统

一、用了matlab 实现DMC 容量迭代的算法如下：

1. 初始化信源分布：.0det a 10,1,0,1)(>>=?==，选置，，k r i r

P k i 一般我选deta=0.000001。

2. }{,)()()()(k ij i

ji

k i ji k i k ij t p p p p t 得反向转移概率矩阵由式∑=。 3.

()()()()(){}111]log exp[]

log exp[+++==∑∑∑k i k i j ij k ji j

ij k ji k i p P t p

t p p 计算由式。

4.

()()()()()()。C t p t P I C k r i s j k ij ji k k k 10011log exp log ,+==++????????????????==∑∑计算由式

5. 若a C C C k k k det )1()

()1(>-++，则k=k+1，转第2步

6.输出迭代次数k 和()1+k C 和1+k P ，终止。

二、了解了信道容量的定义和DMC 信道容量迭代计算方法，我用了matlab 来进行编程进行迭代计算得出信道容量。不足之处在于每迭代一次就输出一次迭代次数直到最后一次迭代。

1)输入：输入信源个数、信宿个数和信道容量的精度，程序能任意生成随机的

信道转移概率矩阵，也可以自己输入信道转移矩阵。

2)输出：输出最佳信源分布和信道容量。将附件里的dmc.m 文件直接run 运行

可以自主输入信道转移概率矩阵，按照程序中提示将那两句代替判

断输入矩阵是否正确的那部分，dmc1.m 运行可以随机生成信道转移

概率矩阵。

三、检验程序

之一：

输入信源个数：2

输入信宿个数：3

输入信道容量的精度： 0.000001

输入信道转移矩阵P ：[0.5000 0.3000 0.2000;0.3000 0.5000 0.2000] 之二：P ：[0.4000 0.3000 0.2000;0.3000 0.5000 0.2000] 之三：P ：[1.1000 -0.3000 0.2000;0.3000 0.5000 0.2000] 之四：P:[0.6 0.4;0.01 0.99]

之五：自动生成信道转移矩阵