Discussion :

[zicos]

Bonjour,
Je souhaite faire des statistiques a partir de fichier composer de '1' ou '0'. exemple pour la premiere boucle je compte le nombre de '1' de 1 a 8 dans une variable et je voudrais compter le nombre de '0' (ou de '1') dans le reste de mon vecteur,

Code :

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

1
2
3
4
5
6
7
8
9
10
for i=1:1472
   StatR(:,i)=load ([DIR,list(i).name]);
   Id=fix(i/9)*8+1;
   Rate(i,1)=size(find(StatR(Id:Id+7,i)),1); % compte le nbre de 1
   Rate(i,2)=(100/8)*Rate(i,1); % calcul le pourcentage
   tmp=size(find(StatR(1:Id,i)),1)+; % cet ici qe j'aimerais compter le nombre de '0' dans tout mon vecteur sauf dans la partie du vecteur ou je viens de 
chercher au dessus
   Rate(i,3)=1472-8-tmp
   Rate(i,4)=(100/(1472-8))*Rate(i,3);  % dans le cas ou on compte le nombre de 0
end

Merci de vos reponses.

[zicos]

Je chercais en fait a savoir si il exitait une option ou command pour chercher partout sauf a un endroit specifier.Voila la solution que j'ai trouver mais j'ai un bug que je ne comprend pas. et un autre pour les huit premieres boucles avec ceci

Code :

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

tmp=size(find(StatR(1:Id,i)),1)+size(find(StatR(Id+7:1464,i)),1);

.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
for i=1:1472
   StatR(:,i)=load ([DIR,list(i).name]);
   Id=fix(i/9)*8+1;
   Rate(1,i)=size(find(StatR(Id:Id+7,i)),1);
   Rate(2,i)=(100/8)*Rate(1,i);
   tmp=size(find(StatR(1:Id,i)),1)+size(find(StatR(Id+7:1464,i)),1);
   Rate(3,i)=1464-tmp;
	Rate(4,i)=(100/1464)*Rate(3,i);
 
end

[Jerome Briot]

En général, dans un vecteur avec uniquement des 0 ou des 1 :

Code :

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

V=[1 1 0 0 1 0 0 0 1]

• pour compter le nombre de 1 :

  Code :

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

  n=sum(V)

• pour compter le nombre de 0 :

  Code :

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

  n=sum(~V)

• pour compter le nombre de 1 dans les k premiers éléments :

  Code :

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

  n=sum(V(1:k))

• pour compter le nombre de 0 dans les k premiers éléments :

  Code :

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

  n=sum(~V(1:k))

• pour compter le nombre de 1 dans les k derniers éléments :

  Code :

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

  n=sum(V(end-k+1:end))

• pour compter le nombre de 0 dans les k derniers éléments :

  Code :

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

  n=sum(~V(end-k+1:end))

[LordPeterPan2]

Pour un vecteur VECT rempli de 1 ou de 0 :

Code :

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

1
2
nnz(VECT)
sum(VECT)

Le premier renvoie le nombre d'éléments non nul de VECT
Le second renvoi la somme des termes de VECT
MatDub à montrer (preuve à l'appui ) que la seconde méthode est plus rapide

[zicos]

Voic la solution que j'ai trouver enfin elle devrait marcher mais j'ai un bug bizarre. j'ai relu plusieurs fois mon code mais je ne vois pas d'erreur qui expliquerais mon bug. Si qq'un pouvait le tester

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
StatR=zeros(1472,1472);
Rate=zeros(4,1472);
for i=1:1472
   tmp=0;
   StatR(:,i)=load ([DIR,list(i).name]);
   Id=fix(i/9)*8+1;
   Rate(1,i)=size(find(StatR(Id:Id+7,i)),1);
   Rate(2,i)=(100/8)*Rate(1,i);
   if (fix(i/9)*8+1 ~=1 )
      tmp=size(find(StatR(1:Id,i)),1)
   end   
   tmp=tmp+size(find(StatR(Id+8:1464,i)),1)
   Rate(3,i)=1464-tmp;
	Rate(4,i)=(100/1464)*Rate(3,i);
 
end

NB: Je suis toujours interesser pour savoir s'il existe un moyen de chercher dans tout un vecteur sauf un endroit specifier.

[zicos]

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
StatR=zeros(1472,1472);
StatR(1:3,1)=1;
StatR(5:8,1)=1;
Rate=zeros(4,1472);
for i=1:1472
   tmp=0;
   %StatR(:,i)=load ([DIR,list(i).name]);
   Id=fix(i/9)*8+1;
   Rate(1,i)=sum(StatR(Id:Id+7,i));
   Rate(2,i)=(100/8)*Rate(1,i);
   if (fix(i/9)*8+1 ~=1 )
      tmp=sum(StatR(1:Id+7,i));
   end   
   tmp=tmp+sum(StatR(Id+8:1464,i));
   Rate(3,i)=1464-tmp;
	Rate(4,i)=((100/1464)*Rate(3,i));
end

le find ou le sum font la meme chose donc toujours le meme bug bizarre apres l'execution du programme ci dessous qui correspond au chargement de mon premier fichier. Apres le calcul Rate(3,i) il ajoute 1.0e+003* dans la premiere case et ensuite toud les chiffres sont diviser par 1000.
Matlab ne peut afficher des milliers???

[Jerome Briot]

le find ou le sum font la meme chose

Euh... NON

Sinon, ton problème est juste un problème de format d'affichage. Essaie ceci :

Code :

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

format long

[zicos]

Euh... NON

Je voulais dire que les deux lignes suivantes donnaient le meme resultat:

Code :

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

1
2
Rate(1,i)=size(find(StatR(Id:Id+7,i)),1);
Rate(1,i)=sum(StatR(Id:Id+7,i));

Le 'format long' ajoute juste des chiffres apres la virgule. Ce n'est qu'une question d'affichage, pas tres important mais bon je trouve ca quand meme bizarre. Une histoire d'option d'Editor sans doute.

Merci pour les reponses. Je resolu