% studio evoluzione libera di un sistema 2D con autovalori complessi
% a parte reale negativa
% varie condizioni iniziali

clear all;close all;

A=[0 1; -2 -2]
B=[1;1];
C=[1 0; 0 1];
D=[0;0];

% autovettori (V) e autovalori (L)
[V,L]=eig(A)
pause

% condizioni iniziali
x0_1=[1;1]
x0_2=[-2;1]
pause

T=0:.01:100;
[X1,Y1]=initial(A,B,C,D,x0_1,T);
[X2,Y2]=initial(A,B,C,D,x0_2,T);
pause

plot(X1(1,1),X1(1,2),'*')
hold on
plot(X1(:,1),X1(:,2))
title('evoluzione libera: traiettoria nello spazio di stato')

pause
hold off

plot(X2(1,1),X2(1,2),'*')
hold on
plot(X2(:,1),X2(:,2))
title('evoluzione libera: traiettoria nello spazio di stato')
