Namespaces
Variants
Actions

Difference between revisions of "Fourier method"

From Encyclopedia of Mathematics
Jump to: navigation, search
m (Convert image to LaTeX)
m (label)
 
(One intermediate revision by the same user not shown)
Line 1: Line 1:
 +
{{TEX|done}}
 
''method of separation of variables''
 
''method of separation of variables''
  
 
A method for finding particular solutions of differential equations
 
A method for finding particular solutions of differential equations
  
$$  \begin{equation}  (Lu)(x,y)=Mu-Nu=0, \quad x \in \mathbf{R}^n, \quad y \in \mathbf{R}^m \tag{1} \end{equation} $$
+
$$  \begin{equation}  (Lu)(x,y)=Mu-Nu=0, \quad x \in \mathbf{R}^n, \quad y \in \mathbf{R}^m, \label{1}\tag{1} \end{equation} $$
  
where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410702.png" /> (respectively, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410703.png" />) are linear differential expressions containing only derivatives with respect to the variables <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410704.png" /> (respectively, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410705.png" />) with coefficients also depending only on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410706.png" /> (respectively, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410707.png" />). A function
+
where $M$ (respectively, $N$) are linear differential expressions containing only derivatives with respect to the variables $x$ (respectively, $y$) with coefficients also depending only on $x$ (respectively, $y$). A function
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410708.png" /></td> <td valign="top" style="width:5%;text-align:right;">(2)</td></tr></table>
+
$$u(x,y)=v(x)w(y)\label{2}\tag{2}$$
  
is a solution of (1) if there is a constant <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f0410709.png" /> such that
+
is a solution of \eqref{1} if there is a constant $\lambda$ such that
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107010.png" /></td> <td valign="top" style="width:5%;text-align:right;">(3)</td></tr></table>
+
$$Mv+\lambda v=0,\quad Nw+\lambda w=0.\label{3}\tag{3}$$
  
 
For example, for the equation of a vibrating string
 
For example, for the equation of a vibrating string
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107011.png" /></td> <td valign="top" style="width:5%;text-align:right;">(4)</td></tr></table>
+
$$u_{xx}-u_{yy}=0,\quad m=n=1,\label{4}\tag{4}$$
  
<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107012.png" />, and the solution (2) takes the form
+
$M\phi=N\phi=\phi''$, and the solution \eqref{2} takes the form
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107013.png" /></td> <td valign="top" style="width:5%;text-align:right;">(5)</td></tr></table>
+
$$u(x,y)=\label{5}\tag{5}$$
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107014.png" /></td> </tr></table>
+
$$=\begin{cases}(c_1x+c_2)(c_3y+c_4),&\lambda=0,\\(c_1e^{\mu x}+c_2e^{-\mu x})(c_3e^{\mu y}+c_4e^{-\mu y}),&\lambda=-\mu^2,\\(c_1\cos\mu x+c_2\sin\mu x)(c_3\cos\mu y+c_4\sin\mu y),&\lambda=\mu^2,\end{cases}$$
  
where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107015.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107016.png" />, are arbitrary constants and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107017.png" />. The solutions (5) satisfy the boundary conditions
+
where $c_i$, $1\leq i\leq4$, are arbitrary constants and $\mu>0$. The solutions \eqref{5} satisfy the boundary conditions
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107018.png" /></td> <td valign="top" style="width:5%;text-align:right;">(6)</td></tr></table>
+
$$u(0,y)=u(\pi,y)=0,\quad y\geq0,\label{6}\tag{6}$$
  
when <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107019.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107020.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107021.png" />. The so-called Fourier series composed of these functions,
+
when $c_1=0$, $\lambda=\mu^2$, $\mu=1,2,\ldots$. The so-called Fourier series composed of these functions,
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107022.png" /></td> <td valign="top" style="width:5%;text-align:right;">(7)</td></tr></table>
+
$$U(x,y)=\sum_{n=1}^\infty\sin nx(a_n\cos ny+b_n\sin ny),\label{7}\tag{7}$$
  
gives the solution of the initial-boundary value problem (4), (6) and initial condition
+
gives the solution of the initial-boundary value problem \eqref{4}, \eqref{6} and initial condition
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107023.png" /></td> <td valign="top" style="width:5%;text-align:right;">(8)</td></tr></table>
+
$$U(x,0)=\phi(x),\quad \frac{\partial U}{\partial y}(x,0)=\psi(x),\quad0\leq x\leq\pi,\label{8}\tag{8}$$
  
if <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107024.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107025.png" /> are the Fourier coefficients
+
if $a_n$ and $nb_n$ are the Fourier coefficients
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107026.png" /></td> </tr></table>
+
$$a_n=\frac2\pi\int\limits_0^\pi\phi(x)\sin nxdx,$$
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107027.png" /></td> </tr></table>
+
$$nb_n=\frac2\pi\int\limits_0^\pi\psi(x)\sin nxdx$$
  
of the data <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107028.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107029.png" /> which are assumed to be sufficiently smooth. In a similar way solutions of initial-boundary value problems can be obtained for more general classes of equations (1), where the role of the theory of Fourier series associated with the expansion (7) is played by the spectral theory of linear operators.
+
of the data $\phi$ and $\psi$ which are assumed to be sufficiently smooth. In a similar way solutions of initial-boundary value problems can be obtained for more general classes of equations \eqref{1}, where the role of the theory of Fourier series associated with the expansion \eqref{7} is played by the spectral theory of linear operators.
  
The Fourier method is closely connected with the special functions solving equations (3) for particular cases of operators <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107030.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107031.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107032.png" />. Many of these functions originally arose in this way. For example, when applying this method to the Helmholtz equation
+
The Fourier method is closely connected with the special functions solving equations \eqref{3} for particular cases of operators $M$ and $N$, $m=n=1$. Many of these functions originally arose in this way. For example, when applying this method to the Helmholtz equation
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107033.png" /></td> </tr></table>
+
$$u_{xx}+u_{yy}+u=0,$$
  
 
one obtains after transition to polar coordinates
 
one obtains after transition to polar coordinates
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/f/f041/f041070/f04107034.png" /></td> </tr></table>
+
$$M=r^2\frac{\partial^2}{\partial r^2}+r\frac{\partial}{\partial r}+r^2,\quad N=-\frac{\partial^2}{\partial\theta^2},$$
  
and the first equation in (3) becomes the Bessel differential equation.
+
and the first equation in \eqref{3} becomes the Bessel differential equation.
  
Generally speaking, a differential equation has a whole family of coordinate systems in which it admits separation of variables, that is, reduces to the form (1). The problem of finding such coordinate systems is closely connected with group properties of differential equations. Applying methods from the theory of Lie groups enables one to describe all the separable solutions of many classical equations in mathematical physics (those of Laplace, Helmholtz and Schrödinger, the wave equation, etc.). In this way also a whole series of relations is obtained in the theory of special functions.
+
Generally speaking, a differential equation has a whole family of coordinate systems in which it admits separation of variables, that is, reduces to the form \eqref{1}. The problem of finding such coordinate systems is closely connected with group properties of differential equations. Applying methods from the theory of Lie groups enables one to describe all the separable solutions of many classical equations in mathematical physics (those of Laplace, Helmholtz and Schrödinger, the wave equation, etc.). In this way also a whole series of relations is obtained in the theory of special functions.
  
 
The method of separation of variables was suggested by J. d'Alembert (1749) for solving the wave equation, the method was developed fairly thoroughly at the beginning of the 19th century by J. Fourier and was formulated in complete generality by M.V. Ostrogradski in 1828.
 
The method of separation of variables was suggested by J. d'Alembert (1749) for solving the wave equation, the method was developed fairly thoroughly at the beginning of the 19th century by J. Fourier and was formulated in complete generality by M.V. Ostrogradski in 1828.

Latest revision as of 15:15, 14 February 2020

method of separation of variables

A method for finding particular solutions of differential equations

$$ \begin{equation} (Lu)(x,y)=Mu-Nu=0, \quad x \in \mathbf{R}^n, \quad y \in \mathbf{R}^m, \label{1}\tag{1} \end{equation} $$

where $M$ (respectively, $N$) are linear differential expressions containing only derivatives with respect to the variables $x$ (respectively, $y$) with coefficients also depending only on $x$ (respectively, $y$). A function

$$u(x,y)=v(x)w(y)\label{2}\tag{2}$$

is a solution of \eqref{1} if there is a constant $\lambda$ such that

$$Mv+\lambda v=0,\quad Nw+\lambda w=0.\label{3}\tag{3}$$

For example, for the equation of a vibrating string

$$u_{xx}-u_{yy}=0,\quad m=n=1,\label{4}\tag{4}$$

$M\phi=N\phi=\phi''$, and the solution \eqref{2} takes the form

$$u(x,y)=\label{5}\tag{5}$$

$$=\begin{cases}(c_1x+c_2)(c_3y+c_4),&\lambda=0,\\(c_1e^{\mu x}+c_2e^{-\mu x})(c_3e^{\mu y}+c_4e^{-\mu y}),&\lambda=-\mu^2,\\(c_1\cos\mu x+c_2\sin\mu x)(c_3\cos\mu y+c_4\sin\mu y),&\lambda=\mu^2,\end{cases}$$

where $c_i$, $1\leq i\leq4$, are arbitrary constants and $\mu>0$. The solutions \eqref{5} satisfy the boundary conditions

$$u(0,y)=u(\pi,y)=0,\quad y\geq0,\label{6}\tag{6}$$

when $c_1=0$, $\lambda=\mu^2$, $\mu=1,2,\ldots$. The so-called Fourier series composed of these functions,

$$U(x,y)=\sum_{n=1}^\infty\sin nx(a_n\cos ny+b_n\sin ny),\label{7}\tag{7}$$

gives the solution of the initial-boundary value problem \eqref{4}, \eqref{6} and initial condition

$$U(x,0)=\phi(x),\quad \frac{\partial U}{\partial y}(x,0)=\psi(x),\quad0\leq x\leq\pi,\label{8}\tag{8}$$

if $a_n$ and $nb_n$ are the Fourier coefficients

$$a_n=\frac2\pi\int\limits_0^\pi\phi(x)\sin nxdx,$$

$$nb_n=\frac2\pi\int\limits_0^\pi\psi(x)\sin nxdx$$

of the data $\phi$ and $\psi$ which are assumed to be sufficiently smooth. In a similar way solutions of initial-boundary value problems can be obtained for more general classes of equations \eqref{1}, where the role of the theory of Fourier series associated with the expansion \eqref{7} is played by the spectral theory of linear operators.

The Fourier method is closely connected with the special functions solving equations \eqref{3} for particular cases of operators $M$ and $N$, $m=n=1$. Many of these functions originally arose in this way. For example, when applying this method to the Helmholtz equation

$$u_{xx}+u_{yy}+u=0,$$

one obtains after transition to polar coordinates

$$M=r^2\frac{\partial^2}{\partial r^2}+r\frac{\partial}{\partial r}+r^2,\quad N=-\frac{\partial^2}{\partial\theta^2},$$

and the first equation in \eqref{3} becomes the Bessel differential equation.

Generally speaking, a differential equation has a whole family of coordinate systems in which it admits separation of variables, that is, reduces to the form \eqref{1}. The problem of finding such coordinate systems is closely connected with group properties of differential equations. Applying methods from the theory of Lie groups enables one to describe all the separable solutions of many classical equations in mathematical physics (those of Laplace, Helmholtz and Schrödinger, the wave equation, etc.). In this way also a whole series of relations is obtained in the theory of special functions.

The method of separation of variables was suggested by J. d'Alembert (1749) for solving the wave equation, the method was developed fairly thoroughly at the beginning of the 19th century by J. Fourier and was formulated in complete generality by M.V. Ostrogradski in 1828.

References

[1] A.V. Bitsadze, "Equations of mathematical physics" , MIR (1980) (Translated from Russian)
[2] U. Miller, "Symmetry and separation of variables" , Addison-Wesley (1977)


Comments

References

[a1] G.W. Bluman, G.D. Cole, "Similarity methods for differential equations" , Springer (1974)
How to Cite This Entry:
Fourier method. Encyclopedia of Mathematics. URL: http://encyclopediaofmath.org/index.php?title=Fourier_method&oldid=32198
This article was adapted from an original article by A.P. Soldatov (originator), which appeared in Encyclopedia of Mathematics - ISBN 1402006098. See original article