离散无记忆信道容量数值模拟

二维信道

clear all

close all

H_a=zeros(1,100);i=1

for p1_0=0.01:0.01:0.04    

    for a=0.01:0.01:1%p0 p1

       H_a(i)=-a*log2(a)-(1-a)*log2(1-a)-(-p1_0*log2(p1_0)-(1-p1_0)*log2(1-p1_0)); 

        i=i+1;

    

    

    end

    plot(H_a)

   

    hold on

    i=1;

end

 legend('0.01','0.02','0.03','0.04')

clear all

close all

H_a2=zeros(10,10);i=1;j=1;

ii=1:1:98;jj=1:1:98;

for p1_0=0.01    

    for a=0.01:0.01:1%p0 p1

        for b=0.01:0.01:1

            if a+b<1

       H_a2(i,j)=-a*log2(a)-b*log2(b)-(1-b-a)*log2(1-b-a)-a*3*p1_0*log2(p1_0)-b*3*p1_0*log2(p1_0)-(1-a-b)*3*p1_0*log2(p1_0); 

       j=j+1 ;

            end

        end

       j=1; i=i+1;

    end

%plot3(ii,jj,H_a2)

   surf(H_a2)

    hold on

    i=1;

end

 legend('0.01') 

三维信道