RAPPORT DEVOIR DE MATLAB : MARGE DE STABILITE D’UN SYSTEME
AVEC LE DIAGRAMME DE BODE
07/05/2024
Objectif :
Tracer et superposer les différentes réponses indicielles du système pour les 3 valeurs de Kp, voir et
commenter l’influence de Kp sur la marge de stabilité du système.
EN BOUCLE FERMEE
Schéma du système :
Fonction de transfert :
Kp
H ( p )= 3 2
avec Kp=0.15; 1 ; 3
0 , 5 p +2 p + p+ Kp
Code Matlab
Réponse indicielle
kp=[0.15 1 3]
num1=[kp(1)]
den1=[0.5 2 1 kp(1)]
H1=tf(num1,den1)
step(H1)
hold on
num2=[kp(2)]
den2=[0.5 2 1 kp(2)]
H2=tf(num2,den2)
step(H2)
num3=[kp(3)]
den3=[0.5 2 1 kp(3)]
H3=tf(num3,den3)
Diagrammes de Bode
kp=[0.15 1 3]
for j=1:3
subplot(2,2,j)
num=[kp(j)]
den=[0.5 2 1 kp(j)]
H=tf(num,den)
bode(H)
end
Résultats du tracé
Réponse indicielle
Diagrammes de Bode
Légende :
Courbe orange : réponse pour Kp = 0.15
Courbe violette : réponse pour Kp = 1
Courbe verte : réponse pour Kp = 3
Step reponse pour chaque système
kp = 3 1
kp = 0.15 1 Commentaire
Kp = 1 :
Courbe orange atteint la stabilité au bout de 30 secondes ;
Courbe orange atteint la stabilité au bout de 40 secondes ;
Courbe orange atteint la stabilité après plus de 100 secondes ;
On pourrait donc conclure que plus la valeur de Kp est petite, plus le système atteint rapidement la
stabilité. Cela se voit notamment dans les résultats des stepinfo
EN BOUCLE OUVERTE
Schéma du système
Fonction de transfert :
Kp
H ( p )= 3 2
0 , 5 p +2 p + p
Code Matlab :
kp=[0.15 1 3]
num1 = kp(1);
den1 = [0.5 2 1 0]
Ho1=tf(num1,den1)
step(Ho1)
hold on
num2 = kp(2);
den2 = [0.5 2 1 0]
Ho2=tf(num2,den2)
step(Ho2)
num3 = kp(3);
den3 = [0.5 2 1 0]
Ho3=tf(num3,den3)
step(Ho3)
Courbe obtenue : Réponse indicielle
Légende :
Courbe orange : réponse pour Kp = 1
Courbe violette : réponse pour Kp = 3
Courbe rouge : réponse pour Kp = 0.15
Diagramme de Bode
Commentaire :
Quelle qu’en soit la valeur de Kp, le système n’atteint jamais la stabilité