Discussion :

[clem3900]

Bonjour,

je souhaiterais réaliser une boucle if comme cela:

Code :

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

1
2
3
4
5
6
7
 
 m=[1:1:10]
 if m > 4
    a=1;
 else
    a=0;
 end

Et je voudrais ensuite créer un tableau avec les différentes valeurs de a.
Ici ce serait donc un tableau de 10 lignes, avec 0 dans les 4 premières lignes, et 1 dans les 6 lignes restantes.


Si quelqu'un pouvait m'aider.

Merci

[shaiHulud]

Bonjour,

tout d'abord, 1:1:10 peut être remplacé par 1:10. Ensuite:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
% avec une boucle
m=[1:1:10];
a= []; % NB: on peut aussi préallouer a via "a= nan(size(m));"
 if m > 4
    a(end+1)=1;
 else
    a(end+1)=0;
 end
 
% version vectorisée, avec indexation booléenne
m=[1:1:10];
a= zeros(size(m)); % initialise à 0
a(m>4)= 1: % mets des 1 aux emplacements ou m>4

[Invité]

Bonjour,

Les crochets m = [1:1:10]; sont inutiles, tu n'effectues pas de concaténation.

On peut faire directement a = m>4;.

[shaiHulud]

Attention aux types:

Code :

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

1
2
3
 
a= m>4; % vecteur de booléens
a= zeros(size(m)); a(m>4)= 1: % vecteur de floats

[Boumbata]

Ou alors:

Code :

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

a = +(m>4); % vecteur de flottants

[clem3900]

Heu, c'est quoi la différence entre booléens, floats et flottants?
La taille des variables stockées?

Et en fait je me suis trompé dans mon code...
Le code que je veux faire c'est quelque chose comme ça:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
l=[1 2 3 4 5 6 7 8 9 10]
 
for m=1:10
 
  a(m>4)=  m./l; 
  a(m<=4)=0;
  ma(:,m)=[a]
 
end

C'est pour ça que je parlais de tableau, parce que je veux récupérer les valeurs de 'a' EN DEHORS de ma boucle for

[Boumbata]

Heu, c'est quoi la différence entre booléens, floats et flottants?
La taille des variables stockées?

Et en fait je me suis trompé dans mon code...
Le code que je veux faire c'est:

Code :

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

1
2
3
4
5
 
for m=1:10
  aa(m>4)= 1; 
  aa(m<=4)=0;
end

C'est pour ça que je parlais de tableau, parce que je veux récupérer les valeurs de 'a' EN DEHORS de ma boucle for

Pour les booléens/flottants/double, je te recommande de lire ce petit passage dans la FAQ Comment déterminer approximativement la quantité de mémoire nécessaire pour le stockage d'une variable ?

Si je comprends bien, tu vas comparer chacun des éléments compris entre 1 et 10 avec 4, et si le i-ième élément est inférieur ou égal à 4, alors a(i) = 0, sinon, a(i) = 1 ?

Dans ce cas, pas besoin de faire de boucles, tu peux vectoriser le calcul en transformant:

Code :

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

1
2
3
4
5
 
aa = zeros(1,10);
for m=1:10
  aa(m)= m>4; 
end

en:

Code :

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

1
2
 
aa = (1:10)>4;

La différence sera que dans la deuxième méthode, aa sera un vecteur de booléens, soit 10 octets pour tout le vecteur, tandis que dans l'autre cas, ce sera un vecteur de Double, soit 10*8 = 80 octets pour le vecteur.

[Invité]

Pour information, on peut aussi obtenir un vecteur de booléens avec une boucle :

Code :

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

1
2
3
4
aa = false(1,10);
for m=1:10
  aa(m)= m>4; 
end

[clem3900]

Boumbata désolé mais j'ai changé le code après que tu ais posté...

Ce que je veux faire c'est si m=5, alors a(5)=m./l, ce qui me donne un vecteur 1x10 (à cause de 'l' qui est de taille 1x10) que j'enregistre dans la 1ère ligne d'un tableau 'ma'. De même pour m=6,7,8,9,10. Et quand m<4, dans mon tableau j'enregistre des lignes de 0.
A la fin j'aurai donc un tableau 10x10 dont 4 lignes seront uniquement des 0.

J'espère que c'est plus clair...

[Invité]

Tu as déjà pas mal d'éléments en main. Je t'invite à effectuer une recherche sur le forum avec les mots clé "indexage logique", ce genre de manipulation ayant déjà été effectué plus d'une fois.

[Boumbata]

a lol oui ça change un peu

Code :

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

1
2
3
4
 
m = 1:10;
l = rand(1,10); % je ne connais pas ton l, donc je prends random
ma = repmat((m.*(m>4)).',size(m))./repmat(l,size(m.'));

La fonction repmat permet de répliquer une matrice n*m fois. Cette méthode permet d'éviter les boucles for.

[shaiHulud]

Outre le stockage en mémoire, le simple exemple suivant plante

Code :

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

1
2
3
 
m = rand(10,1)> 4; %type bool
m==1; % erreur de type

[Invité]

@Boumbata : On peut aussi utiliser bsxfun qui est généralement plus rapide.

@VV33D : m possédant des valeurs logiques, quel est l'intérêt de faire m==1; ?