信道容量的迭代算法的C语言程序.docxVIP

信道容量的迭代算法的C语言程序 信道容量的迭代算法的C语言程序 PAGE 信道容量的迭代算法的C语言程序 C语言程序： #include #include #include int main() { int r,s,i,j,k=0; double p[20]; %存放输入信源概率矩阵 double z[20]; double q[20][20]; %存放信道转移概率矩阵 double F[20][20]; %存放的概率分布矩阵的转置 double x,y,a; double epsilon=1e-5; %门限 double C=; %取初始迭代时的信道容量为一个较大的负数 printf(\n请输入信源符号个数: ); scanf(%d,r); printf(\n请输入信宿符号个数： ); scanf(%d,s); printf(\n请输入信道转移概率矩阵：\n\n); for(i=0;ir;i++) { for(j=0;js;j++) scanf(%lf,q[i][j]); printf(\n); } for(i=0;ir;i++) p[i]=(double)(double)r); %设初始信源分布为等概分布 do { k++; a=C; for(j=0;js;j++) { x=; for(i=0;ir;i++) x=x+(p[i])*(q[i][j]); % x 为的分母部分 if(x0) for(i=0;ir;i++) F[i][j]=(p[i])*(q[i][j])/x; % F为的概率分布矩阵的转置 else for(i=0;ir;i++) F[i][j]=; %的分母部分为0时，令F=0 } y=; for(i=0;ir;i++) { z[i]=; for(j=0;js;j++) { if(F[i][j]0) z[i]=z[i]+(q[i][j]*(log(F[i][j])/log)); %z[i]为的分子部分 } z[i]=(pow,z[i])); y=y+z[i]; %z[i]为的分母部分 } for(i=0;ir;i++) { p[i]=z[i]/y; %更新输入信源概率矩阵 } C=(log(y)/log); %求信道容量单位为“bit” } while(fabs((C-a)/C)epsilon); printf(\n迭代次数为;k=%d\n,k); %输出迭代次数 printf(\n最佳信源分布为：\n\n); for(i=0;ir;i++) { printf(%.3lf ,p[i]); %输出信源的最佳分布，保留3位小数 } printf(\n); printf(\n信道容量为:C=%.3lf bit\n\n,C); %输出信道容量，保留3位小数 } 问题1的执行结果： 问题2 的执行结果：