Ejemplo en Matlab: modos transitorios

% DEMO_CISOIDES
% EJEMPLO CONCEPTUAL: 
% MODOS TRANSITORIOS EN FUNCIÓN DE SU POSICIÓN EN EL PLANO COMPLEJO

figure(1); 
clf; 

% Generamos vector de tiempos (de 0 a 20 segundos)
t = linspace(0,20,1000); 
color = {'r','g','b','y','c','m','k'}; 
for i=1:7, 
	figure(1); 
	
	% Elegimos un punto "s" del plano complejo
	h = subplot(131); 
	hold on; 
	axes(h); 
	axis([-3 1 -2 2]); 
	grid on; 
	plot([-3,1],[0,0],'k','linewidth',2); 
	plot([0,0],[-2,2],'k','linewidth',2);  
	[x,y]=ginput(1); 
	axis([-3 1 -2 2]); 
	grid on; 
    x = round(x*5)/5;
    y = round(y*5)/5;
	s = x + j*y; 
	plot(real(s),imag(s),'x','color',color{i},'markersize',30,'linewidth',5); 
	hold on; 
	
	% Dibujamos e^st
	subplot(132);
    c = exp(s*t);
	plot(real(c),imag(c),color{i},'linewidth',5); 
	axis equal; 
	axis([-1 1 -1 1]); 
	hold on; 
	
	% Dibujamos la respuesta típica real añadiendo su conjugada
	subplot(133); 
	plot(t,(c+conj(c))/2,color{i},'linewidth',3); 
    hold on;
	plot([0,t(end)],[0,0],'k','linewidth',2); 
	plot([0,0],[-5,5],'k','linewidth',2);  
    hold off;
	axis([-0.5,t(end),-2,2]); 
	grid on; 
	hold on; 
end;