stepsize.m

%stepsize源程序
function q=stepsize(n,p)
%通过分解创造步骤大小顺序的一个子带量化和子带(水平的、垂直的、对角线的和近似值子带的最后分解)

if length(p)==2%不明确量化
    q=[];
    qn=2^(8-p(2)+n)*(1+p(1)/2^11);
    for k=1:n
        qk=2^-k*qn;
        q=[q (2*qk) (2*qk)(4*qk)];
    end
    q=[q qk];
else%明确量化
    q=p;
end

q=round(q*100)/100;%大约1/100分配
if any(100*q>65535)
    error('The quantizing steps are not UNIT16 representable.')
end
if any(q==0)
    error('A quantizing step of 0 is not allowed.');
end