计算信源熵的 MATLAB源程序
function H=entropy(P,r)
if (length(find(P<=0))~=0)
error('Not a prob.vector,negative component');
end
if(abs(sum(P)-1)>10e-10)
error('Not a prob.vector,component do not add up to 1');
end
H=(sum(-P.*log2(P)))/(log2(r)+eps);
香农编码的 MATLAB源程序
function [s,L,q]=shannon(p)
%if(length(find(p<=0))~=0)
% error('Not a prob.vector,negative component');
%end
%if(abs(sum(p)-1)>10e-10)
% error('Not a prob.vector,component do not add up to 1')
%end
n=length(p);
x=1:n;
[p,x]=array(p,x);
l=ceil(-log2(p));
P(1)=0;
n=length(p);
for i=2:n
P(i)=P(i-1)+p(i-1);
end
for i=1:n
for j=1:l(i)
temp(i,j)=floor(P(i)*2);
P(i)=P(i)*2-temp(i,j);
end
end
s=[];
for i=1:n
for j=1:l(i)
t=temp(i,j);
评论0
最新资源