Discussion :

[lecteur1001]

Bonjour,

je possède des résultats d'expérience dans un fichier.mat
Ces données sont quelque peu bruitées, je souhaite donc les filtrer à l'aide d'un filtre du 1er ordre avec une constante de temps de 10s.
Sous Simulink, rien de plus bête : j'utilise la fonction 'from Worskpace' avec comme paramètre

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

[time, donnees]

Je relie ce bloc à un bloc 'transfer fcn' dont les arguments sont '1' pour le numérateur et '[10 1]' pour le dénominateur. Et en sortie j'obtiens le signal filtré.
Ma question est bête mais existe-t-il une fonction en code MATLAB qui soit équivalente à ce modèle Simulink ? Je ne m'en sors pas avec la function 'tf'. Suis-je obligé de passer par une intégration avec une fonction ode par exemple ?

Merci de vos réponses

[duf42]

Bonjour,

Je sais pas si ca réponds à ta question mais il y a un exemple de filtre simple dans l'aide de la fonction FILTER.

Duf

[lecteur1001]

Si ça pourrait mais le problème est que la fonction filter définit un filtre discret. Moi je recherche quelque chose de continu.

J'aimerai vraiment savoir quel code se cache derrière le bloc 'Transfer fcn'...

[lecteur1001]

Bon ben... j'ai opté pour la solution de facilité "pas très joli" : j'ai créé un modèle Simulink auquel je fais appel grâce à la fonction 'sim'.
Si quelqu'un a une solution plus élégante je suis preneur

[tubaas]

bonjour
je veux bien essayé de t'aider mais je ne connais pas Simulink.

Ces données sont quelque peu bruitées, je souhaite donc les filtrer à l'aide d'un filtre du 1er ordre avec une constante de temps de 10s.

Je suppose que tu veux appliquer un filtre passe-bas du premier ordre. Cependant 10s me paraît énorme comme constante de temps.

Si ça pourrait mais le problème est que la fonction filter définit un filtre discret. Moi je recherche quelque chose de continu.

Quelque chose de continu... c'est à dire un filtre analogique ? Parce que Simulink applique nécessairement un filtre numérique, et tes données sont discrétisées si je ne m'abuse.
En fait il te faut la fonction de transfert de ton filtre (des vecteurs avec les coefficients du numérateurs et du dénominateur, a et b dans la doc de filter).
Pour les construire si tu ne les connais pas tu as par exemple les fonctions CHEBY1, CHEBY2, BUTTER,... suivant le type de filtre que tu souhaites utilisé mais il faut avoir la Signal Processing Toolbox.
Sinon ça doit pouvoir se trouver si tu précises bien le type de filtre, son ordre, sa fréquence de coupure, éventuellement gain et facteur de qualité...

[lecteur1001]

bonjour
je veux bien essayé de t'aider mais je ne connais pas Simulink.

Je suppose que tu veux appliquer un filtre passe-bas du premier ordre. Cependant 10s me paraît énorme comme constante de temps.

Ben tout dépend de la dynamique des données que j'ai à filtrer. Dans le cas présent, 10s est tout à fait adapté.

Quelque chose de continu... c'est à dire un filtre analogique ? Parce que Simulink applique nécessairement un filtre numérique.

Oui et non... C'est juste qu'avec un filtre discret, tu spécifies la période d'échantillonnage dans le filtre alors qu'avec un filtre "continue", la période d'échantillonnage est celle appliquée au modèle Simulink.

En fait il te faut la fonction de transfert de ton filtre (des vecteurs avec les coefficients du numérateurs et du dénominateur, a et b dans la doc de filter).

Je ne les connais pas... En revanche, je les connais dans le domaine de Laplace. D'où ma volonté initiale de passer par la fonction 'tf'. Du coup, je ne comprends pas quel est l'intérêt de cette fonction finalement.

Pour les construire si tu ne les connais pas tu as par exemple les fonctions CHEBY1, CHEBY2, BUTTER,... suivant le type de filtre que tu souhaites utilisé mais il faut avoir la Signal Processing Toolbox.
Sinon ça doit pouvoir se trouver si tu précises bien le type de filtre, son ordre, sa fréquence de coupure, éventuellement gain et facteur de qualité...

Sans faire preuve de mauvaise volonté, je t'assure que j'ai déjà cherché une solution dans ce sens et je n'y ai rien trouvé. Après, j'ai peut-être mal cherché...
Ce que je ne comprends pas finalement, c'est pourquoi je ne peux pas faire un truc aussi simple via Simulink en code ? Dans Simulink, tu mets ton signal à l'entrée, ton filtre du 1er ordre (dans le domaine de Laplace) et tu obtiens le signal de sortie dans l'intervalle de temps désiré. Là, j'ai l'impression d'être obligé de passer par une fonction de transfert discrète, de chercher les coefficients correspondants...