Probleme champ de vecteurs
Bonjour,
je suis nouveau et etudiant en classe Prépa.
Je désire modéliser un champ de vecteur vitesse d'une equation issue de la mécanique des fluides a l'aide de Python(x,y) et de quiver.
J'ai les équations mais mon résultats n'est pas bon... je n'arrive pas a obtenir le bon champ...Voici le code :
Code:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
| import numpy as np
import matplotlib.pyplot as plt
a=1
N=50
u=np.linspace(-5,5,N)
Vx=np.zeros([N,N])
Vy=np.zeros([N,N])
for kx,x in enumerate(u):
for ky,y in enumerate(u):
Vx[kx,ky]=5*a*(5**(a-1))*(np.sin((1-a)*3*np.pi)/2)
Vy[kx,ky]=5*a*(5**(a-1))*(-np.cos((a-1)*3*np.pi)/2)
plt.quiver(Vx,Vy)
plt.show() |
j'ai pris 5 pour la constante Ko qui y a devant et 5 pour r apres changement d'origine (passage d'un r,teta a x,y)
Si vous pouvez m'aider merci beaucoup^!
Paul
