Discussion :

[zicos]

Salut,

Petit casse tete a resoudre. Realiser le code ci dessous sans boucle for. Je cherche en vain mais trouve pas toujours des erreurs avec des nombres de lignes qui correspondent pas. Donc voici la boucle:

Code :

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

1
2
3
4
5
6
7
8
Nm=1;
for lc=1:Nb_R
   Ed_C(Nm:Nb_C+Nm-1,:)=Ed_Ctmp(:,:);
   Dra_C(Nm:Nb_C+Nm-1,:)=Dra_Ctmp(:,:);
   Or_C(Nm:Nb_C+Nm-1,:)=Or_Ctmp(:,:);
   Rc_C(Nm:Nb_C+Nm-1,:)=Rc_Ctmp(:,:);
   Nm=Nm+Nb_C;
end

Je sais pas si c'est possible, juste pour voir, moi j'abandonne.
Bonne Chance Et Merci d'avance de vos recherches et reponses

[coolzy]

il ne s'agit pas que de coller la meme matrice plusieurs fois dans une plus grande matrice? ca doit pouvoir se faire avec des operations par blocs. Maintenant jsais pas comment lui dire qu'il faut repeter la petite matrice NB_R fois dans la grande...

[zicos]

je bloque la dessus aussi. on peut pas concatener des truc du style 1:1:Nb.
Enfin je pense je sais pas justement.
T'es suite ou en norvege si c'est pa indiscret.

[Caro-Line]

Ce serait mieux si tu nous donnais une idée de la taille de chaque matrice.
Ou sur un petit exemple...
Et peut-être les erreurs générées.

[zicos]

25 colonnes et c'est fixe.
Pour Nb_R et Nb_C, ils sont variables mais pour mes tests valent respectivement 31 et 37.
La version avec la boucle for marche sans erreur, tous mes essais sans boucle se sont solder par une ereur de compile.
A noter que cette boucle est deja a l'interieur d'une autre boucle voici le programme complet pour les interresser.

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
for l=1:Nb_C
   i=1;
   for c=1:5
      Ed_Ctmp(l,i:i+4)=TOC_C(l,1:4:17);
      Dra_Ctmp(l,i:i+4)=TOC_C(l,2:4:18);
      Or_Ctmp(l,i:i+4)=TOC_C(l,3:4:19);
      Rc_Ctmp(l,i:i+4)=TOC_C(l,4:4:20);
      i=i+5;		
   end
end
Nm=1;
for lc=1:Nb_R
   Ed_C(Nm:Nb_C+Nm-1,:)=Ed_Ctmp(:,:);
   Dra_C(Nm:Nb_C+Nm-1,:)=Dra_Ctmp(:,:);
   Or_C(Nm:Nb_C+Nm-1,:)=Or_Ctmp(:,:);
   Rc_C(Nm:Nb_C+Nm-1,:)=Rc_Ctmp(:,:);
   Nm=Nm+Nb_C;
end

J'aimerais bien gagner du temps dans mon programme en eliminant cette boucle, si c'est possible.

[rostomus]

Salut,

Code :

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

1
2
3
4
5
6
7
8
Nm=1;
for lc=1:Nb_R
   Ed_C(Nm:Nb_C+Nm-1,:)=Ed_Ctmp(:,:);
   Dra_C(Nm:Nb_C+Nm-1,:)=Dra_Ctmp(:,:);
   Or_C(Nm:Nb_C+Nm-1,:)=Or_Ctmp(:,:);
   Rc_C(Nm:Nb_C+Nm-1,:)=Rc_Ctmp(:,:);
   Nm=Nm+Nb_C;
end

tu peux la remplacer par:

Code :

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

1
2
3
4
5
 
Ed_C=kron(ones(1,Nb_R),Ed_Ctmp);
Dra_C=kron(ones(1,Nb_R),Dra_Ctmp);
...
...

[Jerome Briot]

J'aimerais bien gagner du temps dans mon programme en eliminant cette boucle, si c'est possible.

As-tu lu la FAQ à propos de la préallocation des matrices ?

[coolzy]

jsuis a Bergen actuellement.