信道容量編碼

《信道容量編碼》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫(kù)。

1、#include#include#includevoidmain(){registerinti,j;registerintk;intr,s;/*信道矩陣的行數(shù)和列數(shù)*/float*p=NULL;float*q=NULL;float*d=NULL;float*a=NULL;float**p_ji=NULL;/*二維數(shù)組p_yx存放信道矩陣*/floatC,u,max,t,IU;/*C為信道容量*/printf("矩陣行數(shù):r=");/*輸入行數(shù)和列數(shù)*/scanf("%d",

2、&r);printf("矩陣列數(shù):s=");scanf("%d",&s);p=(float*)calloc(r,sizeof(float));/*初始化需要用到的數(shù)組*/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);/*計(jì)算IU*/k=1;while((IU-C)>=10e-6)/*迭代計(jì)算*/{/*計(jì)算最佳輸入分布p[i]*/for(i=0;i

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

9、;/*計(jì)算IL*/max=a[0];for(i=0;i=max)max=a[i];}IU=log(max)/log(2.0);/*計(jì)算IU*/k=k+1;}printf("信道矩陣:");/*輸出信道矩陣*/for(i=0;i

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