Extraction de lignes à partir de matrices successives

**tutty** · 13/05/2008, 18h37

salut tout le monde;
j'ai un gros problème surtout que je suis en train d'utiliser un programme qui me donne de fausses resultats sans se rendre compte que maintenant, te la je dois refaire tout mon travil dans un temps limite;
alors j'ai en fait des fichiers tri horires (8 fichiers par jour) chacun comporte une matrice de 9 lignes et 12 colonnes, et comme sortie je dois avoir 9 fichiers qui comporte les donnée sur toute la periode; ce qui reviens a extraire la première ligne de chaque fichier tri horaire et la mettre dans le fichier de sortire1, et ainsi de suite;
a la fin j'aurais 9 fichiers, chacun comporte tout les données ;

pour cela j'ai fais ainsi

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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
 
function store_all
 
 
cd 'donnees_ALADIN';
 
enter = input('Si vous voulez stocker les données du point le plus proche de la station de Tun C, TAPEZ 1\n Si vous voulez stocker les données du point le pp de Mo, TAPEZ 2\n  Si vous voulez stocker les données du point le pp de Je, tapez 3\n   Si vous voulez stocker les données du pp du projet Rawed, TAPEZ 4\n    Si vous stocker les données du pp de halk el menzel, TAPEZ 5\n     Si vous stocker les données du pp de Achtart, TAPEZ 6\n      Si vous stocker les données du pp PT_1, TAPEZ 7\n       Si vous stocker les données du pp PT_2, TAPEZ 8\n        Si vous stocker les données du pp PT_3, TAPEZ 9\n   Si vous voulez quitter l exécution du programme, TAPEZ 10\n') ; % entrer le choix
 
if (enter==1)
      outt = num2str('TU-C');
end    
if (enter==2)
     outt = num2str('MO');
end  
if (enter==3) 
     outt = num2str('JE');
end  
if (enter==4)
     outt = num2str('Achtart');
end  
if (enter==5)   
     outt = num2str('Halkelmenzel');
end  
if (enter==6)   
     outt = num2str('Rawed');
end
if (enter==7)   
     outt = num2str('PT-1');
end
if (enter==8)   
     outt = num2str('PT-2');
end
if (enter==9)   
     outt = num2str('PT-3');
end
if (enter==10) return;
end
 
 
mounth = input ('Si vous voulez travailler pour les données de janvier, tapez 1\n si vous voulez travailler pour les données de février, tapez 2\n  Si vous voulez travailler pour les données de mars, tapez 3\n   Si vous voulez quitter l exécution du programme, tapez 4\n') ; % entrer le choix
if (mounth ==1);
    mois=num2str('ALADIN-JAN-');
    date = num2str('ASCII_200701%02d');
    N =31;
end
if (mounth ==2);
    mois=num2str('ALADIN_-FEB-');
    date = num2str('ASCII_200702%02d');
    N =28;
end
if (mounth ==3);
    mois=num2str('ALADIN-MAR-');
    date = num2str('ASCII_200703%02d');
    N =31;
end
if (mounth ==4) return;
end
 
 
 
 
for j=1:N;
for i=0:3:21;
    k=(j-1)*8+(i/3);
    nom_pt=num2str(j,date);
    heure= num2str(i,'_%02d');
    nom=[nom_pt heure];
    L = load (nom);
    h = (L);
    [m,n] = size(h);
    S1(k+1,:)= h(1,:);
    S2(k+1,:)= h(2,:);
    S3(k+1,:)= h(3,:);
    S4(k+1,:)= h(4,:);
    S5(k+1,:)= h(5,:);
    S6(k+1,:)= h(6,:);
    S7(k+1,:)= h(7,:);
    S8(k+1,:)= h(8,:);
    S9(k+1,:)= h(9,:);
end
end
 
cd '..';
  cd 'ALADIN_POINT';    
  % save data as an ASCII file
 
 
     sortie1= [mois outt];
 
     sortie2= [mois outt];
 
     sortie3= [mois outt];
 
     sortie4= [mois outt];
 
     sortie5= [mois outt];
 
     sortie6= [mois outt];
 
     sortie7= [mois outt];
 
     sortie8= [mois outt];
 
     sortie9= [mois outt];
 
 
save(sortie1, 'S1','-ascii'); 
save(sortie2, 'S2','-ascii'); 
save(sortie3, 'S3','-ascii'); 
save(sortie4, 'S4','-ascii'); 
save(sortie5, 'S5','-ascii');
save(sortie6, 'S6','-ascii');
save(sortie7, 'S7','-ascii');
save(sortie8, 'S8','-ascii');
save(sortie9, 'S9','-ascii');
 
 
cd '..';
 
return

c'est un peux long je sais mais l'important c'est la boucle, maintenant le problème c'est que j'ai toujours des fichiers semblables contenant les données de la 9ème ligne; donc y'a une faute de stockage, aparement ca efface les données précédent a chaque fois et c pour ca que juste la dernière ligne que je retrouve pour tut les fichiers

s'il vous plai aidez moi c vraiment urgent

**Caro-Line** · 13/05/2008, 19h44

Tous tes fichiers sortie1...sortie9 portent le même nom donc forcément ça ne peut pas fonctionner.
Comment souhaites-tu les appeler ?

**Caro-Line** · 13/05/2008, 19h52

Quelques améliorations sinon à ton programme avec quelques commentaires :

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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
function store_all
cd 'donnees_ALADIN';
 
choice=menu('Choix du stockage des données du point le plus proche',...
    {'de la station de Tun C'...
    'de Mo'...
    'de Je'...
    'du projet Rawed'...
    'de halk el menzel'...
    'de Achtart, TAPEZ 6\n'...
    'de PT_1'...
    'de PT_2'...
    'de PT_3'...
    'Quitter l exécution du programme'}) ; % entrer le choix
 
switch choice
    case 1
        outt = 'TU-C';
    case 2
        outt = 'MO';
    case 3
        outt = 'JE';
    case 4
        outt = 'Achtart';
    case 5
        outt = 'Halkelmenzel';
    case 6
        outt = 'Rawed';
    case 7
        outt = 'PT-1';
    case 8
        outt = 'PT-2';
    case 9
        outt = 'PT-3';
    case 10
        return;
end
%Choix du mois
choice=menu('Choix du mois',{'Janvier','Février','Mars','Quitter l exécution du programme'});
switch choice
    case 1;
        mois='ALADIN-JAN-';
        date = 'ASCII_200701%02d';
        N =31;
    case 2
        mois='ALADIN_-FEB-';
        date = 'ASCII_200702%02d';
        N =28;
    case 3
        mois='ALADIN-MAR-';
        date = 'ASCII_200703%02d';
        N =31;
    case 4
        return
end
 
S=cell(1,9);%initialise la cellule S
 
for j=1:N;
    for i=0:3:21;
        k=(j-1)*8+(i/3);
        nom_pt=num2str(j,date);
        heure= num2str(i,'_%02d');
        nom = [nom_pt heure];
        L = load (nom);
        h = (L);%ça sert à quelque chose ça ?
        [m,n] = size(h);%ne sert à rien ??
        for nC=1:9
            S{nC}(k+1,:)= h(nC,:);
        end
    end
end
 
cd '..';
cd 'ALADIN_POINT';
% save data as an ASCII file
for k = 1:9
    sortie{k} = [mois outt]; %a modifier
    %si vraiment les variables doivent s'appeler S1, ...S9
    eval(['S',num2str(k),' = S{',num2str(k),'};']);
    save(sortie{k}, sprintf('S%d',k),'-ascii');
end
cd '..';

Note : je n'ai pas testé puisque je n'ai pas les fichiers qui vont bien.

**tutty** · 14/05/2008, 08h52

merci Caro; je vais tester et je vous tiens au courant

**tutty** · 14/05/2008, 09h10

bein non caro

, j'ai les fichiers qui porte des nom différents mais à l'interieur ce sont toujours les données du dernier point (de la derniere ligne de chaque matrice); je vous envois des fichiers de deux jours, dans ce cas N=2 et les fichiers contiendrons 16 lignes chacun.

merci

**tutty** · 14/05/2008, 09h12

pour les noms de variables; c'est pas important tant que j'aurais le resultat final; mais la je comprends que j'ai un problème de stockage et j'arrive toujours pas à le corriger

**tutty** · 14/05/2008, 09h22

autre chose

pour eclaircir ce que j'ai envoyé, en fait ce sont les données du 1er et du 2 mars 07;

**tutty** · 14/05/2008, 12h16

mon dieu, j'arrive pas a avoir mon résultat, j'ai changé le programme à plusieurs reprise et j'obtient toujours un resultat faux;et le programme ne fournis que les donées des dernières lignes pour les neufs fichiers, juste le nom change;

vous avez essayer de tester Caro??? avec les fichiers que je vous ai envoyé???

merci, j'ai vraiment besoin de votre aide, si vous avez une idée pour que je puisse corriger ou une indication sur la faute

**Caro-Line** · 14/05/2008, 13h21

Peux-tu montrer ce que tu obtiens pour la variable L à la 1ère itération ?

Note : je n'ai absolument pas le temps de tester et pas MATLAB à ma disposition actuellement.

**tutty** · 14/05/2008, 13h56

voilà Caro; ce sont trois fichiers des neufs que j'obtients pour le mois de mars et ils sont tous pareils; ils contiennent les données de la neuvièmes ligne à chaque itération.

sorry si je te dérange, je sais que c'est pas évident n'ayant pas Matlab à disposition, mais c juste que j'arrive pas à avancer et le deadline est troooops proche.

**Jerome Briot** · 15/05/2008, 15h22

Peux-tu au moins nous expliquer ce que sont censés faire les lignes suivantes ?

Code :

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

outt = num2str('TU-C');

Code :

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

1
2
3
4
5
if (mounth ==3);
    mois=num2str('ALADIN-MAR-');
    date = num2str('ASCII_200703%02d');
    N =31;
end

Je ne vois pas pourquoi tu utilises NUM2STR ici

**tutty** · 15/05/2008, 15h52

ah oui!!!! je pense que ce que vous dites est vrai;
en fait je suis pas très forte coté programmation et surtout avec matlab, c un peux nouveau pour moi et c pour ca que je me suis plutot intéréssé au resultat final sans trops faire attention; bon je vai les enlever, tant que ca marche sans num2str, c'etai juste pour nommé les fichiers

**Caro-Line** · 15/05/2008, 21h25

Envoyé par tutty

bein non caro

, j'ai les fichiers qui porte des nom différents mais à l'interieur ce sont toujours les données du dernier point (de la derniere ligne de chaque matrice);

As-tu vraiment des noms de fichiers différents pour sortie1, sortie2. Je ne crois pas. Ils s'appellent tous les 9 [nom outt]. c'est ça qu'il faut que tu changes.
En reprenant la fin du code que je t'ai donné :

Code :

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

1
2
3
4
5
6
for k = 1:9
    sortie{k} = [mois outt '_',num2str(k)]; %ajoute aussi un n°
    %si vraiment les variables doivent s'appeler S1, ...S9
    eval(['S',num2str(k),' = S{',num2str(k),'};']);
    save(sortie{k}, sprintf('S%d',k),'-ascii');
end

**roilait** · 16/05/2008, 10h13

Bonjour,
Je lis ton problème vite fais, mais j'ai vu que tu parles qu'en tes valeur tu ne retrouve que les dernière lignes de tes donnée.

Je te propose d'initialiser ton tableau de stockage, je pense que sa va aller.
par exemple

Code :

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

1
2
3
4
5
6
 
x=donnee
mois_janvier=[];
for i=1:.....
mois_janvier=[mois_janvier x]%ou " mois_janvier=[mois_janvier;x]"
end

Extraction de lignes à partir de matrices successives

MATLAB

Discussions similaires

Partager

Partager