Je veux réalisé un programme qui représente la trajectoire d'un satellite avec la Terre représentée par un cercle bleu avec solve_ivp :

import numpy as np
from math import sqrt
from pylab import plot,show,linspace,Circle,gcf,figure,gca
from scipy import integrate

G, M_T, R_T,Cx = 6.67*10**(-11), 5.972e24, 6.4e6,0.47

circle = Circle((0, 0), R_T, color='b')
figure(figsize=(7,7))
ax = gca()
ax.cla()
ax.add_artist(circle)
temps = np.linspace(0, 8500,10000)

def P(h):
    """
    argument : altitude en km : renvoie la masse volumique de l'air correspondant à l'altitude donnée en argument
    """
    val=np.array([[0,50,100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950,1000],[1.29e-3,8.261e-7,5.173e-10,1.561e-12,1.583e-13,2.644e-14,5.63e-15,1.382e-15,3.868e-16,1.308e-16,5.704e-17,3.163e-17,2.045e-17,1.428e-17,1.034e-17,7.625e-18,5.691e-18,4.289e-18,3.263e-18,2.507e-18,1.945e-18]])
    m,M=[],0
    for k in range(0,21):
        m.append(abs(h-val[0][k])) #compare l'altitude donnée en argument avec la valeur la plus proche
    M=np.min(m)
    for k in range(0,21):
        if M==m[k]:
            return val[1][k]

def f(temps,vecteur):
    """
    argument : une matrice colonne qui contient les CI de la position et de la vitesse en x et y et la fonction P : renvoie la dérivée et l'accélération selon x et y
    """
    Fx=-(G*M_T*vecteur[0])/((vecteur[0]**2+vecteur[1]**2)**(3/2)) #force de gravitation en x
    Fy=-(G*M_T*vecteur[1])/((vecteur[0]**2+vecteur[1]**2)**(3/2)) #force de gravitation en y
    Frx=((P(sqrt(vecteur[0]**2+vecteur[1]**2)))*0.58*Cx*sqrt(vecteur[2]**2+vecteur[3]**2)*vecteur[2])/2 #force de trainée en x
    Fry=((P(sqrt(vecteur[0]**2+vecteur[1]**2)))*0.58*Cx*sqrt(vecteur[2]**2+vecteur[3]**2)*vecteur[3])/2 #force de trainée en y
    return vecteur[2],vecteur[3],Fx+Frx/(43.6),Fy+Fry/(43.6)

x = integrate.solve_ivp(f,temps,(400e3+R_T, 0,0,7.7e3))

#plot(x,y)
#show()

Il y a un problème à la ligne 38 : too many values to unpack (expected 2).

Je voudrais récupérer la liste des résultats pour représenter la trajectoire mais je ne sais pas comment faire.

-
Edité par Runhelm 8 mars 2022 à 10:50:40

Tu peux mettre le message Traceback complet; mais je pense que le paramètre temps que tu passes à la fonction solve_ivp; la fonction attend 2 valeurs  (les bornes de l'intervalle d'intégration) et non pas les 10000 valeurs de temps (https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.solve_ivp.html )

x = integrate.solve_ivp(f,(0,8500),(400e3+R_T, 0,0,7.7e3))
ou pour être plus générique
x = integrate.solve_ivp(f,(temps[0],temps[-1]),(400e3+R_T, 0,0,7.7e3))

Edit: je crois que j'ai laissé involontairement un t devant (0,8500) sur ma 1ère ligne, corrigé

-
Edité par umfred 21 mars 2022 à 11:58:41

Bonjour,

x,y = sol.y

tu es censé récupérer 2 valeurs : x et y

mais sol.y c'est une seule valeur, non ?

sol=integrate.solve_ivp(.....)

donc le sol de ta vidéo correspond au x de ton code initial

Sinon, tu fournis 4 valeurs initiales (et f renvoie 4 valeurs) donc sol.y contient 4 tableaux (un pour chaque variable) donc tu dois recevoir 4 variables qui doivent correspondre à x,y (position) et Vx,Vy (vitesse)

Avec 

import numpy as np
from math import sqrt
from pylab import plot,show,linspace,Circle,gcf,figure,gca
from scipy import integrate

G, M_T, R_T,Cx = 6.67*10**(-11), 5.972e24, 6.4e6,0.47

circle = Circle((0, 0), R_T, color='b')
figure(figsize=(7,7))
ax = gca()
ax.cla()
ax.add_artist(circle)
temps = np.linspace(0, 8500,10000)

def P(h):
    """
    argument : altitude en km : renvoie la masse volumique de l'air correspondant à l'altitude donnée en argument
    """
    val=np.array([[0,50,100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950,1000],[1.29e-3,8.261e-7,5.173e-10,1.561e-12,1.583e-13,2.644e-14,5.63e-15,1.382e-15,3.868e-16,1.308e-16,5.704e-17,3.163e-17,2.045e-17,1.428e-17,1.034e-17,7.625e-18,5.691e-18,4.289e-18,3.263e-18,2.507e-18,1.945e-18]])
    m,M=[],0
    for k in range(0,21):
        m.append(abs(h-val[0][k])) #compare l'altitude donnée en argument avec la valeur la plus proche
    M=np.min(m)
    for k in range(0,21):
        if M==m[k]:
            return val[1][k]

def f(temps,vecteur):
    """
    argument : une matrice colonne qui contient les CI de la position et de la vitesse en x et y et la fonction P : renvoie la dérivée et l'accélération selon x et y
    """
    Fx=-(G*M_T*vecteur[0])/((vecteur[0]**2+vecteur[1]**2)**(3/2)) #force de gravitation en x
    Fy=-(G*M_T*vecteur[1])/((vecteur[0]**2+vecteur[1]**2)**(3/2)) #force de gravitation en y
    Frx=((P(sqrt(vecteur[0]**2+vecteur[1]**2)))*0.58*Cx*sqrt(vecteur[2]**2+vecteur[3]**2)*vecteur[2])/2 #force de trainée en x
    Fry=((P(sqrt(vecteur[0]**2+vecteur[1]**2)))*0.58*Cx*sqrt(vecteur[2]**2+vecteur[3]**2)*vecteur[3])/2 #force de trainée en y
    return vecteur[2],vecteur[3],Fx+Frx/(43.6),Fy+Fry/(43.6)

x,y,dx,dy = integrate.solve_ivp(f,(temps[0],temps[-1]),(400e3+R_T, 0,0,7.7e3))

plot(x,y)
show()

J'ai ValueError: too many values to unpack (expected 4).

Mais donc si j'écris

x = integrate.solve_ivp(f,(temps[0],temps[-1]),(400e3+R_T, 0,0,7.7e3))

x prend x,y,vx,vy comme valeur mais comment récupérer y pour l'afficher ?

Si tu utilises une variable nommé x alors que je l'ai nommé sol et que j'accède aux y de la solution avec sol.y, alors tu devrais pouvoir y accéder de la même façon à savoir x.y

J'ai fait 

import numpy as np
from math import sqrt
from pylab import plot,show,linspace,Circle,gcf,figure,gca
from scipy import integrate
 
G, M_T, R_T,Cx = 6.67*10**(-11), 5.972e24, 6.4e6,0.47
 
circle = Circle((0, 0), R_T, color='b')
figure(figsize=(7,7))
ax = gca()
ax.cla()
ax.add_artist(circle)
temps = np.linspace(0, 8500,10000)
 
def P(h):
    """
    argument : altitude en km : renvoie la masse volumique de l'air correspondant à l'altitude donnée en argument
    """
    val=np.array([[0,50,100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950,1000],[1.29e-3,8.261e-7,5.173e-10,1.561e-12,1.583e-13,2.644e-14,5.63e-15,1.382e-15,3.868e-16,1.308e-16,5.704e-17,3.163e-17,2.045e-17,1.428e-17,1.034e-17,7.625e-18,5.691e-18,4.289e-18,3.263e-18,2.507e-18,1.945e-18]])
    m,M=[],0
    for k in range(0,21):
        m.append(abs(h-val[0][k])) #compare l'altitude donnée en argument avec la valeur la plus proche
    M=np.min(m)
    for k in range(0,21):
        if M==m[k]:
            return val[1][k]
 
def f(temps,vecteur):
    """
    argument : une matrice colonne qui contient les CI de la position et de la vitesse en x et y et la fonction P : renvoie la dérivée et l'accélération selon x et y
    """
    Fx=-(G*M_T*vecteur[0])/((vecteur[0]**2+vecteur[1]**2)**(3/2)) #force de gravitation en x
    Fy=-(G*M_T*vecteur[1])/((vecteur[0]**2+vecteur[1]**2)**(3/2)) #force de gravitation en y
    Frx=((P(sqrt(vecteur[0]**2+vecteur[1]**2)))*0.58*Cx*sqrt(vecteur[2]**2+vecteur[3]**2)*vecteur[2])/2 #force de trainée en x
    Fry=((P(sqrt(vecteur[0]**2+vecteur[1]**2)))*0.58*Cx*sqrt(vecteur[2]**2+vecteur[3]**2)*vecteur[3])/2 #force de trainée en y
    return vecteur[2],vecteur[3],Fx+Frx/(43.6),Fy+Fry/(43.6)
 
sol = integrate.solve_ivp(f,(temps[0],temps[-1]),(400e3+R_T, 0,0,7.7e3))
 
plot(sol.x,sol.y)
show()

Mais ça me dit AttributeError: x

-
Edité par Runhelm 22 mai 2022 à 11:57:37