irefinement.m

来自「SPIHT的一个Matlab程序」· M 代码 · 共 23 行

function dect=irefinement(dect,T,sig_ele)

global iii
stream=[];
ii=1;
while ii<=length(sig_ele)
    stream=[stream input_bit];
    ii=ii+1;
end
%dect(sig_ele), used to fix bug
%stream , To check
for rr=1:length(stream)
    if (dect(sig_ele(rr)) >0),s=0;
    else s=1;
    end
    if stream(rr)=='1'
        dect(sig_ele(rr))=abs(dect(sig_ele(rr)))+T/2; % ??? % T/4 if refinement added first pass
    elseif stream(rr)=='0'
        dect(sig_ele(rr))=abs(dect(sig_ele(rr)))-T/2; % ??? % T/4 if refinement added first pass 
    end
    if s==1,dect(sig_ele(rr))=-dect(sig_ele(rr));end
end