h.m

实现（2

%  H= -plogp-(1-p)log(1-p)

x=1/100:1/100:(1-1/1000);

Z=-x'*log2(x)-(1-x)'*log2(1-x);
m=-x.*log2(x)-(1-x).*log2(1-x);

y=m.^2;
[X,Y]=meshgrid(x,y);
surf(Y,Z,X);