Fourier sine and cosine series

Contents

The cosine series of x

Let f(x) = x. First I calculate the Fourier cosine series of f on the interval [0,1].

syms x k n
evalin(symengine,'assume(k,Type::Integer)');
f = x
 
f =
 
x
 

The following commands compute the partial sum of the Fourier cosine series of f.

a = @(k) 2*int(x*cos(k*pi*x),x,0,1);
fourier_cosine_partial_sum = @(x,n) a(0)/2 + ...
   symsum(a(k)*cos(k*pi*x),k,1,n);

Here are plots of the partial sums for n=2,5,10.

ezplot(f,0,1), hold on
ezplot(fourier_cosine_partial_sum(x,2),0,1), hold off
title('x and partial sum of its cosine series with n=2')
ezplot(f,0,1), hold on
ezplot(fourier_cosine_partial_sum(x,5),0,1), hold off
title('x and partial sum of its cosine series with n=5')
ezplot(f,0,1), hold on
ezplot(fourier_cosine_partial_sum(x,10),0,1), hold off
title('x and partial sum of its cosine series with n=10')

Recall that the cosine series is the Fourier series of the even extension of f. The even extension to [-1,1] is just abs(x). Here are plots of abs(x) and the Fourier cosine series of x on [-1,1].

ezplot(abs(x),-1,1), hold on
ezplot(fourier_cosine_partial_sum(x,2),-1,1), hold off
title('|x| and partial sum of its Fourier series with n=2')
ezplot(abs(x),-1,1), hold on
ezplot(fourier_cosine_partial_sum(x,5),-1,1), hold off
title('|x| and partial sum of its Fourier series with n=5')

The sine series of x

The following commands calculate the nth partial sum of the Fourier sine series of f.

b = @(k) 2*int(x*sin(k*pi*x),x,0,1);
fourier_sine_partial_sum = @(x,n) symsum(b(k)*sin(k*pi*x),k,1,n);

Here are plots of the partial sums for n = 2,5,10.

ezplot(f,0,1), hold on
ezplot(fourier_sine_partial_sum(x,2),0,1), hold off
title('x and partial sum of its sine series with n=2')
ezplot(f,0,1), hold on
ezplot(fourier_sine_partial_sum(x,5),0,1), hold off
title('x and partial sum of its sine series with n=5')
ezplot(f,0,1), hold on
ezplot(fourier_sine_partial_sum(x,10),0,1), hold off
title('x and partial sum of its sine series with n=10')

Recall that the Fourier sine series is the Fourier series of the odd extension of f. The odd extension of f to [-1,1] is x.

ezplot(x,-1,1), hold on
ezplot(fourier_sine_partial_sum(x,10),-1,1), hold off
title('x and partial sum of its Fourier series with n=10')