Discussion :

[abel413]

Bonjour,

A présent j'ai d'autres questions toujours en rapport avec cette discussion

Ce qui m'intéresse c'est donc de stocker l'ensemble des valeurs Imp dans un un vecteur, c'est ce que j'ai fait dans mon vecteur valeurs_Imp_global ; j'ai ensuite besoin de toutes sortes d'informations : d'abord réaliser un histogramme des valeurs, c'est que je pense avoir fait ci-dessous mais je ne suis pas très satisfait. Je stocke respectivement dans nelements et xcenter le nombre d'éléments atteignant une certaine valeur et le centre des différentes branches de l'histogramme. Je vous laisse récupérer mon code ci-dessous et le fichier en PJ pour tester dans un premier temps ou sinon regarder après le code.

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
pathname = uigetdir('C:\MATLAB701\work\IMP\DATA\', 'Pick a Directory');
nomfichier0=fullfile(pathname,'PhData.int');
nomfichier1=fullfile(pathname,'Imp1.int');
nomfichier2=fullfile(pathname,'Imp2.int');
nomfichier3=fullfile(pathname,'Imp3.int');
nomfichier4=fullfile(pathname,'Imp4.int');
nomfichier5=fullfile(pathname,'Imp5.int');
nomfichier6=fullfile(pathname,'Imp6.int');
fid0=fopen(nomfichier0);
fid1=fopen(nomfichier1);
fid2=fopen(nomfichier2);
fid3=fopen(nomfichier3);
fid4=fopen(nomfichier4);
fid5=fopen(nomfichier5);
fid6=fopen(nomfichier6);
valeurs_Imp1=fread(fid1,'short');
valeurs_Imp2=fread(fid2,'short');
valeurs_Imp3=fread(fid3,'short');
valeurs_Imp4=fread(fid4,'short');
valeurs_Imp5=fread(fid5,'short');
valeurs_Imp6=fread(fid6,'short');
valeurs_Imp_global=[valeurs_Imp1;valeurs_Imp2;valeurs_Imp3;valeurs_Imp4;valeurs_Imp5;valeurs_Imp6];
 
figure(1); %Makes a new figure
hist(valeurs_Imp_global); %Makes a histogram using the data specified in the first input. The second input is the number of bins.
[nelements,xcenter]=hist(valeurs_Imp_global);
pourcentage=(nelements/sum(nelements))*100
title('Histogram');
% filename=fopen('HIST.xls');
% xlswrite(filename,nelements,xcenter)

Je n'arrive pas bien à interpréter l’histogramme et les valeurs obtenues :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
nelements =
 
     1    13    17    27    22    23     7     4     1     5
 
 
xcenter =
 
 -853.7000 -705.1000 -556.5000 -407.9000 -259.3000 -110.7000   37.9000  186.5000  335.1000  483.7000
 
pourcentage =
 
    0.8333   10.8333   14.1667   22.5000   18.3333   19.1667    5.8333    3.3333    0.8333    4.1667

Pour les xcenter, d'accord c'est le centre des différents rectangles. Pour les nelements, il est indiqué combien il y a d'éléments qui indiquent une valeur donnée mais laquelle ?

D'autre part, est-il judicieux de faire un histogramme de la sorte ou est-ce plus simple autrement ?

Enfin j'aurais besoin d'effectuer un calcul de pourcentage : tel valeur est atteinte un certain nombre de fois : j'ai donc envie de diviser chaque élément du vecteur éléments par le nombre d’éléments total et de multiplier par 100 mais à nouveau je ne comprend pas quel valeur est atteinte.

[le fab]

salut

nelement est le nombre d'éléments "dans une barre", ou autour du "xcenter"
dans ton cas tes xcenter sont espacés de 148.6, une barre correspond donc à toutes les valeurs telles que
xcenter(i) - 148.6/2<= x <= xcenter(i) + 148.6/2
par exemple pour la 3ème barre :
il y a 17 élément compris entre -630.8 et -482.2 et cela représente 14.1667 % de ton échantillon

dit autrement, la valeur 556.5 (+- 74.3) a une probabilité de 14.1667 %
par contre faire un calcul de probabilité sur une valeur précise (556.5 exactement) n'a pas beaucoup de sens me semble t'il

D'autre part, est-il judicieux de faire un histogramme de la sorte ou est-ce plus simple autrement ?

oui totalement
cela permet de bien voir la répartition de ton échantillon, et en l'occurrence, qu'il ne semble pas suivre une répartition gaussienne

[abel413]

Merci pour ces précisions le_fab !

D'autre part, j'aurais besoin de pouvoir écrire ces données (nelements, xcenter et pourcentage) dans un fichier Excell. Le problème c'est que je voudrais le faire au fur et à mesure mais je ne sais pas comment le faire à partir de ma fonction. En gros je voudrais que dans l'Excell il s'inscrive ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
nelements 
 
           1           0           0           0           0           0           1         137      163895      346686
 
 
xcenter 
 
   10.0235   30.0706   50.1177   70.1647   90.2118  110.2589  130.3060  150.3530  170.4001  190.4472
 
 
pourcentage 
 
    0.0002         0         0         0         0         0    0.0002    0.0268   32.0910   67.8818

Et au fur et à mesure que je traite les différent sous-répertoires c'est-à-dire qu'en dessous de nelements j'ai donc 10 valeurs, et en dessous encore 10 et ainsi de suite jusqu'à avoir traité mes répertoires, pareil pour xcenter et pourcentage. Evidemment, cela m'oblige à ouvrir tous ces répertories manuellement. En fait je préférais une automatisation à nouveau : tous les sous-répertoires de C:\MATLAB701\work\IMP\DATA\ sont examinés et les données écrites dans le Excell.

[le fab]

salut

regarde xlswrite

[abel413]

Oui merci le_fab, c'est que je fais en fait mais je n'arrive pas à "automatiser" le processus. Je peux écrire une ligne mais à chaque fois et pour chaque fichier seulement et pas l'un après l'autre.

[le fab]

tu parcours tous tes fichiers (dans tous tes repertoires) avec la fonction dir et des boucles for
à chaque calculs tu met tes résultats dans une grosse matrice (ou tableau de cellules les résultats n'ont pas la même tailles à chaque calcul)
à la fin tu écris cette grosse matrice via xlswrite

[abel413]

Oui c'est ce qui m'a d'ailleurs été demandé de faire mais je n'arrive pas à comprendre le fonctionnement de dir. Pourrais-tu donner un exemple, pas forcément dans le cas que je veux traiter, mais disons qu'on a un répertoire C:/work et qu'il contient les sous répertoire C:/work/s1 et C:/work/s2 lesquels contiennent respectivement les fichiers 1.txt et 2.txt ; je voudrais donc parcourir ces deux sous-répertoires automatiquement, faire des calculs avec les données de 1 et 2 et stocker les résultats dans une cellule ou matrice et écrire l'ensemble dans un fichier excell. Evidemment si tu as un autre exemple plus "simple" je suis preneur.

[le fab]

un debut d'exemple si toujours le même nombre de barres
(sinon tableau de cellules à la place de la matrice)

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
34
35
36
37
38
39
40
 
% initialisation matrice résultat
mResults = [];
 
% repertoire de tête
sDir = dir('C:/work') 
 
% supression de tout ce qui n'est pas dossier
sDir(~[sDir.isdir]) = [];
 
% on enlève les '.' et '..'
sDir(strcmp({sDir(:).name},'.')) = [];
sDir(strcmp({sDir(:).name},'..')) = [];
 
% compteur de fichier lus
nbFiles = 0;
 
% on parcours tous les dossiers retenus
% a la recherche de .txt
for ii=1:numel(sDir)
  sFichiers = dir(fullfile(sDir(ii),'*.txt'));
 
  % on parcours tous les fichier
  for jj = 1:numel(sFichier)
 
    % chargement du fichier sFichier{ii}
    % calculs
 
    % increment
    nbFiles = nbFiles +1;
 
    % stockage dans mResults
    mResults(3*(nbFiles -1)+1,:) = nbelements;
    mResults(3*(nbFiles -1)+2,:) = xcenter;
    mResults(3*(nbFiles -1)+3,:) = pourcentage;
  next
next
 
% ecriture dans un fichier excel
xlswrite('resultats.xls',mResults);

[abel413]

Bonjour et merci de cette réponse le_fab. Je suis en train d'essayer d'adapter ce code à mon problème et surtout de le comprendre mais j'ai beaucoup de mal ! D'après ce que je comprend, ce que je dois écrire se trouve après le deuxième for ; à partir de là, ce que je voudrais réaliser, c'est ce que j'ai entamé dans le code de mon premier post, à savoir obtenir le vecteur valeurs_Imp_global=[valeurs_Imp1;valeurs_Imp2;valeurs_Imp3;valeurs_Imp4;valeurs_Imp5;valeurs_Imp6]; ; avec sFichiers je voudrais donc parcourir l'ensemble des .int à l'exception de PhData.int, voilà pourquoi je voudrais mettre un while :

Code :

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

1
2
3
 while (sFichiers ~= 'PhData.int') 
 
    valeurs_Imp_global= ?

Cependant je suis bloqué car je pense qu'il faudrait mettre des indices mais je n'arrive pas à bien le voir. bref c'est ça le problème, après je pourrais faire mon histogramme et les différents calculs nbelements ect.

EDIT:

Pour clarifier encore plus, je me trouve dans C:\MATLAB701\work\IMP\DATA\ et dedans j'ai différents sous-répertoire : 12-2-1, 29-4-1 ect et dans chacun d'eux j'ai les différents .int qui m'intéresse, comme tu peux le voir dans cette photo :



EDIT2

Ce qui est embêtant c'est rien qu'en testant ce code :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
sDir = dir('C:\MATLAB701\work\IMP\DATA\') 
 
% supression de tout ce qui n'est pas dossier
sDir(~[sDir.isdir]) = [];
 
% on enlève les '.' et '..'
sDir(strcmp({sDir(:).name},'.')) = [];
sDir(strcmp({sDir(:).name},'..')) = [];
 
for ii=1:numel(sDir)
  sFichiers = dir(fullfile(sDir(ii),'*.int'));
end

J'obtiens déjà le message d'erreur suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
Warning: Inputs must be character arrays or cell arrays of strings.
> In fullfile at 31
  In test_auto at 11
??? Error using ==> strrep
Conversion to double from struct is not possible.
 
Error in ==> fullfile at 31
   f = strrep(f,'/','\');
 
Error in ==> test_auto at 11
  sFichiers = dir(fullfile(sDir(ii),'*.int'));

[Invité]

Bonjour,

sDir de la ligne :

Code :

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

sDir = dir('C:\MATLAB701\work\IMP\DATA\')

est un tableau de structures contenant les champs

    name
    date 
    bytes
    isdir
    datenum

donc sDir(ii) est une structure.
Toi tu veux juste récupérer le nom sDir(ii).name.
Attention de plus que ces noms se réfèrent au dossier 'C:\MATLAB701\work\IMP\DATA\', il te faut donc le rajouter devant :

Code :

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

fullfile('C:\MATLAB701\work\IMP\DATA\', sDir(ii).name,'*.int'));

[abel413]

Bonjour et merci pour toutes cette aide, j'ai finalement trouvé une solution à mon problème mais légèrement différente. Je la poste si cela peut être utile :

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
34
35
36
37
38
function valeurs_Imp_global=ecriture_1ph_log(repertoire)
nomfichier0=strcat('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\',repertoire,'\','PhData.int');
nomfichier1=strcat('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\',repertoire,'\','Imp1.int');
nomfichier2=strcat('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\',repertoire,'\','Imp2.int');
nomfichier3=strcat('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\',repertoire,'\','Imp3.int');
nomfichier4=strcat('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\',repertoire,'\','Imp4.int');
nomfichier5=strcat('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\',repertoire,'\','Imp5.int');
nomfichier6=strcat('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\',repertoire,'\','Imp6.int');
fid0=fopen(nomfichier0);
fid1=fopen(nomfichier1);
fid2=fopen(nomfichier2);
fid3=fopen(nomfichier3);
fid4=fopen(nomfichier4);
fid5=fopen(nomfichier5);
fid6=fopen(nomfichier6);
valeurs_Imp1=fread(fid1,'short');
valeurs_Imp2=fread(fid2,'short');
valeurs_Imp3=fread(fid3,'short');
valeurs_Imp4=fread(fid4,'short');
valeurs_Imp5=fread(fid5,'short');
valeurs_Imp6=fread(fid6,'short');
 
%l'ensemble des données d'impédance sont stockées dans un unique vecteur
%valeurs_Imp_global
 
valeurs_Imp_global=[valeurs_Imp1;valeurs_Imp2;valeurs_Imp3;valeurs_Imp4;valeurs_Imp5;valeurs_Imp6];
 
%Détermination de l'histogramme et de nelements, xcenter en recourant à une
%analyse logarithmique ; pour ce faire, les valeurs à considérer doivent
%être nécéssairement positives : on retranche au vecteur d'impédance la 
%valeur minimale d'impédance et on augmente le résultat de 1 avant de
%passer au log.
 
[nelements,xcenter]=hist(20*log(valeurs_Imp_global- min(valeurs_Imp_global) + 1));
 
%calcul du pourcentage des nelements parmis l'ensemble des élements
 
pourcentage=(nelements/sum(nelements))*100;

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
34
35
36
37
38
39
40
41
42
43
%Ce script permet d'écrire de façon ordonnées les différentes grandeurs
%nelements, xcenter et pourcentage dans un fichier Excell pour les
%différentes données d'impédances receuillies.
 
 
clear all
mResults = [];
ALPHA='ABCDEFGHIJKLMNOPQRSTUVWXYZ';
sDir = dir('C:\MATLAB701\work\IMP\CODE_1PH\DONNEES_1PH\') 
sDir(~[sDir.isdir]) = [];
sDir(strcmp({sDir(:).name},'.')) = [];
sDir(strcmp({sDir(:).name},'..')) = [];
i=1;
j=1;
PP=[]; T=[];
while(j>0)
   try
        repertoire=sDir(i,1).name;
        i=i+1;
    catch
 
        j=0;
    end
 
end
nfichier=i-1;
for i=1:nfichier
    repertoire=sDir(i,1).name;
 
    valeurs_Imp_global=ecriture_1pH_NON_log(repertoire);
    [nelements,xcenter]=hist(valeurs_Imp_global);
    pourcentage=(nelements/sum(nelements))*100;
    P=[nelements;xcenter;pourcentage];
    PP=[PP;P];
    T1=strcat(cellstr(repertoire),'-N');
    T2=strcat(cellstr(repertoire),'-XC');
    T3=strcat(cellstr(repertoire),'-%');
    T=[T; T1;T2;T3];
end
rangeT=strcat('A1:A',num2str(size(T,1)));
rangeP=strcat('B1:',ALPHA(size(PP,2)+1),num2str(size(PP,1)));
xlswrite('resultats_NON_log_1ph', T,'P',rangeT);
xlswrite('resultats_NON_log_1ph', PP,'P', rangeP);

L’utilisation combinée de ces deux fonctions permet de stocker dans un fichier Excell les différentes données dont j'ai besoin l'une à la suite des autres avec en plus le nom du fichier et le type de données regarder (nelements, xcenter et pourcentage).

Vous trouverez un exemple d'Excell généré en PJ.

[abel413]

Je me pose une nouvelle question. Admettons que j'ai des données et que je veuille en faire l'histogramme comme ceci :

Code :

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

[nelements,xcenter]=hist(valeurs_Imp_global);

Dans nelements je stocke donc les éléments atteignant une certaine valeur et xcenter le centre des rectangles constitutifs de l'histogramme.

J'ai un doute pour le calcul des pourcentages des éléments parmi l'ensemble des éléments :

Code :

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

pourcentage=(nelements/sum(nelements))*100;

ou pour un nombre différent :

Code :

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

1
2
[nelements_1024,xcenter_1024]=hist(valeurs_Imp_global,1024)
pourcentage_1024=1024*(nelements/sum(nelements))*100;

En fait j'ai du mal à comprendre la signification de ce deuxième paramètre 1024.

[le fab]

J'ai un doute pour le calcul des pourcentages des éléments parmi l'ensemble des éléments :

Code :

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

pourcentage=(nelements/sum(nelements))*100;

c'est ca

Code :

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

1
2
[nelements_1024,xcenter_1024]=hist(valeurs_Imp_global,1024)
pourcentage_1024=1024*(nelements/sum(nelements))*100;

non, ne pas multiplier par 1024

En fait j'ai du mal à comprendre la signification de ce deuxième paramètre 1024.

c'est simplement le découpage que tu souhaites (aka le nombre de barres de ton histogramme)
tu prend un nombre faible pour avoir une tendance générale, on nombre grand pour voir plus finement les "discontinuités"
par défaut il te fait 10 barres (ou il regroupes les données dans 10 sous ensembles)

[abel413]

Merci le_fab, effectivement c'est plus clair.