Matlab绘制二进制熵函数曲线实验

一：实验目的

• 熟悉工作环境及Matlab软件
• 掌握绘图函数的运用
• 理解熵函数表达式及其性质

 二：实验内容

1.用 Matlab 软件（matlab2010a版本）绘制二进熵函数曲线。 

 

 

• 二元信源

 

2．绘制三元信源的熵

 三元信源

 

 

3．绘制平均互信息量图形 

 

对于二元对称信道的输入概率空间为

 

 

 

三：实验代码 

1.

p=0.00001:0.001:1;
h=-p.*log2(p)-(1-p).*log2(1-p);
plot(p,h);
title('二进熵函数曲线');
ylabel('H(p,1-p)')

 2.

p=linspace(eps,1-eps,100);

q=linspace(eps,1-eps,100);

[P,Q]=meshgrid(p,q);

P_Q=P+Q;

for n=1:100

    for m=1:100

        if P_Q(n,m)>=1

            Q(n,m)=nan;

        end

    end

end

H=-P.*log2(P)-Q.*log2(Q)-(1-P-Q).*log2(1-P-Q);

mesh(P,Q,H)

title('三维熵函数的图形')

3.

[w,p]=meshgrid(0.00001:0.001:1);
h=-(w.*(1-p)+(1-w).*p).*log2(w.*(1-p)+(1-w).*p)-(w.*p+(1-w).*(1-p)).*log2(w.*p+(1-w).*(1-p))+(p.*log2(p)+(1-p).*log2(1-p))
meshz(w,p,h)
title('平均互信息量')
ylabel('H')

四：实验结果