%-------------------------------Système linéaire---------------------------
clc
clf
close all
clear all
disp('						FONCTION DE TRANSFERT ')
disp(' ');
disp('	Maintenant nous allons analyser le comportement d''un système');
disp('  Dans les domaines fréquentiel et temporel en utilisant ');
disp('   les représentations graphiques suivantes :')
disp('  indicielle; BODE; NYQUIST;BLACK;etc...  ')
disp(' ')
disp(' pour cela la fonction de transfert du système est nécessaire pour l''étude ')
disp(' ');
disp(' ');
disp('	Taper une touche pour continuer');
pause
clc
disp('	EXEMPLE d''UN SYSTEME DU SECOND ORDRE ')
disp('  Edition et Affichage d''une fonction de Transfert:')
disp(' ')
disp('	une FT est définie comme étant le rapport de deux polynômes')
disp('	 édités suivant les puissances décroissantes')
disp(' ');
disp(' ');
disp(' soit un système du second ordre défini comme suite:')
disp('  ')
disp('  ')
disp('  de pulsation propre : 2rd/s')
disp(' de gain statique :1')
disp(' de coefficient d''amortissement 0,2')
disp(' ')
disp('                1')
disp('soit G(p)=---------------')
disp('           p^2+0,8.p+4')
disp('  ')
disp(' tapez une touche pour continuer')
pause
%keyboard
%clc
disp('  ')
disp('  ')
disp('	Entrer le numérateur et le dénominateur de la FT en entrant ')
disp('entre crochets les coefficients du polynôme')
disp('  soit :')
disp('  ')
disp('       num=[1]')
disp('       den=[1 0.8 4] ')
disp('  ')
disp('ensuite tapez return en toute lettre pour continuer')
keyboard
disp(' avec le numérateur (num) et le dénominateur (den) que vous venez de créer')
disp('  ')

disp(' Définir la Fonction de transfert nommée G avec la commande suivante:')
disp('  ')
disp(' G=tf(num,den) ')
disp('  ')
disp('A vous maintenant')
disp('ensuite return pour continuer')
keyboard
disp('  ')
disp(' pensez à vérifier la présence des variables que vous venez de créer')
disp('  avec la commande whos')
disp('   ')
disp('ensuite  return pour continuer')
keyboard
clc
disp('votre fonction de transfert est bien ?')
disp('  ')
disp('   ')
G
disp('sinon recommencer')
disp('  ')
disp('si oui tapez return')
keyboard
clc

 disp('                  ANALYSE TEMPORELLE ')
 disp(' ')
 disp(' Tracer la réponse indicielle du système avec la commande suivante:')
 disp('  ')
disp(' step(G) ')
disp('  ')
disp(' pour plus d''info taper help step ')
disp(' ')
disp('return pour continuer'	)
keyboard

disp('pour analyser les caractéristiques de cette réponse:')
disp(' ')
disp('faites un click droit sur la figure obtenue')
disp(' et sélectionner "grid" ')
disp('  ')
disp('ensuite  "characteristics"')
disp(' peak response : dépassement')
disp(' setting time: temps de réponse')
disp('rise time : temps de montée')
disp('steady state : valeur finale')
disp(' ')
disp(' ')
disp('essayer également la commande "properties"')
disp('qui vous permettra d''ajuster les paramètres de l''affichage')
disp(' ')
disp('  ')
disp(' ensuite return pour continuer')
keyboard
clc
disp('                       ANALYSE FREQUENTIELLE  ')
disp('')
disp('             Tracer du diagramme de BODE du système:')
disp(' ')
disp('si vous souhaitez garder la courbe précédente, affichée sur la figure1')
disp('activer tout d''abord la figure2 avec la commande suivante:')
disp(' figure(2)')
disp('  ')
disp(' ensuite tracer le diagramme de bode avec la commande:')
disp('  ')
disp('bode(G)')
%figure(2);
disp(' ')
disp('return pour continuer')
keyboard

disp('un click droit sur la figure vous permettra d''afficher les caractéristiques')
disp(' ')
disp('peak response: gain max')
disp(' minimum stability margins: marge de phase mini')
disp(' ')
disp(' ')
disp('return pour continuer')
keyboard
clc
disp('Tracer le diagramme de NICHOLS avec les commandes suivantes :')
disp('  ') 
disp('tout d''abord la commande :')
disp('  ')
disp(' figure(3) (si vous désirez changer de figure)')
disp('  ')
disp('ensuite la commande:')
disp('  ')
disp('nichols(G)')
disp('  ')
%figure(3)
disp('de même afficher les caractéristiques')
disp('  ')
disp('  ')
disp('ensuite return pour continuer')
keyboard

disp('   Tracer le diagramme de NYQUIST avec les commandes suivantes  :')
disp('  ')
disp(' tout d''abord la commande : ')
disp('  ')
disp('    figure(4) (si vous désirez changer de figure)')
disp('  ')
disp(' ensuite la commande :')
disp('  ')
disp(' nyquist(G)')
disp('  ')
%figure(4)
disp('afficher les caractéristiques')
disp('  ')
disp('  ')
disp(' return pour continuer ')
keyboard
clc



%-------------------------------------FIN----------------------------------
%----
clc
disp(' ')
disp(' ')
disp(' ')
disp('	ATTENTION  LE PROGRAMME INIMAT EST TERMINE: ')
disp(' ')
disp('	vous allez rentrer sous le fonctionnement autonome de MATLAB:   ')
disp(' ')
disp(' ')
disp('	Vous gardez acces à toutes les variables crées dans inimat, ')
disp('	avec les même commandes que vous venez d''expérimenter ')
disp(' ')
disp(' ')
disp(' ')
disp('         Maintenant nous allons procéder différemment ')
disp(' à l''analyse fréquentielle et temporelle d''un système linéaire')
disp(' en utilisant un outil graphique appelé LTI VIEWER GUI ')
disp('           voir fiche de guidance LTI VIEWER GUI')
disp('  ')
disp(' .......Frappez une touche quelquonque pour retourner à MATLAB....')
disp('  ')
disp('         ensuite taper la commande suivante:   "ltiview"')
pause
