## Diagramme de Bode

import numpy as np
import matplotlib.pyplot as plt

plt.close('all')

Q = 1
f0 = 1500
w0 = 2*np.pi*f0

N=1000

f = np.linspace(f0/10,10*f0,N)
w = 2*np.pi*f
H = 1/ np.sqrt(1+ (Q*(w/w0-w0/w))**2 ) 

GdB = 20*np.log(H) / np.log(10)


plt.figure(figsize=(8,6))
plt.semilogx(f,GdB)
plt.xticks(fontsize=14)
plt.yticks(fontsize=14)
plt.xlabel(" $f$ (Hz)",fontsize = 16)
plt.ylabel(" $G_{dB}$",fontsize = 16)
plt.title("Diagramme de Bode",fontsize = 18)
plt.xscale('log')
plt.grid(True,which="both")
plt.show()