%%%%%%%%%%%%%%%%%%%%%%%%
%%% Programme rebi.m %%%
%%%%%%%%%%%%%%%%%%%%%%%%

close all; clear all

% Définition d'une représentation d'état du système
A=input('matrice A du système?')     
B=input('matrice B du système?')
C=input('matrice C du système?')
[n,m]=size(B); [p,n]=size(C);

% Placement de pôles 
vitesse=input('partie réelle maximale des pôles en BF (entrer une valeur négative) ?')
dzeta=input('amortissement minimal des pôles en BF?')
Pbf=placement_poles_rebi(vitesse,dzeta,A,n);

% Représentation d'état du système augmenté
Aa=[A zeros(n,1);C 0];
Ba=[B;0];

% Calcul de Ka 
Ka=acker(Aa,Ba,Pbf);
Kp=Ka(1:n);
Ki=Ka(n+1);

% Tracé des lieux de Bode et calcul des marges de module
[Mmi,Mmo,Mmci,Mmco]=calcul_marges_1(A,B,C,Kp,Ki,n,m,p)