Influence de pôles complexes

%------------------------------------
% comparaison : influence de pôles complexes
% réponse indicielle 
% (c) JPaul Molina
%------------------------------------
clear all;                            % réinitialise
echo on
%---------------------------------
% réponse indicielle 2éme ordre
% influence de pôles complexes
%-----------------------------------
echo off

tmax=15; % 15s 
t=(0:0.01:tmax)';
K=3.06;                    
a=1;  % coeff de p²
b=1;  % coeff de p
denom=[a b 1];                         % dénominateur ap²+ bp + 1
	% on définit la fonction à étudier
denom1=conv( denom , [1  0.55 K] );  % effectue (ap²+ bp + 1)(p²+0.5p+K)
denom2=[1  0.55 K]; % p²+0.5p+K
	% calcule la réponse indicielle y(T)   
[y1,X,T]=step(K,denom1,t); 
[y2,X,T]=step(K,denom2,t); 
clf;                            		% efface la figure
plot(t,y1,t,y2);                	% tracé des 2 courbes

ymax1=max(y1);
ymax2=max(y2);

% mise en place des légendes
title('comparaison (réponse indicielle) ');
xlabel('temps en s');    		% légende abscisse
ylabel('y(t)');          		% légende ordonnée  
	%prépare l'affichage
str1=sprintf('                  K\nF1(p)= ------------------------\n            (ap²+bp+1)(p²+0.5p+K)');
str2=sprintf('                  K\nF2(p)= ----------------\n            p²+0.5p+K');

% écrit les textes relatifs à chaque courbe
text( 3.5, ymax1, str1 );
text( 1, ymax2, str2 );

grid;                    % dessine les axes