#Importation des bibliothèques
from scipy import signal
import numpy as np     
import matplotlib.pyplot as plt
import matplotlib.ticker as tck
#Fonction permettant de récupérer les coefficients du numérateur b et denominateur a

b, a = signal.butter(??, ???, '???', analog=True, output='ba')

#Fonction calculant pour chaque fréquence , la transmittance H

w, h = signal.freqs(b, a, worN=np.logspace(1,5,5000))
GdB=20*np.log10(abs(h))

#Tracé du diagramme de Bode en gain
plt.figure(0)
plt.plot(w, GdB,color='black')
plt.xscale('log')
plt.yticks(np.arange(-650, 20, 50))
ax=plt.gca()
ax.yaxis.set_minor_locator(tck.MultipleLocator(base=10))
plt.title('Système passe-bas de Butterworth - ordre   ')

plt.xlabel('Fréquence [Hz]')

plt.ylabel('Gain [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()