Discussion :

[noone2006]

bonjour tout le monde, voila j'ai un petit souci avec mon programme pas pour le fonctionnement mais pour la durée d'exécution, il met 20 pour terminer l'exécution .

remarque:
voila qq données complementair pour mieu comprendre mon code :

resol=2;
vect{:}
ans =
detail: [1x93312 double]
ans =
detail: [1x23328 double]

code{:}
ans =
book: [1x29 double]
ans =
book: [1x50 double]

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
etiq = cell(resol,1);
load vect
load code
 
ii=resol;
 
while ii>0
 
ccc = vect{ii}.detail;
codeb = code{ii}.book;
[mc nc] = size(ccc);
[mc nkb] = size(codeb);
etiq_det(1:nc) =0;
 
  for i = 1:nc
          d=100000000;
          for k = 1:nkb
              d1 = ( ccc(i) - codeb(k) )^2;
 
              if(d1<d)
                 l = k;
                 d = d1;
              end
          end
 
          etiq_det(i) =l;
 
          etiq{ii} = struct('detail',etiq_det);
 
 
  end
  ii=ii-1;
end

voila mon code si qqun pouvait m'aider avec une idée ou une astuce je serai ravi , Merci d'avance et je vous signal que c pas mon pc qui est lent ( j'ai un core duo avec 1 giga de ram ).
Merci d'avance

[noone2006]

bonjour, j'ai trouvé la solution a mon problem, il fallait fair sortir le terme :etiq{ii} = struct('detail',etiq_det);
de la boucle, donc comme ca il n'accède a la structure une fois pour chaque résolution, merci comme même.
voila le nouveau code ( au fait c'est un code pour la générer l'étiquette après quantification ).

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
etiq = cell(resol,1);
load vect
load code
 
ii=resol;
 
while ii>0
 
ccc = vect{ii}.detail;
codeb = code{ii}.book;
[mc nc] = size(ccc);
[mc nkb] = size(codeb);
etiq_det(1:nc) =0;
 
  for i = 1:nc
          d=100000000;
          for k = 1:nkb
              d1 = ( ccc(i) - codeb(k) )^2;
 
              if(d1<d)
                 l = k;
                 d = d1;
              end
          end
 
          etiq_det(i) =l;
 
  end
 
  etiq{ii} = struct('detail',etiq_det);
  ii=ii-1;
end

[Jerome Briot]

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
 
  for i = 1:nc
          d=100000000;
          for k = 1:nkb
              d1 = ( ccc(i) - codeb(k) )^2;
 
              if(d1<d)
                 l = k;
                 d = d1;
              end
          end
 
          etiq_det(i) =l;
 
  end

La boucle sur "k" n'est pas nécessaire ici :

Code :

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

1
2
3
4
5
 
  for i = 1:nc
          d1 = (ccc(i) - codeb).^2;
          [d,etiq_det(i)] = min(d1);                
  end

A tester bien entendu