Discussion :

[oreldu71]

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
41
42
43
44
45
46
47
48
49
50
from pylab import*


M=loadtxt('data.txt')
print M
taille =shape (M)
nb= size (M)
print "taille de la matrice", taille
print "nombre d'element", nb

def moyg (x): 
        N=len(x)
        print N
        M=zeros((N-2,))

        for i in range(0,N-2):
            a=array([1,2,1])
            l=(x[i,1])*a[0]
            s=(x[(i+1),1])*a[1]
            s1=(x[(i+2),1])*a[2]
            S=s+s1+l
            M[i]=S/3
          
        return M
        

Tmoy=moyg(M)
print Tmoy
t=M[0:141,0]
x=t
y=Tmoy
xlabel("Temps ( minutes)")
ylabel("Temperature (C)")
title( " Moyenne Ponderee ")  
plot(x,y)
show ()


f=M[:,1]
f1=arange(0,len(f))
h=[1,2,1]

c=zeros (M.shape)
for i in range (0,len(f)-2):
    h=array ([1,2,1])
    h1[i]=h[0]*(x[i])
    h2[i]=h[1]*(x[i,1])
    h3[i]=h[2]*(x[i,2])
    H=(h1+h2+h3)/3
print H

[lg_53]

Ca fait quasiment la meme chose que la fonction moyg. Donc si tu comprends ce que fait moyg, tu dois aussi comprendre ce que fait le bout de code que tu as mis en rouge.

Note quand même que quand tu écris :

Code :

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

1
2
3
4
for i in range (0,len(f)-2):
    ### blabla
    H=(h1+h2+h3)/3
print H

et bien là la boucle for ne sert à rien pour H. Car H aura forcément la valeur calculée lors du dernier passage dans la boucle (à chaque itération tu écrase la valeur de H, pour y stocker la nouvelle valeur)

Ensuite là tu utilises des fonctionalités numpy, qui est construit pour éviter les boucles for grace au slice. Du style

Code :

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

1
2
for i in range(0,N-2):
            M[i]=( x[(i+1),1]+x[(i+2),1] )/3

va plutot s'écrire directement

Code :

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

M=( x[1:(N-1),1]+x[2:N,1] )/3

C'est un fondamental de numpy ca. Tu devrais te documenter un peu dessus, ca répondrait surement à plusieurs de tes interrogations...

[oreldu71]

Merci j'ai bien compris