clear all
close all
clc

% parametri del sistema
y0 = 0; % deve essere a 0 per avere la linearità
b = 0.5;
m = 0.6;
load('dati_lab3_es3') % vettori x e x2
n = 0:50;

% faccio girare il sistema
y = [y0, popolazione(y0, m, b, x)];
% e lo rappresento
figure
stem(n,y,'b-')
grid on
xlabel('tempo n')
ylabel('numerosità della popolazione y(n)')
title('y(0)=0, b=0.5, m=0.6, ingresso x')


% verifico la linearità
a1 = 2;
a2 = 0.5;
y = [y0, popolazione(y0, m, b, a1*x+a2*x2)];
y1 = [y0, popolazione(y0, m, b, x)];
y2 = [y0, popolazione(y0, m, b, x2)];
% e la rappresento
figure
stem(n,y,'b-x')
hold on;
stem(n,a1*y1+a2*y2,'r-')
hold off
grid on
xlabel('tempo n')
ylabel('numerosità della popolazione y(n)')
title('verifica della linearità del sistema')