当前位置:文档之家› 数学物理方程第七章_变分法及其应用

数学物理方程第七章_变分法及其应用

第 7 章 变分法及其应用
在数学物理中能够精确求解的边值问题或固有值问题并不多。因此,在实际工作中,经 常采用各种近似的方法来解决具体的问题,变分法就是其中最有力的方法之一。所谓变分法 就是求函数极值的方法,下面先介绍什么是泛函以及泛函极值,然后再简要介绍求泛函极值 的方法以及它的一些应用。
7.1 泛函和泛函极值
∫ J[ y(x)] = x1 F (x, y, y′)dx x0
(7.1.3)
式中,被积函数 F (x, y, y′) 称为核。
在实际工作中,为了完成某项任务,我们首先要分析实际问题特殊现象与一般规律 之间的关系,然后建立数学上的表达式。如求连接两个定点的曲线段中弧长最短的曲线
方程 y(x) ,即求自变量 y = y* ,使泛函
1
=C
2gy 1+ y'2
y(1 +
y′2) =
1 2gC 2
= 2r
引入变量代换 x = x(θ ) ,并设
则由式(7.1.14)可得
y′ = cot θ 2
(7.1.12) (7.1.13) (7.1.14)
上式对θ 求导,得

所以
y = 2r sin 2 θ = r(1− cosθ ) 2
y′ dx = r sinθ dθ
J[ y(x) + εϕ(x)] ε =0
=
0
由于
J[
y(x)
+
εϕ ( x)]
=
J[ y(x)]
+
∫b ∂F [ a ∂y
εϕ ( x)
+
∂F ∂y′
εϕ ′(x)]dx
则有
∫b ∂F [
εϕ ( x)
+
∂F
εϕ ′(x)]dx
=
0
a ∂y
∂y ′
(7.1.10)
以 ε 乘式(7.1.10),且
δy(x) = εϕ(x)
r(1

cosθ
)
7.2 变分法在固有值问题中的应用
本节我们介绍利用变分法解决固有值问题。在第 2 章里我们学过,所谓固有值问题就是在一 定的边界条件下,求解含有参数的微分方程。为了表达上的方便,将需要求解的常微分方程 写为
L[ y(x)] = λy(x)
(7.2.1)
其中 L 是一个作用在函数 y(x) 上的线性微分算符,假定所讨论的固有值问题是劝函数
δy(x) = εϕ(x)
式中, ε 是任意小的实数;ϕ(x) 是充分光滑的任意函数,并且满足条件
ϕ(a) = 0, ϕ(b) = 0
这样,函数
y(x) + εϕ(x)
满足边界条件式(7.1.5)。因此,泛函
J[ y(x) + εϕ(x)]
当 ε = 0 时取最小值 J[ y(x)] ,从而有
d dε
y(x0 ) = y0 , y(x1) = y1
的一切可微函数 y(x) 的集合,这里的每一个元素对应着 xy 平面上由点 P0 (x0 , y0 ) 到点
P1(x1, y1) 的一条光滑曲线 y = y(x) 。用 L 表示曲线上 P0 P1 段的弧长,则
∫ L = x1 1+ y′2 dx x0
大(小)值、取极值的必要条件是 dy dx
x= x0
= 0 。下面我们仿照函数微分的概念来定义泛
函的变分概念,进而导出泛函极值存在的必要条件。
设 y, y0 是集合 C 的元素,称 δy = y − y0 为函数 y 在 y0 处的变分。
这里的 δy 是 x 的函数,它与 ∆y 的区别在于:变分 δy 反映的是整个函数的改变,
式(7.1.11)称为欧拉-拉格朗日方程,简记为 E-L 方程,这就是泛函 J[ y(x)] 有极限
的必要条件,也就是说, y = y(x) 使泛函式(7.1.6)取极小值,则 y = y(x) 一定使欧
拉-拉格朗日方程式(7.1.11)满足边界条件式(7.1.5)的解。
我们把满足 E-L 方程边值问题的解称为驻留函数,对应的积分曲线称为驻留曲线。 严格地讲,E-L 方程边值问题的解满足变分问题的必要条件,因此它是否是极值函数, 还需作进一步的判别。在实际问题中,极值的存在性通常给出问题时己经肯定了,这样,

y(0) = 0, y( p) = q
的所有连续函数 y(x) 中,求出一个函数 y* 使泛函式(7.1.4)取最小值。
对泛函求极值的问题称为变分问题,使泛函取极值的函数称为变分问题的解,也称 为极值函数。
在微分学中,求函数 y = y(x) 的 极值是求自变量 x 的值,当 x 取这些值时, y 取极
而 ∆y 表示的是同一个函数 y(x) 因 x 的不同值而产生的差异。在本书,我们总是假定
y(x) 和 F (x, y, y′) 都是充分光滑的,且 y(x) 在两个端点处固定,即
y(a) = y1, y(b) = y2
式中, y1, y2 是两个常数。
下面我们考虑泛函
∫ J[ y(x)] = b F (x, y, y′)dx a
L[ϕn (x)] = λnϕn (x)
这样,将线性微分算符 L ,同时作用在式(7.2.5)两侧,则
∑ ∑ L[ y(x)] = L[ Cnϕn (x)] = L[Cnϕn (x)]
n
n
以 y(x) 乘上式,同时作积分并记为 J[ y(x)] ,即
(7.2.7) (7.2.8) (7.2.9)
J[ y(x)] = ∫ y(x)L[ y(x)]dx
m
∫ ∑ ∑ J[ y(x)] = Cmϕm (x) ⋅ L[ Cnϕn (x)] n
∑∑ ∫ =
CnCm ϕn (x)L[ϕn (x)]dx
mn
∑∑ ∫ =
CnCmλn ϕm (x)ϕn (x)dx
我们以前研究的函数是指这样一种现象,对于数集 A 中的任一个元素 z ,数集 B 中存在
一个元素 w 与之对应,我们就说 w 是 z 的一个函数,记为 w(z) 。在自然现象中,不仅存在
这样的数与数的对应,还存在着其他种种性质不同的对应关系。我们看下面的问题。
设 C 为区间[x0 , x1] 上满足条件
y2 x2
− −
y1 x1
(x − x1) +
y1
∫ J[ y(x)] = p 1+ y′2 dx 0 2gy

y(0) = 0, y( p) = q
这样
F (x, y, y') = F( y, y') = 1+ y′2 2gy
其 E-L 方程为
由于 所以有
∂F ∂y

d dx
⎜⎜⎝⎛
∂F ∂y′
⎟⎟⎠⎞
当函数 y(x) 有微小改变且变为 y(x) + δy(x) 时,利用
(7.1.5) (7.1.6)
上式可推出
F
(
x,
y
+
δy,
y′
+
δy′)
=
F
( x,
y,
y′)
+
∂F ∂y
δy
+
∂F ∂y′
δy′
∫ J ( y + ∆y) − J ( y) =
b
[
∂F
δy
+
∂F
δy′]dx
a ∂y ∂y′
上式称为 J ( y) 的变分,记为δJ ( y) ,即
J[ y(x)] 是 y(x) 的泛函,这样我们就引出一系列重要结论。
(7.2.10)
引理 1 泛函式(7.2.10)的极小值等于相应的固有值问题的最小固有值 λ1 ,而使
J[ y(x)] 取这一极小值的极值函数就是相应于固有值 λ1 所对应的固有函数ϕ1(x) 。
证明 将展开式(7.2.5)代入泛函式(7.2.10),则
当一个实际现象已知其有唯一的极值存在,而这时也只得到一个驻留函数,则可以判定
这个驻留函数就是极值函数。
下面我们来解决本章开始部分的两个例题。
例 1 最短距离问题

∫ J[ y(x)] = x1 1+ y′2 dx x0
因为 F = 1− y′2 ,所以
E-L 方程为 则有
∂F = 0, ∂F = y′
∫ δJ ( y) =
b
[
∂F
δy
+
∂F
δy′]dx
a ∂y ∂y′
(7.1.7)
下面我们证明,泛函 J ( y) 取极值的必要条件是
δJ ( y) = 0
或者
∂F ∂y

d dx
⎜⎜⎝⎛
∂F ∂y′
⎟⎟⎠⎞
=
0
设 y = y(x) 使泛函 J ( y) 取极值,取函数 y(x) 变分的特殊形式为
(7.1.8) (7.1.9)
=
0
d [F ( y, y′) − y′ ∂F ]
dx
∂y′
=
∂F ∂y
y′ +
∂F ∂y′
y′′

y′′
∂F ∂y′

y′
d dx
⎜⎜⎝⎛ −
∂F ∂y′
⎟⎟⎠⎞
=
0
F(
y,
y′)

∂F ∂y′
=
C
将(7.1.2)代入式(7.1.13)
1+ y′2 − y'
y'
=C
2gy
2gy 1+ y'2
由此得
素 J 与之对应,称 J 是 y(x) 的泛函数,记作
J = J[ y(x)]
(7.1.2)
相关主题