# -*- coding: utf-8 -*-
# Created on Sat Sep 31 08:15:00 2024
# @author: mhebding
# modelX le script pour confronter modele et experience

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

# Vos mesures experimentales ici
mesures_x = [100, 190, 230, 250, 260, 300, 310, 380, 500, 1000, 2000]
mesures_y = [0.002, 0.01, 0.015, 0.018, 0.019, 0.02, 0.017, 0.011, 0.007, 0.003, 0.002]

# Votre modele physique ici
def modele(x):
    " renvoie I(f) "
    E = 4
    R = 200
    L = 0.3
    C = 10**-6
    I = E/(R**2+(L*2*pi*x-1/(C*2*pi*x))**2)**0.5
    return I

# Trace et sauvegarde du graphe. Ne pas modifier
xmin, xmax = min(mesures_x), max(mesures_x)
modele_x = list(np.linspace(xmin, xmax, 500)) # modele avec 500 points
modele_y = [modele(x) for x in modele_x]
plt.clf ()
plt.plot(mesures_x, mesures_y, "+", label='points experimentaux')
plt.plot(modele_x, modele_y, label='modele physique')
plt.title ("Experience et modele")
plt.legend()
plt.savefig('modelX')
plt.show()