Discussion :

[jmbonnaz]

Bonjour,

J'ai le code suivant

Code :

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

1
2
3
 for trial =1:10
X([2*trial-1,2*trial],1:6)=matX(A,B,trial);
end

.

La fonction matX me renvoie une matrice de taille 2*6 que je veux ranger dans la variable X en fonction de la variable trial.

C'est-àdire que quand trial=1, je voudrais que la valeur que me renvoie la fonction aille dans les 2 premières lignes de X, quand trial=2, la valeur aille dans la 3eme et la 4eme ligne ...

Ca j'y arrive mais le probleme est que lorsque trial égal 2 par exemple, les valeurs que j'ai calculé pour trial=1 (les 2 premières lignes) disparaissent. Ainsi je n'arrive pas à avoir une matrice X complète (seulement 2 lignes suivant la valeur trial).

Comment puis je faire pour ne pas perdre ce que j'ai calculé avant ?


A titre indicatif le code de ma fonction est le suivant :

Code :

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

1
2
3
4
5
6
7
function X=matX(A,B,trial)
 
sprintf('Calcul de X')
 
for j=1:6
    X(2*trial-1:2*trial,j)=inv([A(:,1) A(:,j+1)]'*[A(:,1) A(:,j+1)])*[A(:,1) A(:,j+1)]'*B(:,j);
end

.

A est de taille [150000*2].
B est de taille [150000*1

[Jerome Briot]

La variable X dans le premier code et la variable X dans la fonction ne sont pas identiques. Elles sont indépendantes.

La fonction matX devient donc :

Code :

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

1
2
3
4
5
6
7
function X=matX(A,B)
 
sprintf('Calcul de X')
 
for j=1:6
    X(1,j)=inv([A(:,1) A(:,j+1)]'*[A(:,1) A(:,j+1)])*[A(:,1) A(:,j+1)]'*B(:,j);
end

et le code:

Code :

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

1
2
3
for trial =1:10
   X([2*trial-1,2*trial],1:6)=matX(A,B);
end

Note : si A contient des nombres réels, utilise .' à la place de ' pour la transposé. Voir la à ce sujet ici

[jmbonnaz]

merci