我的信道容量迭代.asv

来自「信息理论与编码——关于信道容量迭代的matlab仿真程序」· ASV 代码 · 共 29 行

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);