import matplotlib.pyplot as plt
import numpy as np
from math import *

# Importation des données
x=[0,0.75,1.65,2.6,3.75,5.05,6.4,7.8,9.3,10.9,12.55,14.2,15.8,17.5,19.2,20.9,22.3,23.7,25.1]
y=[4.15,3.7,3.3,2.8,2.4,1.9,1.5,1.1,0.7,0.5,0.3,0.2,0.25,0.35,0.45,0.7,1.1,1.5,1.85]
Vx=[]
Vy=[]

min_x = min(x)
max_x = max(x)
min_y = min(y)
max_y = max(y)

# Paramètres de la chronophotographie
N = len(x)
dt =   # Intervalle de temps entre deux positions

for i in range(1, N-1):
    vy=
    Vy.append(vy)          # Le calcul est effectué ici entre deux points entourant la position considérée
    vx=
    Vx.append(vx)
    v= (sqrt(vy**2+vx**2))
    print("à t=",round(dt*i,2),"s, la vitesse est de",round(v,2),'m/s')
    echelle=0.02
    vy=vy*echelle
    vx=vx*echelle

# Tracé des vecteurs vitesse
    plt.quiver(x[i],y[i],vx,vy,color="red",scale=0.4,scale_units='xy')

N2 = len(Vx)
for i in range(1, N2-1):
    dVy=(Vy[i+1]-Vy[i-1])          # Le calcul de la variation de vitesse
    dVx=(Vx[i+1]-Vx[i-1])
    dV=(sqrt(dVy**2+dVx**2))

    echelle=0.4
    dVy=dVy*echelle
    dVx=dVx*echelle
# Tracé des vecteurs variation de vitesse
    plt.quiver(x[i+1],y[i+1],dVx,dVy,color="blue",scale=0.5,scale_units='xy')

#Configuration du graphique
axes=plt.gca()
plt.axis('equal')
axes.set_xlim(0.8*min_x, 1.2*max_x)
axes.set_ylim(0.8*min_y, 1.2*max_y)
plt.xlabel("x (en m)")
plt.ylabel("y (en m)")
plt.title('Trajectoire d''un mobile autoporté accroché à un fil ',fontsize=8)
plt.plot(x, y, "go", markersize=8)
plt.grid()
plt.show()