%用迭代法計算信道容量;
%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