Discussion :

[new_matlab]

Salut tout le monde,

Bon j'ai une question assez technique en fait. J'utilise la toolbox System Identification Tool pour faire du fitting sur des courbes. J'arrive à faire marcher la toolbox dans ce sens là pas de soucis.

Etant donné que je manipule des modèle autoregressif à moyenne ajusté (ARMA), je sais que y(t) est calculé en fonction des échantillons précédents (le nombre dépendant des différents ordres du filtre) . La toolbox me sort les coefficients du filtre. Moi ce que j'aimerais avoir c'est la valeur de y(t) calculé à partir des échantillons précédents.

Vous pensez qu'on puisse faire ca à partir des coefs du filtre ARMA sorti?

Merci par avance de votre aide

[phryte]

Bonjour.

Moi ce que j'aimerais avoir c'est la valeur de y(t) calculé à partir des échantillons précédents.

Pour avoir les y(nT) il faut programmer l'équation récurrente :
Exemple pour un ordre 3
Calcul des coef ARMA:

Code :

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

1
2
3
4
5
6
7
8
9
[b,a]=butter(4,0.5,'low')
[h,w]=freqz(b,a,64)
figure(1)
plot(w,abs(h))
hold on
[bb,aa]=invfreqz(h,w,3,3)
[hh,ww]=freqz(bb,aa)
plot(ww,abs(hh),':')
grid

Equation récurrente :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
y(1)=x(1);%INITIALISATION des valeur passées du signal d'entré
y(2)=x(2);
y(3)=x(3);
for n=4:N
    y(n)=bb(1)*x(n)+bb(2)*x(n-1)+bb(3)*x(n-2)+bb(4)*x(n-3)-aa(2)*y(n-1)-aa(3)*y(n-2)-aa(4)*y(n-3);%Equation récurrente
    y(n-3)=y(n-2);
    y(n-2)=y(n-1);
    y(n-1)=y(n);
end
figure(4)
plot(t,y)
grid

[new_matlab]

Merci beaucoup pour ce code...
Néanmoins , il y a qqchose que je ne comprends pas. Pourquoi met on à jour les valeurs de y dans la boucle consacrée à l'équation réccurrente?

[phryte]

Bonjour.

les valeurs de y dans la boucle consacrée à l'équation réccurrente?

Effectivement cela n'est pas nécessaire lorsqu'on travaille avec un tableau. Le faire systématiquement est une déformation du temps réel lorsqu'on travaille sur des variables.

[new_matlab]

oui c'est ce que je pensais effectivement...
Dans le même genre de question, dites moi ce que vous en pensez : Je fais des prédictions à base d'ARMA par segments (je segmente mon signal). Je générè les coefficients de mon ARMA pour chaque segment à l'aide de la fonction ARMAX de Matlab.
Je sais que plus le segment en question est long et plus les coefficients permettront une prédiction meilleure.

Voilà ma question : est-il judicieux de calculer mes coefficients ARMA sur un segment disons de 40 points et de prendre des ordres de 3 et 3 ? Ce qui fait que quand je réinjecterais mes coefficients ARMA dans mon équation réccurente, j'aurais une précision bien plus importante que si j'avais calculé mes ARMA seulement sur 3 points.

Qu'en pensez vous ?

Merci

[mihaispr]

Salut!

Oui. Tu auras une precision plus importante quand le nombre des points est 40.

La precision c'est plus bonne si tu as de nombreux points!

C'est comme dans la regression polynomiale.

Bonne chance,

Michel