% geração de funções degrau unitário e delta de Dirac,
% a partir das funções d índice k e f índice k

clc
clearvars
close all

% variáveis
ls=10;
k=5; % 5 (início), 0.5 e 0.05 (final)
s=linspace(-ls,ls,500);

% função dk
dk=stepfun(s,k/2);
dk=dk+((s/k)+0.5).*(stepfun(s,-k/2)-stepfun(s,k/2));

% função fk
fk=zeros(1,length(s));
fk=fk+(1/k)*(stepfun(s,-k/2)-stepfun(s,k/2));

% visualização
figure(1)
subplot(2,1,1)
plot(s,dk,'b')
hold on
plot([0 0],[0 1*1.1],'k-')
hold off
axis([-ls,ls,0,1*1.1])
grid on
title('Função d_k')
grid on
subplot(2,1,2)
plot(s,fk,'r')
hold on
plot([0 0],[0 (1/k)*1.1],'k-')
hold off
axis([-ls,ls,0,(1/k)*1.1])
grid on
title('Função f_k')
grid on