Как найти оператор лапласа

Текущая версия страницы пока не проверялась опытными участниками и может значительно отличаться от версии, проверенной 19 марта 2022 года; проверки требуют 2 правки.

Опера́тор Лапла́са (лапласиа́н, оператор дельта) — дифференциальный оператор, действующий в линейном пространстве гладких функций и обозначаемый символом  Delta. Функции F он ставит в соответствие функцию

{displaystyle Delta F={partial ^{2}F over partial x_{1}^{2}}+{partial ^{2}F over partial x_{2}^{2}}+ldots +{partial ^{2}F over partial x_{n}^{2}}}

в n-мерном пространстве.

Оператор Лапласа эквивалентен последовательному взятию операций градиента и дивергенции: Delta=operatorname{div},operatorname{grad}, таким образом, значение оператора Лапласа в точке может быть истолковано как плотность источников (стоков) потенциального векторного поля  operatorname{grad}F в этой точке. В декартовой системе координат оператор Лапласа часто обозначается следующим образом Delta=nablacdotnabla=nabla^2[1], то есть в виде скалярного произведения оператора набла на себя. Оператор Лапласа симметричен.

Оператор Лапласа для вектора {displaystyle mathbf {A} =A_{x}mathbf {i} +A_{y}mathbf {j} +A_{z}mathbf {k} }:

{displaystyle Delta mathbf {A} ={Delta }A_{x}mathbf {i} +{Delta }A_{y}mathbf {j} +{Delta }A_{z}mathbf {k} ={biggl (}{frac {partial ^{2}A_{x}}{partial x^{2}}}+{frac {partial ^{2}A_{x}}{partial y^{2}}}+{frac {partial ^{2}A_{x}}{partial z^{2}}}{Biggr )}mathbf {i} +{biggl (}{frac {partial ^{2}A_{y}}{partial x^{2}}}+{frac {partial ^{2}A_{y}}{partial y^{2}}}+{frac {partial ^{2}A_{y}}{partial z^{2}}}{Biggr )}mathbf {j} +{biggl (}{frac {partial ^{2}A_{z}}{partial x^{2}}}+{frac {partial ^{2}A_{z}}{partial y^{2}}}+{frac {partial ^{2}A_{z}}{partial z^{2}}}{Biggr )}mathbf {k} }[2]

Лапласиан вектора – тоже вектор.

Другое определение оператора Лапласа[править | править код]

Оператор Лапласа является естественным обобщением на функции нескольких переменных обычной второй производной функции одной переменной. В самом деле, если функция  f (x) имеет в окрестности точки  x_0 непрерывную вторую производную  f''(x), то, как это следует из формулы Тейлора

 f(x_0+r)=f(x_0)+rf'(x_0)+frac{r^2}{2}f''(x_0)+o(r^2), при rto 0,,
 f(x_0-r)=f(x_0)-rf'(x_0)+frac{r^2}{2}f''(x_0)+o(r^2), при rto 0,

вторая производная есть предел

 f''(x_0)=limlimits_{r to 0} frac{2}{r^2} left{ frac{f(x_0+r)+f(x_0-r)}{2}-f(x_0) right}.

Если, переходя к функции  F от  k переменных, поступить таким же образом, то есть для заданной точки  M_0(x_1^0,x_2^0, ... ,x_k^0) рассматривать её  k -мерную шаровую окрестность  Q_r радиуса  r и разность между средним арифметическим

 frac{1}{sigma(S_r)}intlimits_{S_r}Fdsigma

функции  F на границе  S_r такой окрестности с площадью границы  sigma(S_r) и значением  F(M_0) в центре этой окрестности  M_0, то в случае непрерывности вторых частных производных функции  F в окрестности точки  M_0 значение лапласиана  Delta F в этой точке есть предел

 Delta F(M_0)=limlimits_{r to 0} frac{2k}{r^2} left{frac{1}{sigma(S_r)}intlimits_{S_r}F(M)dsigma -F(M_0) right}.

Одновременно с предыдущим представлением для оператора Лапласа функции  F, имеющей непрерывные вторые производные, справедлива формула

 Delta F(M_0)=limlimits_{r to 0} frac{2(k+2)}{r^2} left{frac{1}{omega(Q_r)}intlimits_{Q_r}F(M)domega -F(M_0) right}, где  omega(Q_r) — объём окрестности  Q_r.

Эта формула выражает непосредственную связь лапласиана функции с её объёмным средним в окрестности данной точки.

Доказательство этих формул можно найти, например, в[3].

Вышеизложенные пределы, во всех случаях, когда они существуют, могут служить определением оператора Лапласа функции  F. Такое определение предпочтительнее обычного определения лапласиана, предполагающего существование вторых производных рассматриваемых функций, и совпадает с обычным определением в случае непрерывности этих производных.

Выражения для оператора Лапласа в различных криволинейных системах координат[править | править код]

В произвольных ортогональных криволинейных координатах в трёхмерном пространстве q_1, q_2, q_3:

Delta f (q_1, q_2, q_3) = operatorname{div},operatorname{grad},f(q_1, q_2, q_3) =
=frac{1}{H_1H_2H_3}left[ frac{partial}{partial q_1}left( frac{H_2H_3}{H_1}frac{partial f}{partial q_1} right) + frac{partial}{partial q_2}left( frac{H_1H_3}{H_2}frac{partial f}{partial q_2} right) +  frac{partial}{partial q_3}left( frac{H_1H_2}{H_3}frac{partial f}{partial q_3} right)right],
где H_i  — коэффициенты Ламе.

Цилиндрические координаты[править | править код]

В цилиндрических координатах вне прямой  r=0:

 Delta f 
= {1 over r} {partial over partial r}
  left( r {partial f over partial r} right) 
+ {partial^2f over partial z^2}
+ {1 over r^2} {partial^2 f over partial varphi^2}

Сферические координаты[править | править код]

В сферических координатах вне начала отсчёта (в трёхмерном пространстве):

 Delta f 
= {1 over r^2} {partial over partial r}
  left( r^2 {partial f over partial r} right) 
+ {1 over r^2 sin theta} {partial over partial theta}
  left( sin theta {partial f over partial theta} right) 
+ {1 over r^2sin^2 theta} {partial^2 f over partial varphi^2}

или

 Delta f 
= {1 over r} {partial^2 over partial r^2}
  left( rf right) 
+ {1 over r^2 sin theta} {partial over partial theta}
  left( sin theta {partial f over partial theta} right) 
+ {1 over r^2 sin^2 theta} {partial^2 f over partial varphi^2}.

В случае если  f=f(r) в n-мерном пространстве:

 Delta f =  {d^2 fover dr^2} + {n-1 over r } {dfover dr}.

Параболические координаты[править | править код]

В параболических координатах (в трёхмерном пространстве) вне начала отсчёта:


Delta f= frac{1}{sigma^{2} + tau^{2}} 
left[
frac{1}{sigma} frac{partial }{partial sigma} 
left( sigma frac{partial f}{partial sigma} right) +
frac{1}{tau} frac{partial }{partial tau} 
left( tau frac{partial f}{partial tau} right)right] +
frac{1}{sigma^2tau^2}frac{partial^2 f}{partial varphi^2}

Цилиндрические параболические координаты[править | править код]

В координатах параболического цилиндра вне начала отсчёта:

Delta F(u,v,z)={frac  {1}{c^{2}(u^{2}+v^{2})}}left[{frac  {partial ^{2}F}{partial u^{2}}}+{frac  {partial ^{2}F}{partial v^{2}}}right]+{frac  {partial ^{2}F}{partial z^{2}}}.

Общие криволинейные координаты и римановы пространства[править | править код]

Пусть на гладком многообразии X задана локальная система координат и g_{ij} — риманов метрический тензор на X, то есть метрика имеет вид

ds^2 =sum^n_{i,j=1}g_{ij} dx^idx^j .

Обозначим через g^{ij} элементы матрицы (g_{ij})^{-1} и

g = operatorname{det} g_{ij} = (operatorname{det} g^{ij})^{-1}.

Дивергенция векторного поля F, заданного координатами F^i (и представляющего дифференциальный оператор первого порядка sum_i F^ifrac{partial}{partial x^i}) на многообразии X вычисляется по формуле

operatorname{div} F = frac{1}{sqrt{g}}sum^n_{i=1}frac{partial}{partial x^i}(sqrt{g}F^i),

а компоненты градиента функции f — по формуле

(nabla f)^{j}=sum _{{i=1}}^{n}g^{{ij}}{frac  {partial f}{partial x^{i}}}.

Оператор Лапласа — Бельтрами на X:

Delta f=operatorname {div}(nabla f)={frac  {1}{{sqrt  {g}}}}sum _{{i=1}}^{n}{frac  {partial }{partial x^{i}}}{Big (}{sqrt  {g}}sum _{{k=1}}^{n}g^{{ik}}{frac  {partial f}{partial x^{k}}}{Big )}.

Значение Delta f является скаляром, то есть не изменяется при преобразовании координат.

Применение[править | править код]

С помощью данного оператора удобно записывать уравнения Лапласа, Пуассона и волновое уравнение. В физике оператор Лапласа применим в электростатике и электродинамике, квантовой механике, во многих уравнениях физики сплошных сред, а также при изучении равновесия мембран, плёнок или поверхностей раздела фаз с поверхностным натяжением (см. Лапласово давление), в стационарных задачах диффузии и теплопроводности, которые сводятся, в непрерывном пределе, к обычным уравнениям Лапласа или Пуассона или к некоторым их обобщениям.

Вариации[править | править код]

  • Оператор Д’Аламбера — обобщение оператора Лапласа для гиперболических уравнений. Включает в себя вторую производную по времени.
  • Векторный оператор Лапласа — обобщение оператора Лапласа на случай векторного аргумента.

См. также[править | править код]

  • Оператор набла
  • Оператор Д’Аламбера
  • Уравнение Лапласа
  • Гармоническая функция
  • Матрица Кирхгофа
  • В.Г.Воднев, А.Ф.Наумович, Н.Ф.Наумович “Математический словарь высшей школы”. Издательство МПИ 1984.

Примечания[править | править код]

  1. Стоит избегать обозначения для оператора Лапласа в виде квадрата оператора набла, поскольку из такой записи непонятно, скалярное или векторное произведение подразумевается под возведением в квадрат.
  2. В.Г.Воднев, А.Ф.Наумович, Н.Ф.Наумович “Математический словарь высшей школы”. Издательство МПИ 1984. Статья “Оператор Лапласа” и “Ротор векторного поля”.
  3. Тиман А. Ф., Трофимов В. Н. Введение в теорию гармонических функций. М. Наука. 1968 г. 208с.

Ссылки[править | править код]

  • MathWorld description of Laplacian Архивная копия от 17 июня 2020 на Wayback Machine

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols nablacdotnabla, nabla ^{2} (where nabla is the nabla operator), or Delta . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).

The Laplace operator is named after the French mathematician Pierre-Simon de Laplace (1749–1827), who first applied the operator to the study of celestial mechanics: the Laplacian of the gravitational potential due to a given mass density distribution is a constant multiple of that density distribution. Solutions of Laplace’s equation Δf = 0 are called harmonic functions and represent the possible gravitational potentials in regions of vacuum.

The Laplacian occurs in many differential equations describing physical phenomena. Poisson’s equation describes electric and gravitational potentials; the diffusion equation describes heat and fluid flow, the wave equation describes wave propagation, and the Schrödinger equation describes the wave function in quantum mechanics. In image processing and computer vision, the Laplacian operator has been used for various tasks, such as blob and edge detection. The Laplacian is the simplest elliptic operator and is at the core of Hodge theory as well as the results of de Rham cohomology.

Definition[edit]

The Laplace operator is a second-order differential operator in the n-dimensional Euclidean space, defined as the divergence (nabla cdot) of the gradient (nabla f). Thus if f is a twice-differentiable real-valued function, then the Laplacian of f is the real-valued function defined by:

{displaystyle Delta f=nabla ^{2}f=nabla cdot nabla f}

(1)

where the latter notations derive from formally writing:

{displaystyle nabla =left({frac {partial }{partial x_{1}}},ldots ,{frac {partial }{partial x_{n}}}right).}

Explicitly, the Laplacian of f is thus the sum of all the unmixed second partial derivatives in the Cartesian coordinates xi:

{displaystyle Delta f=sum _{i=1}^{n}{frac {partial ^{2}f}{partial x_{i}^{2}}}}

(2)

As a second-order differential operator, the Laplace operator maps Ck functions to Ck−2 functions for k ≥ 2. It is a linear operator Δ : Ck(Rn) → Ck−2(Rn), or more generally, an operator Δ : Ck(Ω) → Ck−2(Ω) for any open set Ω ⊆ Rn.

Motivation[edit]

Diffusion[edit]

In the physical theory of diffusion, the Laplace operator arises naturally in the mathematical description of equilibrium.[1] Specifically, if u is the density at equilibrium of some quantity such as a chemical concentration, then the net flux of u through the boundary V of any smooth region V is zero, provided there is no source or sink within V:

{displaystyle int _{partial V}nabla ucdot mathbf {n} ,dS=0,}

where n is the outward unit normal to the boundary of V. By the divergence theorem,

{displaystyle int _{V}operatorname {div} nabla u,dV=int _{partial V}nabla ucdot mathbf {n} ,dS=0.}

Since this holds for all smooth regions V, one can show that it implies:

{displaystyle operatorname {div} nabla u=Delta u=0.}

The left-hand side of this equation is the Laplace operator, and the entire equation Δu = 0 is known as Laplace’s equation. Solutions of the Laplace equation, i.e. functions whose Laplacian is identically zero, thus represent possible equilibrium densities under diffusion.

The Laplace operator itself has a physical interpretation for non-equilibrium diffusion as the extent to which a point represents a source or sink of chemical concentration, in a sense made precise by the diffusion equation. This interpretation of the Laplacian is also explained by the following fact about averages.

Averages[edit]

Given a twice continuously differentiable function {displaystyle f:mathbb {R} ^{n}to mathbb {R} }, a point {displaystyle pin mathbb {R} ^{n}} and a real number {displaystyle h>0}, we let {displaystyle {overline {f}}_{B}(p,h)} be the average value of f over the ball with radius h centered at p, and {displaystyle {overline {f}}_{S}(p,h)} be the average value of f over the sphere (the boundary of a ball) with radius h centered at p. Then we have:[2]

{displaystyle {overline {f}}_{B}(p,h)=f(p)+{frac {Delta f(p)}{2(n+2)}}h^{2}+o(h^{2})quad {text{for}};;hto 0}

and

{displaystyle {overline {f}}_{S}(p,h)=f(p)+{frac {Delta f(p)}{2n}}h^{2}+o(h^{2})quad {text{for}};;hto 0.}

Density associated with a potential[edit]

If φ denotes the electrostatic potential associated to a charge distribution q, then the charge distribution itself is given by the negative of the Laplacian of φ:

{displaystyle q=-varepsilon _{0}Delta varphi ,}

where ε0 is the electric constant.

This is a consequence of Gauss’s law. Indeed, if V is any smooth region with boundary V, then by Gauss’s law the flux of the electrostatic field E across the boundary is proportional to the charge enclosed:

{displaystyle int _{partial V}mathbf {E} cdot mathbf {n} ,dS=int _{V}operatorname {div} mathbf {E} ,dV={frac {1}{varepsilon _{0}}}int _{V}q,dV.}

where the first equality is due to the divergence theorem. Since the electrostatic field is the (negative) gradient of the potential, this gives:

{displaystyle -int _{V}operatorname {div} (operatorname {grad} varphi ),dV={frac {1}{varepsilon _{0}}}int _{V}q,dV.}

Since this holds for all regions V, we must have

{displaystyle operatorname {div} (operatorname {grad} varphi )=-{frac {1}{varepsilon _{0}}}q}

The same approach implies that the negative of the Laplacian of the gravitational potential is the mass distribution. Often the charge (or mass) distribution are given, and the associated potential is unknown. Finding the potential function subject to suitable boundary conditions is equivalent to solving Poisson’s equation.

Energy minimization[edit]

Another motivation for the Laplacian appearing in physics is that solutions to Δf = 0 in a region U are functions that make the Dirichlet energy functional stationary:

{displaystyle E(f)={frac {1}{2}}int _{U}lVert nabla frVert ^{2},dx.}

To see this, suppose f : UR is a function, and u : UR is a function that vanishes on the boundary of U. Then:

{displaystyle left.{frac {d}{dvarepsilon }}right|_{varepsilon =0}E(f+varepsilon u)=int _{U}nabla fcdot nabla u,dx=-int _{U}u,Delta f,dx}

where the last equality follows using Green’s first identity. This calculation shows that if Δf = 0, then E is stationary around f. Conversely, if E is stationary around f, then Δf = 0 by the fundamental lemma of calculus of variations.

Coordinate expressions[edit]

Two dimensions[edit]

The Laplace operator in two dimensions is given by:

In Cartesian coordinates,

{displaystyle Delta f={frac {partial ^{2}f}{partial x^{2}}}+{frac {partial ^{2}f}{partial y^{2}}}}

where x and y are the standard Cartesian coordinates of the xy-plane.

In polar coordinates,

{displaystyle {begin{aligned}Delta f&={frac {1}{r}}{frac {partial }{partial r}}left(r{frac {partial f}{partial r}}right)+{frac {1}{r^{2}}}{frac {partial ^{2}f}{partial theta ^{2}}}\&={frac {partial ^{2}f}{partial r^{2}}}+{frac {1}{r}}{frac {partial f}{partial r}}+{frac {1}{r^{2}}}{frac {partial ^{2}f}{partial theta ^{2}}},end{aligned}}}

where r represents the radial distance and θ the angle.

Three dimensions[edit]

In three dimensions, it is common to work with the Laplacian in a variety of different coordinate systems.

In Cartesian coordinates,

{displaystyle Delta f={frac {partial ^{2}f}{partial x^{2}}}+{frac {partial ^{2}f}{partial y^{2}}}+{frac {partial ^{2}f}{partial z^{2}}}.}

In cylindrical coordinates,

{displaystyle Delta f={frac {1}{rho }}{frac {partial }{partial rho }}left(rho {frac {partial f}{partial rho }}right)+{frac {1}{rho ^{2}}}{frac {partial ^{2}f}{partial varphi ^{2}}}+{frac {partial ^{2}f}{partial z^{2}}},}

where rho represents the radial distance, φ the azimuth angle and z the height.

In spherical coordinates:

{displaystyle Delta f={frac {1}{r^{2}}}{frac {partial }{partial r}}left(r^{2}{frac {partial f}{partial r}}right)+{frac {1}{r^{2}sin theta }}{frac {partial }{partial theta }}left(sin theta {frac {partial f}{partial theta }}right)+{frac {1}{r^{2}sin ^{2}theta }}{frac {partial ^{2}f}{partial varphi ^{2}}},}

or

{displaystyle Delta f={frac {1}{r}}{frac {partial ^{2}}{partial r^{2}}}(rf)+{frac {1}{r^{2}sin theta }}{frac {partial }{partial theta }}left(sin theta {frac {partial f}{partial theta }}right)+{frac {1}{r^{2}sin ^{2}theta }}{frac {partial ^{2}f}{partial varphi ^{2}}},}

where φ represents the azimuthal angle and θ the zenith angle or co-latitude.

In general curvilinear coordinates (ξ1, ξ2, ξ3):

{displaystyle Delta =nabla xi ^{m}cdot nabla xi ^{n}{frac {partial ^{2}}{partial xi ^{m},partial xi ^{n}}}+nabla ^{2}xi ^{m}{frac {partial }{partial xi ^{m}}}=g^{mn}left({frac {partial ^{2}}{partial xi ^{m},partial xi ^{n}}}-Gamma _{mn}^{l}{frac {partial }{partial xi ^{l}}}right),}

where summation over the repeated indices is implied,
gmn is the inverse metric tensor and Γl mn are the Christoffel symbols for the selected coordinates.

N dimensions[edit]

In arbitrary curvilinear coordinates in N dimensions (ξ1, …, ξN), we can write the Laplacian in terms of the inverse metric tensor, {displaystyle g^{ij}}:

{displaystyle Delta ={frac {1}{sqrt {det g}}}{frac {partial }{partial xi ^{i}}}left({sqrt {det g}}g^{ij}{frac {partial }{partial xi ^{j}}}right),}

from the Voss-Weyl formula[3] for the divergence.

In spherical coordinates in N dimensions, with the parametrization x = RN with r representing a positive real radius and θ an element of the unit sphere SN−1,

{displaystyle Delta f={frac {partial ^{2}f}{partial r^{2}}}+{frac {N-1}{r}}{frac {partial f}{partial r}}+{frac {1}{r^{2}}}Delta _{S^{N-1}}f}

where ΔSN−1 is the Laplace–Beltrami operator on the (N − 1)-sphere, known as the spherical Laplacian. The two radial derivative terms can be equivalently rewritten as:

{displaystyle {frac {1}{r^{N-1}}}{frac {partial }{partial r}}left(r^{N-1}{frac {partial f}{partial r}}right).}

As a consequence, the spherical Laplacian of a function defined on SN−1RN can be computed as the ordinary Laplacian of the function extended to RN∖{0} so that it is constant along rays, i.e., homogeneous of degree zero.

Euclidean invariance[edit]

The Laplacian is invariant under all Euclidean transformations: rotations and translations. In two dimensions, for example, this means that:

{displaystyle Delta (f(xcos theta -ysin theta +a,xsin theta +ycos theta +b))=(Delta f)(xcos theta -ysin theta +a,xsin theta +ycos theta +b)}

for all θ, a, and b. In arbitrary dimensions,

{displaystyle Delta (fcirc rho )=(Delta f)circ rho }

whenever ρ is a rotation, and likewise:

{displaystyle Delta (fcirc tau )=(Delta f)circ tau }

whenever τ is a translation. (More generally, this remains true when ρ is an orthogonal transformation such as a reflection.)

In fact, the algebra of all scalar linear differential operators, with constant coefficients, that commute with all Euclidean transformations, is the polynomial algebra generated by the Laplace operator.

Spectral theory[edit]

The spectrum of the Laplace operator consists of all eigenvalues λ for which there is a corresponding eigenfunction f with:

{displaystyle -Delta f=lambda f.}

This is known as the Helmholtz equation.

If Ω is a bounded domain in Rn, then the eigenfunctions of the Laplacian are an orthonormal basis for the Hilbert space L2(Ω). This result essentially follows from the spectral theorem on compact self-adjoint operators, applied to the inverse of the Laplacian (which is compact, by the Poincaré inequality and the Rellich–Kondrachov theorem).[4] It can also be shown that the eigenfunctions are infinitely differentiable functions.[5] More generally, these results hold for the Laplace–Beltrami operator on any compact Riemannian manifold with boundary, or indeed for the Dirichlet eigenvalue problem of any elliptic operator with smooth coefficients on a bounded domain. When Ω is the n-sphere, the eigenfunctions of the Laplacian are the spherical harmonics.

Vector Laplacian[edit]

The vector Laplace operator, also denoted by nabla ^{2}, is a differential operator defined over a vector field.[6] The vector Laplacian is similar to the scalar Laplacian; whereas the scalar Laplacian applies to a scalar field and returns a scalar quantity, the vector Laplacian applies to a vector field, returning a vector quantity. When computed in orthonormal Cartesian coordinates, the returned vector field is equal to the vector field of the scalar Laplacian applied to each vector component.

The vector Laplacian of a vector field mathbf {A} is defined as

{displaystyle nabla ^{2}mathbf {A} =nabla (nabla cdot mathbf {A} )-nabla times (nabla times mathbf {A} ).}

In Cartesian coordinates, this reduces to the much simpler form as

{displaystyle nabla ^{2}mathbf {A} =(nabla ^{2}A_{x},nabla ^{2}A_{y},nabla ^{2}A_{z}),}

where A_x, A_y, and A_z are the components of the vector field mathbf {A} , and nabla ^{2} just on the left of each vector field component is the (scalar) Laplace operator. This can be seen to be a special case of Lagrange’s formula; see Vector triple product.

For expressions of the vector Laplacian in other coordinate systems see Del in cylindrical and spherical coordinates.

Generalization[edit]

The Laplacian of any tensor field mathbf {T} (“tensor” includes scalar and vector) is defined as the divergence of the gradient of the tensor:

{displaystyle nabla ^{2}mathbf {T} =(nabla cdot nabla )mathbf {T} .}

For the special case where mathbf {T} is a scalar (a tensor of degree zero), the Laplacian takes on the familiar form.

If mathbf {T} is a vector (a tensor of first degree), the gradient is a covariant derivative which results in a tensor of second degree, and the divergence of this is again a vector. The formula for the vector Laplacian above may be used to avoid tensor math and may be shown to be equivalent to the divergence of the Jacobian matrix shown below for the gradient of a vector:

{displaystyle nabla mathbf {T} =(nabla T_{x},nabla T_{y},nabla T_{z})={begin{bmatrix}T_{xx}&T_{xy}&T_{xz}\T_{yx}&T_{yy}&T_{yz}\T_{zx}&T_{zy}&T_{zz}end{bmatrix}},{text{ where }}T_{uv}equiv {frac {partial T_{u}}{partial v}}.}

And, in the same manner, a dot product, which evaluates to a vector, of a vector by the gradient of another vector (a tensor of 2nd degree) can be seen as a product of matrices:

{displaystyle mathbf {A} cdot nabla mathbf {B} ={begin{bmatrix}A_{x}&A_{y}&A_{z}end{bmatrix}}nabla mathbf {B} ={begin{bmatrix}mathbf {A} cdot nabla B_{x}&mathbf {A} cdot nabla B_{y}&mathbf {A} cdot nabla B_{z}end{bmatrix}}.}

This identity is a coordinate dependent result, and is not general.

Use in physics[edit]

An example of the usage of the vector Laplacian is the Navier-Stokes equations for a Newtonian incompressible flow:

{displaystyle rho left({frac {partial mathbf {v} }{partial t}}+(mathbf {v} cdot nabla )mathbf {v} right)=rho mathbf {f} -nabla p+mu left(nabla ^{2}mathbf {v} right),}

where the term with the vector Laplacian of the velocity field muleft(nabla ^2 mathbf{v}right) represents the viscous stresses in the fluid.

Another example is the wave equation for the electric field that can be derived from Maxwell’s equations in the absence of charges and currents:

{displaystyle nabla ^{2}mathbf {E} -mu _{0}epsilon _{0}{frac {partial ^{2}mathbf {E} }{partial t^{2}}}=0.}

This equation can also be written as:

{displaystyle Box ,mathbf {E} =0,}

where

{displaystyle Box equiv {frac {1}{c^{2}}}{frac {partial ^{2}}{partial t^{2}}}-nabla ^{2},}

is the D’Alembertian, used in the Klein–Gordon equation.

Generalizations[edit]

A version of the Laplacian can be defined wherever the Dirichlet energy functional makes sense, which is the theory of Dirichlet forms. For spaces with additional structure, one can give more explicit descriptions of the Laplacian, as follows.

Laplace–Beltrami operator[edit]

The Laplacian also can be generalized to an elliptic operator called the Laplace–Beltrami operator defined on a Riemannian manifold. The Laplace–Beltrami operator, when applied to a function, is the trace (tr) of the function’s Hessian:

{displaystyle Delta f=operatorname {tr} {big (}H(f){big )}}

where the trace is taken with respect to the inverse of the metric tensor. The Laplace–Beltrami operator also can be generalized to an operator (also called the Laplace–Beltrami operator) which operates on tensor fields, by a similar formula.

Another generalization of the Laplace operator that is available on pseudo-Riemannian manifolds uses the exterior derivative, in terms of which the “geometer’s Laplacian” is expressed as

{displaystyle Delta f=delta df.}

Here δ is the codifferential, which can also be expressed in terms of the Hodge star and the exterior derivative. This operator differs in sign from the “analyst’s Laplacian” defined above. More generally, the “Hodge” Laplacian is defined on differential forms α by

{displaystyle Delta alpha =delta dalpha +ddelta alpha .}

This is known as the Laplace–de Rham operator, which is related to the Laplace–Beltrami operator by the Weitzenböck identity.

D’Alembertian[edit]

The Laplacian can be generalized in certain ways to non-Euclidean spaces, where it may be elliptic, hyperbolic, or ultrahyperbolic.

In Minkowski space the Laplace–Beltrami operator becomes the D’Alembert operator Box or D’Alembertian:

{displaystyle square ={frac {1}{c^{2}}}{frac {partial ^{2}}{partial t^{2}}}-{frac {partial ^{2}}{partial x^{2}}}-{frac {partial ^{2}}{partial y^{2}}}-{frac {partial ^{2}}{partial z^{2}}}.}

It is the generalization of the Laplace operator in the sense that it is the differential operator which is invariant under the isometry group of the underlying space and it reduces to the Laplace operator if restricted to time-independent functions. The overall sign of the metric here is chosen such that the spatial parts of the operator admit a negative sign, which is the usual convention in high-energy particle physics. The D’Alembert operator is also known as the wave operator because it is the differential operator appearing in the wave equations, and it is also part of the Klein–Gordon equation, which reduces to the wave equation in the massless case.

The additional factor of c in the metric is needed in physics if space and time are measured in different units; a similar factor would be required if, for example, the x direction were measured in meters while the y direction were measured in centimeters. Indeed, theoretical physicists usually work in units such that c = 1 in order to simplify the equation.

The d’Alembert operator generalizes to a hyperbolic operator on pseudo-Riemannian manifolds.

See also[edit]

  • Laplace–Beltrami operator, generalization to submanifolds in Euclidean space and Riemannian and pseudo-Riemannian manifold.
  • The vector Laplacian operator, a generalization of the Laplacian to vector fields.
  • The Laplacian in differential geometry.
  • The discrete Laplace operator is a finite-difference analog of the continuous Laplacian, defined on graphs and grids.
  • The Laplacian is a common operator in image processing and computer vision (see the Laplacian of Gaussian, blob detector, and scale space).
  • The list of formulas in Riemannian geometry contains expressions for the Laplacian in terms of Christoffel symbols.
  • Weyl’s lemma (Laplace equation).
  • Earnshaw’s theorem which shows that stable static gravitational, electrostatic or magnetic suspension is impossible.
  • Del in cylindrical and spherical coordinates.
  • Other situations in which a Laplacian is defined are: analysis on fractals, time scale calculus and discrete exterior calculus.

Notes[edit]

  1. ^ Evans 1998, §2.2
  2. ^ Ovall, Jeffrey S. (2016-03-01). “The Laplacian and Mean and Extreme Values” (PDF). The American Mathematical Monthly. 123 (3): 287–291. doi:10.4169/amer.math.monthly.123.3.287. S2CID 124943537.
  3. ^ Archived at Ghostarchive and the Wayback Machine: Grinfeld, Pavel. “The Voss-Weyl Formula”. YouTube. Retrieved 9 January 2018.
  4. ^ Gilbarg & Trudinger 2001, Theorem 8.6
  5. ^ Gilbarg & Trudinger 2001, Corollary 8.11
  6. ^ MathWorld. “Vector Laplacian”.

References[edit]

  • Evans, L. (1998), Partial Differential Equations, American Mathematical Society, ISBN 978-0-8218-0772-9
  • The Feynman Lectures on Physics Vol. II Ch. 12: Electrostatic Analogs
  • Gilbarg, D.; Trudinger, N. (2001), Elliptic Partial Differential Equations of Second Order, Springer, ISBN 978-3-540-41160-4.
  • Schey, H. M. (1996), Div, Grad, Curl, and All That, W. W. Norton, ISBN 978-0-393-96997-9.

Further reading[edit]

  • The Laplacian – Richard Fitzpatrick 2006

External links[edit]

  • “Laplace operator”, Encyclopedia of Mathematics, EMS Press, 2001 [1994]
  • Weisstein, Eric W. “Laplacian”. MathWorld.
  • Laplacian in polar coordinates derivation
  • equations on the fractal cubes and Casimir effect

Опера́тор Лапла́са (Лапласиан) — дифференциальный оператор, действующий в линейном пространстве гладких функций и обозначаемый символом {displaystyle Delta }. Функции {displaystyle F } он ставит в соответствие функцию {displaystyle left({partial ^{2} over partial x_{1}^{2}}+{partial ^{2} over partial x_{2}^{2}}+ldots +{partial ^{2} over partial x_{n}^{2}}right)F}.

В произвольных ортогональных криволинейных координатах {displaystyle q_{1}, q_{2}, q_{3}}:

{displaystyle Delta f(q_{1}, q_{2}, q_{3})=operatorname {div} ,operatorname {grad} ,f(q_{1}, q_{2}, q_{3})={frac {1}{H_{1}H_{2}H_{3}}}left[{frac {partial }{partial q_{1}}}left({frac {H_{2}H_{3}}{H_{1}}}{frac {partial f}{partial q_{1}}}right)+{frac {partial }{partial q_{2}}}left({frac {H_{1}H_{3}}{H_{2}}}{frac {partial f}{partial q_{2}}}right)+{frac {partial }{partial q_{3}}}left({frac {H_{1}H_{2}}{H_{3}}}{frac {partial f}{partial q_{3}}}right)right]},

где {displaystyle H_{i} } – коэффициенты Ламе.

В цилиндрических координатах:

{displaystyle Delta f={1 over r}{partial  over partial r}left(r{partial f over partial r}right)+{partial ^{2}f over partial z^{2}}+{1 over r^{2}}{partial ^{2}f over partial phi ^{2}}}

В сферических координатах:


{displaystyle Delta f={1 over r^{2}}{partial  over partial r}left(r^{2}{partial f over partial r}right)+{1 over r^{2}sin theta }{partial  over partial theta }left(sin theta {partial f over partial theta }right)+{1 over r^{2}sin ^{2}theta }{partial ^{2}f over partial phi ^{2}}}

или

{displaystyle Delta f={1 over r}{partial ^{2} over partial r^{2}}left(rfright)+{1 over r^{2}sin theta }{partial  over partial theta }left(sin theta {partial f over partial theta }right)+{1 over r^{2}sin ^{2}theta }{partial ^{2}f over partial phi ^{2}}.}

Оператор Лапласа часто записывается следующим образом {displaystyle nabla ^{2}}, то есть скалярное произведение оператора набла на себя.

Применение

Через данный оператор удобно записывать уравнения Лапласа, Пуассона и волновое уравнение, хотя наиболее простой вид последнее принимает с использованием оператора Даламбера (Даламбертиана).

См. также

  • Оператор набла

cs:Laplaceův operátor
nl:Laplace-operator
pl:Operator Laplace’a
sk:Laplaceov operátor
sl:Laplaceov operator
sr:Лапласов оператор
sv:Laplaceoperatorn

“Del Squared” перенаправляется сюда. Для использования в других целях, см Del Squared (значения) .

В математике , то оператор Лапласа или лапласиан является дифференциальный оператор задается дивергенции от градиента о наличии скалярной функции на евклидовом пространстве . Обычно обозначается символами , (где – оператор набла ), или . В декартовой системе координат лапласиан задается суммой вторых частных производных функции по каждой независимой переменной . В других системах координат , таких как цилиндрические и сферические координаты , лапласиан также имеет полезную форму. Неформально, лапласиан Δ f ( p ) функции f в точке p измеряет, насколько среднее значение f по маленьким сферам или шарам с центром в p отклоняется от f ( p ) .
 набла  cdot  набла набла ^ {2} набла  Дельта

Оператор Лапласа назван в честь французского математика Пьера-Симона де Лапласа (1749–1827), который первым применил оператор к изучению небесной механики , где оператор дает постоянную, кратную плотности массы, когда он применяется к гравитационному исследованию. потенциал из-за распределения массы с данной плотностью. Решения уравнения Δ f = 0 , теперь называемого уравнением Лапласа , являются так называемыми гармоническими функциями и представляют возможные гравитационные поля в областях вакуума .

Лапласиан встречается в дифференциальных уравнениях, которые описывают многие физические явления, такие как электрические и гравитационные потенциалы , уравнение диффузии для потока тепла и жидкости , распространение волн и квантовую механику . Лапласиан представляет собой плотность потока от потока градиента функции. Например, чистая скорость, с которой химическое вещество, растворенное в жидкости, движется к некоторой точке или от нее, пропорциональна лапласиану химической концентрации в этой точке; выраженное символически, результирующее уравнение является уравнением диффузии. По этим причинам он широко используется в науке для моделирования различных физических явлений. Лапласиан – простейший эллиптический оператор, он лежит в основе теории Ходжа, а также результатов когомологий де Рама . В обработке изображений и компьютерном зрении оператор Лапласа использовался для различных задач, таких как обнаружение пятен и краев .

Определение

Оператор Лапласа – это дифференциальный оператор второго порядка в n- мерном евклидовом пространстве , определяемый как дивергенция ( ) градиента ( ). Таким образом, если – дважды дифференцируемая вещественнозначная функция , то лапласиан функции определяется следующим образом:
 набла  cdot  nabla fж ж

 Delta f =  nabla ^ 2 f =  nabla  cdot  nabla f

( 1 )

где последние обозначения происходят от формального написания:

{ displaystyle  nabla f =  left ({ frac { partial f} { partial x_ {1}}},  ldots, { frac { partial f} { partial x_ {n}}}  right ).}

Эквивалентно, лапласиан f является суммой всех несмешанных вторых частных производных в декартовых координатах x i :

 Delta f =  sum_ {i = 1} ^ n  frac { partial ^ 2 f} { partial x ^ 2_i}

( 2 )

Как дифференциальный оператор второго порядка, оператор Лапласа отображает функции C k в функции C k −2 при k ≥ 2 . Выражение ( 1 ) (или, что эквивалентно ( 2 )) определяет оператор Δ: C k ( R n ) → C k −2 ( R n ) или, в более общем смысле, оператор Δ: C k (Ω) → C k – 2 (Ω) для любого открытого множества Ω .

Мотивация

Диффузия

В физической теории диффузии оператор Лапласа (через уравнение Лапласа ) естественным образом возникает при математическом описании равновесия . В частности, если у есть плотность в равновесии некоторой величины , такие как химические концентрации, то результирующий поток из U через границу любой гладкой области V равен нуль, при условии , что нет источника или раковин в пределах V :

 int _ { partial V}  nabla и  cdot  mathbf {n} , dS = 0,

где п есть внешнее единичный вектор нормали к границе V . По теореме о дивергенции ,

 int_V  operatorname {div}  nabla u , dV =  int _ { partial V}  nabla u  cdot  mathbf {n} , dS = 0.

Поскольку это верно для всех гладких областей V , можно показать, что это влечет:

 operatorname {div}  nabla u =  Delta u = 0.

Левая часть этого уравнения – оператор Лапласа. Сам оператор Лапласа имеет физическую интерпретацию неравновесной диффузии как степень, в которой точка представляет источник или сток химической концентрации, в смысле, уточненном уравнением диффузии .

Средние

Учитывая дважды непрерывно дифференцируемую функцию , точку и действительное число , мы позволяем быть средним значением по шару с радиусом с центром в , и быть средним значением по сфере с радиусом с центром в . Тогда у нас есть:
{ displaystyle f:  mathbb {R} ^ {n}  to  mathbb {R}}{ displaystyle p  in  mathbb {R} ^ {n}}{ displaystyle h> 0}{ displaystyle { overline {f}} _ {B} (p, h)}жчасп{ displaystyle { overline {f}} _ {S} (p, h)}жчасп

{ displaystyle { overline {f}} _ {B} (p, h) = f (p) + { frac { Delta f (p)} {2 (n + 2)}} h ^ {2} + o (h ^ {2})  quad { mbox {for}} ; ; h  to 0}

а также

{ displaystyle { overline {f}} _ {S} (p, h) = f (p) + { frac { Delta f (p)} {2n}} h ^ {2} + o (h ^ {2})  quad { mbox {for}} ; ; h  to 0.}

Плотность, связанная с потенциалом

Если φ обозначает электростатический потенциал, связанный с распределением заряда q , то само распределение заряда задается отрицанием лапласиана φ :

{ displaystyle q = -  varepsilon _ {0}  Delta  varphi,}

где ε 0 – электрическая постоянная .

Это следствие закона Гаусса . Действительно, если V – любая гладкая область, то по закону Гаусса поток электростатического поля E пропорционален приложенному заряду:

{ displaystyle  int _ { partial V}  mathbf {E}  cdot  mathbf {n} , dS =  int _ {V}  operatorname {div}  mathbf {E} , dV = { frac {1} { varepsilon _ {0}}}  int _ {V} q , dV.}

где первое равенство следует из теоремы о расходимости . Поскольку электростатическое поле представляет собой (отрицательный) градиент потенциала, теперь это дает:

{ displaystyle -  int _ {V}  operatorname {div} ( operatorname {grad}  varphi) , dV = { frac {1} { varepsilon _ {0}}}  int _ {V} q , dV.}

Итак, поскольку это верно для всех областей V , мы должны иметь

{ displaystyle  operatorname {div} ( operatorname {grad}  varphi) = - { frac {1} { varepsilon _ {0}}} q}

Тот же подход подразумевает, что отрицательным элементом лапласиана гравитационного потенциала является распределение масс . Часто приводится распределение заряда (или массы), а связанный с ним потенциал неизвестен. Нахождение потенциальной функции с подходящими граничными условиями эквивалентно решению уравнения Пуассона .

Минимизация энергии

Другая причина появления лапласиана в физике состоит в том, что решения для Δ f = 0 в области U являются функциями, которые делают функционал энергии Дирихле стационарным :

{ displaystyle E (f) = { frac {1} {2}}  int _ {U}  lVert  nabla f  rVert ^ {2} , dx.}

Чтобы убедиться в этом, предположим , что F  : UR является функцией, и у  : UR является функцией , которая обращается в нуль на границе U . Потом:

{ displaystyle  left. { frac {d} {d  varepsilon}}  right | _ { varepsilon = 0} E (f +  varepsilon u) =  int _ {U}  nabla f  cdot  nabla u , dx = -  int _ {U} u ,  Delta f , dx}

где последнее равенство следует из первого тождества Грина . Этот расчет показывает, что если Δ f = 0 , то E неподвижен вокруг f . Наоборот, если E стационарно вокруг f , то Δ f = 0 по основной лемме вариационного исчисления .

Координатные выражения

Два измерения

Оператор Лапласа в двух измерениях определяется выражением:

В декартовой системе координат ,

{ displaystyle  Delta f = { frac { partial ^ {2} f} { partial x ^ {2}}} + { frac { partial ^ {2} f} { partial y ^ {2} }}}

где х и у являются стандартными декартовы координаты по ху плоскости.

В полярных координатах ,

{ displaystyle { begin {align}  Delta f & = { frac {1} {r}} { frac { partial} { partial r}}  left (r { frac { partial f} { partial r}}  right) + { frac {1} {r ^ {2}}} { frac { partial ^ {2} f} { partial  theta ^ {2}}} \ & = {  frac { partial ^ {2} f} { partial r ^ {2}}} + { frac {1} {r}} { frac { partial f} { partial r}} + { frac {1} {r ^ {2}}} { frac { partial ^ {2} f} { partial  theta ^ {2}}},  end {align}}}

где r представляет собой радиальное расстояние, а θ – угол.

Три измерения

В трех измерениях обычно работают с лапласианом в различных системах координат.

В декартовой системе координат ,

{ displaystyle  Delta f = { frac { partial ^ {2} f} { partial x ^ {2}}} + { frac { partial ^ {2} f} { partial y ^ {2} }} + { frac { partial ^ {2} f} { partial z ^ {2}}}.}

В цилиндрических координатах ,

{ displaystyle  Delta f = { frac {1} { rho}} { frac { partial} { partial  rho}}  left ( rho { frac { partial f} { partial  rho }}  right) + { frac {1} { rho ^ {2}}} { frac { partial ^ {2} f} { partial  varphi ^ {2}}} + { frac { частичный ^ {2} f} { partial z ^ {2}}},}

где представляет собой радиальное расстояние, φ – азимутальный угол, а z – высоту.
 rho

В сферических координатах :

{ displaystyle { begin {align}  Delta f & = { frac {1} {r ^ {2}}} { frac { partial} { partial r}}  left (r ^ {2} { frac { partial f} { partial r}}  right) + { frac {1} {r ^ {2}  sin  theta}} { frac { partial} { partial  theta}}  left ( sin  theta { frac { partial f} { partial  theta}}  right) + { frac {1} {r ^ {2}  sin ^ {2}  theta}} { frac {  partial ^ {2} f} { partial  varphi ^ {2}}},  end {выровнено}}}

где φ представляет собой азимутальный угол, а θ – зенитный угол или со-широту .

В общих криволинейных координатах ( ξ 1 , ξ 2 , ξ 3 ):

{ displaystyle  Delta =  nabla  xi ^ {m}  cdot  nabla  xi ^ {n} { frac { partial ^ {2}} { partial  xi ^ {m}  partial  xi ^ { n}}} +  nabla ^ {2}  xi ^ {m} { frac { partial} { partial  xi ^ {m}}} = g ^ {mn}  left ({ frac { partial ^ {2}} { partial  xi ^ {m}  partial  xi ^ {n}}} -  Gamma _ {mn} ^ {l} { frac { partial} { partial  xi ^ {l }}}Правильно),}

где подразумевается суммирование по повторяющимся индексам ,
g mn – обратный метрический тензор, а Γ l mn – символы Кристоффеля для выбранных координат.

Размеры N

В произвольных координатах криволинейных в N измерениях ( £ , 1 , …, ξ N ), то можно записать лапласиан в терминах обратного метрического тензора , :
{ displaystyle g ^ {ij}}

{ displaystyle  Delta = { frac {1} { sqrt { det g}}} { frac { partial} { partial  xi ^ {i}}}  left ({ sqrt { det g }} g ^ {ij} { frac { partial} { partial  xi ^ {j}}}  right),}

из формулы Фосса – Вейля для расходимости .

В сферических координатах в N измерениях , с параметризацией x = R N, где r представляет положительный действительный радиус, а θ – элемент единичной сферы S N −1 ,

{ displaystyle  Delta f = { frac { partial ^ {2} f} { partial r ^ {2}}} + { frac {N-1} {r}} { frac { partial f} { partial r}} + { frac {1} {r ^ {2}}}  Delta _ {S ^ {N-1}} f}

где Δ S N −1 – оператор Лапласа – Бельтрами на ( N – 1) -сфере, известный как сферический лапласиан. Два члена с радиальной производной могут быть эквивалентно переписаны как:

{ displaystyle { frac {1} {r ^ {N-1}}} { frac { partial} { partial r}}  left (r ^ {N-1} { frac { partial f}) { partial r}}  right).}

Как следствие, сферический лапласиан функции , определенный на S N -1R N может быть вычислен как обычный лапласиан функции распространяется на R N ∖ {0} , так что он является постоянным вдоль лучей, то есть, однородные степеней нуль.

Евклидова инвариантность

Лапласиан инвариантен относительно всех евклидовых преобразований : вращений и сдвигов . Например, в двух измерениях это означает, что:

{ Displaystyle  Дельта (е (Икс  соз  тета -у  грех  тета + а, х  грех  тета + у  соз  тета + Ь)) = ( Дельта е) (х  соз  тета - y  sin  theta + a, x  sin  theta + y  cos  theta + b)}

для всех θ , a и b . В произвольных размерах,

{ Displaystyle  Delta (е  circ  rho) = ( Delta f)  circ  rho}

всякий раз, когда ρ – вращение, и аналогично:

{ Displaystyle  Delta (е  circ  tau) = ( Delta f)  circ  tau}

всякий раз, когда τ – перевод. (В более общем смысле это остается верным, когда ρ – ортогональное преобразование, такое как отражение .)

Фактически, алгебра всех скалярных линейных дифференциальных операторов с постоянными коэффициентами, коммутирующих со всеми евклидовыми преобразованиями, является алгеброй полиномов, порожденной оператором Лапласа.

Спектральная теория

Спектр оператора Лапласа состоит из всех собственных значений А , для которых есть соответствующий собственная функция F с:

-  Delta f =  lambda f.

Это известно как уравнение Гельмгольца .

Если Ω является ограниченная область в R п , то собственные функции лапласиана являются ортонормированный базис для гильбертова пространства L 2 (Q) . Этот результат по существу следует из спектральной теоремы о компактных самосопряженных операторах , примененной к обратному лапласиану (компактному в силу неравенства Пуанкаре и теоремы Реллиха – Кондрахова ). Также можно показать, что собственные функции являются бесконечно дифференцируемыми функциями. В более общем смысле, эти результаты верны для оператора Лапласа – Бельтрами на любом компактном римановом многообразии с краем или даже для проблемы собственных значений Дирихле любого эллиптического оператора с гладкими коэффициентами в ограниченной области. Когда Ω является n- сферой , собственными функциями лапласиана являются сферические гармоники .

Векторный лапласиан

Оператор Лапласа вектор , также обозначается , является дифференциальным оператором , определенным над векторным полем . Векторный лапласиан аналогичен скалярному лапласиану; в то время как скалярный лапласиан применяется к скалярному полю и возвращает скалярную величину, векторный лапласиан применяется к векторному полю , возвращая векторную величину. При вычислении в ортонормированных декартовых координатах возвращаемое векторное поле равно векторному полю скалярного лапласиана, примененного к каждому компоненту вектора.
 набла ^ {2}

Вектор лапласиана из векторного поля определяется как
 mathbf {A}

 nabla ^ {2} { mathbf {A}} =  nabla ( nabla  cdot { mathbf {A}}) -  nabla  times ( nabla  times { mathbf {A}}).

В декартовых координатах это сводится к гораздо более простой форме, как

 nabla ^ {2} { mathbf {A}} = ( nabla ^ {2} A_ {x},  nabla ^ {2} A_ {y},  nabla ^ {2} A_ {z}),

где , и – компоненты векторного поля , а слева от каждой компоненты векторного поля находится (скалярный) оператор Лапласа. Это можно рассматривать как частный случай формулы Лагранжа; см. тройное произведение вектора .
A_xA_yA_z mathbf {A}  набла ^ {2}

Для выражения векторного лапласиана в других системах координат см. Del в цилиндрических и сферических координатах .

Обобщение

Лапласиан любого тензорного поля ( «тензор» включает в себя скалярное и векторное) определяется как дивергенции от градиента тензора:
 mathbf {T}

{ displaystyle  nabla ^ {2}  mathbf {T} = ( nabla  cdot  nabla)  mathbf {T}.}

Для частного случая, когда – скаляр (тензор нулевой степени), лапласиан принимает знакомый вид.
 mathbf {T}

Если – вектор (тензор первой степени), градиент представляет собой ковариантную производную, которая приводит к тензору второй степени, и дивергенция этого снова является вектором. Формула для векторного лапласиана, приведенная выше, может использоваться, чтобы избежать тензорной математики, и может быть показана как эквивалентная дивергенции матрицы Якоби, показанной ниже для градиента вектора:
 mathbf {T}

{ displaystyle  nabla  mathbf {T} = ( nabla T_ {x},  nabla T_ {y},  nabla T_ {z}) = { begin {bmatrix} T_ {xx} & T_ {xy} & T_ { xz} \ T_ {yx} & T_ {yy} & T_ {yz} \ T_ {zx} & T_ {zy} & T_ {zz}  end {bmatrix}}, { text {where}} T_ {uv}  Equiv { frac { partial T_ {u}} { partial v}}.}

И таким же образом скалярное произведение, которое оценивается как вектор, вектора градиентом другого вектора (тензор 2-й степени) можно рассматривать как произведение матриц:

{ displaystyle  mathbf {A}  cdot  nabla  mathbf {B} = { begin {bmatrix} A_ {x} & A_ {y} & A_ {z}  end {bmatrix}}  nabla  mathbf {B} = { begin {bmatrix}  mathbf {A}  cdot  nabla B_ {x} &  mathbf {A}  cdot  nabla B_ {y} &  mathbf {A}  cdot  nabla B_ {z}  end { bmatrix}}.}

Эта идентичность является результатом, зависящим от координат, и не является общим.

Использование в физике

Примером использования векторного лапласиана являются уравнения Навье-Стокса для ньютоновского потока несжимаемой жидкости :

 rho  left ({ frac { partial { mathbf {v}}} { partial t}} + ({ mathbf {v}}  cdot  nabla) { mathbf {v}}  right) =  rho { mathbf {f}} -  nabla p +  mu  left ( nabla ^ {2} { mathbf {v}}  right),

где член с вектором лапласианом скорости поля представляет собой вязкие напряжения в жидкости.
 mu  left ( nabla ^ 2  mathbf {v}  right)

Другой пример – волновое уравнение для электрического поля, которое может быть получено из уравнений Максвелла в отсутствие зарядов и токов:

 nabla ^ 2  mathbf {E} -  mu_0  epsilon_0  frac { partial ^ 2  mathbf {E}} { partial t ^ 2} = 0.

Предыдущее уравнение также можно записать как:

 Box ,  mathbf {E} = 0,

куда

 Box  Equiv  frac {1} {c ^ 2}  frac { partial ^ 2} { partial t ^ 2} -  nabla ^ 2,

– даламбертиан , используемый в уравнении Клейна – Гордона .

Обобщения

Версия лапласиана может быть определена везде, где имеет смысл функционал энергии Дирихле , а именно теория форм Дирихле . Для пространств с дополнительной структурой можно дать более явное описание лапласиана следующим образом.

Оператор Лапласа – Бельтрами

Лапласиан также может быть обобщен до эллиптического оператора, называемого оператором Лапласа – Бельтрами, определенным на римановом многообразии . Оператор Даламбера обобщается до гиперболического оператора на псевдоримановых многообразиях . Оператор Лапласа – Бельтрами, примененный к функции, представляет собой след ( tr ) гессиана функции :

{ Displaystyle  Delta f =  OperatorName {tr} { big (} H (f) { big)}}

где след берется относительно обратного к метрическому тензору . Оператор Лапласа – Бельтрами также может быть обобщен на оператор (также называемый оператором Лапласа – Бельтрами), который работает с тензорными полями , по аналогичной формуле.

Другое обобщение оператора Лапласа, доступное на псевдоримановых многообразиях, использует внешнюю производную , в терминах которой «лапласиан геометра» выражается как

{ displaystyle  Delta f =  delta df.}

Здесь δ – кодифференциал , который также может быть выражен через звезду Ходжа и внешнюю производную. Этот оператор отличается по знаку от «лапласиана аналитика», определенного выше. В более общем смысле лапласиан «Ходжа» определяется на дифференциальных формах α формулой

{ displaystyle  Delta  alpha =  delta d  alpha + d  delta  alpha.}

Это известно как оператор Лапласа – де Рама , который связан с оператором Лапласа – Бельтрами тождеством Вайтценбека .

Даламбертиан

Лапласиан может быть определенным образом обобщен на неевклидовы пространства, где он может быть эллиптическим , гиперболическим или ультрагиперболическим .

В пространстве Минковского оператор Лапласа-Бельтрами становится оператором Даламбера или даламбертиана:
Коробка

{ displaystyle  square = { frac {1} {c ^ {2}}} { frac { partial ^ {2}} { partial t ^ {2}}} - { frac { partial ^ { 2}} { partial x ^ {2}}} - { frac { partial ^ {2}} { partial y ^ {2}}} - { frac { partial ^ {2}} { partial z ^ {2}}}.}

Это обобщение оператора Лапласа в том смысле, что это дифференциальный оператор, инвариантный относительно группы изометрий основного пространства, и он сводится к оператору Лапласа, если ограничиваться функциями, не зависящими от времени. Общий знак метрики здесь выбран таким, чтобы пространственные части оператора допускали отрицательный знак, что является обычным соглашением в физике частиц высоких энергий . Оператор Даламбера также известен как волновой оператор, потому что это дифференциальный оператор, фигурирующий в волновых уравнениях , а также часть уравнения Клейна – Гордона , которое сводится к волновому уравнению в безмассовом случае.

Дополнительный коэффициент c в метрике необходим в физике, если пространство и время измеряются в разных единицах; аналогичный коэффициент потребовался бы, если, например, направление x измерялось в метрах, а направление y – в сантиметрах. Действительно, физики-теоретики обычно работают с такими единицами, что c = 1 , чтобы упростить уравнение.

Смотрите также

  • Оператор Лапласа – Бельтрами , обобщение на подмногообразия в евклидовом пространстве, риманово и псевдориманово многообразие.
  • Вектор Лапласа оператор, обобщение лапласиана в векторных полей .
  • Лапласиан в дифференциальной геометрии .
  • Дискретный оператор Лапласа является конечно-разностным аналогом непрерывного лапласиана, определенный на графиках и сетках.
  • Лапласиан – распространенный оператор в обработке изображений и компьютерном зрении (см. Лапласиан гаусса , детектор капель и масштабное пространство ).
  • Список формул в римановой геометрии содержит выражение для лапласиана в терминах символов Кристоффеля.
  • Лемма Вейля (уравнение Лапласа) .
  • Теорема Ирншоу, которая показывает, что стабильная статическая гравитационная, электростатическая или магнитная подвеска невозможна.
  • Del в цилиндрических и сферических координатах .
  • Другими ситуациями, в которых определяется лапласиан, являются: анализ фракталов , исчисление шкалы времени и дискретное внешнее исчисление .

Примечания

использованная литература

  • Эванс, Л. (1998), Уравнения в частных производных , Американское математическое общество, ISBN 978-0-8218-0772-9
  • Feynman, R .; Leighton, R; Сэндс, М. (1970), «Глава 12: Электростатические аналоги», Лекции Фейнмана по физике , 2 , Аддисон-Уэсли-Лонгман
  • Gilbarg, D .; Трудингер Н. (2001), Эллиптические дифференциальные уравнения с частными производными второго порядка , Springer, ISBN 978-3-540-41160-4.
  • Schey, HM (1996), Div, Grad, Curl, и все такое , WW Norton, ISBN 978-0-393-96997-9.

дальнейшее чтение

  • http://farside.ph.utexas.edu/teaching/em/lectures/node23.html

внешние ссылки

  • «Оператор Лапласа» , Математическая энциклопедия , EMS Press , 2001 [1994]
  • Вайсштейн, Эрик В. «Лапласиан» . MathWorld .
  • Вывод лапласиана в полярных координатах

Опера́тор Лапла́са (лапласиа́н) — дифференциальный оператор, действующий в линейном пространстве гладких функций и обозначаемый символом  Delta. Функции F он ставит в соответствие функцию left({partial^2 over partial x_1^2} + {partial^2 over partial x_2^2} + ldots  + {partial^2 over partial x_n^2}right)F.

Оператор Лапласа часто обозначается следующим образом nabla^2, то есть в виде скалярного произведения оператора Набла на себя. Оператор Лапласа эквивалентен также последовательному взятию операций градиента и дивергенции: Delta=nablacdotnabla=operatorname{div},operatorname{grad}, таким образом значение оператора Лапласа в точке может быть истолковано как плотность источников (стоков) потенциального векторного поля  operatorname{grad}F в этой точке.

Содержание

  • 1 Другое определение оператора Лапласа
  • 2 Выражения для оператора Лапласа в различных криволинейных системах координат
    • 2.1 Цилиндрические координаты
    • 2.2 Сферические координаты
    • 2.3 Параболические координаты
    • 2.4 Цилиндрические параболические координаты
  • 3 Применение
  • 4 Вариации и обобщения
  • 5 См. также
  • 6 Литература
  • 7 Внешние ссылки

Другое определение оператора Лапласа

Оператор Лапласа является естественным обобщением на функции нескольких переменных обычной второй производной функции одного переменного. В самом деле, если функция  f (x) имеет в окрестности точки  x_0 непрерывную вторую производную  f''(x), то, как это следует из формулы Тейлора

 f(x_0+r)=f(x_0)+rf'(x_0)+frac{r^2}{2}f''(x_0)+o(r^2), при rto 0,,
 f(x_0-r)=f(x_0)-rf'(x_0)+frac{r^2}{2}f''(x_0)+o(r^2), при rto 0,

вторая производная есть предел

 f''(x_0)=limlimits_{r to 0} frac{2}{r^2} left{ frac{f(x_0+r)+f(x_0-r)}{2}-f(x_0) right}.

Eсли, переходя к функции  F от  k переменных, поступить таким же образом, т.е. для заданной точки  M_0(x_1^0,x_2^0, ... ,x_k^0) рассматривать её  k -мерную шаровую окрестность  Q_r радиуса  r и разность между средним арифметическим

 frac{1}{sigma(S_r)}intlimits_{S_r}Fdsigma

функции  F на границе  S_r такой окрестности с площадью границы  sigma(S_r) и значением  F(M_0) в центре этой окрестности  M_0, то в случае непрерывности вторых частных производных функции  F в окрестности точки  M_0 значение лапласиана  Delta F в этой точке есть предел

 Delta F(M_0)=limlimits_{r to 0} frac{2k}{r^2} left{frac{1}{sigma(S_r)}intlimits_{S_r}F(M)dsigma -F(M_0) right}.

Одновременно с предыдущим представлением для оператора Лапласа функции  F, имеющей непрерывные вторые производные, справедлива формула

 Delta F(M_0)=limlimits_{r to 0} frac{2(k+2)}{r^2} left{frac{1}{omega(Q_r)}intlimits_{Q_r}F(M)domega -F(M_0) right}, где  omega(Q_r) – объём окресности  Q_r.

Эта формула выражает непосредственную связь лапласиана функции с её объёмным средним в окрестности данной точки.

Доказательство этих формул можно найти, например, в [1].

Вышеизложенные пределы, во всех случаях, когда они существуют, могут служить определением оператора Лапласа функции  F. Такое определение предпочтительнее обычного определения лапласиана, предполагающего существование вторых производных рассматриваемых функций, и совпадает с обычным определением в случае непрерывности этих производных.

Выражения для оператора Лапласа в различных криволинейных системах координат

В произвольных ортогональных криволинейных координатах в трехмерном пространстве q_1, q_2, q_3:

Delta f (q_1, q_2, q_3) = operatorname{div},operatorname{grad},f(q_1, q_2, q_3) =
=frac{1}{H_1H_2H_3}left[ frac{partial}{partial q_1}left( frac{H_2H_3}{H_1}frac{partial f}{partial q_1} right) + frac{partial}{partial q_2}left( frac{H_1H_3}{H_2}frac{partial f}{partial q_2} right) +  frac{partial}{partial q_3}left( frac{H_1H_2}{H_3}frac{partial f}{partial q_3} right)right],
где H_i — коэффициенты Ламе.

Цилиндрические координаты

В цилиндрических координатах вне прямой  r=0:

 Delta f 
= {1 over r} {partial over partial r}
  left( r {partial f over partial r} right) 
+ {partial^2f over partial z^2}
+ {1 over r^2} {partial^2 f over partial varphi^2}

Сферические координаты

В сферических координатах вне начала отсчёта:

 Delta f 
= {1 over r^2} {partial over partial r}
  left( r^2 {partial f over partial r} right) 
+ {1 over r^2 sin theta} {partial over partial theta}
  left( sin theta {partial f over partial theta} right) 
+ {1 over r^2 sin^2 theta} {partial^2 f over partial varphi^2}

или

 Delta f 
= {1 over r} {partial^2 over partial r^2}
  left( rf right) 
+ {1 over r^2 sin theta} {partial over partial theta}
  left( sin theta {partial f over partial theta} right) 
+ {1 over r^2 sin^2 theta} {partial^2 f over partial varphi^2}.

В случае если  f=f(r)

 Delta f =  {d^2 fover dr^2} + {n-1 over r } {dfover dr}.

Параболические координаты

В параболических координатах (в трёхмерном пространстве) вне начала отсчёта:


Delta f= frac{1}{sigma^{2} + tau^{2}} 
left[
frac{1}{sigma} frac{partial }{partial sigma} 
left( sigma frac{partial f}{partial sigma} right) +
frac{1}{tau} frac{partial }{partial tau} 
left( tau frac{partial f}{partial tau} right)right] +
frac{1}{sigma^2tau^2}frac{partial^2 f}{partial varphi^2}

Цилиндрические параболические координаты

В координатах параболического цилиндра вне начала отсчёта:

Delta F(u,v,z) = frac{1}{c^2(u^2+v^2)} left[ frac{partial^2 F }{partial u^2}+ frac{partial^2 F }{partial v^2}right] + frac{partial^2 F }{partial z^2}

Применение

С помощью данного оператора удобно записывать уравнения Лапласа, Пуассона и волновое уравнение, хотя наиболее простой вид последнее принимает с использованием оператора Д’Aламбера (Даламбертиана). Впрочем, последний представляет собой не более, чем оператор Лапласа в пространстве Минковского (формально пространство Минковского можно ввести для любого поля, подчиняющегося волновому уравнению, хотя, конечно, параметр c может быть в каждом конкретном случае своим, например, скорость звука).

В физике оператор Лапласа применим в электростатике и электродинамике, во многих уравнениях физики сплошных сред, а также при изучении равновесия мембран, пленок или поверхностей раздела фаз с поверхностным натяжением (см. Лапласово давление), стационарных задач диффузии и теплопроводности, которые сводятся в непрерывном пределе к обычным уравнениям Лапласа или Пуассона или к некоторым их обобщениям.

Вариации и обобщения

  • Оператор Д’Аламбера

См. также

  • Оператор набла
  • Уравнение Лапласа
  • Гармоническая функция

Литература

  1. Тиман А.Ф., Трофимов В.Н. Введение в теорию гармонических функций. М. Наука. 1968г. 208с.

Внешние ссылки

  • MathWorld description of Laplacian

Wikimedia Foundation.
2010.

Добавить комментарий