F = 5;
Fd = 1000*F;
Td = 1/Fd;
t = 0:Td:(10/F);

s = zeros(1, length(t));

for i=1:length(t)
    s(i)=sin(2*pi*F*t(i));
end

d = 1/1000;
T = 0:d:(10/F);
res = zeros(1, length(T));

for i=1:length(T)
    for n=1:length(s)
        res(i) = res(i) + s(n) * sinc((T(i)/Td) - (n-1));
    end
end

figure();
hold on;
plot(t, s);
plot(T, res, "r--");
xlabel("t");
ylabel("u(t)");
legend("Исходный сигнал", "Восстановленный сигнал");