#Importation bibliothèques
import numpy as np
import matplotlib.pyplot as plt
import scipy as signal
import cmath as cm
import matplotlib.ticker as tck

####
Npt=1000
f=np.logspace(1,5,Npt)
H0=????    # valeur de l'amplification statique
f0=????    # valeur de la fréquence propre en Hz
Q=?????    # valeur du facteur de qualité du système

H=H0/(1-f**2/f0**2+1j*f/(Q*f0))  # transmittance isochrone complexe pour un passe-bas d'ordre 2
H6=??????                        # transmittance isochrone complexe pour un passe-bas d'ordre 6
GdB=????????????????????         # expression littérale du gain



######Graphe

plt.figure(0)

plt.plot(f,GdB,color='black')

plt.title('Diagramme de Bode en gain du passe-bas - ordre 6 ')
plt.yticks(np.arange(int(min(GdB))-2, int(max(GdB))+10, 10))
ax=plt.gca()
ax.yaxis.set_minor_locator(tck.MultipleLocator(base=5))
plt.xscale('log')

plt.xlabel('f en Hz')
plt.ylabel('G en dB')
plt.grid(b=True, which='major',axis="both", color='black', linestyle='-')
plt.grid(b=True, which='minor',axis='both', color='grey', linestyle='--')

plt.show()