《实验二---一般信道容量迭代算法(共3页).doc》由会员分享，可在线阅读，更多相关《实验二---一般信道容量迭代算法(共3页).doc（3页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上实验二 一般信道容量迭代算法1 实验目的掌握一般离散信道的迭代运算方法。2 实验要求1） 理解和掌握信道容量的概念和物理意义2） 理解一般离散信道容量的迭代算法3） 采用Matlab编程实现迭代算法4） 认真填写试验报告3算法步骤初始化信源分布（一般初始化为均匀分布），置迭代计数器k=0，设信道容量相对误差门限为，0，可设； 如果，转向；置迭代序号,转向；输出和的结果；停止。4代码P=input(转移概率矩阵P=)e=input(迭代精度e=)r,s=size(P);n=0;C=0;C_k=0;C_k1=0;X=ones(1,r)/r;A=zeros(1,r);B=z
2、eros(r,s);%初始化各变量while(1) n=n+1; for i=1:r for j=1:s B(i,j)=log(P(i,j)/(X*P(:,j)+eps); if P(i,j)=0 B(i,j)=0; else end end A(1,i)=exp(P(i,:)*B(i,:); end C_k=log2(X*A); C_k1=log2(max(A); if (abs(C_0-C_1)e) C=C_0; %检验是否满足迭代精度 fprintf(迭代次数n=n,n)%满足则输出相关的被求量 fprintf(信道容量:C=n,C) disp(最佳分布P(x)=);disp(X); b
3、reak; else %不满足则继续迭代 X=(X.*A)/(X*A); continue; endend5.计算下列信道的信道容量例一：P = 0.9800 0.0200 0.0500 0.950e = 1.0000e-005迭代次数: n=5信道容量: C=0.比特/符号最佳分布P(x)= 0.5129 0.4871例二：P = 0.6000 0.4000 0.0100 0.9900e =1.0000e-005迭代次数: n=17信道容量: C=0.比特/符号最佳分布P(x)= 0.4238 0.5762例三：P = 0.7900 0.1600 0.0500 0.0500 0.1500 0.8000e =1.0000e-005迭代次数: n=6信道容量: C=0.比特/符号最佳分布P(x)=0.5009 0.4991思考题：迭代精度指的是什么？它对计算结果的影响？迭代精度小数点的位数与迭代次数相同，迭代精度越小，信道容量越大。专心-专注-专业