サンプリング定理

clear; close all; clc;

% fs = 3;
time = linspace(0,1,100);

plot(time,SF(time,3),'DisplayName','Partial Sum (fs=3)');
hold on ;
plot(time,SF(time,4),'DisplayName','Partial Sum (fs=4)');
plot(time,SF(time,5),'DisplayName','Partial Sum (fs=5)');

plot(time,sin(2*pi*time),'DisplayName','x(t)');
xlabel('time');
ylabel('signal & partial sum');
legend();
big;

function ps=SF(time,fs)
  ps = 0;
  t = linspace(0,1,fs+1)
  for k=1:length(t)
      ps = ps + sin(2*pi*t(k))*sinc(fs*(time-t(k)));
  end  
end