import numpy as np
import matplotlib.pyplot as plt


# Saisie des valeurs expérimentales dans 3 tableaux
F=np.array([])
Us=np.array([])
Ue=np.array([])

# Définition du Gain et linéarisation de l'expression
G = Us/Ue
X= A COMPLETER
Y= A COMPLETER

# Ajustement linéaire des points expérimentaux
a1, b1 = np.polyfit(X,Y,1)

# Tracé du nuage de points expérimentaux linéarisé et de la droite modèle.

plt.title("Tracé de 1/G^2 en fonction de f^2")
# Tracé des points expérimentaux
plt.plot(X,Y,'o',label="points expérimentaux")
plt.xlabel('f^2 en Hz^2')
plt.ylabel('1/G^2')
# Tracé de la droite modèle
plt.plot(X,(a1*X+b1),'r',label="droite modèle")
plt.legend()
plt.grid()
plt.show()

# Estimation de la fréquence de coupure à partir du coéfficient directeur et affichage de la valeur.
FC= A COMPLETER
print('Fréquence de coupure fc =',FC,'Hz')