Discussion :

[david_Montreal]

Bonjour,

J'ai un petit souci pour calculer les éléments d'une matrice dans une boucle "for" . Je souhaiterai pouvoir calculer la matrice M pour n'importe quelle dimension de x. Mon code (qui ne fonctionne pas) est le suivant:

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
x=1:0.1:1.5
 
for n=1:1:length(x)
 
        g(n)=tan(pi*x(n));
        t(n)=tan(3*pi*x(n));
 
        A(n)=1-t(n).*g(n)*8;
        B(n)=t(n)+8*g(n);
        C(n)=t(n);        
        D(n)=1;
 
        M(:,n)=[A(n) B(n) ; C(n) D(n)];
 
end

En vous remerciant d'avance de me remettre dans le droit chemin!

David

[Jerome Briot]

Quelle doit être la taille finale de M ?

[Caro-Line]

Juste un truc :
là :

Code :

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

M(:,n)=[A(n) B(n) ; C(n) D(n)];

Tu essayes de mettre une matrice 2*2 ([A B ;C D]) dans un vecteur colonne.
Ca ne va pas rentrer...

[galaadb]

Bonjour,

j'ai eu le même problème récemment.
En fait si je comprends bien, tu cherches à réaliser une "suite" de matrices non?
Ta matrice M est une matrice 2x2 dont les éléments varient en fonction de x, c'est ça?
si c'est le cas, il faut ajouter une 3e dimension à ta matrice je crois:

Code :

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

1
2
 
M(1:2,1:2,n)=[A(n), B(n) ; C(n), D(n)];

Enfin, dans mon cas ça avait marché. J'espère que j'ai bien compris ton problème...

[david_Montreal]

Merci pour vos réponses en temps réel!

Oui c'est ca, la matrice M reste une matrice 2x2 mais la longueur des éléments A, B,C et D varie si je change le nombre de points pour x. Je viens d'esseyer ta suggestion. Cela fonctionne mais le hic c'est que je voudrai pouvoir stocker toutes les valeurs dans la matrices M car je dois ensuite faire un produit matricielle avec une autre matrice de dimension identique.

Resultats désiré

M=
[A1 A2 A3 | B1 B2 B3
---------|---------
C1 C2 C3 | D1 D2 D3]

Resultats obtenu

x =
1 1.25 1.5


M(:,:,1) =
1 -1.3471e-015
-3.6739e-016 1

M(:,:,2) =
9 7
-1 1

M(:,:,3) =
-7.9025e+031 4.5365e+016
1.8146e+015 1

Bonjour,

j'ai eu le même problème récemment.
En fait si je comprends bien, tu cherches à réaliser une "suite" de matrices non?
Ta matrice M est une matrice 2x2 dont les éléments varient en fonction de x, c'est ça?
si c'est le cas, il faut ajouter une 3e dimension à ta matrice je crois:

Code :

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

1
2
 
M(1:2,1:2,n)=[A(n), B(n) ; C(n), D(n)];

Enfin, dans mon cas ça avait marché. J'espère que j'ai bien compris ton problème...

[Jerome Briot]

Sans chercher à trop optimiser ton code, tu peux concaténer A, B, C et D après la boucle FOR-END :

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
x=1:0.1:1.5
 
for n=1:1:length(x)
 
        g(n)=tan(pi*x(n));
        t(n)=tan(3*pi*x(n));
 
        A(1,n)=1-t(n).*g(n)*8;
        B(1,n)=t(n)+8*g(n);
        C(1,n)=t(n);        
        D(1,n)=1;
 
end 
 
M =[A B ; C D];

[david_Montreal]

Cette ultime concaténation aura eu raison de ma matrice récalcitrante!

Problème résolu, encore merci à tous et la prochaine fois je ferai un meilleur usage des balises et autres...

Sans chercher à trop optimiser ton code, tu peux concaténer A, B, C et D après la boucle FOR-END :

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
x=1:0.1:1.5
 
for n=1:1:length(x)
 
        g(n)=tan(pi*x(n));
        t(n)=tan(3*pi*x(n));
 
        A(1,n)=1-t(n).*g(n)*8;
        B(1,n)=t(n)+8*g(n);
        C(1,n)=t(n);        
        D(1,n)=1;
 
end 
 
M =[A B ; C D];

[Jerome Briot]

Et donc maintenant, sans la boucle FOR-END :

Code :

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

1
2
3
x=1:0.1:1.5
 
M = [1-tan(3*pi*x).*tan(pi*x)*8 tan(3*pi*x)+8*tan(pi*x) ; tan(3*pi*x) ones(1,length(x))]

Voila du code Matlab propre qui s'optimise encore comme ceci :

Code :

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

1
2
3
4
5
6
x=1:0.1:1.5
 
t=tan(3*pi*x);
tt=tan(pi*x);
 
M = [1-t.*tt*8 t+8*tt ; t ones(1,length(x))]

Ensuite, on va s'arrêter la parce que ça deviendrait trop indigeste pour le moment