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'} )
Dernière Modification : Lun 19 Mars 2007 13:28
Copyright © 1999-2010 Jean-Paul Molina Tous droits réservés.