Discussion :

[belton]

Bonjour à tous,

Je m'initie à Matlab (R2008b)depuis quelques semaines et je fais maintenant face à un problème d'ordre syntaxique (probablement, du moins je l'espère...) assez énervant. Je cherche en effet à créer un matrice à N dimensions en la remplissant au fur et à mesure par des matrices incrémentées de N-1 dimensions (pardon pour la clarté... ^^). Plus simplement (prenons N=3) je cherche à construire une matrice cubique constituée d'un ''empilement'' de matrices n*p. Le code exposé ci-dessous éclaircira peut-être mes dires...

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
B = [];
C = [];
 
for i=1:5             
    for j=1:3         
        for k=1:3     
 
        B(j,k) = blablabla;     % ''blablabla'' simule dans cet exemple l'endroit où sont sockées mes valeurs...
        end
    end
 
    B(:,:)       % à ce stade B(j,k) est bon pour chaque '' i ''
 
    C(i,:,:) = C(i,B(:,:))    % Cette syntaxe ne marche évidement pas, sniff...
 
end 
 
C(4,:,:)    % J'aimerais donc pouvoir lire la ''matrice'' correspondant par exemple à B lorsque i=4.

Voilà, j'ai beaucoup de mal à trouver la syntaxe pour ce problème, dont l'algorithme paraît pourtant simple.... Il est peut être faux...je ne sais pas...Si quelqu'un pouvait m'aider je lui en serais grandement reconnaissant ^^ !!!!!!!
Merci d'avance!!

[Jerome Briot]

Je ne suis pas sûr de tout bien comprendre mais je ferais plutôt ceci :

Code :

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

C(:,:,i) = B;

[belton]

Merci de ta réponse.
En fait ce type de syntaxe ne passe pas car C et B ne sont pas de la même dim:

dim(C) = dim(B)+1

Pour être plus explicite, dans le cas de matrices une dimension, soit des vecteurs, je cherche à les ranger un à un dans une matrice d'ordre 2 (un tableau). ex:

B1 = [1 2 3 4 5]
B2 = [4 8 6 9 67]

et je voudrais les ranger dans C(:,:) = [B1 B2].

Pour cet exemple pas de problème, sauf que dans mon cas les B(i) sont d'ordres supérieurs (typiquement d'ordre 2) et sont calculés au fur et à mesure dans une boucle. J'aimerais donc les ranger au fur et à mesure dans une matrice C(:,:,:) d'ordre 3. C(1,:,:) correspondrait donc au premier ''tableau'' des B calculés, C(2,:,:) au deuxième, etc.....

J'espère avoir été plus clair ^^!
Merci

[kmaniche]

Essaye de faire ceci :

Code :

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

1
2
 
C = [B1 B2] ; % On appelle concaténation horizontale

mais celà ne marche pas :

Code :

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

1
2
 
C = [B1 ; B2] ; % On appelle concaténation verticale

[belton]

Salut kmaniche,

en fait étant dans une boucle for (qui décompte le nombre '''d'étages"" de ma matrice cubique), le C doit se remplir au fur et à mesure que B est calculé, je ne peux donc pas utiliser C = [B1 ...... Bn] ....

[kmaniche]

Et la solution que je t'ai proposée :

Code :

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

1
2
3
4
 
for i=1:n
    C = [C B] ;
end

[Jean Dumoncel]

Salut,

je n'ai pas bien compris pourquoi la solution de Dut ne convenait pas? Tu cherches bien à mettre des tableaux 2D les uns à la suite des autres pour obtenir un tableau 3D, non? Est-ce que B a toujours les mêmes dimensions à chaque itération de i?

Pourrais-tu donner un exemple simple de la tête que dois avoir ton tableau final? Quelle doivent être ces dimensions?

[kmaniche]

J'ai pensé que tu cherché à concaténer tes matrices horizontalement. Explique d'avantage. Ton problème n'a pas l'aire d'être trop compliqué

[belton]

En réalité je n'ai accès à Matlab qu'à mon travail, et je n'ai donc pas pu tester le code que me propose Dut....
Le vrai problème est en effet que mes matrices B n'ont pas le même nombre d'éléments, ce qui complique la création de la matrice. Typiquement, je cherche à obtenir un ''immeuble'': rdc je construit 10 appartements sous la forme 2*5. Ensuite je fais des calculs et je construit au premier étage 18 apparts (genre 3*6). Puis un deuxièmes étages avec 32 apparts, un troisième avec 8, etc....
PS: désolé pour l'exemple, je ne suis pas très fort en explications... (heureusement pour nos jolies têtes blondes je ne suis pas prof...)

[Caro-Line]

Donc ce n'est tout simplement pas possible (si j'ai bien compris).

Tu peux changer de façon de voir les choses et utiliser un tableau de cellules.

[belton]

typiquement

B1 ferais 180000 par 4
B2 ferais 31000 par 4
B3 fearis 64000 par 4 (jmax varie mais pas kmax fixé à 4)
.
.
.
B20 ..........

Donc créer une matrice C de taille 180000*4*20. Le tableau B ayant une taille variable, il faudrait remplir de 0 les ''cases'' sans nombre...

[belton]

Salut Caro_Line,

ça paraît être en effet une meilleure approche!!
Tu pense que je vais pouvoir remplir le tableau de cellules au fur et à mesure?

[kmaniche]

typiquement

B1 ferais 180000 par 4
B2 ferais 31000 par 4
B3 fearis 64000 par 4 (jmax varie mais pas kmax fixé à 4)
.
.
.
B20 ..........

Donc créer une matrice C de taille 180000*4*20. Le tableau B ayant une taille variable, il faudrait remplir de 0 les ''cases'' sans nombre...

Pourquoi un nombre aussi important de cellules
B1 ferais 180000 par 4 ou bien 18000

[belton]

... c'est bien 180 000....
C'est des données d'enregistrements de capteurs hautes fréquences (de l'ordre de la dizaine de kHz)

[kmaniche]

Et la solution que je t'ai proposée :

Code :

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

1
2
3
4
 
for i=1:n
    C = [C B] ;
end

[belton]

C'est vrai que c'est à tester. Je pense par contre que sans passer par le système de tableau de cellules je ne pourrais pas m'en sortir, comme m'a conseiller Caro_Line. La variabilité de jmax va me poser des pb...

[Caro-Line]

Là j'ai un peu plus de temps, et oui le tableau de cellules c'est l'idéal dans ton cas je pense :

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
%B = [];ça ça sert à rien puisque la taille de B va changer pour chaque i
tailleC = 3;
C = cells(1,tailleC);%C sera de taille 1*tailleC
 
for i=1:tailleC           
    %je suppose que c'est ici que sont définis ligneB et colB
    B = zeros(ligneB,colB); %preallocation
    for j=1:ligneB         
        for k=1:colB   
 
        B(j,k) = blablabla;     % ''blablabla'' simule dans cet exemple l'endroit où sont sockées mes valeurs...
        end
    end
 
    C{i} = B;    % Tout simplement
 
end

[Jerome Briot]

Ce qui donnerait sans B :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
tailleC = 3;
C = cells(1,tailleC);%C sera de taille 1*tailleC
 
for i=1:tailleC           
    %je suppose que c'est ici que sont définis ligneB et colB
    C{i} = zeros(ligneB,colB); %preallocation
    for j=1:ligneB         
        for k=1:colB   
 
        C{i}(j,k) = blablabla;     % ''blablabla'' simule dans cet exemple l'endroit où sont sockées mes valeurs...
        end
    end
end