close all;

N = 1000;
x = linspace(-5,5,N);

amp = exp(-abs(x));
plot(x,amp,'--'); hold on;
plot(x,-amp,'--');
y = amp.*sin(3/2*pi*x);
plot(x,y);

x1 = -5:0.5:5;
y1 = exp(-abs(x1)).*sin(3/2*pi*x1);
stem(x1,y1,'r'); hold on;

x2 = -4.75:0.5:4.75;
y2 = exp(-abs(x2)).*sin(3/2*pi*x2);
stem(x2,y2,'k');
fs = 14;
title('Illustration of spin echo aliasing','fontsize',fs);
ylabel('Spin echo (arb. scale)','fontsize',fs);
xlabel('Time (arb. scale)','fontsize',fs);