Discussion :

[stephane2000]

Bonjour,
Je voulais simuler un AR(1), j'ai fais ce code mais ça marche pas

Code :

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

1
2
3
4
5
6
7
8
9
%AR1 :  yt=a+b*yt-1+et
T=3000 ;
y(0)=10;
a=3;
b=0.7 ;  
et=randn(T,1);
for k=1:T
y(k)=a+b* y(k-1)+ e(k) ;
end ;

[soft001]

Dans matlab on pas un élément du tableau d'ordre 0;

Code :

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

1
2
y(0) n'existe pas
y(1)=..

En plus pour la ligne n° 8 si k=0 tu vas retomber dans le cas expliquant auparavant.
donc tu doit commencer ton boucle avec 2

Code :

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

1
2
3
4
5
6
7
8
9
10
%AR1 :  yt=a+b*yt-1+et
T=3000 ;
y(1)=10;
a=3;
b=0.7 ;  
et=randn(T,1);
for k=2:T
y(k)=a+b* y(k-1)+ et(k) ;
end ;
plot( y)

[Jerome Briot]

mais ça marche pas

"ça marche pas" ne veut rien dire. Il faut être plus précis.

Les résultats ne sont pas ceux attendus ?
MATLAB renvoie un message d'erreur ?

Heureusement que dans ce cas les erreurs était assez flagrantes et que donc soft001 a pu te donner une bonne réponse.

Sinon tu peux supprimer le ; inutile après le end de la boucle FOR-END.

[stephane2000]

Peut-être j'ai fait une faute de frappe
Ce que a répondu soft001 marche bien.
Pour simuler un AR(1), à mon avis jusqu'ici le code est incomplet
Comment peut-on générer 1000 fois le modèle et doit-on travailler avec a et b estimés?