forecastupdate.m

此为基于二代提升小波的阈值压缩算法源码

%**************************************************************************
% 函数功能：计算提升方案中的一个更新或预测环节。
% 输入参数：x:输入向量；F：劳伦多项式系数向量；DF：该多项式单项最高次数；S：要输出向量y的长度；Lstype非空进行整数提升
% 返回值：预测或更新环节得到的向量y
%**************************************************************************
function y = forecastupdate(x,F,DF,S,LStype)

lF = length(F);
sx = size(x);
y  = zeros(sx);
lx = length(x);
for j=1:lF
    t = F(j)*x;
    k = DF-j+1; 
    if     k>0 , t = t(1+k:end); t(end+k)= 0;
    elseif k<0 , t(1-k:end) = t(1:end+k); t(1:-k) = 0;
    end
    y = y + t(1:lx);
end
d = S-lx;
if d>0 , y(end+d) = 0; elseif d<0 , y = y(1:S); end
if ~isempty(LStype) , y = fix(y); end