Discussion :

[jenaCa]

Salut,
voici mon problème

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
n=10;
 
for iexample=1:3
    if iexample==1
        phi=0.6;
    elseif iexample==2
        phi=0.7;
    else
        phi=0.8;
    end
 
sig_x=1; 
sig_e=sig_x*sqrt(1-phi.^2);
e=normrnd(0,sig_e,n,1);
x(1)=rand();
    for t=2:n
      x(t)=phi*x(t-1)+e(t);
    end    
end

Je souhaite que les résultats soient rassemblés dans une matrice de taille n*3 sachant que x est un vecteur colonne de longueur n.
Merci par avance de votre aide.

[dionisis]

Salut jenaCa,

Le code peut etre simplifie si tu vectorises tes operations (ie en evitant les boucles).

Voici une version qui devrait faire le travail (je n ai vectorise qu une partie de ton code, je pense qu on peu vectoriser la totalite, je le laisse comme ca pour que ce soit plus lisible).

Code :

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

1
2
3
4
5
6
7
8
9
10
 
nr = 10;
phi = 0.6:0.1:0.8;
nc = numel(phi);
sig_e = sqrt(1-phi.^2);
e = repmat(sig_e,nr,1).*rand(nr,nc);
x = [rand(1,nc); zeros(nr-1,nc)];
for t = 2:nr
	x(t,:) = phi.*x(t-1,:) + e(t,:);
end

Cordialement

[jenaCa]

Salut Dionisis,

ça marche avec la vectorisation ,Merci beaucoup.

Amicalement,