%% Solução do exemplo 3.8 do Lathi

clear all; close all; clc

%% Dados do problema

% Equação diferença:
%   y[n] - 0.5 y[n-1] = x[n]
%   y[n] = x[n] + 0.5 y[n-1]
%
% Sinal de entrada:
%   x[n] = n.^2 .* u[n]
%
% Condição inicial:
%   y[-1] = 16;

%% Solução

% Vetor de amostras
n = -1:20;

% Sinal de entrada
x = n.^2 .* u(n);

% Sinal de saída
y = 16;

% Solução interativa
for i = 0:20
    
    % Adequando o índice à notação vetorial do Matlab
    j = i+2;
    
    y(j) = x(j) + 0.5*y(j-1);
    
end

%% Exibição dos resultados

% Plotando o sinal de entrada
figure
stem(n,x)
grid on
title('Sinal de entrada')

% Plotando o sinal de saída
figure
stem(n,y)
grid on
title('Sinal de saída')














