clear all
close all
clc

% parametri del sistema
y0 = 0; % deve essere a 0 per avere la linearità
n = 1:50;
b = 0.5+0.5*sin(n*pi/10); % provate a vedere cosa cambia se il sin viene 
                          % moltiplicato per 0 invece che per 0.5 !!!
m = 0.6;
x = zeros(size(n));
x(1) = 50;

% evoluzione #1
y = popolazione_tv(y0, m, b, x);

% traslo l'ingresso
x2 = zeros(size(n));
x2(5) = 50;

% evoluzione #2
y2 = popolazione_tv(y0, m, b, x2);

% mostro i risultati
figure
stem(n,y,'b-x')
hold on
stem(n,y2,'r-o')
hold off
grid on
xlabel('tempo n')
ylabel('numerosità della popolazione y(n)')
title('verifica della tempo-varianza con b(n)=0.5+0.5*sin(n*pi/10)')
legend('x(n) = 50 delta(n-1)','x(n) = 50 delta(n-5)')