Examen 2021-2022 D'identification Et Modélisation
Examen 2021-2022 D'identification Et Modélisation
Examen 2021-2022 D'identification Et Modélisation
I MSEI / GEM
ENSAM Identification/ Modélisation des
Casablanca EXAMEN 2021/ 2022 Prof:A. TOUATI
Questions préliminaires:
1. Citer les types de boucles de régulation industrielles que vous connaissez?
2. Quelle est la différence entre un régulateur dit à sens d'action direct et un régulateur a
sensd'actioninverse ?
3. Le régulateur PID que vous avez 'habitude de dimensionner est-il à sens d'action direct
ou inverse ?
Problème N°1
Soit le système définit par le schéma synoptique suivant
1 yft
b g + +bo042 1+a0 q-2
u(t)
etT
1. De quel type de modèle de bruit s'agit-il ? Exprimer la relation récurrente qui relie y{k)
1) ?
A partir des données entrées/sorties du système, la méthode des moindres carrés MMC
sera appliquée pour identifier les paramètres du modèle. La réponse indicielle du système
a donné pour les tro is premiers points d'échantillonnage (T=0.2s) les valeurs suivantes
yl1)-1.12; y(2)= 0.95; y(3)=1. On prendra pour les valeurs initiales : u(0) = 1 et y{0) =
0.
Problème N°2:
On consideère un procédé d'entrée u et de sotie y qui sera modélisé soit par sa présentation
d'état
Les grandeurs xi, Xa et xa du vecteur d'état sont toutes mesurables, de même que la sortie y.
Calculer dans ces conditions les gains du retour d'état [Ka Ka Kal permettant de placer les
Problème N°3:
On considère le contrôle du systême illustrê ci-dessous. Le vecteur d'état est: = [X1 *21",
Le système en boucle ouverte est instable, par conséquent une commande Optimale par
plant
r:0
k
Schéma synoptique delacommande Optimale par retour d'état.
1. Apartir du schéma bloc du système a Plant », identifier les relations entre les variables
d'état x1, X2 avec l'entrée de commande u(t) et la sortie y(t) ? En déduire alors matrices
--asonooz-. F i n du sujet -
Annexe
|e(k)
**********************t**********
H()
Expression du prédicteur à un pas
Gq) (k)
Se-1)H)+(1 yo) u(k)
GC)
y(k)
* * * * * * - * * * * * - - - * -
b P ++b,P + bo () =|
Y(P) 0 0 0 1
x(0)+ |:|u()
U(p) p" + an-1Pn-1 + +a,P +ao -d0 -01 -n-1
yc)=[bo bm 0 0 x(0)
Matrice de commandabilité Q [B AB AB AT-1B]
Fonction de transfert en BF. Ggr(P) = C.(pl,-A + BK)-1 B
-
1- 4
+6
-2
com(B) -10 2 2
-2
Indice de performance
dt=-xPxl8°=x(0)Pz(0)|
Commande Optimale: Détermination de (HTP+PH) = -In
la matrice de pondération P. Avec H = A- BK