%% Script per il confronto tra risposta forzata (calcolata con lsim)
%% risposta a regime permanente (calcolata dalla Teoria usando la risposta in frequenza)
%% calcolo della componente di transitorio della risposta
%% trovata per differenza


clear all
close all

dT = 0.01;
H = tf(1,[1 2 2]);

t=[0:dT:100];
omega=1;
u = sin(omega*t)-cos(omega*t);

Homega = 1/((j*omega)^2 + 2*j*omega + 2);

yrp = abs(Homega)*sin(omega*t+phase(Homega))-abs(Homega)*cos(omega*t+phase(Homega));

yf = lsim(H,u,t)';


ytr = yf-yrp;

figure
%plot(t,yf)
%hold on
%plot(t,yf2,':r')
plot(t,yf,'b')
hold on
plot(t,yrp,'r:')
plot(t,ytr,'g')

title('Risposta Forzata')
