Discussion :

[lolita2017]

Bonjour à tous,
Comme je suis débutante en python, j'ai rencontré des problèmes à résoudre mon problème qui réside dans le fait de construire une matrice B en fonction de deux vecteurs H et a.
Le problème c'est que la formulation de B est complexe et j'ai pas su comment gérer les listes, array et matrice.
Donc voici mon code:

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
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])

merci.

[raph_prog]

Bonjour,

Je pense savoir pourquoi ton programme ne fonctionne pas : tout d'abord, il faut savoir que un "vecteur" en math n'est pas un type "vecteur" en Python (c'est vrai pour beaucoup d'autre langage), et que une "matrice" en math, c'est un "tableau" en Python (valable pour d'autres langage aussi). Par contre un "vecteur" en math, est un "tableau" en Python. Le "tableau" en python est directement lier à un array.
Pour illustrer voici un exemple :

vecteur_math=(1 2 3) --> tableau_python[3]= [1,2,3] (aussi égal à une list[3])
matrice_math=((1 2 3) (1 2 3)) --> tableaux_python[3][2]=[[1,2,3];[1,2,3]]

Voila avec cet exemple on peut voir que ton vecteur est un array (ou list) unidimensionnel, et ta matrice un array multidimensionnel.
Donc rien qu'avec la classe array (si tu ne sait pas ce qu'est une classe dit toi que c'est une grosse bibliothèque) tu peut résoudre ton problème (lire la doc est un bon exercice).

Si tu veux y arriver autrement, tu peut mais ce sera plus long et plus compliqué à déboguer.

[lolita2017]

Bonjour Raph,
merci pour l'explication. Cependant, j'ai pas pu écrire correctement mon code.
Pour le résultat, au lieu d'avoir une matrice B, j'ai une liste d'array au nombre (Nfirst+Ndown+Nup) et chaque array est de longueur Nep.
Que dois-je faire pour corriger le problème?

[lg_53]

C'est pas la 3eme fois que tu poste ce même problème sur ce forum ? ...

Fais tout avec des arrays, comme tu l'as fait pour Man, Mirr, etc ...
Là delta est encore une liste, H aussi, a également.
Débarasse toi de TOUS les append et de TOUTES les listes. S'il en reste c'est que c'est pas comme ca qu'il faut faire.
Quand on a des problèmes de dimension de tableaux, c'est plus facile de les résoudre si on fixe la taille du tableau dès le début...

[lolita2017]

Je pense que j'ai réussi à résoudre le problème.
J'ai transformé le tout en matrice.
tous les vecteurs sont initialisés à des matrices de zéros.

merci