%	test sul calcolo dei coefficienti di Fourier

%	Giuseppe Rodriguez,  16-07-10
%	Email: rodriguez@unica.it

%fun = @(x) x;
fun = @(x) x.^2;
%fun = @(x) sin(2*x)+.5*cos(5*x);

n = 8;
a = -pi;
b = pi;
m = 128;

t = [0:n-1]'/n * (b-a) + a;
f = fun(t);
fs = fftshift(f);
c = 1/n * fft(fs);
c = fftshift(c);
k = [-n/2:n/2-1]';

figure(1)
clf
stem(k,abs(c),'b*')
title('Modulo dei coefficienti di Fourier')

%cex = [zeros((m-n)/2,1); c; zeros((m-n)/2,1)];
cex = fftshift(c);
cex = [c(1:n/2); zeros(m-n,1); c(n/2+1:n)];
sf = m*ifft(cex);
sf = fftshift(sf);
x = [0:m-1]'/m * (b-a) + a;

figure(2)
clf
plot(x,[fun(x) sf])
title('Serie di Fourier')
