matrice en python avec boucles imbriquées

Version imprimable

Bonjour à tous, comme je suis débutante sous python, j'ai rencontré de gros problèmes.
pour construire ma matrice, voici mon code:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
>>>mu0 = 4 * np.pi * 1e-7
>>>alpha=1.4982*10**(-4)
>>>c=0.03097619
>>>k=498.13655
>>>Ms=3.62970e+05
>>>delta =1
>>>H = [0]
>>>Man = [0]
>>>dMirrdH = [0]
>>>Mirr = [0]
>>>M = [0]
>>>B= [0]
>>>for i in range(0,154,1):
              for j in range ( 0,5000,10):
                       Man.append(Ms * (1 / np.tanh((H[j + 1] + alpha * M[j]) / (((119+650*i)/(1+i)))) - ((119+650*i)/(1+i)) / (H[j + 1] + alpha * M[j])))
                       dMirrdH.append((Man[j+1] - M[j]) / (k * delta[j+1] - alpha * (Man[j + 1] - M[j])))
                       Mirr.append(Mirr[j] + dMirrdH[j + 1] * (H[j+1] - H[j]))
                       M.append(c * Man[j + 1] + (1 - c) * Mirr[j + 1])
                       B.append(mu0 * (H[j+1]) + mu0 * M[j+1])
                       X=np.array([[B for j in range(0,5000,10)] for i in range(0,154,1)])
                       print(X)

les problèmes:
1/ je dois avoir une matrice de taille (m,n), alors que en ecrivant X.shape: le résultat est (m,n,n+1)??
2/pour l'incrément j ça avance normalement, mais pour i le résultat est toujours pour le dernier i=153. c'est à dire pour i colonnes, le calcul se fait pour i=153.

Avez vous une idée comment résoudre ce problème.

19/10/2017, 23h57
lg_53
Quand tu utilises des matrices comme tu dis, il faut utiliser la bibliothèque dédié à cela, à savoir numpy. Tu t'embeteras beaucoup moins.

Ensuite tu n'auras qu'à écrire quelque chose du genre :
Code:

1 2 3 4 5 import numpy as np A=np.zeros((1500,2000)) ### avec la forme qui te plait for i in range(1500): for j in range(2000): A[i,j]= i+j ### Mettre ici ce qu'il faut dans ton cas

J'ai essayé avec ce code, mais toujours des erreurs!

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
 
from math import*
import numpy as np
from matplotlib import pyplot as plt
 
a=[]
Nep=10
for j in range(0,Nep):
        a.append((119+650*j)/(1+j))
a=np.array(a)
 
Nfirst = 11 
Ndown =21
Nup=21
 
c=0.03097619
k=498.13655
Ms=3.62970e+05
mu0 = 4 * np.pi * 1e-7 
alpha=1.4982*10**(-4)
DeltaH = 10
 
delta = [0]
Man = np.zeros((54,10))
dMirrdH = []
Mirr = np.zeros((10,54))
M = [0]
B=np.zeros((54,10))
 
H = [0]
 
for j in range(Nfirst):
    H.append(H[j] + DeltaH)
for j in range(Ndown):
    H.append(H[-1] - DeltaH)
for j in range(Nup):
    H.append(H[-1] + DeltaH)
 
H=np.array(H)
 
delta = [0]
for j in range(len(H) - 1):
    if H[j + 1] > H[j]:
        delta.append(1)
    else:
        delta.append(-1)
 
 
for j in range (Nfirst + Ndown + Nup) :
        for i in range (0,Nep):
 
                Man[j,i]= (Ms * (1 / np.tanh((H[j + 1] + alpha * H[j]) / a[i]) - a[i] / (H[j + 1] + alpha * M[j])))
                dMirrdH[j]=((Man[j+1,i] - M[j]) / (k * delta[j+1] - alpha * (Man[j + 1] - M[j])))
                Mirr[j]=(Mirr[j] + dMirrdH[j + 1] * (H[j+1] - H[j]))
                M[j+1]=(c * Man[j + 1,i] + (1 - c) * Mirr[j + 1])
                B[j+1,i]=(mu0 * (H[j+1]) + mu0 * M[j+1])

22/10/2017, 13h41
lg_53

Tu ne construit pas tes tableaux comme il est bon de le faire.

Ce que tu fais : Tu créer une liste, que tu fais grossir petit à petit (avec la fonction append). Enfin tu termines par convertir la grosse liste obtenu en un array.

Ce qu'il faut faire : Oublier que les listes existent, et ne pas utiliser la fonction append. Dès le début déclarer des arrays, qui ont la bonne dimension et que l'on rempli de zéros par défaut. Un fois que tu as ces tableaux, là tu les parcours pour les remplir.

Réécrit ton code sans liste et sans append. Ensuite essaie avec des matrices toutes petites (du genre 5 par 5) pour pouvoir en faire des print, et voir ce qui se passe sur un petit cas. Quand le petit cas fonctionneras, tu pourras essayer avec des dimensions plus grandes.

Ensuite plonge toi un peu dans un cours ou un tuto, car visiblement, quand tu écris

Code:

X=np.array([[B for j in range(0,5000,10)] for i in range(0,154,1)])

tu ne comprends pas ce que tu fais...
24/10/2017, 10h06
lolita2017

merci.