#Importation bibliothèques
import numpy as np
import matplotlib.pyplot as plt
import scipy as signal
import cmath as cm
import pylab
import matplotlib.ticker as tck
pi=np.pi
####
Npt=5000 #nombre de points à tracer
f=np.linspace(0,5,Npt)
H0=1 #valeur théorique de l'amplification statique
f0=1.5 #valeur de la fréquence propre du système

#Premier système
Q= ???? #valeur du facteur de qualité
H=??? #formule littérale de la transmittance du système
modH= ???  #formule littérale du gain du système


#deuxième système
Q1=??? #valeur du facteur de qualité
H1=??? #formule littérale de la transmittance du système
modH1=???  #formule littérale du gain du système


#troisième système
Q2=??? #valeur du facteur de qualité
H2=??? #formule littérale de la transmittance du système
modH2=???  #formule littérale du gain du système


#quatrième système
Q3=??? #valeur du facteur de qualité
H3=??? #formule littérale de la transmittance du système
modH3=??  #formule littérale du gain du système


######Graphe

plt.figure(0)

plt.plot(f,modH,'k-',label='Q=2.50')
plt.plot(f,modH1,'r-', label='Q=1.5')
plt.plot(f,modH2,'g-', label='Q=0.707')
plt.plot(f,modH3,'b-', label='Q=0.500')
plt.legend()
plt.title('Diagramme théorique du module de H - système passe bas - ordre 2')
ax = plt.gca()

ax.xaxis.set_major_locator(tck.MultipleLocator(base=1))
ax.xaxis.set_minor_locator(tck.MultipleLocator(base=0.2))
ax.yaxis.set_major_locator(tck.MultipleLocator(base=0.5))
ax.yaxis.set_minor_locator(tck.MultipleLocator(base=0.1))


plt.xlabel('f en Hz')
plt.ylabel('(module de H)=Xm/Xa,  sans unité')
plt.grid(b=True, which='major',axis='both', color='grey', linestyle='-')
plt.grid(b=True, which='minor',axis='both', color='grey', linestyle='--')




plt.show()