% fft_demo
n=64;
x1=linspace(0,2*pi,n+1)';
x=x1(1:n);
h1=ones(n+1,1);
h1(1:n/2+1)=x1(1:n/2+1);
%%h1(n/2+1:n+1)=-ones(n/2+1,1);
%%h1(n/2+1)=0;
h1(n/2+2:n+1)=2*pi-x1(n/2+2:n+1);
figure(1)
%%plot(x1,h1,'o',[0 pi],[1 1],[pi,2*pi],[-1 -1]);
%%axis([0 7 -2 2]);
plot(x1,h1,'.',[0 pi],[0 pi],[pi,2*pi],[pi 0]);
axis([0 7 0 7]);
hold on;
h=h1(1:n);
c=fft(h);
x4=linspace(0,2*pi,4*n+1);
c4=zeros(4*n,1);
c4(1:n/2)=c(1:n/2);
c4(3*n+n/2+1:4*n)=c(n/2+1:n);
h4=4*ifft(c4);
h41=[h4;h4(1)];
plot(x4,h41);
figure(2)
%hold off;
loglog(1:n/2,abs(c(1:n/2)/n),'.')
axis([1 n/2 10^(-5) 1]);