信道容量编码

《信道容量编码》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、#include#include#includevoidmain(){registerinti,j;registerintk;intr,s;/*信道矩阵的行数和列数*/float*p=NULL;float*q=NULL;float*d=NULL;float*a=NULL;float**p_ji=NULL;/*二维数组p_yx存放信道矩阵*/floatC,u,max,t,IU;/*C为信道容量*/printf("矩阵行数:r=");/*输入行数和列数*/scanf("%d",

2、&r);printf("矩阵列数:s=");scanf("%d",&s);p=(float*)calloc(r,sizeof(float));/*初始化需要用到的数组*/q=(float*)calloc(r,sizeof(float));d=(float*)calloc(r,sizeof(float));a=(float*)calloc(r,sizeof(float));p_ji=(float**)calloc(r,sizeof(float));for(i=0;i

3、sizeof(float));printf("Pleasematrixofchannel:");/*输入信道矩阵*/for(i=0;i

4、

5、t>=1.000001){pri

6、ntf("Error!");exit(-1);}}for(i=0;i

7、j]/q[j]);}a[i]=exp(d[i]);}u=0.0;for(i=0;i=max)max=a[i];}IU=log(max)/log(2.0);/*计算IU*/k=1;while((IU-C)>=10e-6)/*迭代计算*/{/*计算最佳输入分布p[i]*/for(i=0;i

8、出信源的概率分布q[j]*/{q[j]=0;for(i=0;i

9、;/*计算IL*/max=a[0];for(i=0;i=max)max=a[i];}IU=log(max)/log(2.0);/*计算IU*/k=k+1;}printf("信道矩阵:");/*输出信道矩阵*/for(i=0;i

10、]*/for(i=0;i