clf
clear all
set(0,'defaultaxesfontsize',22);
t=-6:.001:6;
t=t+0.0000001;
%
x=sin(pi*t);
%
plot(t,x,'Linewidth',3)
Ts=0.5;
ts=-6.5:Ts:6.5;
%
ts=ts-0.25;
xs=sin(pi*ts);
%
xr=zeros(size(t));
for k=1:length(ts),
xr=xr+Ts*xs(k)*(sin((pi/Ts)*(t-ts(k)))./(pi*(t-ts(k))));
end
plot(t,x,'Linewidth',3)
grid
axis([min(t) max(t) -2 2])
set(gca,'XTick',[min(t):0.5:max(t)],'YTick',[-1:5])
hold on
plot([0 0],[-2 2],'k','Linewidth',3)
plot([-6 6],[0 0],'k','Linewidth',2)
plot(ts,xs,'rx','MarkerSize',18,'Linewidth',4)
plot(t,xr,'r','Linewidth',3)
for kk=1:length(ts),
Y(kk,:)=Ts*xs(kk)*(sin((pi/Ts)*(t-ts(kk)))./(pi*(t-ts(kk))));
end
ko=7;
plot(t,Y(4+ko,:),'--',t,Y(5+ko,:),'--',t,Y(6+ko,:),'--',t,Y(7+ko,:),'--', ...
     t,Y(8+ko,:),'--',t,Y(9+ko,:),'--',t,Y(10+ko,:),'--','Linewidth',2)
title('Sampling Theory');
hold off

