close all
clear all
clc

Time=0:0.2:100;  % Asse di tempi

u=0; % Per la risposta libera l'ingresso è nullo

x0=[1 -1]'; % Vettore delle condizioni iniziali

OPTIONS=odeset('RelTol',1e-3,'AbsTol',1e-6);
[t,x_lin]=ode23('lin',Time,x0,OPTIONS,u); %Integrazione numerica

% L'uscita, nei sistemi lineari è una cominazione lineare degli stati

y=2*x_lin(:,1)+x_lin(:,2);

%Grafici delle traiettorie degli stati e delle uscite

plot(t,y);
grid
xlabel('t')
ylabel('y')

figure

plot(t,x_lin(:,1))
grid
xlabel('t')
ylabel('x1')

figure

plot(t,x_lin(:,2))
grid
xlabel('t')
ylabel('x2')


%Grafici dei movimenti degli stati

figure

plot3(t,x_lin(:,1),x_lin(:,2),'r')

grid
xlabel('t')
ylabel('x1')
zlabel('x2')