2entropy.m

现代通信系统（Proakis）第二版

function h=bientrop(p)
% h=bientrop(p) Returns the binary entropy function of the components 
%		of the vevtor p.
n=length(p);
for i=1:n
p1=[p(i) 1-p(i)];
h(i)=entropy(p1);
end