Discussion :

[diego45]

Bonjour à tous.

J'ai un vecteur colonne Beta de taille n;
J'ai une matrice de covariance notée C de taille [n,n];

Je dois déterminer cette formule :

Code :

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

1
2
3
4
5
 for i=1:n 
        for j = 1:n 
            D=sqrt((sum(Beta(i,1)*Beta(j,1)*C(i,j))
        end 
    end

Pensez-vous que cette formule fasse bien la chose ?

[Jerome Briot]

Plutôt comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
D2 = 0;
 
for i = 1:n
    for j = 1:n
        D2 = D2 + Beta(i)*Beta(j)*C(i,j);
    end
end
 
D = sqrt(D2);

[diego45]

Merci pour ta réponse, en fait C est réestimée tout les 250 jours alors j'ai une boucle for dessus.

Donc je pensais adapter ton code ainsi :

Code :

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

1
2
3
4
5
6
7
8
9
D2(i,1) = 0;
 
for k = 1:n
    for j = 1:n
        D2(i,1) = D2(i,1) + Beta(k)*Beta(j)*C(k,j);
    end
end
 
D(i,1) = sqrt(D2(i,1));

Mais l'ordre de grandeur est absolument pas ce que je dois avoir !

[Dombrai]

sinon il me semble que ta formule s’écrit sans boucle

Code :

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

1
2
 
D2 = Beta'*C*Beta;