fichier csv et python

Version imprimable

Bonjour
Mon problème est simple lorsque je lance ce programme:

Code:

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
import matplotlib.pyplot as plt ; import numpy as n


#### Coordonnées des points en fonction du temps  ####
fic = open("Mouvement_NH.csv",'r')                #
donnee=fic.readlines()                               #
x=[] ; y=[] ; dt = 10*86400 # 10 jours entre 2 points#
for point in donnee :                                #
    point=point.split() # chaque point est un triplet#
    if point[0][0].isnumeric():  #
        x.append(float(point[1].replace(',','.')))   #
        y.append(float(point[2].replace(',','.')))   #
######################################################
###  Calcul de vx, vy et v pour le vecteur vitesse ###
vx=[] ; vy=[] ; v=[]                               #
for i in range(1,len(x)-1):                          #
    vx.append((x[i+1]-x[i-1])/(2*dt))                #
    vy.append((y[i+1]-y[i-1])/(2*dt))                #
    v.append((vx[-1]**2+vy[-1]**2)**0.5)             #
######################################################
##  Calcul de ax et ay pour le vecteur accélération ##
ax=[] ; ay=[]  ; a=[]                                #
for i in range(1,len(vx)-1):                         #
    ax.append((vx[i+1]-vx[i-1])/(2*dt))              #
    ay.append((vy[i+1]-vy[i-1])/(2*dt))              #
    a.append((ax[-1]**2+ay[-1]**2)**0.5)             #
######################################################


########  Tracé des champs de vecteurs:  #############
vec_v=plt.quiver(x[1:-1],y[1:-1],vx,vy,color="green")#
plt.quiverkey(vec_v,X=0.7, Y=1.05, U=5e4,label='V')  #
                                                     #
vec_a=plt.quiver(x[2:-2],y[2:-2],ax,ay,color="blue") #
plt.quiverkey(vec_a,X=0.2, Y=1.05, U=1e-2,label='a') #
                                                     #
plt.plot(x,y,"ro") ; plt.axis('equal')               #
plt.show()                                           #
######################################################
###########     Traitement statistique     ###########
v = n.array(v) ; a = n.array(a) # données sous numpy #
v_moy=n.mean(v) ; a_moy = n.mean(a) # Moyennes       #
k=a_moy*149e9/v_moy**2                               #
####  Calcul des incertitudes U sur a, v et k : ######
U_v=n.std(v,ddof=1)/(len(vx)**0.5)                   #
U_a=n.std(a,ddof=1)/(len(a)**0.5)                    #
U_k = (k*((2*U_v/v_moy)**2+(U_a/a_moy)**2)**0.5)     #
#### Affichage des résultats statistiques :    #######
plt.text(-0.5e11,0,
    "V = ("+str(int(v_moy/1e2)/10) +chr(177)
    +str(int(U_v/1e2)/10)+")km/s\n"+


    "a = ("+ str(int(a_moy*1e5)/1e2)+chr(177)
    +str(int(U_a*1e5)/1e2)+')mm/s²\n'+


    "k = ("+ str(int(k*100)/100)+chr(177)
    +str(int(U_k*100)/100)+")")

Python ne semble pas reconnaître mes variables x et y (ligne en rouge ci-dessus). j'ai pourtant nommé dans excel (et regressi) les colonnes comme ils se doivent mais quelques chose (de bête, certainement) doit m'échapper. Je précise que que je suis novice en la matière.
Merci par avance

04/08/2020, 07h42
tyrtamos

Bonjour,

Il faudrait avoir un exemple de la variable 'point'. Si elle est une liste de 3 valeurs numériques, la 1ère est atteinte par 'point[0]' et non 'point[0][0]'. Mais si cette 1ère valeur est numérique, le 'point[0][0]' devrait générer une erreur. Si c'est une chaîne, on obtient son 1er caractère qui n'est pas numérique. Il faudrait donc obtenir une valeur numérique par conversion en entier ou flottant, et modifier la ligne:

Code:

if point[0].isnumeric():

Bonjour et merci pour l'aide.
j'ai donc remplacé comme suggéré mais maintenant, rien ne s'affiche dans le graphe. Du coup je me demande si mon fichier excel converti en csv est bien fait? J'ai trois colonnes, toutes numériques x, y et dt:

x y dt

m m s

0,001057785 -0,00035327 0

-0,042642379 -0,04996529 0,04

-0,058937355 -0,077362973 0,08

-0,071528927 -0,10698209 0,12

-0,078935734 -0,13882264 0,16

-0,082639138 -0,170663191 0,2

-0,083379819 -0,203984697 0,24

-0,077454373 -0,234344291 0,28

-0,070047566 -0,266184842 0,32

-0,05449327 -0,296544436 0,36

-0,039679656 -0,322461163 0,4

-0,019681276 -0,346156457 0,44

0,003279827 -0,369111272 0,48

0,029203652 -0,38762322 0,52

0,058090201 -0,403173256 0,56

0,086236069 -0,412799469 0,6

0,115863298 -0,420204248 0,64

0,14771257 -0,422425682 0,68

0,179561841 -0,420944726 0,72

0,20844839 -0,415761381 0,76

0,238075619 -0,406875646 0,8

0,265480806 -0,392806566 0,84

0,291404632 -0,377997007 0,88

0,314365734 -0,359485059 0,92

0,335845475 -0,338751678 0,96

0,352881132 -0,314315906 1

0,368435428 -0,289139657 1,04

0,378064277 -0,261001496 1,08

0,386952446 -0,23138238 1,12

0,39139653 -0,202503741 1,16

0,39139653 -0,172884624 1,2

0,388433807 -0,143265508 1,24

0,383249042 -0,115867825 1,28

0,372138831 -0,088470142 1,32

0,359547259 -0,062553415 1,36

0,343252283 -0,038858121 1,4

0,324735265 -0,017384262 1,44

0,302514843 0,00334912 1,48

0,278813059 0,020380112 1,52

0,253629914 0,033708714 1,56

0,226965408 0,043334927 1,6

0,198078859 0,051480185 1,64

0,16845163 0,055923052 1,68

0,141046443 0,055923052 1,72

merci

04/08/2020, 19h59
wiztricks

Salut,

Citation:

Envoyé par fmarioli

mais maintenant, rien ne s'affiche dans le graphe. Du coup je me demande si mon fichier excel converti en csv est bien fait? J'ai trois colonnes, toutes numériques x, y et dt:

Pour vérifier que vous retrouvez bien dans les listes x et y ce que vous attendez, vous pouvez faire print(x), print(y) et comparer.

- W
15/08/2020, 22h15
Archimède

1 pièce(s) jointe(s)
Salut, ta première colonne est apparemment celle des x et la seconde, celle des y...
donc il faudrait faire plutôt :
Code:

1 2 3 x.append(float(point[0].replace(',','.'))) # y.append(float(point[1].replace(',','.'))) #
Avec tes valeurs :
Pièce jointe 577339

Tes vecteurs accélération sont presque centripètes...

x	y	dt
m	m	s
0,001057785	-0,00035327	0
-0,042642379	-0,04996529	0,04
-0,058937355	-0,077362973	0,08
-0,071528927	-0,10698209	0,12
-0,078935734	-0,13882264	0,16
-0,082639138	-0,170663191	0,2
-0,083379819	-0,203984697	0,24
-0,077454373	-0,234344291	0,28
-0,070047566	-0,266184842	0,32
-0,05449327	-0,296544436	0,36
-0,039679656	-0,322461163	0,4
-0,019681276	-0,346156457	0,44
0,003279827	-0,369111272	0,48
0,029203652	-0,38762322	0,52
0,058090201	-0,403173256	0,56
0,086236069	-0,412799469	0,6
0,115863298	-0,420204248	0,64
0,14771257	-0,422425682	0,68
0,179561841	-0,420944726	0,72
0,20844839	-0,415761381	0,76
0,238075619	-0,406875646	0,8
0,265480806	-0,392806566	0,84
0,291404632	-0,377997007	0,88
0,314365734	-0,359485059	0,92
0,335845475	-0,338751678	0,96
0,352881132	-0,314315906	1
0,368435428	-0,289139657	1,04
0,378064277	-0,261001496	1,08
0,386952446	-0,23138238	1,12
0,39139653	-0,202503741	1,16
0,39139653	-0,172884624	1,2
0,388433807	-0,143265508	1,24
0,383249042	-0,115867825	1,28
0,372138831	-0,088470142	1,32
0,359547259	-0,062553415	1,36
0,343252283	-0,038858121	1,4
0,324735265	-0,017384262	1,44
0,302514843	0,00334912	1,48
0,278813059	0,020380112	1,52
0,253629914	0,033708714	1,56
0,226965408	0,043334927	1,6
0,198078859	0,051480185	1,64
0,16845163	0,055923052	1,68
0,141046443	0,055923052	1,72

Pour ton problème d'affichage :

Code:

1
2
3
4
5
6
7
8
9
10
11
 
fic = open("orbite.csv",'r')
donnee=fic.readlines()
dt = 10*86400 # 10 jours entre 2 points#
x,y=[],[]
 
for point in donnee :
    point=point.replace('\n','').replace('\r','')
    point=point.split(';')
    x.append(float(point[0].replace(',','.')))
    y.append(float(point[1].replace(',','.')))

Ceci dit la méthode présentée dans le Bordas n'est pas la plus compacte !

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
 
import matplotlib.pyplot as plt ; import numpy as n
import csv, os, sys
 
os.chdir(sys.path[0])# pour pointer dans le dossier du fichier source
#### Coordonnées des points en fonction du temps  ####
fic = open("orbite.csv",'r')                #
x=[] ; y=[] ; dt = 10*86400 # 10 jours entre 2 points#
 
reader = csv.reader(fic,delimiter=';')
donnees = list( zip(*reader) )
x= [ float( donnees[0][i].replace(',','.') ) for i in range( len(donnees[0])) ] #on garde la colonne n°1(indice 0) des x
y= [ float( donnees[1][i].replace(',','.') ) for i in range(len(donnees[1])) ] #on garde la colonne n°2(indice 1) des y

Bonne deuxième loi de Newton et bon courage pour Kepler !

Rq: tu as un problème à la fin de ton code

Code:

1
2
3
4
5
6
7
8
9
10
11
12
 
'''plt.text(-0.5e11,0,
    "V = ("+str(int(v_moy/1e2)/10) +chr(177)
    +str(int(U_v/1e2)/10)+")km/s\n"+
 
 
    "a = ("+ str(int(a_moy*1e5)/1e2)+chr(177)
    +str(int(U_a*1e5)/1e2)+')mm/s²\n'+
 
 
    "k = ("+ str(int(k*100)/100)+chr(177)
    +str(int(U_k*100)/100)+")")'''

Je n'ai pas eu le temps de regarder...