Discussion :

[jeaneudes75]

Bonjour je suis entrain de faire un programme pour tracer la phase d'un diagramme de bode avec ce programme:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
def phase(FT,omega):
    argument=angle(FT(omega))*180/pi
    return argument

Je calcul la phase pour plusieurs omega, cependant lorsque je multiplie un 1er ordre avec un deuxième ordre et que j'affiche la phase, à partir de -180° la courbe de phase augment au lieu de descendre jusqu'à -270°.( Pour un 1er et 2eme ordre ça marche bien).
Je pense qu'il faut mettre une condition dans la fonction phase pour corriger celà mais je ne vois pas comment faire, si quelqu'un pourrait m'aider.
Cordialement.

[marco056]

Il faudrait peut-être le programme en entier pour comprendre le pb.

[jeaneudes75]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
 
def H3(omega):
    return H1(omega)*H2(omega) ## H1 est une foncion de transfert du 1er ordre et H2 du 2eme ordre.
 
def Graphe_Bode2(FT,omega_min,omega_max,n):
    fig,ax=plt.subplots(nrows=2,ncols=1,sharex=True)
    ax[0].semilogx(Liste_Bode(FT,omega_min,omega_max,n)[2],Liste_Bode(FT,omega_min,omega_max,n)[0]) ##Liste_Bode[2] correspond à la liste des valeurs de omega et Liste_Bode[0] la liste des valeurs du gain.
    ax[1].semilogx(Liste_Bode(FT,omega_min,omega_max,n)[2],Liste_Bode(FT,omega_min,omega_max,n)[1]) ##Liste_Bode[1] correspond à la liste des valeurs de la phase.
    plt.show()
 
Graphe_Bode2(H3,0,200,2000)

Voila ce que ça me fait sur le graphe pour la phase de la fonction de transfert du 3eme ordre:



Il faudrait modifier la fonction phase pour attribuer la bonne valeur de la phase après -180° mais je ne sais pas comment faire.

[marco056]

Ce ne serait pas un pb de physique ?
C'est quoi H1(w) et H2(w) ?
Cela devient lourd de devoir vous extraire des informations ?
Votre programme est confidentiel ou pas ?

[Chelmy88]

Et paf, ça fait des chocapics !