MATLAB迭代法计算信道容量

%用迭代法计算信道容量;
%inputp最佳输入分布（不唯一）;
%k迭代步数;
%c2信道容量;
function [c2 inputp  k]=C（Peps）
if nargin~=2
    eps=1e-10;
end
[rs]=size（P）;
back=[];%后验概率矩阵;
inputp=ones（r1）/r;%初始化输入概率分布;
c1=1;%信道容量初始化;
%..........................................................................
for k=0:1000
    for i=1:r
        for j=1:s
            lamda=sum（inputp.*P（:j））;
            back（ji）=inputp（i）*P（ij