Discussion :

[myaxto]

Bonjour, j'aurais besoin d'un peu d'aide sur un problème sans doute un peu nul.
Il y a deux choses:
- la première c'est que je n'arrive pas à réaliser un calcul à partir d'une valeur d'une liste (ligne 39) déterminée avant
- la deuxième chose est que je voudrais stocker les données calculées Pcum et stockage dans un fichier externe mais n'ai aucune idée de comment faire.
Voila mon programme:

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
 
import csv
import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import quad
 
energie = [] #creation des listes
angle = []
section = []
Pcum = []
stockage = []
 
fichier = csv.reader(open("donnees-energie5.csv","r"))
fichier2 = open("stockage2.txt","a")
 
for row in fichier: #transferts des donnees dans les listes
    energie.append(row[0])
    angle.append(row[1])
    section.append(row[2])
 
plt.plot(angle,section)
 
i = 0
C = 0
 
angle = [float(x) for x in angle] #transformation en reels
section = [float(y) for y in section]
 
def fonction(x,y): #fonction a integrer
    return np.sin(x)*y
 
x = angle[i]
y = section [i]
 
while C<10: #integration pour toutes les valeurs
    I = quad(fonction,0,angle[i],args=(y))
    stockage.insert(i,I)
    C = C+1
    i = i+1
    Pcum[i] = stockage[i]*2*np.pi

[VinsS]

Salut,

Je suppose que tu as une erreur "out of range".

Le nom i est simplement un index et tu l'incrémentes avant d'utiliser la valeur à cet index, qui n'existe pas.

Code :

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

1
2
3
4
5
    stockage.insert(i,I)
    C = C+1
    i = i+1
    Pcum[i] = stockage[i]*2*np.pi

Ton code devrait se résumer à ça:

Code :

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

1
2
3
4
5
6
 
while C < 10:
    I = quad(fonction,0,angle[i],args=(y))
    stockage.append(I)
    Pcum.append(stockage[-1]*2*np.pi)
    C = C+1

[lg_53]

Pour écrire un tableau numpy dans un fichier tu peux te servir de np.savetxt(filename)

https://docs.scipy.org/doc/numpy-1.1...y.savetxt.html

et il est réimportable avec la commande np.loadtxt(filename)

[myaxto]

Merci beaucoup pour ta réponse, c'était bien ça mais impossible de trouver l'erreur tout seul.
Pourrais juste tu m'aider pour le second point qui est l'écriture dans un autre fichier.
Je voudrais copier la liste stockage dans un fichier txt mais je n'y arrive pas, python me répond que il ne peux pas transférer une liste en écriture (c'est peut-être pas très claire je peux me réexpliquer).

[myaxto]

En fait j'ai toujours le même problème : can't multiply sequence by non-int of type 'float' pour le calcul de Pcum. Et le coup de l'index qui ne marche pas c'est parce que j'utilisais l'index après l'avoir changé c'est ca ?
Tu aurais pas un exemple pour la fonction savetesxt parce que j'avoue que pour un débutant je ne comprend pas ce que dit la doc.

[lg_53]

En fait j'ai toujours le même problème : can't multiply sequence by non-int of type 'float' pour le calcul de Pcum. Et le coup de l'index qui ne marche pas c'est parce que j'utilisais l'index après l'avoir changé c'est ca ?
Tu aurais pas un exemple pour la fonction savetesxt parce que j'avoue que pour un débutant je ne comprend pas ce que dit la doc.

Il y a un exemple au bas du lien que je t'ai donné ...

[myaxto]

Il y a un exemple au bas du lien que je t'ai donné ...

J'ai bien vu l'exemple mais j'ai toujours le problçme qu'il me faudrait les écrire en colonne et que je n'arrive seulement à le faire en ligne.
Voila le petit test de programme que j'ai fait:

Code :

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

1
2
3
4
import numpy as np
M = np.array([[4,7,2,9,1,3,2,7,1]])
print(M)
np.savetxt('test.txt',M, fmt='%10.5f',newline='\n',)

Et avec ce programme:

Code :

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

1
2
3
4
5
6
while C<10: #integration pour toutes les valeurs
    I = quad(fonction,0,angle[i],args=(y))
    stockage.append(I)
    Pcum.append(stockage[i]*2*np.pi)
    C = C+1
    i = i+1

J'obtient toujours la même erreur. De plus pourquoi ce n'est pas stockage[i] dans ton changement ?
Merci pour ta patience et très réponses précédentes.