import numpy as np
import matplotlib.pyplot as plt

f = np.array([1e2 ,1e3, 1.6e3, 2e3, 3e3, 4e3, 5e3, 1e4, 1e5])
Us = np.array([])
phi = np.array([])

Ue =   # Amplitude du signal d'entrée
G= Us/Ue
GdB = 20*np.log(G)

plt.figure(figsize=(18, 5))

# ---- Gain G(w) ----
plt.subplot(1, 3, 1)
plt.plot(f, G, marker='o', label="G(w)", linewidth=2)
plt.xscale("log")
plt.xlabel("Fréquence (Hz)")
plt.ylabel("Gain")
plt.title("Diagramme de Bode en gain")
plt.grid(which='both')

# ---- Gain dB GdB(w) ----
plt.subplot(1, 3, 2)
plt.plot(f, GdB, marker='x', label="GdB", linewidth=2)
plt.xscale("log")
plt.xlabel("Fréquence (Hz)")
plt.ylabel("GdB")
plt.title("Diagramme de Bode en Gain (dB)")
plt.grid(which='both')
plt.show()

# ---- Phase φ(w) ----
plt.subplot(1, 3, 3)
plt.plot(f, -phi, marker='o', label="φ(w)", linewidth=2)
plt.xscale("log")
plt.xlabel("Fréquence (Hz)")
plt.ylabel("φ (°)")
plt.title("Diagramme de Bode en phase")
plt.grid(which='both')