Influence d'un zéro

%------------------------------------
% ajout d'un zéro sur un 2ème ordre
% réponse indicielle 
% (c) JPaul Molina
%------------------------------------
clear all;                            	% réinitialise
echo on
%---------------------------------
% réponse indicielle 
% ajout d'un zéro sur 2éme ordre
%-----------------------------------
echo off
tmax=input('représentation : temps maxi en s ? ');% 20 par exemple 
t=(0:0.01:tmax)';
K=input('Entrer K : ');                   		% 1 par exemple  
a=input('Entrer le coefficient de p² : ');		% 0.1 par exemple  
b=input('Entrer le coefficient de p  : ');		% 0.2 par exemple 
denom=[a b 1];                         	% dénominateur ap²+ bp + 1
	% numérateurs
num1=[1/10 1]; % tau = 0.1
num2=[  4 1];  % tau = 4
num3=[10  1];  % tau = 10

% calcule la réponse indicielle y(T)   
[y1,X,T]=step(K*num1,denom,t); 
[y2,X,T]=step(K*num2,denom,t); 
[y3,X,T]=step(K*num3,denom,t); 
clf;                                 		% efface la figure
plot(t,y1,t,y2,t,y3);                		% tracé
	% calcul  des maxi
ymax1=max(y1);
ymax2=max(y2);
ymax3=max(y3);
	% mise en place des légendes
title('Influence de l''ajout d''un zéro sur un 2ème ordre(réponse indicielle) ');
xlabel('temps en s');    			% légende abscisse
ylabel('y(t)');          			% légende ordonnée  
	% écrit les textes au niveau du 1er maxi
text( t(ymax1), ymax1, 'tau = 0.1' );
text( t(ymax2), ymax2, 'tau = 4' );
text( t(ymax3), ymax3, 'tau = 10' );

grid;                    			% dessine les axes

echo on
% placer la légende en cliquant à l'endroit souhaité
echo off

%le texte qui suit sera affiché au point cliqué
gtext( {'            K(1+tau p)','F(p)= ----------------','            ap²+bp+1'} )