Discussion :

[franchouze]

Bonjour,

Quelqu'un aurait il un code (simple) permettant d'effectuer l'estimation des coefficients AR du filtre correspondant à une DSP ?

Je n'ai pas la signal processing toolbox de matlab...

Merci pour votre aide.

[mihaispr]

Bonjour,

Quelqu'un aurait il un code (simple) permettant d'effectuer l'estimation des coefficients AR du filtre correspondant à une DSP ?

Je n'ai pas la signal processing toolbox de matlab...

Merci pour votre aide.

Ok.

Premierement tu dois chercher sur google des liens qui puisse t'aider:

http://www.lmgtfy.com/?q=ar+coeffici...+matlab+source

Apres ca tu dois installer signal processing toolbox pour matlab.

http://www.mathworks.com/products/signal/tryit.html (tu peux l'utiliser seulement pour 15 jours)


Tu le mets dans le dossier toolbox du logiciel matlab.

C'est tout ca.



Si tu peux m'expliquer que representent les coefficients AR , que tu desires faire en detail ca serait mieux.


Bonne chance,

Michel

[franchouze]

Bonjour mihaispr,

C'est vrai que je n'étais pas très clair...
AR signifie autorégressive.

Mon but est de calculer un bruit rose à partir d'une DSP. J'ai trouvé des méthodes rapides fréquentielles mais la DSP du bruit obtenue est très bruité.
Pour ceux que ca interesse :
http://www.dsprelated.com/showmessage/104231/1.php (post de Martin Gerhold).

La méthode AR est une méthode permettant de faire cela (entre autre j'imagine) : à partir uniquelment d'une densité spectrale de puissance (et donc d'une autocorrelation de mon signal), elle permet de calculer un bruit ayant la même DSP que celle d'origine. Il faut pour cela calculer un certain nombre de coefficients.Pour les calculer le lien ci dessous parle de l'algorithme de Durbin. Après quelques lectures de diverses documents... je ne comprends pas vraiment comment il marche et comment à partir de quelques coefficients on arrive à reproduire un signal avec les bonnes propriétés...

http://publication.lal.in2p3.fr/2002...L/node258.html

j'aurais donc juste voulu une routine qui permettait de calculer ces coefficients et qui n'utilise pas la signal processing toolbox. 15 jours d'utilisation de cette toolbox, cela risque d'être juste pour ce que j'ai à faire.

[mihaispr]

J'ai vois la formule de la densite spectrale pour un signal mais je n'ai pas compris une chose.

http://en.wikipedia.org/wiki/Spectral_density

Tu dis:

calculer un bruit rose

Le terme rose a quoi fais reference?

Si tu veux m'expliquer un peu.

Je crois que tu as trouve deja des liens dont on parle de l'algorithme Durbin en cherchant sur google.

http://ccrma.stanford.edu/~jos/latti...algorithm.html

Bonne chance,

Michel

[franchouze]

Le bruit le plus simple est appelé bruit blanc : la DSP est constante sur toutes les fréquences : dans ton signal, toutes les fréquences apparaissent de façon equiprobable. Donc sur ton spectre, les fréquences allant de 0 à FréquenceEchantillonnage/2 ont la même amplitude.

Un bruit rose est un bruit qui n'est pas à DSP constante. Certaines fréquence sont majoritaires. Par exemple on peut imaginer vouloir générer un bruit allant de FréquenceEchantillonnage/4 à FréquenceEchantillonnage/3.

Effectivement j'avais déjà vu le lien que tu fournis. Il est à l'origine de mon mal de crâne d'hier ;-).

Ja vais approfondir la chose.

[mihaispr]

Je suis desole que je ne peux t'aider plus.

Mais essaye a installer Signal Processing Toolbox meme si tu l'installe pour 15 jours.


Bonne chance,

Michel

[robotron]

Pour répondre à la question originale, il te suffit d'estimer quelques coefficients de l'auto corrélation de ton signal et après tu peux utiliser l'algo de Levinson qui résoud les équations de Yule-Walker pour estimer tes coefs du modèle AR.

Apparemment ces deux fonctions sont dans la signal processing tb mais tu dois probablement pouvoir les trouver sur le net et sinon cela peut se programmer à la main... (je l'ai déjà fait pour la corrélation et pour levinson).

[franchouze]

merci robotron

je vais essayer de trouver ces codes.