R Gulation 1
R Gulation 1
R Gulation 1
système automatise
2.
a) Les paramètres du système :
On peut déterminer graphiquement :
Le gain statique : K = 3
La constante de temps :
Elle correspond a y(t) = 0,63.k=1,9 =2
Le temps de reponse
2) Réponse fréquentielle :
a) Le Diagramme de Bode :
Pour tracer le diagramme de Bode de la fonction G(p), on utilise le script suivant sous
MATLAB :
N=[3];
D=[21];
w= logspace(-1, 1, 100);
[gain,phase]=bode(N,D,w) ;
subplot(2, 1, 1);
semilogx(w, 20.*log10(gain),'r');
title('Diagramme de Bode en GAIN');grid;
subplot(2, 1, 2);
semilogx(w,phase);
title('Diagramme de Bode en PHASE');grid ;
Pulsation de coupure 0 :
b) Le diagramme de Nyquist :
Pour tracer le diagramme de Nyquist de la fonction G(p), on utilise le script
suivant sous MATLAB :
N = [3];
D = [2 1];
w = logspace(-3,3,100);
[reel, imaginaire] = nyquist(N,D,w);
plot(reel,imaginaire,'r');
title('Diagramme de Nyquist');grid;
Avec b0 = a0 =1 , a1 = 7 , a2 =10
1) Etude préliminaire :
Y ( p) 1 1
G1 ( p)
U ( p) 10 p 7 p a1 (5 p 1)(2 p 1)
2
Constantes de temps : 1 =5 et 2 =2
Gain statique : K=1.
Pulsation propre :
0 = a0 / a2 ===> 0 = 0,31 rad/s
Facteur d’amortissement :
a1
2 a0 a2 ===>
= 1,1
Le régime du système :
2) Réponse indicielle :
a) Simulation de la réponse indicielle :
N=[1];
D=[10 7 1];
w = logspace(-2,1,100);
[gain,phase] = bode(N,D,w);
subplot(2,1,1);
semilogx(w,20.*log10(gain),'r');
title('Diagramme de Bode en
GAIN');grid;
subplot(2,1,2);
semilogx(w,phase);
title('Diagramme de Bode en
PHASE');grid;
Commentaares :
Le diagramme de Bode en gain a deux asymptotes, la premiere en basses
fréquences ( << 0 )d'équation Gdb = 20log(k) = 0dB, et une autre en hautes
fréquences ( >> 0 )d'une équation de pente -40dB/décade. La phase de notre
système varie de 0 au -180 et passe par les environnements de -90 la
pulsation non amortie 0 .
Pulsation de coupure 0 :
Daprés le diagramme de Bode en gain, on peut déterminer la pulsation de
coupure 0 , pour Iaquelle |G|dB= Gmax -3dB= -3 dB
Graphiquement 0 = 2. 101 ===> 0 = 0.2 rad/s.
b) Le diagramme de Nyquist :
Pour tracer le diagramme de Nyquist de la fonction G1(p) , on utilise le script
suivant sous MATLAB
N = [1];
D =[10 7 1];
w= logspace(-3, 3, 100) ;
[reel,imaginaire]=nyquist(N,D,w) ;
plot(reel, imaginaire, 'r');
title('Diagramme de Nyquist'); grid ;
xlabel('Partie Réelle') ;
ylabel('Partie Imaginaire') ;
1) Etude preliminaire :
1 1
Y ( p) G1 ( p) Y ( p) 10 p 7 p 1
2
Y ( p) 1 Y ( p) 2
X ( p ) 1 G1 ( p) X ( p) 1 1 X ( p) 10 p 2 7 p 1 X ( p) 5 p 2 7 p 1
10 p 2 7 p 1 2
Gain statique : K = 0,5
Pulsation propre :
1
0 = ==> 0 = 0,44 rad/s
5
Facteur d’amortissement :
= 7 0 /4 ==> = 0,78
t D1 = ==> t D1 11,41s
1 2
2) Reponse indicielle :
a) Simulation de la reponse indicielle :
Commentawes :
Nous remarquons que le systeme en boucle fermée sans correcteur est
devenu plus rapide, mais il présente une erreur importante.
b) Les parametres du systeme .
Conclusion partielle .
En faisant une comparaison les résultats théoriques de K, et t D% et t D1 et
celles trouvés graphiquement, on constate qu'ils sont presque les memes.
Et on petit conclure que le bouclage d'un système sans un correcteur
n'est suffisant d'où la nécessité d'introduire un correcteur Pl pour arriver a
réduire le temps de réponse et l'erreur statique.
V. Simulation de boucle avec correction :
Soit le système du second ordre suivant en boucle fermée avec correcteur
proportionnelle
1) Etude preliminaire :
Kc Kc
Y ( p) Kc.G1 ( p) 10 p 7 p 1
2
Kc 1 Kc
X ( p ) 1 Kc.G1 ( p) 1 Kc 10 p 7 p 1 Kc
2
10
p2
7
p 1
10 p 7 p 1
2
1 Kc 1 Kc
Kc
Gain statique : K=
1 Kc
Pulsation propre :
1 10 1 Kc
0
0 1 Kc
2
10
Facteur d’amortissement :
2. 7 49
0 1 Kc 40.(1 Kc)
Le traçage de K, 0 et en fonction de Kc ;
Kc= 0:0.1:50;
K=(Kc./(1+Kc));
w=sqrt((1+Kc)/10);
z=sqrt(49./(40*(1+Kc)));
subplot(2, 3, 1);
plot(Kc,K,'r');
title ('Variations de K=f(Kc)');
subplot(2, 3, 2);
plot(Kc,w,'g');
title (' Variations de w=flKc)');
subplot(2, 3, 3);
plot(Kc,z,'b');
title(' Variations de z=f(Kc)de ')
Conclusion :
D'après les courbes ci-dessus, on
remarque qu'une petite variation
pour les
petites valeurs de Kc influence
sur le gain statique K beaucoup
plus que pour
des valeurs de Kc proche de 50.
Ainsi que la pulsation propre
croit avec
l'augmentation de Kc. Par contre,
le coefficient d'amortissement diminue en
donnant lieu des oscillations au niveau de la sortie du système.
2) Réponse indicielle :
a) Simulation de la réponse indicielle :
Commentaires :
Nous remarquons que l'augmentation de Kc entraine
La rapidité du systeme.
L'augmentation du gain statique K
L'augmentation de la pulsation propre 0
La diminution du coefficient d'amortissement , car la réponse
devient de plus en plus oscillante.
p Kc
Et : X( )= lim p 0 pX(p) = lim p 0 =1 Donc = 0,99 ==> Kc=99
p 1 Kc
49
Facteur d’amortissement = ==> = 0,11
40.(1 Kc)
1 Kc
Pulsation propre : 0 = ==> 0 = Kc .3,16 ras/s
10