function [f] = sinc_reconstruct(t, t_sampled, F_sampled, Ns, dts)
% Function to calculate a sinc resconstruction filter
  f = 0;
  for i = 1 : Ns+1
    f = f + F_sampled(i)*sinc((t-t_sampled(i))/dts);
  end
end