%%%%%%%%%%%%%%%%%%%%%%%
%%% Programme pii.m %%%
%%%%%%%%%%%%%%%%%%%%%%%

clear all; close all;

% Définition du système à corriger
s=tf('s');
P=input('Fonction de transfert du système à corriger?')

% Définition de la pulsation de coupure et de la marge de phase souhaitées
w0=input('Pulsation de coupure souhaitée pour le système corrigé?')
Mp=input('Marge de phase souhaitée (en radians)?')

% Caractéristiques du PI à appliquer
[gain_P_w0,phase_P_w0]=bode(P,w0); 
phase_P_w0=phase_P_w0/180*pi; 
gain_K_PI=1/gain_P_w0; phase_K_PI=-pi+Mp-phase_P_w0;

% Calcul du régulateur PI
K_PI=calcul_pi(gain_K_PI,phase_K_PI,w0);

% Diagrammes de Bode
L=P*K_PI;
figure; bode(L,K_PI,'-.',P,'--'); grid

% NOTA BENE: même remarque que pour le programme pdf.m