信道容量迭代算法验证
一、实验目的
1.熟悉信道容量的迭代算法。
2.学习如何将复杂的公式转化成程序。
3.熟悉c++编译码过程
二、实验要求
1,学习vs2012软件编程和调试方法;
2,输入:任意一个转移概率矩阵,包括信源符号个数,信宿符号个 数,信道转移概率,在程序运行时从键盘输入。
3,输出信道容量C 迭代次数,最佳信源分布。
三、实验仪器,设备
笔记本电脑,vs2012软件
四、实验原理
已知:信源符号个数r,新宿符号份额数s、信道转移概率矩阵p;
输入:任意的一个信道转移概率矩阵。r、s、p在运行时从键盘输入。
输出:最佳信源分布P‘,信道容量C。
4.1信道容量的含义
一个信道,若给定输入分布Q(x) ,和转移概率分布P(y|x) ,则输出概率分布为:Θ(y)=∑xQ(x)P(y|x)
对DMC信道有Θ(y)=∑xQ(x)P(y|x)
当研究联合 空间XY时,则信道输出和输入之间的平均互信息量为:
=H(Y)−H(Y|X)
=H(X)−H(X|Y)
从上述公式中,当信道输入为X,接收端得到的平均信息量 I(X;Y)是信道输入概率和信道转移概率的函数。一般干扰越大,传输的信息量就越小。现在研究信道哦给定,即 {P(y|x)}给定的情况下,如何选择输入分布 QkQk使传输的信息量 I(X;Y)达到最大。
信道容量的定义:
C=maxI(X;Y)
C=maxI(X;Y)
即C为改变输入分布时,使每个符号所能含有的平均互信息量的最大值。相应的输入分布称为最佳输入分布信道容量的迭代算法。
五、实验方案与步骤
5.1信道容量的迭代算法的步骤如下:
- 取出时分布 r(0),置k=0。
- 由公式
计算 。
4.由公式
计算。
5. 若 |C(k+1)−C(k)|≤δ,进入第7步。
6. 令 k=k+1, 转向第二步。
7. 程序结束输出最佳输入分布及C(k+1)。
5.2程序设计流程图:
5.3关键程序代码
#include<stdio.h>
#include<math.h>
#include<stdlib.h>int main()
{
int r, s, i, j, k = 0;
double p[20]; //存放输入信源概率矩阵
double z[20];
double q[20][20]; //存放信道转移概率矩阵
double F[20][20]; //存放(|)iipxy的概率分布矩阵的转置
double x, y, a;
double epsilon = 1e-5; //门限
double C = -1000.0; //取初始迭代时的信道容量为一个较大的负数
printMessage();
printf("!-------请输入信源符号个数:---------!\n");
scanf("%d", &r);
printf("!-------请输入信宿符号个数:---------!\n");
scanf("%d", &s);
printf("!-------请输入信道转移概率矩阵---------!\n\n");
for (i = 0; i < r; i++)
{
for (j = 0; j < s; j++)
scanf("%lf", &q[i][j]);
printf("\n");
}
for (i = 0; i < r; i++)
p[i] = (double)(1.0 / (double)r);
//设初始信源分布为等概分布
do
{
k++;
a = C;
for (j = 0; j < s; j++)
{
x = 0.0;
for (i = 0; i < r; i++)
x = x + (p[i])*(q[i][j]);// x 为(|)iipxy的分母部分
if (x > 0)
for (i = 0; i < r; i++)
F[i][j] = (p[i])*(q[i][j]) / x;// F为(|)iipxy的概率分布矩阵的转置
else
for (i = 0; i < r; i++)
F[i][j] = 0.0;//(|)iipxy的分母部分为0时,令F=0
}
y = 0.0;
for (i = 0; i < r; i++)
{
z[i] = 0.0;
for (j = 0; j < s; j++)
{
if (F[i][j] > 0)
z[i] = z[i] + (q[i][j] * (log(F[i][j]) / log(2.0)));//z[i]为ip的分子部分
}
z[i] = (pow(2.0, z[i]));
y = y + z[i];//z[i]为ip的分母部分
}
for (i = 0; i < r; i++)
{
p[i] = z[i] / y;//更新输入信源概率矩阵
}
C = (log(y) / log(2.0));//求信道容量单位为“bit”
} while (fabs((C - a) / C) > epsilon);
printf("!----------迭代次数为:k=%d --------!\n", k);//输出迭代次数
printf("!----------最佳信源分布为:--------!\n\n");
for (i = 0; i < r; i++)
{
printf("%.3lf", p[i]); //输出信源的最佳分布,保留3位小数
}
printf("\n");
printf("!----------信道容量为:C=%.3lf bit --------!\n\n", C); //输出信道容量,保留3位小数
getchar();
getchar();
}
六、实验数据与结果分析
6.1实验结果图,如下图所示:
6.2实验结果分析:
根据图片的结果信息得到了信源迭代次数的最佳信源分布,并且计算了出了最优的信道容量。经过验算得到结果正确,代码编译成功。
6.3实验过程中的问题及其解决
1,首先,不了解信道迭代算法,最后回到书中,复习了下才明白怎么回事。
2,然后在设计算法的时候,设计的流程结构有问题,一直算不出来最佳信源分布,最后经过多次修改才得到结果。
3,运行的时候,输入完函数结果没出来就自动结束进程,一直以为是代码的原因,最后发现是因为没有延时函数,其实代码结果早就出来了。最后加了一个延时函数,得到结果。
七、实验总结
7.1、实验结论;
验证了信道容量迭代算法的正确性,解决了将公式变成代码的问题。
7.2、心得体会;
1,根据本次实验,我学习到vs2012软件编程和调试方法;为今后的实验打下了坚实的基础。
2,了解的信道容量迭代函数的原理,更加明白信道容量的具体意义。
3,知道如何将公式变成具体的代码,为以后的工作,学习打下了坚实的基础。
4.懂得了在面对困难的时候不放弃不抛弃,总会看见成功的。
7.3、对本实验方法、手段及过程等的改进建议。
对本次实验是建立在迭代的基础上的,首先要深刻理解信道容量的迭代算法,其次是,将复杂的公式转化成为程序。通过本次实验我了解信道容量的迭代算法的基本知识,学会的如何将复杂的公式用语言描述出来,掌握了高级语言数值计算过程的程序设计和调试的相关知识,对于将实际问题公式转换成程序的能力有所提高。