Rampe

%------------------------------------
% réponse à une rampe 2ème ordre
% (c) JPaul Molina
%------------------------------------

clear all;                            % réinitialise
echo on
%---------------------------------
% réponse rampe 2éme ordre
%-----------------------------------
echo off
tmax=input('représentation : temps maxi en s ? ');    % 5 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

% calcule la réponse à une rampe y(T)   
[y,x]=lsim(K,denom,t,t); 
clf;                                 			% efface la figure
plot(t,y,t,K*t);                           			% tracé
	% mise en place des légendes
title('réponse à une rampe 2ème ordre');
text(0.1, 0.5*K,'Tracé de Kt'); 
xlabel('temps en s');    				% légende abscisse
ylabel('y(t)');          				% légende ordonnée  
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','F(p)= ---------------','          ap²+bp+1'} )