《信道容量的matlab程序》由会员分享，可在线阅读，更多相关《信道容量的matlab程序（2页珍藏版）》请在金锄头文库上搜索。
1、代码：functionmain()clc;p=1/21/401/4010000101/401/41/2;channel_cap(p,0.001)%Matlab实现离散信道容量的迭代算法%功能：利用迭代算法计算离散信道的容量%参数解释%C：信道容量%P：转移概率矩阵%B：中间变量矩阵%e:信道容限%X：输入概率分布%n：迭代次数functionchannel_cap(P,e)n=0;C=0;C_0=0;C_1=0;r,s=size(P);fori=1:rif(sum(P(i,:)=l)%检测概率转移矩阵是否行和为1.error(概率转移矩阵输入有误！)return;endforj=1:sif(P
2、(i,j)l)%检测概率转移矩阵是否负值或大于1error(概率转移矩阵输入有误！)return;endendendX=ones(1,r)/r;A=zeros(1,r);B=zeros(r,s);while(1)n=n+1;fori=1:rforj=1:sB(i,j)=log(P(i,j)/(X*P(:,j)+eps);endA(1,i)=exp(P(i,:)*B(i,:);endC_0=log2(X*A);C_1=log2(max(A);if(abs(C_0-C_1)e)%满足迭代终止条件停止迭代C=C_0;fprintf(迭代次数：n=%dn,n)fprintf(信道容量：C=%f比特/符号n,C)break;%满足后输出结果并退出elseX=(X.*A)/(X*A);continue;endend结果：迭代次数：n=6信道容量：C=1.321928比特/符号
