Perturbation

%----------------------------------------------------------
% influence perturbation (réponse indicielle 2ème ordre)
%----------------------------------------------------------
clear all;                            % réinitialise
echo on
%------------------------------------------------------
% influence perturbation (réponse indicielle 2ème ordre)
%--------------------------------------------------------
echo off
tmax=5;      		% 5s 
t=(0:0.01:tmax)';
num1=[16 0];               % 16p
num2=[2 16];               % 2p+16
denom=[0.2 3 16];        % dénominateur 0.2p²+ 3p + 16

% calcule la réponse indicielle y(T)   
[y1,X,T]=step(num1,denom,t); 
[y2,X,T]=step(num2,denom,t);

clf;                                 	% efface la figure
plot(t,y1,t,y2);                           % tracé
% mise en place des légendes
title('influence d''une perturbation (réponse indicielle 2ème ordre)');
xlabel('temps en s');    		% légende abscisse
ylabel('y(t)');          		% légende ordonnée

posX=round( length(t)/5 );
	%prépare l'affichage
str1=sprintf('                  16p\nF1(p)= ------------------\n            0.2p²+3p+16');
str2=sprintf('                  2p+16\nF2(p)= ----------------\n            0.2p²+3p+16');

% écrit les textes relatifs à chaque courbe
text( t(posX,:), y1(posX,:), str1 );
text( t(posX,:), y2(posX,:), str2 );

grid;                    % dessine les axes