statmoments.m

来自「Matlab數位影音處理 影像方面的處理」· M 代码 · 共 24 行

function [v,unv]=statmoments(p,n)
lp=length(p);
if (lp~=256) & (lp~=65536)
    error('P must be a 256- or 65536- element vector.');
end
g=lp-1;
p=p/sum(p);
p=p(:);
z=0:g;
z=z./g;
m=z*p;
z=z-m;
v=zeros(1,n);
v(1)=m;
for j=2:n
    v(j)=(z.^j)*p;
end
if nargout >1
    unv=zeros(1,n);
    unv(1)=m.*g;
    for j=2:n
        unv(j)=((z*g).^j)*p;
    end
end