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;