信息量及信道容量的计算.doc

下载文档

141
0
约1.8千字
约 3页
2018-06-25 发布于河南
举报
版权申诉
保障服务

信息量及信道容量的计算.doc

1、本文档共3页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

信息量及信道容量的计算

#include iostream #include cmath #include string using namespace std; int main() { int i,j,k,m,n; char r; char A=Y,B=N; double x[20],p[12][12],q[12][12],y[20]; cout输入信源x的个数N=; cinn; cout输入所需信源概率:endl; for(i=1;i=n;i++) { coutx(i)=; cinx[i-1]; } cout输入信道矩阵:endl; for(i=1;i=n;i++) { for(k=1;k=n;k++) { cout请输入P(yk/xi)=; cinp[i-1][k-1]; } } //for(;;); // { cout1信源熵 2条件熵 3联合熵 4交互熵endl; cout选择所要计算的熵值:endl; cinm; if(m==1) { double H=0,h; for(int j=1;j=n;j++) { h=-x[j-1]*log10(x[j-1])/log10(2); H=H+h; } cout信源熵为:Hendl; } else if(m==2) { double H1=0,h1=0 ,H2=0,h2=0; for(i=1;i=n;i++) { for(j=1;j=n;j++) { q[i-1][j-1]=p[i-1][j-1]*x[i-1]; //cout联合概率yixjq[i-1][j-1] endl; h1=-q[i-1][j-1]*log10(p[i-1][j-1])/log10(2); H1=H1+h1; } } y[j-1]=y[j-1]+q[i-1][j-1]; cout信源yiy[j-1] endl; h2=-q[i-1][j-1]*log10(q[i-1][j-1]/y[j-1])/log10(2); H2=H2+h2; } cout条件熵H(Y/X)为:H1endl; cout条件熵H(x/y)为:H2endl; } else if(m==3) { double H3=0,h3=0; for(i=1;i=n;i++) { for(j=1;j=n;j++) { q[i-1][j-1]=p[i-1][j-1]*x[i-1]; H3=-q[i-1][j-1]*log10(q[i-1][j-1])/log10(2); H3=h3+H3; } } cout联合熵H(xy)为:H3endl; } else if(m==4) { double H4=0,h4=0; for(i=1;i=n;i++) { for(j=1;j=n;j++) { q[i-1][j-1]=p[i-1][j-1]*x[i-1]; H4=-q[i-1][j-1]*log10(q[i-1][j-1]/x[i-1])/log10(2); H4=h4+H4; } } cout交互熵I(x;y)=I(y;x)为:H4endl; } // cout是否继续计算，继续请输入Y，退出请输入N!endl; // cout请输入Y或者N，进行选择: ; // cinr; // if (r==A) // {continue;} // if (r==B) // {break;} //} return 0; }