clear all; close all; clc

%% Eixo do tempo para as operações

n = -10:10;

%% Impulso unitário

figure
stem(n,d(n)); xlabel('n'); grid on; title('\delta[n]')

%% Degrau

figure
stem(n,u(n)); xlabel('n'); grid on; title('u[n]')

%% Exponencial com base real

n = -50:50;

% Lambda maior ou menor que zero
lambda = [0.5 -0.5 0];

figure
for i = 1:length(lambda)
    
    % Cálculo do exponencial
    x = exp(lambda(i)*n);
    
    % Exibindo os resultados
    subplot(3,1,i); stem(n,x); xlabel('n'); grid on; ylabel('e^{\lambda n}'); title(['\lambda = ' num2str(lambda(i))])
    
end

%% Exponencial com base complexa

n = -10:10;

% Parte real do lambda maior ou menor que zero
a = [0.5 -0.5 0];

lambda = a + 1i*0.7;

figure
for i = 1:length(lambda)
    
    % Cálculo do exponencial
    x = exp(lambda(i)*n);
    
    % Exibindo os resultados
    subplot(3,1,i); stem(n,x); xlabel('n'); grid on; ylabel('e^{\lambda n}'); title(['\lambda = ' num2str(lambda(i))])
    
end








