calcul function de transfert

Bonjour a tous

Je souhaite calculer la fonction de transfert de la formule suivante:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
 %s= opérateur de Laplace
%j’espère avoir vecteur K de 1x4 fonction de transfert
 
K=Kc*inv(s*eye(4,4)-Ar+Br*Kc+K0*Cr+K0*Dr*Kc)*K0;
 
            Ar=[-89.6431   -0.4713   -0.1287         0
   -0.6725   -0.0119   -0.0776         0
  309.8703    0.6060    0.4269         0
    1.0000         0         0         0];
 
 
   K0 =[ 0.0788   -0.0020   -0.2616   -0.0008
   -0.0020    0.9520   -0.2463   -0.0010
   -0.2616   -0.2463    2.2022    0.0039
   -0.0008   -0.0010    0.0039    1.0000];
 
 Kc =[  421.8509 -134.5831  128.4216  127.1500];
 
Cr =[ 1     0     0     0
     0     1     0     0
     0     0     1     0
     0     0     0     1];
 
    Dr =[  0
           0
           0
           0]
 
Br = [0.0225
    0.0003
    0.0073
         0];