function d=D4(s);

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

N=length(s);
G=zeros(1,8);
G(1)= -.230378;
G(2)= .714846; 
G(3)= -.630881;
G(4)= -.027984;
G(5)= .187035; 
G(6)= .030841;
G(7)= -.032883; 
G(8)= -.010597; 
S=conv(s,G);
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
	d=u;
else
	d=[u 0];
end;