clc
clear all
n=0:9999;
x = sin(pi/5*n);
spike = repmat([1 zeros(1,9)],1,1000);
x1 = x + spike;
plot(abs(fftshift(fft(x1,1000))))