function d=D2(s);

%                         d=D2(s);
%
% Algorithme de filtrage-décimation correspondant à l'analyse 
% multirésolution d'Ingrid Daubechies à 4 paramètres (operateur R_\varphi)
% agissant sur un signal digital ligne s.

N=length(s);
G=zeros(1,4);
G(1)=-.48296;
G(2)=.83652;
G(3)=-.22414;
G(4)=-.12941;
S=conv(G,s);
m=fix((N+3)/2);
u=ones(1,m); 
for l=1:m, 
    u(l)=S(2*l-1);
end
if mod(m,2) == 0
    d=u;
else
    d=[u,0]; 
end;            