function r=R2(s);

%                        r=R2(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);
H=zeros(1,4);
H(1)=-.12941;
H(2)=.22414;
H(3)=.83652;
H(4)=.48296;
S=conv(H,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
	r=u;
else
	r=[u 0];
end;