entropy.m

来自「覆盖了全书的信号与线性系统、随机过程、模拟调制、模拟数字转换、基带数字传输、带限」· M 代码 · 共 11 行

function h=entropy(p)
%		H=ENTROPY(P) returns the entropy function of 
%		the probability vector p.
if length(find(p<0))~=0,
  error('Not a prob. vector, negative component(s)')
end
if abs(sum(p)-1)>10e-10,
  error('Not a prob. vector, components do not add up to 1')
end
h=sum(-p.*log2(p));