信息论实验二报告(6)

for i=1:r

Pb(j,1 )=Pb(j, 1)+Pa(i)*Pba(i,j,1 );

%求信源的输出概率矩阵，即Pb(j,1)

end

end

for i=1:r

temp(i)=0;

for j=1:s

temp(i)=temp(i) + Pb(j,1)*exp(S(m) *d(i,j));

% temp为临时项，求Pha(i,j,2)时表达式的分母