嫡计算， function h = entropys(p) % H = ENTROPY(P) returns the entropy function of % the probability vector p; if length(find(p10e-10, error('ERROR Not a prob. vector, component do not