clear all
clf
N=10;
L=2^N;
m=L/8;
n=[1:L];
x=[ones(1,m)/sqrt(m),zeros(1,L-m)];
figure(1)
plot(x)
figure(2)
y=fft(x)/sqrt(L);
plot(real([y y]))
freq=2*pi/L;
%sincx=sin(m*n*freq)./(m*2*n*freq);
%hold on
%plot(sincx,'r')