[traitement de nombres] Comment pourrais je faire ça sous Linux ?

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : [traitement de nombres] Comment pourrais je faire ça sous Linux ?

Mjules

Parle dans le vide

Bonsoir,

voilà le pb :
j'ai plusieurs fichiers texte (codage ASCII) contenant 2 colonnes de nombres. Ce sont des fichiers décrivant des courbes de diffraction X.
Ma première colonne contient des angles en degrées.
Je cherche à appliquer cette formule sur ces angles (ça me permettra de comparer des diagrammes fait dans différents appareils) :

angle_final (en degré) = arcsin{(lambda'/lambda)*sin(angle_initial)}

j'ai commencé par faire ça avec un tableur mais je suis vite limité (entre 3500 et 6000 valeurs à chaque fois, ça devient vite lourd).
Je cherches donc un moyen d'automatiser ces opérations avec un script/programme.

L'ennui c'est que je ne sais pas quoi utiliser ; j'ai quelques notions de BASH (je suis sous linux) et j'avais pensé à awk, mais celui ci ne semble pas faire les arcsin.

je m'en remets donc à vous pour me donner des pistes ou des noms de softs ou bien des liens. Enfin bref tout ce qui pourrait m'aider.

Merci d'avance.

---------------
Celui qui pose une question est idiot 5 minutes. Celui qui n'en pose pas le reste toute sa vie. | Membre du grand complot pharmaceutico-médico-scientifico-judéo-maçonnique.

Publicité

uriel

blood pt.2

bon je dirais pour des calculs mathematique le Fortran est ce qu'il y a de mieux, mais bon pas top , avec du Perl tu devrais t'en sortir :jap:

---------------
IVG en france

Mjules

Parle dans le vide

et c'est pas trop compliqué à apprendre le perl ? tu aurais des bons liens ?

HS : au fait, tu t'en es sortie avec ton DHCP ?

Message édité par Mjules le 03-02-2004 à 19:41:26

---------------
Celui qui pose une question est idiot 5 minutes. Celui qui n'en pose pas le reste toute sa vie. | Membre du grand complot pharmaceutico-médico-scientifico-judéo-maçonnique.

uriel

blood pt.2

Mjules a écrit :

et c'est pas trop compliqué à apprendre le perl ? tu aurais des bons liens ?

HS : au fait, tu t'en es sortie avec ton DHCP ?

Pour Perl : http://forum.hardware.fr/forum2.ph [...] 097&cat=10

compliqué... la syntaxe est un peu special mais c'est pas trop compliqué, maintenant la plupart de langages te permette de faire ce que tu veux, je prends Perl parce que je le connais un peu [:spamafote]

HS: :jap: pas encore verifier vu que j'ai pas la machine sous la main, mais je teste et je réup le sujet

---------------
IVG en france

Evadream -jbd-

Salut, je comprends pas trop
Tu as des fichiers texte contenant 2 colonnes de nombres, la première colonne réprésente des angles en dégré, et la deuxième, on s'en fiche ?

C'est "angle_initial" que l'on doit remplacer dans ta formule ?

@+

Evadream -jbd-

Bon, je viens de me mettre au python y'a 10 minutes [:ddr555]

Pour un fichier essai.tmp de cette forme :

1.2 5.8
3.9 2.8
1.2 5.1

Et un script python qui ressemble à çà (indulgence demandé)

Code :

#!/usr/bin/python
# -*- coding: utf-8 -*-
import math
# angle_final (en degré) = arcsin{(lambda'/lambda)*sin(angle_initial)}
fichier = open('essai.tmp','r')
a = 1.9
b = 2.47
# on lit tout le fichier qu'on place dans lignes, ligne par ligne
lignes = fichier.readlines()
# on lit ligne par ligne
for ligne in lignes :
# split sépare les élements de la ligne (par défaut, l'espace est l'élément séparateur
donnees = ligne.split()
if len(donnees) == 2 :
s = "Angle initial:" + donnees[0] + " "
p = "Angle final:%f" % math.asin((a/b)*math.sin(float(donnees[0])))
print s+p
fichier.close()

J'obtiens la sortie suivante :

$ ./files.py
Angle initial:1.2 Angle final:0.799422
Angle initial:3.9 Angle final:-0.557482
Angle initial:1.2 Angle final:0.799422
$

Ca à l'air pas trop mal, à moins que je réponde à coté de la plaque

Enfin tout ca pour te dire que ca semble pas extremement complexe, d'autant plus que la doc python est super bien faite !

@++

uriel

blood pt.2

c'est vrai que Python est peut être plus lisible et facile d'acces maintenant que tu en parles

---------------
IVG en france

Mjules

Parle dans le vide

Evadream -jbd- a écrit :

oui exactement :jap:

---------------
Celui qui pose une question est idiot 5 minutes. Celui qui n'en pose pas le reste toute sa vie. | Membre du grand complot pharmaceutico-médico-scientifico-judéo-maçonnique.

Mjules

Parle dans le vide

Evadream -jbd- a écrit :

Bon, je viens de me mettre au python y'a 10 minutes [:ddr555]

Pour un fichier essai.tmp de cette forme :

1.2 5.8
3.9 2.8
1.2 5.1

Et un script python qui ressemble à çà (indulgence demandé)

Code :

#!/usr/bin/python
# -*- coding: utf-8 -*-
import math
# angle_final (en degré) = arcsin{(lambda'/lambda)*sin(angle_initial)}
fichier = open('essai.tmp','r')
a = 1.9
b = 2.47
# on lit tout le fichier qu'on place dans lignes, ligne par ligne
lignes = fichier.readlines()
# on lit ligne par ligne
for ligne in lignes :
# split sépare les élements de la ligne (par défaut, l'espace est l'élément séparateur
donnees = ligne.split()
if len(donnees) == 2 :
s = "Angle initial:" + donnees[0] + " "
p = "Angle final:%f" % math.asin((a/b)*math.sin(float(donnees[0])))
print s+p
fichier.close()

J'obtiens la sortie suivante :

$ ./files.py
Angle initial:1.2 Angle final:0.799422
Angle initial:3.9 Angle final:-0.557482
Angle initial:1.2 Angle final:0.799422
$

ouaouh !
c'est ça qu'il me faut :love: sauf que moi je veux avoir en sortie un fichier avec colonne1 modifié et colonne2 intact.

bon je vais regarder python
Merci

---------------
Celui qui pose une question est idiot 5 minutes. Celui qui n'en pose pas le reste toute sa vie. | Membre du grand complot pharmaceutico-médico-scientifico-judéo-maçonnique.

chaica

Be gorgeous, use python.

Publicité

Evadream -jbd-

Mjules a écrit :

ouaouh !
c'est ça qu'il me faut :love: sauf que moi je veux avoir en sortie un fichier avec colonne1 modifié et colonne2 intact.

bon je vais regarder python
Merci

C'est cool
Si tu veux pas te prendre la tete, tu remplaces :

if len(donnees) == 2 :
s = "Angle initial:" + donnees[0] + " "
p = "Angle final:%f" % math.asin((a/b)*math.sin(float(donnees[0])))
print s+p

Par

if len(donnees) == 2 :
s = "%f " % math.asin((a/b)*math.sin(float(donnees[0]))) + donnees[1]
print s

Et tu fais :

./nomduscript > nouveaufichier

et hop

Message édité par Evadream -jbd- le 03-02-2004 à 21:49:42

Mjules

Parle dans le vide

Evadream -jbd- a écrit :

C'est cool
Si tu veux pas te prendre la tete, tu remplaces :

if len(donnees) == 2 :
s = "Angle initial:" + donnees[0] + " "
p = "Angle final:%f" % math.asin((a/b)*math.sin(float(donnees[0])))
print s+p

Par

if len(donnees) == 2 :
s = "%f " % math.asin((a/b)*math.sin(float(donnees[0]))) + donnees[1]
print s

Et tu fais :

./nomduscript > nouveaufichier

et hop

La grande classe :wahoo:

Merci 1000 fois :jap:

---------------
Celui qui pose une question est idiot 5 minutes. Celui qui n'en pose pas le reste toute sa vie. | Membre du grand complot pharmaceutico-médico-scientifico-judéo-maçonnique.

Taz

bisounours-codeur

fichier = open('essai.tmp','r')

a = 1.9
b = 2.47

# on lit tout le fichier qu'on place dans lignes, ligne par ligne
lignes = fichier.readlines()

# on lit ligne par ligne

t'emmerdes pas

for ligne in open('fichier'):
# ...

Evadream -jbd-

Ah ok ! Très bien ! Merci !

Evadream -jbd-

Mjules a écrit :

La grande classe :wahoo:
Merci 1000 fois :jap:

Pas de soucis, bon courage !

@+

Taz

bisounours-codeur

juste une réflexion

le % prends en argument droit un tuple. python fait automatiquement de packing, unpacking
comme dans

f, = sequence
ou
pair = first, second

donc "..." % arg fonctionne
seulement, je pense que ce n'est pas une bonne habitude, et ça peut se révéler problématique si arg est lui même un tuple

donc je recommande d'écrire

"..." % (arg,)

Evadream -jbd-

Je note ! Merci pour cette remarque

Message édité par Evadream -jbd- le 03-02-2004 à 22:55:13

Taz

bisounours-codeur

eih j'espère que t'as surtout compris que pas besoin de fermer explicitement les fichiers en python, ça se fait tout seul à la destruction (plus de références sur les descripteur de fichiers)

Evadream -jbd-

Ta remarque s'applique même lorsque qu'on fait un fichier = open('monnom'), ou bien c'était pour appuyer ta remarque précédente où le open se faisait dans le for ?

Message édité par Evadream -jbd- le 03-02-2004 à 23:02:26

Taz

bisounours-codeur

oui. en règle générale, ne prends pas la peine de fermer un fichier. le seul cas ou tu le fermes explicitement, c'est quand tu veux faire une autre opération (réouverture dans un autre mode, suppression, déplacement, renommage, etc) afin d'assurer la cohérence

Evadream -jbd-

Très bien, encore merci

chaica

taz a écrit :

Rapidement :

http://www.python.org/doc/current/tut/node9.html

String and Unicode objects have one unique built-in operation: the % operator (modulo). This is also known as the string formatting or interpolation operator. Given format % values (where format is a string or Unicode object), % conversion specifications in format are replaced with zero or more elements of values. The effect is similar to the using sprintf() in the C language. If format is a Unicode object, or if any of the objects being converted using the %s conversion are Unicode objects, the result will also be a Unicode object.

If format requires a single argument, values may be a single non-tuple object. 2.8 Otherwise, values must be a tuple with exactly the number of items specified by the format string, or a single mapping object (for example, a dictionary).

Vraiment pour l'habitude alors

Message édité par chaica le 03-02-2004 à 23:21:22

Mjules

Parle dans le vide

Bonjour,

alors voilà, après 2 jours d'apprentissage, j'ai finalement réussi à pondre ça :

Code :

#!/usr/bin/python
# -*- coding: utf-8 -*-
# Ce petit programme permet de convertir des angles obtenu par diffraction X avec une anticathode donnée
# en les angles qu'on obtiendrais avec une autre anticathode
# pour me contacter : mjules chez ifrance point com
# v 0.1
#
#NOTE : le fichier d'entrée doit être composé de 2 colonnes angles et intensités, dans cet ordre et séparé par un ou plusieurs espaces
import math
# angle_final (en degr?) = arcsin{(lambda'/lambda)*sin(angle_initial)}
fichier_anal = raw_input('Entrer le nom du fichier à analyser :') # on entre le fichier à analyser
try: # on vérifie que ce fichier existe
f = open(fichier_anal, 'r')
f.close()
except:
print "le fichier", fichier_anal, "est introuvable"
cath1 = raw_input('Entrer l\'élément utilisé pour l\'analyse (Co, Cr, Cu) : ') # on entre le métal de l'anticathode d'analyse
if cath1 == 'Co' : # on assigne la valeur de la raie à la variable utilisée dans l'équation
anticath_anal = 1.788
elif cath1 == 'Cu' :
anticath_anal = 1.541
elif cath1 == 'Cr' :
anticath_anal = 2.2798
else :
print "Les autres éléments ne sont pas encore supporté, désolé" # si l'élément est inconnu, on le dit
cath2 = raw_input('Entrer l\'élément vers lequel les angles doivent être convertis (Co, Cr, Cu) : ') #idem pour l'anticathode vers laquelle on converti
if cath2 == 'Co' :
anticath_result = 1.788
elif cath2 == 'Cu' :
anticath_result = 1.541
elif cath2 == 'Cr' :
anticath_result = 2.2798
else :
print "Les autres éléments ne sont pas encore supporté, désolé"
# on lit tout le fichier qu'on place dans lignes, ligne par ligne
fichier_result = fichier_anal + cath1 + "vers" + cath2 # on donne un nom au fichier de résultat
# fichier original suivi des 2 métaux utilisés
result = open(fichier_result, 'w') # on crée le fichier de résultat
# on lit ligne par ligne
for ligne in open(fichier_anal, 'r') :
# split s?pare les ?lements de la ligne (par d?faut, l'espace est l'?l?ment s?parateur
donnees = ligne.split()
if len(donnees) == 2 :
pre_resultat = (anticath_result/anticath_anal)*math.sin(float(donnees[0])*(math.pi/360)) # si le résultat du sinus de l'angle multiplié par le rapport
if (pre_resultat > 1) : # est supérieur à 1, l'arcsinus risque de faire la tête
# alors on met la valeur 0
print "certains angles sont trop grands pour la conversion, ils seront remplacés par des 0"
s = "0 " + donnees[1] + "\n"
else :
s = "%f " %((360/math.pi)*math.asin(float(pre_resultat)),) + donnees[1] + "\n" # on écrit le résultat dans le fichier de sortie
result.write(str(s))

ça marche pas trop mal mais ya encore du boulot dessus à faire.

Merci beaucoup pour m'avoir donné les pistes et à bientôt sur la subcat python parce que je sens que je vais pas tarder à avoir besoin d'un peu d'aide

A+

---------------
Celui qui pose une question est idiot 5 minutes. Celui qui n'en pose pas le reste toute sa vie. | Membre du grand complot pharmaceutico-médico-scientifico-judéo-maçonnique.

Publicité

FORUM HardWare.fr

Programmation

Divers

[traitement de nombres] Comment pourrais je faire ça sous Linux ?

Sujets relatifs
Adaptations entre Windows et Linux	[VB .NET] Traitement vba trop lent ??
configurer MySQL et PHP sur Linux	[ ASM INTEL LINUX ] - Besoin d'aide pour comprendre qques trucs Svp
portage vers linux	Un IDE C# sous linux qui soit bien?
Détecter si l'OS est Win ou linux...	les nombres complexes en pascal
Remplir un tableau de nombres flottants avec une certaine valeur	Windows : timer - coordonnées / Borland Windows <==> Linux
Plus de sujets relatifs à : [traitement de nombres] Comment pourrais je faire ça sous Linux ?

Page générée en 0.154 secondes