function r=R4(s);

%                        r=R4(s); 
%
% Algorithme de filtrage-décimation correspondant à l'analyse 
% multirésolution de Daubechies à 8 paramètres (opérateur R_\varphi)
% agissant sur un signal digital ligne s. 

N=length(s);
H=zeros(1,8);
H(1)= -.010597;
H(2)= .032883; 
H(3)= .030841;
H(4)= -.187035;
H(5)= -.027984; 
H(6)= .630881; 
H(7)= .714846; 
H(8)= .230378; 
S=conv(s,H);
m=fix((N+7)/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;