Filtre Médian : pblm sur valeurs initiales

Fermé
OReLo - 18 mai 2010 à 11:50
Bonjour à tous,

Je cherche à créer un filtre médian 50 sous matlab afin de récupérer seulement l'enveloppe d'un signal bruité. jusque là tout va bien...

Mais en y regardant dans le détail, l'enveloppe que j'obtiens est correct mais la première valeur ne me convient pas du tout.

Auriez-vous une explication...??!

Un exemple simple :

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%TEST FILTRE MEDIAN

% Génération d'un sinus et ajout d'un bruit blanc Gaussien

fe = 8000; % Fréquence d'échantillonnage
N = 8192; % Nombre de points de la séquence

% Axe des temps
t = (1:N)/fe;

% Génération du sinus
f0 = 1200;
sinus = sin(2*pi*f0*t);

% Génération du bruit
sigma = 0.25; % variance du bruit
moy = 0; % moyenne
bruit = moy + sigma*randn(1,N);

signal = sinus + bruit; % les 2 vecteurs sont de même longueur

figure(1);
plot(t,signal)

signal2 = medfilt1(signal,50)

figure(2);
plot(t,signal2)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Merci beaucoup du petit coup de pouce, pied, boule et main que vous pourrez m'apporter...!