function f=dwvlissee(h2d,s,M);

%            f=dwvlissee(h2d,s,M);
%
%Calcule la transformée de Wigner-Ville discrète lissée par une fenêtre
%h2d. Le signal (échantillonné avec un pas tau) est répété par
%périodisation aux fenêtres latérales. L'entrée est un vecteur ligne.
%L'abscisse est le temps; l'axe des ordonnées est l'axe des fréquences. Le
%domaine des fréquences est [0, pi/tau] périodisé.  

I=dwv(s,M);
a=sum(sum(h2d));
f=(1/a)*conv2(I,h2d,'same');