Discussion :

[mfontan]

Bonjour,

je sollicite encore une fois vos lumières pour m'aider à résoudre mon problème.

je souhaite extraire un ensemble de nombre dans un fichier text dont la position (ligne) peut évoluer mais je sais en revanche qu'il se présente toujours de la même manière.

voici la partie du code (qui commence à la ligne i) qui m'intéresse:

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
%   -------------------
                                               R E A L   E I G E N V A L U E S
                                         (BEFORE AUGMENTATION OF RESIDUAL VECTORS)
   MODE    EXTRACTION      EIGENVALUE            RADIANS             CYCLES            GENERALIZED         GENERALIZED
    NO.       ORDER                                                                       MASS              STIFFNESS
        1         1        6.827670E+00        2.612981E+00        4.158688E-01        1.000000E+00        6.827670E+00
        2         2        3.453286E+01        5.876466E+00        9.352686E-01        1.000000E+00        3.453286E+01
        3         3        1.079509E+02        1.038994E+01        1.653610E+00        1.000000E+00        1.079509E+02
        4         4        2.563271E+02        1.601022E+01        2.548105E+00        1.000000E+00        2.563271E+02
        5         5        5.359315E+02        2.315019E+01        3.684468E+00        1.000000E+00        5.359315E+02
        6         6        9.561032E+02        3.092092E+01        4.921217E+00        1.000000E+00        9.561032E+02
        7         7        1.649111E+03        4.060925E+01        6.463162E+00        1.000000E+00        1.649111E+03
        8         8        2.519466E+03        5.019429E+01        7.988668E+00        1.000000E+00        2.519466E+03
        9         9        3.893863E+03        6.240083E+01        9.931399E+00        1.000000E+00        3.893863E+03
       10        10        5.391129E+03        7.342431E+01        1.168584E+01        1.000000E+00        5.391129E+03
  
%   -------------------

je cherche seulement a extraire les données de la colonne définie par le premier et dernier nombre en gras.


mon code actuel permet de charge le fichier grâce a fopen, j'arrive a trouver l'indice de "R E A L E I G E N V A L U E S" et je sais que les nombres qui m'intéressent sont toujours écrits avec la meme syntaxe, à la même position.

d'où ma question, comment extraire et sauvegarder ces nombres (boucle pour la colonne, ok) ? si vous connaissez les opérateurs, je suis preneur que je poursuive mon script et que je revienne vers vous si besoin.

je vous remercie de votre aide par avance,

Maxime

[Jerome Briot]

Une solution vite faite.

Le code récupère le bloc de ligne compris entre la ligne contenant "R E A L E I G E N V A L U E S" et la première ligne vide suivante.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
X = textread('test.txt','%s','delimiter','\n');
 
idxA = find(strcmp(X,'R E A L   E I G E N V A L U E S'));
 
idxB = find(cellfun('isempty',X));
 
idxC = find((idxB-idxA)>0,1,'first');
 
subX = X(idxA+4:idxB(idxC)-1);
 
E = strread(char(subX).','%*d%*d%f%*f%*f%*f%*f')

[mfontan]

Merci pour le code DUT,

il doit y avoir une petite erreur, j obtiens ce message :

Code :

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

1
2
3
4
5
6
7
 
 
??? Error using ==> minus
Matrix dimensions must agree.
 
Error in ==> oui at 11
idxC = find((idxB-idxA)>0,1,'first') ;

J 'ai essaye de modifier la taille de la matrice idxA et idxB mais ça ne marche pas...

[Invité]

Dans ce cas il existe dans ton fichier plusieurs lignes contenant 'R E A L E I G E N V A L U E S'.
Il te suffit de remplacer cela par une partie unique à l'extrait que tu nous montres dans ton 1er message et d'adapter l'indice:

Code :

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

subX = X(idxA+4:idxB(idxC)-1);

en conséquence.

[mfontan]

Vous êtes des champions!

Merci beaucoup!

A moi de bien apprendre les commandes textes! en tout cas ça me fait bien avancer,

merci encore, c est chouette!