Influence d'un pôle et d'un zéro

%------------------------------------
% influence de pôle et de zéro sur un 2ème ordre
% réponse indicielle 
% (c) JPaul Molina
%------------------------------------
clear all;                            	% réinitialise
echo on
%-------------------------------------------------
% réponse indicielle 
% influence d'ajout de pôle et zéro sur 2éme ordre
%-----------------------------------------------------
echo off
tmax=input('représentation : temps maxi en s ? ');% 10 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
	% on définit la fonction à étudier
num1=K*[1  10];  	% K(p+10)  par exemple au numérateur de la fonction
denom1=conv( denom , [1  5] );  % effectue (ap²+ bp + 1)(p+5)

% calcule la réponse indicielle y(T)   
[y1,X,T]=step(K,denom,t); 	% second ordre classique
[y2,X,T]=step(num1,denom1,t); 
clf;                            		% efface la figure
plot(t,y1,t,y2);                	% tracé des 2 courbes
posX = round(length(t)/6);
	% mise en place des légendes
title('Influence de pôles 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 relatifs à chaque courbe
text( t(posX,:), y1(posX,:), 'classique' );
text( t(posX,:), y2(posX,:), 'modifié' );

grid;                    		% dessine les axes
echo on
% placer le texte du 2ème ordre modifié en cliquant à l'endroit souhaité
echo off
	%le texte qui suit sera affiché au point cliqué
gtext( {'                  K(p+10)','F1(p)= ------------------------','            (ap²+bp+1)(p+5)'} )

echo on
% placer le texte du 2ème ordre CLASSIQUE en cliquant à l'endroit souhaité
echo off
gtext( {'                 K','F(p)= ------------------','            ap²+bp+1'} )