Aide exercice matlab

[invite80a3841b]

bnjr tt le monde,je serai trés reconnaissance si vous m'aidez à résoudre ce probleme:je ne sais pas comment tracer la commande par retour d'etat d'un systeme (dx/dt=Ax+Bu) sachant que u=-kx+kr r.g entamé le programmemais j'ai pas pu le terminer.merci d'avance .
A=[-0.006868 0.01395 0 -9.81;-0.09055 -0.3151 235.91 0;0.0003894 -0.003366 -0.4285 0;0 0 1 0]
B=[0.0000573,2.94;-5.44,0;-1.158,0;0,0];;
C=[1 0 0 0;0 0 0 1];
d=zeros(4,2);
I=[1 0 0 0;0 0 0 1];
[V,D]=eig(A)
r=0;
j=[1 0;0 1]
p=[-0.5+ 0.5i;-0.5-0.5i;-2.5+2.5i;-2.5-2.5i];
k=place(A,B,p)
k1=real(k)
kr=[-c*[[A-B*k1]^-1]*B]^-1*j
h=B*[-k1]
l=B*kr*r;
m=[A+h]
-----

[invited50d225a]

bsr,je te donne un exemple et tu applique la meme choz à tn exercice
tu ouvre une page script tu ecri
le programme qui defini la fonction qui contient le systeme diff et le correcteur K
:
function Xd=essai2(t,X)
uc=10;
p=[-1;-2;-3];
A=[0 1 0;0 0 1;0 0 0];
B=[0;0;1];
C=[1 0 0];
K=place(A,B,p);
u=-K*X+uc;
Xd=A*X+B*u;

et tu ouvre une nouvelle page et tu ecri le prog qui permet de visualiser l allure des variables
t0=0;
tf=20;
X0=[0;0;0]; %valeurs initiales des variables
[t,X]=ode45('essai2',[t0 tf],X0); % fonction prédéfini en matlab qui permet la resolution num des equa diff
figure(1)
plot(t,X(:,1))
figure(2)
plot(t,X(:,2))
figure(3)
plot(t,X(:,3))

bonne chance

[invite80a3841b]

merci bcp gse