How close is my Fourier Sine series to a function f(x)?
Answer: L_2(N) = -\frac{1}{2} a_n^2 + \int_\Omega f^2(x) dx
This question came up recently while discussing PDE (partial differential equations) solution techniques. The final result is quite interesting.
Our sine series approximates the function f(x) over the domain \Omega as
f(x) \approx \sum {{a_n}\sin (n\pi x)} You can determine the coefficients with the formula. a_n = 2\int_\Omega f(x) sin(n\pi x) dx
One good measure of the error between the sine series and the function is the L_2, pronounced "el squared", error. L_2 = \int_{\Omega} \left( u(x) - f(x) \right)^2 dx
Substitute for the sine series to obtain. L_2 = \int_{\Omega} \left({{a_n}\sin (n\pi x)}- f(x) \right)^2 dx
Expand the terms to obtain L_2 = \int_{\Omega} a_n^2 \sin^2(n\pi x) dx -2\int_{\Omega} a_n f(x) \sin(n \pi x) dx + \int_{\Omega} f^2(x) dx
Now, these integrals are particularly interesting. The first integral is a constant 0.5 a_n^2. The second contains the definition of a_n. The third only contains the function f^2(x). This simplifies to L_2(N) = -\frac{1}{2} a_n^2 + \int_\Omega f^2(x) dx
Interestingly, if we let the error become zero, the following is an identity a_n^2 = 2 \int_\Omega f^2(x) dx Knowledge of this identity will allow you to quickly compute integrals of squared trig functions.
Neat! Do you have any interesting Fourier results? Let me know in the comments.
No comments:
Post a Comment