close all
f = 1;
T = 1/f;
t = 0 : T/10 : 10*T;
function1 = sin(2 * pi * f * t); %обычная функция
fd = 2*f; %новая частота по заданию
Td = 1/fd;
td = 0 : Td : 10*T;
function2 = sin(2 * pi * f * td);%дискретизированная ф-ия
plot(t,function1,td, function2);%строим графики
figure(2);

temp = zeros(length(td),length(t));%восстановление функции
for i = 1: length(td)
 temp(i, :) = sinc(((t - td(i)) * fd)) .* function2(i);%ф-ла 1.72 модифицированная
end
yr = sum(temp);%складываем отдельные синусоиды
plot(t,function1, td, function2, t, yr);%выводим все три графика
%поиграйся с частотой fd - больше 2fd сигнал почти восстанавливается