Discussion :

[oomph]

Bonjour,

je cherche un format de nombre equivalent à celui appelé F9 en fortran.
celui-ci donne un nombre ecrit si il est positif:
0.1234E02 (le zero est important)

soit si il est negatif:
-.1234E02

(ici le format en fortran c'est F9.4)

Je n'ai pas trouvé un tel format sous matlab (j'ai essayé de bidouiller avec le format '.4e' mais ça n'a pas donné grand chose...)

Avez vous des idées pour pouvoir récuperer ce format de nombre?

Merci de votre participation

[Jerome Briot]

Je ne vois pas où se situe le problème mais ta question n'est pas très claire

Que comptes-tu faire avec ce format ?

Par exemple, à la lecture :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>> sscanf('0.1234E02','%f')
 
ans =
 
   12.3400
 
>> sscanf('-.1234E02','%f')
 
ans =
 
  -12.3400

[oomph]

ah oui j'ai pas donné le contexte, du coup c'est pas tres claire...

Je veux créer des fichiers contenant des plages de valeurs, à partir d'un autre fichier contenant ses valeurs, mais sous un autre format. Actuellement je récupère mes valeurs comme char dans le premier fichier (format du nombre: -1.593667828E-002), et je veux les transposer dans mon nouveau fichier sous le format: -.159367E-01

Pour cela j'utilise la fonction str2num, puis num2str pour le formater comme je le veux... (et c'est là que ça coince!)

Es ce plus claire? (je suis pas convaincue de mon explication...)
c'est pour preparer un calcul dans un logiciel.

[Jerome Briot]

Je ne connais pas le Fortran mais ne serait-ce pas le format E.4 plutôt que F.4 que tu cherches à convertir ?

Essaie quand même ceci :

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
function str = conv2E4(x)
 
if x==0
   str = sprintf('%#.4E',x);
   return
end
 
p = floor(log10(abs(x)));
 
str = sprintf('%.4E',x);
 
if x>0
 
    str = sprintf('0.%c%sE%c%02d',str(1),str(3:5),44-sign(p),abs(p+1));
 
else
 
    str = sprintf('-.%c%sE%c%02d',str(2),str(4:6),44-sign(p),abs(p+1));
 
end
 
if ispc
   str(end-2) = []; 
end

Exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
x = 1.593667828E-002;
conv2E4(x)
 
x = 12.34;
conv2E4(x)
 
x = -1.593667828E-002;
conv2E4(x)
 
x = -12.34;
conv2E4(x)
 
x = 0;
conv2E4(x)

donne :

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
ans =
 
0.1593E-01
 
 
ans =
 
0.1234E+02
 
 
ans =
 
-.1593E-01
 
 
ans =
 
-.1234E+02
 
 
ans =
 
0.0000E+00

Une autre solution plus élégante serait de passer par un fichier MEX Fortran pour faire la conversion "numérique => chaîne de caractères" avec le vrai format E.4

[oomph]

yes ça marche merci!!!

quelques petites modifs pour l'incorporer dans mon code et c'est nikel!!


Pour info en Fortran mes lignes de format etaient:

Code :

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

1
2
011	FORMAT(F9.6,'	 ',F9.4,'	 ',F9.4,'		 ',F12.6)
012	FORMAT(A18,' Poste',I1)

donc j'imagine que F9.4 et E9.4 doivent correspondre...

Je ne passe pas par un fichier mex, car je dois envoyer des infos au fortran, et il n'est pas ecrit de façon à ce que ce soit possible... (il attend une entrée clavier, et je lui envoi des arguments...)

Bref merci beaucoup pour ton aide!