Re: Integrals of Fourier Series
- Subject: [mg3255] Re: Integrals of Fourier Series
- From: paul at earwax.pd.uwa.edu.au (Paul Abbott)
- Date: 21 Feb 1996 13:09:17 -0600
- Approved: usenet@wri.com
- Distribution: local
- Newsgroups: wri.mathgroup
- Organization: University of Western Australia
- Sender: daemon at wri.com
George Oster wrote:
> Suppose I want to substitue a Fourier Series into an integral:
>
> u[x_] := Sum[A[n] Sin[n Pi x/L], {n, 1, Infinity}]
>
> Integrate[(u''[x])^2, {x, 0, L}]
>
> This has an easy analytical solution that I can't get Mma to find, > because Mma doesn't know that Sum and Integrate commute, and that
> Sin[n Pi] = 0 for all integer n.
>
> How to do this?
Instead of working with the Sum, just focus on the summand:
Remove[u]
u[n_][x_] = Sin[n Pi x/L];
The generic (cross-product) term of the integrand is
u[n]''[x] u[m]''[x]
2 2 4 m Pi x n Pi x
m n Pi Sin[------] Sin[------]
L L
---------------------------------
4
L
and the integral can be evaluated to yield
Integrate[%, {x, 0, L}]
2 2 3
(m n Pi (m Sin[(m - n) Pi] + n Sin[(m - n) Pi] -
m Sin[(m + n) Pi] + n Sin[(m + n) Pi])) /
3 2 2
(2 L (m - n ))
When m!=n,
% /. Sin[n_ Pi] -> 0
0
Taking the limit as m->n:
Limit[%%, m->n] /. Sin[n_ Pi] -> 0
4 4
n Pi
------
3
2 L
one finds that the doubly-infinite summation reduces to
Sum[a[n]^2 %, {n, 1, Infinity}]
4 4 2
n Pi a[n]
Sum[------------, {n, 1, Infinity}]
3
2 L
_________________________________________________________________
Paul Abbott
Department of Physics Phone: +61-9-380-2734
The University of Western Australia Fax: +61-9-380-1014
Nedlands WA 6907 paul at physics.uwa.edu.au
AUSTRALIA http://www.pd.uwa.edu.au/Paul
_________________________________________________________________