problème lors d'un filtrage pass bas d'un signal numérique

[loupou]

Bonjour,

je rencontre un problème lorsque je filtre un signal numérique. Le signal à filtrer est le suivant:
stoch_etat_eq.png

lorsque je le filtre j'obtient le résultat suivant :
stoch_etat_eq_filtre.png

On peut voir que le filtrage filtre supprime une partie de l'information qu'il faut garder (amplitude à première vue) car la courbe bleu devrait alterner entre les 2 courbes en points tillés.
Voici le code que j'ai implémenté pour le filtrage du signal. D'après l'article que nous utilisons pour la modélisation, le filtrage passe bas à une fréquence de coupure de 10000 years. voici les lignes de codes issue d'un TP de traitement du signal:
code_filtre.PNG
où la variable T est le temps, et Y le signal à filtrer.

Peut être quelqu'un pourra m'aider à résoudre ce problème..

Cordialement
-----

[Gwinver]

Bonsoir.

Est-il possible de faire un zoom sur une ou plusieurs zones "bleues" denses.

[loupou]

Oui bien sur, a quelle niveau voulez vous que je zoom sur xlim ?

je vous ai rajouté quelques infos sur le signal que trace:
fonction_code.PNGvariable_code.PNGfonct_temp.PNG

[Gwinver]

Bonsoir.
L'idée est d'avoir une idée plus précise de cette zone par exemple. Les points sont-ils répartis également de part et d'autres du trait en pointillé ou y a t'il une majorité en dessous?

:

P.S. : je ne suis pas assez familier de Python pour avoir une idée précise de l'algorithme utilisé. Mais le graphe "filtré" fait penser à une sorte d'écrêtage: les partis basses des groupes de valeurs au-dessus de la moyenne et les parties hautes des groupes en dessous des valeur moyenne.

[A voir en vidéo sur Futura]

[JPL]

le filtrage passe bas à une fréquence de coupure de 10000 years.

Le correcteur orthographique a encore frappé ? 10 000 Hz je suppose ?

[Gwinver]

Il faut bien regarder les figures proposées

[Henrix]

Supprimé par moi-même

[loupou]

Bonjour,

tout d'abord merci pour vos réponses.

Gwinver

Voici le zoom en question:

Les valeur sont pondérées par une valeur aléatoire de moyenne 0 et d'ecart type 1.

Non c'est bien 10000 years, l'unité est en année non en seconde mais le raisonnement reste le même. Je veux filtrer le bruit avec les fréquences supérieur à 1/Tfiltre où Tfiltre = 10000 pour garder seulement les petites fréquences.

Je ne suis pas familier avec l'écretage...

[Gwinver]

Bonjour.

Le zoom montre bien une oscillation avec une moyenne proche du pointillé.

Il y a aussi peut-être un effet lié aux transitions rapides entre les créneaux du bas et ceux du haut.
En termes de FFT, cette transition peut correspondre à une fréquence élevée qui est "rabotée" par le filtrage.
Serait-il possible de n'effectuer le filtrage que sur la partie zoomée?

[loupou]

Oui biensur,



voici le résultat pour la même fenêtre que précédemment

[Gwinver]

Bonsoir.

J'imagine que la courbe verte est l'ancien pointillé et la courbe bleue le résultat du filtrage.

En regardant à nouveau les courbes, je me rends compte que le signal est polarisé, c'est à dire entre +9°C et +15°C environ.
Autrement dit, ce serait décomposable en signaux alternatifs ajoutés à une valeur de polarisation.
Cet effet est-il compatible avec l'outil utilisé?

Si tel est le cas, il serait possible de retrancher aux échantillons leur valeur moyenne soit environ 12°C.

voici les lignes de codes issue d'un TP de traitement du signal

Comment se présentait le lot de données de ce TP?

[loupou]

bonjour,

Je n'ai pas compris qu'est ce qui doit être compatible avec l'outils utilisé ?

Le lot de données est issue de l'itération temporelle sur une équa diff, le bruit est déterminer par une variable aléatoire, les formules sont présentes sur les lignes de codes que j'ai mis plus haut.

[Gwinver]

Bonsoir.

En règle générale, les signaux analysés sont alternatifs, la valeur moyenne est de zéro.
Ici, il s'agit d'un signal polarisé, la valeur moyenne (estimée d'après le graphique) est de l'ordre de 12°C. La valeur minimale est environ 8°C.
Je me demande donc si la méthode de filtrage utilisée est compatible avec cette "polarisation".
Quelle était la situation lors du TP dans lequel cette méthode a été utilisée?

Il serait possible d'essayer de recentrer les valeurs sur zéro en retranchant leur valeur moyenne.