Le lissage consiste à retirer le bruit d'un signal à variation lente. L'algorithme Savitsky-Gloay (du nom de ses concepteurs) consiste à effectuer une modélisation de type polynomiale (degré 2 essentiellement) sur un certains nombre de points adjacents au point considéré. Ensuite, il suffit de recalculer la variable dépendente (ordonnée la plus part du temps) par rapport au modèle.
Donc en fait, il faut et il suffit de maîtriser la modélisation (et donc les calculs matriciels) et d'écrire une routine itérative. Là je crois pouvoir dire qu'on parle de joli code, et je ne suis pas sûr qu'il soit dispo tel quel...
Bon, l'algoritthme date un peu, et de nouvelles technologies, non plus de lissage mais de débruitage, donnent d'excellents résultats; la méthode par Transformée Rapide de Fourrier est elle même suplantée par les Ondelettes (calculs fractals), mais çà c'est encore une autre histoire ;-)