Equation Laplace SQUELETTE AVEC AIDE
Par Ramy • 30 Septembre 2018 • 1 873 Mots (8 Pages) • 428 Vues
...
# Ordonnée du point de départ dans les tableaux
# Si segment horizontal
# Abcisse point d’arrivée (Rq YF=YD)
# remplit le segment dans C avec la valeur 1
# remplit le segment dans G avec la valeur 0
# remplit le segment dans D avec la valeur 0
# remplit le segment dans H avec la valeur 0
# remplit le segment dans B avec la valeur 0
# remplit le segment dans F avec la valeur v
# Si segment vertical
# Ordonnée point d’arrivée (Rq XF=XD)
# remplit le segment dans C avec la valeur 1
# remplit le segment dans G avec la valeur 0
# remplit le segment dans D avec la valeur 0
# remplit le segment dans H avec la valeur 0
# remplit le segment dans B avec la valeur 0
# remplit le segment dans F avec la valeur v
# Si boite
# Abscisse du point d’arrivée
# Ordonnée du point d’arrivée
# remplit la zone dans C avec la valeur 1
# remplit la zone dans G avec la valeur 0
# remplit la zone dans D avec la valeur 0
# remplit la zone dans H avec la valeur 0
# remplit la zone dans B avec la valeur 0
# remplit la zone dans F avec la valeur v
# Sinon erreur de saisie...
return None
#___________________________________________________________________________
# Calcule et renvoie de la norme de la matrice U
def norme(self):
return ()
# Rq : pour les array de numpy A*A réalise le produit terme à terme...
#___________________________________________________________________________
# Une itération de Gauss-Seidel avec sur-relaxation
def gauss_seidel(self,w,umw):
"""
Effectue une itération de Gauss-Seidel avec sur-relaxation
parcours complet sur a matrice U pour modifer ses valeurs
variables :
w : valeur du paramètre de sur-relaxation w
umw : 1-w (pour ne pas le recalculer à chaque case !)
"""
# Pour i et j parcourant toute la zone de travail, implémenter l’éq. 19
return None
#___________________________________________________________________________
# ni itérations de Gauss-Seidel avec sur-relaxation
def iterations(self,ni,w):
"""
Effectue ni itération de Gauss-Seidel avec sur-relaxation
après chaque itération, la norme de U est affichée avec print.
variables :
ni : nombre d’itérations
w : valeur du paramètre de sur-relaxation w
"""
return None
## Mise en oeuvre
# Initialisation
poisson = Poisson(6)
poisson.laplacien()
poisson.valeur_bords(0)
# Condensateur plan : définition des zones de potentiel fixé
# Mise en oeuvre de la résolution
poisson.iterations(100,poisson.omega)
# Affichage du potentiel en niveau de gris
matshow(poisson.U,cmap=’gray’)
# Affichage de lignes équipotentielles
contour(poisson.U,20)
...