偏微分方程(Partial Differential Equation I)
偏微分方程(Partial Differential Equation II)
偏微分方程(Partial Differential Equation III)
偏微分方程(Partial Differential Equation IV)
参考文献:
《数学物理方程》| 季孝达
《数学物理方法》| 吴崇试
《数学物理方法》| 梁昆淼
MOOC北京大学《数学物理方法》| 吴崇试 、高春媛
变分法初步
泛函的概念
泛函: 泛函是将任意给定的函数映射为一个数,简单的说就是以整个函数为自变量的函数。这个概念可以看成是函数概念的推广。
对于一元函数,以函数集合
M
=
{
y
(
x
)
∣
a
⩽
x
⩽
b
}
M=\{y(x)|a⩽x⩽b\}
M={y(x)∣a⩽x⩽b} 为定义域,几何上表示某一平面曲线的集合。定义域内的函数
y
(
x
)
y(x)
y(x) 与数值存在映射,称为泛函
J
[
y
]
J[y]
J[y] 。
示例 1:(极小曲线问题)设在
O
x
y
Oxy
Oxy 平面上有一簇曲线
y
(
x
)
y(x)
y(x) ,其长度为泛函
J
[
y
]
=
∫
a
b
1
+
y
′
2
d
x
J[y]=\int_a^b\sqrt{1+y'^2}\mathrm dx
J[y]=∫ab1+y′2dx
显然,
J
[
y
]
J[y]
J[y] 的数值依赖于整个函数
y
(
x
)
y(x)
y(x) 的改变而改变。对于函数,给定一个
x
x
x 值,有一个函数值与之对应,对于泛函,则必须给出某一区间上的函数
y
(
x
)
y(x)
y(x) ,才能得到一个泛函值
J
[
y
]
J[y]
J[y] 。
(定义在同一区间上的)函数不同,泛函值当然不同。为了强调泛函值
J
[
y
]
J[y]
J[y] 与函数
y
(
x
)
y(x)
y(x) 之间的依赖关系,常常又把函数
y
(
x
)
y(x)
y(x) 称为自变函数。
泛函的形式可以是多种多样的,本课程中只限于用积分形式定义的泛函。
对于自变函数为一元函数
y
(
x
)
y(x)
y(x) ,则泛函为
J
[
y
]
=
∫
a
b
L
(
x
,
y
,
y
′
)
d
x
(1.1)
J[y]=\int_a^bL(x,y,y')\mathrm dx\tag{1.1}
J[y]=∫abL(x,y,y′)dx(1.1)
其中
L
L
L 是已知函数,具有连续的二阶偏导数。
如果自变函数是二元函数
u
(
x
,
y
)
u(x,y)
u(x,y) ,则泛函为
J
[
u
]
=
∬
S
L
(
x
,
y
,
u
,
u
x
,
u
y
)
d
x
d
y
(1.2)
J[u]=\iint\limits_SL(x,y,u,u_x,u_y)\mathrm dx\mathrm dy\tag{1.2}
J[u]=S∬L(x,y,u,ux,uy)dxdy(1.2)
对于更多个自变量的多元函数,也有类似的定义。
示例 2:(最速下降问题)如图,在重力作用下,一质点从
(
x
0
,
y
0
)
(x_0,y_0)
(x0,y0) 点沿平面曲线
y
(
x
)
y(x)
y(x) 无摩擦自由下滑到
(
x
1
,
y
1
)
(x_1,y_1)
(x1,y1) ,则所需的时间为
y
(
x
)
y(x)
y(x) 的泛函
J
[
y
]
=
∫
x
0
x
1
d
s
2
g
(
y
0
−
y
)
=
∫
x
0
x
1
1
+
y
′
2
2
g
(
y
0
−
y
)
d
x
\begin{aligned} J[y]= & \int_{x_0}^{x_1}\cfrac{\mathrm ds}{\sqrt{2g(y_0-y)}} \\ = & \int_{x_0}^{x_1}\cfrac{\sqrt{1+y'^2}}{\sqrt{2g(y_0-y)}}\mathrm dx \end{aligned}
J[y]==∫x0x12g(y0−y)ds∫x0x12g(y0−y)1+y′2dx
这里,要求
y
(
x
)
y(x)
y(x) 一定通过端点
(
x
0
,
y
0
)
(x_0,y_0)
(x0,y0) 和
(
x
1
,
y
1
)
(x_1,y_1)
(x1,y1)
示例 3:(极小曲面问题)设在空间上的光滑曲面簇
u
(
x
,
y
)
u(x,y)
u(x,y) ,其曲面面积定义了泛函
J
[
y
]
=
∬
S
1
+
u
x
2
+
u
y
2
d
x
d
y
J[y]=\iint\limits_S\sqrt{1+u_x^2+u_y^2}\mathrm dx\mathrm dy
J[y]=S∬1+ux2+uy2dxdy
泛函的极值
变分法基本引理:设
f
(
x
)
f(x)
f(x) 在
[
a
,
b
]
[a,b]
[a,b] 上连续,若对于任意满足边界条件
h
(
a
)
=
h
(
b
)
=
0
h(a)=h(b)=0
h(a)=h(b)=0 的函数
h
(
x
)
h(x)
h(x) 均有
∫
a
b
f
(
x
)
h
(
x
)
d
x
=
0
\int_a^bf(x)h(x)\mathrm dx=0
∫abf(x)h(x)dx=0
则必有
f
(
x
)
=
0
f(x)=0
f(x)=0 。由于
h
(
x
)
h(x)
h(x) 的随意性,此引理可用反证法证明。对于多元函数重积分的情况也有类似的引理。
变分的概念:对于函数
f
(
x
)
f(x)
f(x) 假设自变量
x
x
x 不变,改变函数的形式得到一个与原函数稍有差别的新函数
f
ˉ
(
x
)
=
f
(
x
)
+
δ
g
(
x
)
\bar f(x)=f(x)+δg(x)
fˉ(x)=f(x)+δg(x) 。其中,
g
(
x
)
g(x)
g(x) 是任意连续函数,
δ
δ
δ 是微小系数,即对于
∀
ϵ
>
0
,
∣
δ
y
(
x
)
∣
<
ϵ
\forallϵ>0,|δy(x)|<ϵ
∀ϵ>0,∣δy(x)∣<ϵ 。
对于函数的任意自变量,函数
f
(
x
)
f(x)
f(x) 由于形式上的微小改变而得到的改变量称为该函数的变分。自变函数的变分其实是自变量微分的推广。如下图,函数
y
(
x
)
y(x)
y(x) 的变分
δ
y
(
x
)
δy(x)
δy(x)
以一元函数为例,总结变分的几条运算性质
(1) 由于变分独立于函数自变量,所以变分与微分可以交换
δ
f
′
=
(
δ
f
)
′
δf'=(δf)'
δf′=(δf)′
(2) 线性运算
δ
(
α
f
+
β
g
)
=
α
⋅
δ
f
+
β
⋅
δ
g
δ(\alpha f+\beta g)=\alpha\cdot δf+\beta\cdot δg
δ(αf+βg)=α⋅δf+β⋅δg
(3) 乘积的变分运算
δ
(
f
⋅
g
)
=
(
δ
f
)
⋅
g
+
f
⋅
(
δ
g
)
δ(f\cdot g)=(δf)\cdot g+f\cdot (δg)
δ(f⋅g)=(δf)⋅g+f⋅(δg)
(4) 积分的变分运算(只需把定积分看成级数和即可证明)
δ
∫
a
b
f
d
x
=
∫
a
b
(
δ
f
)
d
x
δ\int_a^bf\mathrm dx=\int_a^b(δf)\mathrm dx
δ∫abfdx=∫ab(δf)dx
(5) 复合函数的变分运算,其法则和微分运算完全相同,例如(注意变分和自变量无关)
δ
F
(
x
,
y
,
y
′
)
=
∂
F
∂
y
δ
y
+
∂
F
∂
y
′
δ
y
′
δF(x,y,y')=\cfrac{∂F}{∂y}δy+\cfrac{∂F}{∂y'}δy'
δF(x,y,y′)=∂y∂Fδy+∂y′∂Fδy′
一元函数泛函的极值:类似于函数的极值问题,若对自变函数为
y
0
(
x
)
y_0(x)
y0(x) 及其附近的自变函数
y
0
(
x
)
+
δ
y
(
x
)
y_0(x)+δy(x)
y0(x)+δy(x) ,恒有
J
[
y
]
⩽
J
[
y
0
+
δ
y
]
J[y]⩽J[y_0+δy]
J[y]⩽J[y0+δy] ,则称泛函
J
[
y
]
J[y]
J[y] 在自变函数为
y
0
(
x
)
y_0(x)
y0(x) 时取得极小值。类似的,可以定义泛函的极大值。极大值和极小值统称为极值。
可以仿照函数的方法,导出泛函取极值的必要条件。对于一元函数
y
(
x
)
y(x)
y(x) 的泛函
J
[
y
]
=
∫
a
b
L
(
x
,
y
,
y
′
)
d
x
J[y]=\int_a^bL(x,y,y')\mathrm dx
J[y]=∫abL(x,y,y′)dx
泛函的差值为
J
[
y
+
δ
y
]
−
J
[
y
]
=
∫
a
b
[
L
(
x
,
y
+
δ
y
,
y
′
+
δ
y
′
)
−
L
(
x
,
y
,
y
′
)
]
d
x
J[y+δy]-J[y]=\int_a^b[L(x,y+δy,y'+δy')-L(x,y,y')]\mathrm dx
J[y+δy]−J[y]=∫ab[L(x,y+δy,y′+δy′)−L(x,y,y′)]dx
考虑到函数的变分
δ
y
(
x
)
δy(x)
δy(x) 足够小,可以将上式被积函数在极值函数附近作泰勒展开
J
[
y
+
δ
y
]
−
J
[
y
]
=
∫
a
b
[
(
δ
y
∂
∂
y
+
∂
∂
y
′
δ
y
′
)
L
+
1
2
!
(
δ
y
∂
∂
y
+
∂
∂
y
′
δ
y
′
)
2
L
+
⋯
]
d
x
=
∫
a
b
(
δ
y
∂
L
∂
y
+
∂
L
∂
y
′
δ
y
′
)
d
x
+
1
2
!
∫
a
b
(
δ
y
∂
∂
y
+
∂
∂
y
′
δ
y
′
)
2
L
d
x
+
⋯
\begin{aligned} J[y+δy]-J[y]= & \int_a^b[(δy\cfrac{∂}{∂y}+\cfrac{∂}{∂y'}δy')L +\cfrac{1}{2!}(δy\cfrac{∂}{∂y}+\cfrac{∂}{∂y'}δy')^2L+\cdots]\mathrm dx \\ = & \int_a^b(δy\cfrac{∂L}{∂y}+\cfrac{∂L}{∂y'}δy')\mathrm dx + \cfrac{1}{2!}\int_a^b(δy\cfrac{∂}{∂y}+\cfrac{∂}{∂y'}δy')^2L\mathrm dx+\cdots \end{aligned}
J[y+δy]−J[y]==∫ab[(δy∂y∂+∂y′∂δy′)L+2!1(δy∂y∂+∂y′∂δy′)2L+⋯]dx∫ab(δy∂y∂L+∂y′∂Lδy′)dx+2!1∫ab(δy∂y∂+∂y′∂δy′)2Ldx+⋯
其中定义
δ
J
[
y
]
=
∫
a
b
(
δ
y
∂
L
∂
y
+
δ
y
′
∂
L
∂
y
′
)
d
x
(2.1)
δJ[y]=\int_a^b(δy\cfrac{∂L}{∂y}+δy'\cfrac{∂L}{∂y'})\mathrm dx\tag{2.1}
δJ[y]=∫ab(δy∂y∂L+δy′∂y′∂L)dx(2.1)
是泛函
J
[
y
]
J[y]
J[y] 的一级变分。
δ
2
J
[
y
]
=
∫
a
b
(
δ
y
∂
∂
y
+
∂
∂
y
′
δ
y
′
)
2
L
d
x
(2.2)
δ^2J[y]=\int_a^b(δy\cfrac{∂}{∂y}+\cfrac{∂}{∂y'}δy')^2L\mathrm dx\tag{2.2}
δ2J[y]=∫ab(δy∂y∂+∂y′∂δy′)2Ldx(2.2)
是泛函
J
[
y
]
J[y]
J[y] 的二级变分,依次可以定义三级及以上高级变分。
函数
f
(
x
)
f(x)
f(x) 取极值的必要条件为
f
′
(
x
)
=
0
f'(x)=0
f′(x)=0 ,微分形式为
d
f
(
x
)
=
f
′
(
x
)
d
x
=
0
\mathrm df(x)=f'(x)\mathrm dx=0
df(x)=f′(x)dx=0 。
注意到
δ
y
(
x
)
δy(x)
δy(x) 和
δ
y
′
(
x
)
δy'(x)
δy′(x) 永远是微量,因此舍弃掉二次项及以上高次项。和函数取极值类似,泛函取极值的必要条件为泛函的一级变分为零
δ
J
[
y
]
=
∫
a
b
(
δ
y
∂
L
∂
y
+
δ
y
′
∂
L
∂
y
′
)
d
x
=
0
(2.3)
δJ[y]=\int_a^b(δy\cfrac{∂L}{∂y}+δy'\cfrac{∂L}{∂y'})\mathrm dx=0\tag{2.3}
δJ[y]=∫ab(δy∂y∂L+δy′∂y′∂L)dx=0(2.3)
将第二项分部积分
δ
J
[
y
]
=
∂
L
∂
y
′
δ
y
∣
a
b
+
∫
a
b
(
∂
L
∂
y
−
d
d
x
∂
L
∂
y
′
)
δ
y
d
x
=
0
δJ[y]=\cfrac{∂L}{∂y'}δy\Big|_a^b+\int_a^b(\cfrac{∂L}{∂y} -\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂L}{∂y'})δy\mathrm dx=0
δJ[y]=∂y′∂Lδy∣∣∣ab+∫ab(∂y∂L−dxd∂y′∂L)δydx=0
根据变分法基本引理,
y
(
x
)
y(x)
y(x) 为极值时必须满足
∂
L
∂
y
−
d
d
x
∂
L
∂
y
′
=
0
(2.4)
\cfrac{∂L}{∂y}-\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂L}{∂y'}=0\tag{2.4}
∂y∂L−dxd∂y′∂L=0(2.4)
和边界条件
∂
L
∂
y
′
∣
b
δ
y
(
b
)
−
∂
L
∂
y
′
∣
a
δ
y
(
a
)
=
0
(2.5)
\cfrac{∂L}{∂y'}\Big|_bδy(b)-\cfrac{∂L}{∂y'}\Big|_aδy(a)=0\tag{2.5}
∂y′∂L∣∣∣bδy(b)−∂y′∂L∣∣∣aδy(a)=0(2.5)
方程 (2.4) 称为 Euler-Lagrange 方程。
如果
L
=
L
(
x
,
y
′
)
L=L(x,y')
L=L(x,y′) ,则 Euler-Lagrange 方程化为
d
d
x
∂
L
∂
y
′
=
0
\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂L}{∂y'}=0
dxd∂y′∂L=0 ,它的首次积分
F
y
′
=
C
1
F_{y'}=C_1
Fy′=C1
如果
L
=
L
(
y
,
y
′
)
L=L(y,y')
L=L(y,y′) ,则 Euler-Lagrange 方程的首次积分
F
−
y
′
F
y
′
=
C
1
F-y'F_{y'}=C_1
F−y′Fy′=C1
如果
L
=
L
(
x
,
y
)
L=L(x,y)
L=L(x,y) , Euler-Lagrange 方程是一个隐函数方程
F
y
=
0
F_y=0
Fy=0
值得指出的是,E-L方程只是泛函有极值的必要条件,并不是充分条件。就是说,当泛函有极值时,E-L方程成立,解可能不止一个,它们只是极值函数的候选者。
(1) 如果边界值
y
(
a
)
=
α
,
y
(
b
)
=
β
y(a)=\alpha,y(b)=\beta
y(a)=α,y(b)=β,即两个端点固定的情形,此时
δ
y
(
a
)
=
δ
y
(
b
)
=
0
δy(a)=δy(b)=0
δy(a)=δy(b)=0
一定满足边界条件 (2.5)
(2) 如果边界值
y
(
a
)
,
y
(
b
)
y(a),y(b)
y(a),y(b) 可以任意取值,此时
δ
y
(
a
)
,
δ
y
(
b
)
δy(a),δy(b)
δy(a),δy(b) 也可以任意取值,所以极值满足
∂
L
∂
y
′
∣
a
=
∂
L
∂
y
′
∣
b
=
0
\cfrac{∂L}{∂y'}\Big|_a=\cfrac{∂L}{∂y'}\Big|_b=0
∂y′∂L∣∣∣a=∂y′∂L∣∣∣b=0
此边界条件称为自然边界条件。
(3) 如果一端固定,如
y
(
a
)
=
α
,
y
(
b
)
y(a)=\alpha,y(b)
y(a)=α,y(b) 自由滑动,所以极值满足
y
(
a
)
=
α
,
∂
L
∂
y
′
∣
b
=
0
y(a)=\alpha,\cfrac{∂L}{∂y'}\Big|_b=0
y(a)=α,∂y′∂L∣∣∣b=0
二元函数泛函的极值: 对于二元函数
u
(
x
,
y
)
u(x,y)
u(x,y) 的泛函
J
[
u
]
=
∬
S
L
(
x
,
y
,
u
,
u
x
,
u
y
)
d
x
d
y
J[u]=\iint\limits_SL(x,y,u,u_x,u_y)\mathrm dx\mathrm dy
J[u]=S∬L(x,y,u,ux,uy)dxdy
二元函数的泛函取极值的必要条件依然为泛函的一级变分为零。
δ
J
[
u
]
=
∬
S
δ
L
(
x
,
y
,
u
,
u
x
,
u
y
)
d
x
d
y
=
∬
S
[
∂
L
∂
u
δ
u
+
∂
L
∂
u
x
δ
u
x
+
∂
L
∂
u
y
δ
u
y
]
d
x
d
y
=
∬
S
[
∂
L
∂
u
δ
u
+
∂
L
∂
u
x
(
δ
u
)
x
+
∂
L
∂
u
y
(
δ
u
)
y
]
d
x
d
y
=
∬
S
[
∂
L
∂
u
−
∂
∂
x
(
∂
L
∂
u
x
)
−
∂
∂
y
(
∂
L
∂
u
y
)
]
δ
u
d
x
d
y
+
∬
S
[
∂
∂
x
(
∂
L
∂
u
x
δ
u
)
+
∂
∂
y
(
∂
L
∂
u
y
δ
u
)
]
d
x
d
y
=
0
\begin{aligned} δJ[u]= & \iint\limits_SδL(x,y,u,u_x,u_y)\mathrm dx\mathrm dy \\ =& \iint\limits_S[\cfrac{∂L}{∂u}δu+\cfrac{∂L}{∂u_x}δu_x+\cfrac{∂L}{∂u_y}δu_y] \mathrm dx\mathrm dy \\ =& \iint\limits_S[\cfrac{∂L}{∂u}δu+\cfrac{∂L}{∂u_x}(δu)_x+\cfrac{∂L}{∂u_y}(δu)_y] \mathrm dx\mathrm dy \\ =& \iint\limits_S\left[\cfrac{∂L}{∂u} -\cfrac{∂}{∂x}\left(\cfrac{∂L}{∂u_x}\right) -\cfrac{∂}{∂y}\left(\cfrac{∂L}{∂u_y}\right)\right] δu\mathrm dx\mathrm dy \\ +& \iint\limits_S\left[\cfrac{∂}{∂x}\left(\cfrac{∂L}{∂u_x}δu\right) +\cfrac{∂}{∂y}\left(\cfrac{∂L}{∂u_y}δu\right)\right]\mathrm dx\mathrm dy \\ = &0 \end{aligned}
δJ[u]====+=S∬δL(x,y,u,ux,uy)dxdyS∬[∂u∂Lδu+∂ux∂Lδux+∂uy∂Lδuy]dxdyS∬[∂u∂Lδu+∂ux∂L(δu)x+∂uy∂L(δu)y]dxdyS∬[∂u∂L−∂x∂(∂ux∂L)−∂y∂(∂uy∂L)]δudxdyS∬[∂x∂(∂ux∂Lδu)+∂y∂(∂uy∂Lδu)]dxdy0
上式中使用格林公式
∬
S
(
∂
Q
∂
x
−
∂
P
∂
y
)
d
x
d
y
=
∫
∂
S
P
d
x
+
Q
d
y
\iint\limits_S(\cfrac{∂Q}{∂x}-\cfrac{∂P}{∂y})\mathrm dx\mathrm dy= \int_{∂S} P\mathrm dx+Q\mathrm dy
S∬(∂x∂Q−∂y∂P)dxdy=∫∂SPdx+Qdy
可得到
∬
S
[
∂
∂
x
(
∂
L
∂
u
x
δ
u
)
+
∂
∂
y
(
∂
L
∂
u
y
δ
u
)
]
d
x
d
y
=
∫
∂
S
[
−
∂
L
∂
u
y
d
x
+
∂
L
∂
u
x
d
y
]
δ
u
\iint\limits_S\left[\cfrac{∂}{∂x}\left(\cfrac{∂L}{∂u_x}δu\right) +\cfrac{∂}{∂y}\left(\cfrac{∂L}{∂u_y}δu\right)\right]\mathrm dx\mathrm dy =\int_{∂S} \left[-\cfrac{∂L}{∂u_y}\mathrm dx+\cfrac{∂L}{∂u_x}\mathrm dy\right]δu
S∬[∂x∂(∂ux∂Lδu)+∂y∂(∂uy∂Lδu)]dxdy=∫∂S[−∂uy∂Ldx+∂ux∂Ldy]δu
因此,二元函数泛函取极值的必要条件(积分形式)为
δ
J
[
u
]
=
∬
S
[
∂
L
∂
u
−
∂
∂
x
(
∂
L
∂
u
x
)
−
∂
∂
y
(
∂
L
∂
u
y
)
]
δ
u
d
x
d
y
+
∫
∂
S
[
−
∂
L
∂
u
y
d
x
+
∂
L
∂
u
x
d
y
]
δ
u
=
0
δJ[u]=\iint\limits_S\left[\cfrac{∂L}{∂u} -\cfrac{∂}{∂x}\left(\cfrac{∂L}{∂u_x}\right) -\cfrac{∂}{∂y}\left(\cfrac{∂L}{∂u_y}\right)\right] δu\mathrm dx\mathrm dy \\ +\int_{∂S} \left[-\cfrac{∂L}{∂u_y}\mathrm dx+\cfrac{∂L}{∂u_x}\mathrm dy\right]δu=0
δJ[u]=S∬[∂u∂L−∂x∂(∂ux∂L)−∂y∂(∂uy∂L)]δudxdy+∫∂S[−∂uy∂Ldx+∂ux∂Ldy]δu=0
根据变分法基本引理,极值满足 Euler-Lagrange 方程(微分形式)
∂
L
∂
u
−
∂
∂
x
(
∂
L
∂
u
x
)
−
∂
∂
y
(
∂
L
∂
u
y
)
=
0
(2.6)
\cfrac{∂L}{∂u} -\cfrac{∂}{∂x}\left(\cfrac{∂L}{∂u_x}\right) -\cfrac{∂}{∂y}\left(\cfrac{∂L}{∂u_y}\right)=0\tag{2.6}
∂u∂L−∂x∂(∂ux∂L)−∂y∂(∂uy∂L)=0(2.6)
和边界条件
∫
∂
S
[
−
∂
L
∂
u
y
d
x
+
∂
L
∂
u
x
d
y
]
δ
u
=
0
(2.7)
\int_{∂S} \left[-\cfrac{∂L}{∂u_y}\mathrm dx+\cfrac{∂L}{∂u_x}\mathrm dy\right]δu=0\tag{2.7}
∫∂S[−∂uy∂Ldx+∂ux∂Ldy]δu=0(2.7)
(1) 若函数满足第一类边界条件
u
∣
∂
S
=
ϕ
(
x
,
y
)
u\Big|_{∂S}=\phi(x,y)
u∣∣∣∂S=ϕ(x,y),即边界数值固定,因此
δ
u
(
x
,
y
)
∣
∂
S
=
0
δu(x,y)\Big|_{∂S}=0
δu(x,y)∣∣∣∂S=0
则一定满足边界条件 (2.7)
(2) 如果允许函数在边界面上自由取值,则极值必须满足边界条件
[
∂
L
∂
u
x
d
y
−
∂
L
∂
u
y
d
x
]
∣
∂
S
=
0
\left[\cfrac{∂L}{∂u_x}\mathrm dy-\cfrac{∂L}{∂u_y}\mathrm dx\right]\Big|_{∂S}=0
[∂ux∂Ldy−∂uy∂Ldx]∣∣∣∂S=0
类似的可以讨论三元及以上函数泛函和泛函极值的情况。
泛函的条件极值
条件极值:和函数求条件极值类似,泛函
J
[
y
,
z
]
=
∫
a
b
F
(
x
,
y
,
z
,
y
′
,
z
′
)
d
x
J[y,z]=\int_a^bF(x,y,z,y',z')\mathrm dx
J[y,z]=∫abF(x,y,z,y′,z′)dx
在约束条件
G
(
x
,
y
,
z
)
=
0
G(x,y,z)=0
G(x,y,z)=0
的条件极值也可用Lagrange 乘数法求解。
设约束条件
G
=
0
G=0
G=0 确定的微分方程为
z
=
ϕ
(
x
,
y
)
z=\phi(x,y)
z=ϕ(x,y),带入泛函可得
J
[
y
]
=
∫
a
b
F
1
(
x
,
y
,
y
′
)
d
x
J[y]=\int_a^bF_1(x,y,y')\mathrm dx
J[y]=∫abF1(x,y,y′)dx
其中
F
1
(
x
,
y
,
y
′
)
=
F
(
x
,
y
,
ϕ
,
y
′
,
ϕ
x
+
ϕ
y
y
′
)
F_1(x,y,y')=F(x,y,\phi,y',\phi_x+\phi_yy')
F1(x,y,y′)=F(x,y,ϕ,y′,ϕx+ϕyy′) ,实际上用消元法把条件极值问题转化为无条件极值,所以极值满足的方程为
∂
F
1
∂
y
−
d
d
x
∂
F
1
∂
y
′
=
0
\cfrac{∂F_1}{∂y}-\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂F_1}{∂y'}=0
∂y∂F1−dxd∂y′∂F1=0
于是
(
∂
F
∂
y
−
d
d
x
∂
F
∂
y
′
)
+
ϕ
z
(
∂
F
∂
z
−
d
d
x
∂
F
∂
z
′
)
=
0
(\cfrac{∂F}{∂y}-\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂F}{∂y'}) +\phi_z(\cfrac{∂F}{∂z}-\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂F}{∂z'})=0
(∂y∂F−dxd∂y′∂F)+ϕz(∂z∂F−dxd∂z′∂F)=0
又因为
G
(
x
,
y
,
z
)
=
0
,
ϕ
z
=
−
G
y
G
z
G(x,y,z)=0,\phi_z=-\frac{G_y}{G_z}
G(x,y,z)=0,ϕz=−GzGy 所以
1
G
y
(
∂
F
∂
y
−
d
d
x
∂
F
∂
y
′
)
=
1
G
z
(
∂
F
∂
z
−
d
d
x
∂
F
∂
z
′
)
\cfrac{1}{G_y}(\cfrac{∂F}{∂y}-\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂F}{∂y'}) =\cfrac{1}{G_z}(\cfrac{∂F}{∂z}-\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂F}{∂z'})
Gy1(∂y∂F−dxd∂y′∂F)=Gz1(∂z∂F−dxd∂z′∂F)
等式两端都是关于
x
x
x 的函数,记为
−
λ
(
x
)
-λ(x)
−λ(x) ,则
{
F
y
−
d
d
x
F
y
′
+
λ
(
x
)
G
y
=
0
F
z
−
d
d
x
F
z
′
+
λ
(
x
)
G
z
=
0
\begin{cases} F_y-\cfrac{\mathrm d}{\mathrm dx}F_{y'}+λ(x)G_y=0 \\ F_z-\cfrac{\mathrm d}{\mathrm dx}F_{z'}+λ(x)G_z=0 \end{cases}
⎩⎪⎨⎪⎧Fy−dxdFy′+λ(x)Gy=0Fz−dxdFz′+λ(x)Gz=0
故得出了泛函条件极值的拉格朗日乘数法。
等周问题:在约束条件
J
1
[
y
]
=
∫
a
b
G
(
x
,
y
,
y
′
)
d
x
=
L
J_1[y]=\int_a^bG(x,y,y')\mathrm dx=L
J1[y]=∫abG(x,y,y′)dx=L
下求泛函
J
[
y
]
=
∫
a
b
F
(
x
,
y
,
y
′
)
d
x
J[y]=\int_a^bF(x,y,y')\mathrm dx
J[y]=∫abF(x,y,y′)dx
的极值问题,称为等周问题。其中,
L
L
L 为常数,边界值固定
δ
y
(
a
)
=
δ
y
(
b
)
=
0
δy(a)=δy(b)=0
δy(a)=δy(b)=0
可以证明等价于泛函
J
0
[
y
]
=
J
[
y
]
−
λ
J
1
[
y
]
J_0[y]=J[y]-λJ_1[y]
J0[y]=J[y]−λJ1[y]
的极值问题。因此必要条件是
(
∂
∂
y
−
d
d
x
∂
∂
y
′
)
(
F
−
λ
G
)
=
0
\left(\cfrac{∂}{∂y}-\cfrac{\mathrm d}{\mathrm dx}\cfrac{∂}{∂y'}\right)(F-λG)=0
(∂y∂−dxd∂y′∂)(F−λG)=0
微分方程的变分法
泛函取极值的必要条件(Euler-Lagrange 方程)是常微分方程或偏微分方程,它和自变函数的定解条件结合起来,就构成一个定解问题。泛函的条件极值问题,其必要条件中出现待定参数,它和齐次边界条件结合起来,就构成本征值问题。
现在研究将微分方程的定解问题或本征值问题转化为泛函的无条件极值或条件极值问题,这称为微分方程的变分法。不难理解,本征函数正好泛函的极值函数,而本征值正好是泛函的极值。
示例 1:考虑二阶线性方程
d
d
x
[
k
(
x
)
y
′
]
+
q
(
x
)
y
−
f
(
x
)
=
0
(
a
<
x
<
b
)
\cfrac{\mathrm d}{\mathrm dx}[k(x)y'] +q(x)y-f(x)=0\quad(a<x<b)
dxd[k(x)y′]+q(x)y−f(x)=0(a<x<b)
根据变分法基本引理,方程可来源于积分
∫
a
b
[
d
d
x
(
k
y
′
)
+
q
y
−
f
]
δ
y
d
x
=
0
\int_a^b\left[\cfrac{\mathrm d}{\mathrm dx}(ky') +qy-f\right]δy\mathrm dx=0
∫ab[dxd(ky′)+qy−f]δydx=0
因为已知函数
q
(
x
)
,
f
(
x
)
q(x),f(x)
q(x),f(x) 与变分
δ
y
δy
δy 无关,因此变分计算中可看做常数,上式中第二三项可直接化为
∫
a
b
q
y
δ
y
d
x
=
1
2
δ
∫
a
b
q
y
2
d
x
∫
a
b
f
δ
y
d
x
=
δ
∫
a
b
f
y
d
x
\int_a^bqyδy\mathrm dx=\cfrac{1}{2}δ\int_a^bqy^2\mathrm dx \\ \int_a^bfδy\mathrm dx=δ\int_a^bfy\mathrm dx
∫abqyδydx=21δ∫abqy2dx∫abfδydx=δ∫abfydx
第一项可用分部积分法化为
∫
a
b
[
d
d
x
(
k
y
′
)
]
δ
y
d
x
=
k
(
x
)
y
′
δ
y
(
x
)
∣
a
b
−
∫
a
b
k
(
x
)
y
′
(
δ
y
)
′
d
x
\int_a^b\left[\cfrac{\mathrm d}{\mathrm dx}(ky')\right]δy\mathrm dx =k(x)y'δy(x)\Big|_a^b-\int_a^bk(x)y'(δy)'\mathrm dx
∫ab[dxd(ky′)]δydx=k(x)y′δy(x)∣∣∣ab−∫abk(x)y′(δy)′dx
若讨论第一边值问题,给定
y
(
a
)
=
α
,
y
(
b
)
=
β
y(a)=\alpha,\quad y(b)=\beta
y(a)=α,y(b)=β ,则
δ
y
(
a
)
=
δ
y
(
b
)
=
0
δy(a)=δy(b)=0
δy(a)=δy(b)=0 ,于是
∫
a
b
[
d
d
x
(
k
y
′
)
]
δ
y
d
x
=
−
∫
a
b
k
(
x
)
y
′
(
δ
y
)
′
d
x
=
−
1
2
δ
∫
a
b
k
(
x
)
y
′
2
d
x
\int_a^b\left[\cfrac{\mathrm d}{\mathrm dx}(ky')\right]δy\mathrm dx =-\int_a^bk(x)y'(δy)'\mathrm dx =-\cfrac{1}{2}δ\int_a^bk(x)y'^2\mathrm dx
∫ab[dxd(ky′)]δydx=−∫abk(x)y′(δy)′dx=−21δ∫abk(x)y′2dx
所以二阶线性方程可转化为泛函
J
[
y
]
=
∫
a
b
[
k
(
x
)
y
′
2
−
q
(
x
)
y
2
−
2
f
(
x
)
y
]
d
x
J[y]=\int_a^b[k(x)y'^2-q(x)y^2-2f(x)y]\mathrm dx
J[y]=∫ab[k(x)y′2−q(x)y2−2f(x)y]dx
取极值的必要条件。
几乎所有的物理和力学的基本规律都可陈述为规定某一泛函极值问题。于此,变分法使许多重要的物理问题及技术问题得以解决。
费马原理:光线永远沿用时最短的路径传播。
对于动力学系统,都遵循 Hamilton原理:在一个动力学系统中,质点系的真实运动满足积分
J
[
u
]
=
∫
t
0
t
1
(
T
−
V
)
d
t
J[u]=\int_{t_0}^{t_1}(T-V)\mathrm dt
J[u]=∫t0t1(T−V)dt
有极值的必要条件,即
δ
J
[
u
]
=
0
δJ[u]=0
δJ[u]=0 。其中
T
,
V
T,V
T,V 分别为
t
t
t 时刻系统的总动能和总势能,函数
L
=
T
−
V
L=T-V
L=T−V 称为 Lagrange 函数。
示例 2:有限长弦的强迫振动问题,
f
(
x
,
t
)
f(x,t)
f(x,t)表示横向力密度,
u
(
x
,
t
)
u(x,t)
u(x,t) 表示横向位移。
弦的总动能和势能分别为
T
=
∫
a
b
1
2
ρ
(
∂
u
∂
t
)
2
d
x
,
V
=
∫
a
b
[
1
2
T
(
∂
u
∂
x
)
2
−
f
(
x
,
t
)
u
]
d
x
T=\int_{a}^{b}\cfrac{1}{2}\rho(\cfrac{∂u}{∂t})^2\mathrm dx,\quad V=\int_{a}^{b}[\cfrac{1}{2}T(\cfrac{∂u}{∂x})^2-f(x,t)u]\mathrm dx
T=∫ab21ρ(∂t∂u)2dx,V=∫ab[21T(∂x∂u)2−f(x,t)u]dx
其中
ρ
\rho
ρ 是弦的线密度,
T
T
T 是张力,弦的端点为
a
,
b
a,b
a,b。由动力学理论知道,
u
(
x
,
t
)
u(x,t)
u(x,t) 满足泛函
J
[
u
]
=
∫
t
0
t
1
d
t
∫
a
b
[
ρ
2
(
∂
u
∂
t
)
2
−
T
2
(
∂
u
∂
x
)
2
+
f
(
x
,
t
)
u
]
d
x
J[u]=\int_{t_0}^{t_1}\mathrm dt\int_{a}^{b} [\cfrac{\rho}{2}(\cfrac{∂u}{∂t})^2-\cfrac{T}{2}(\cfrac{∂u}{∂x})^2+f(x,t)u]\mathrm dx
J[u]=∫t0t1dt∫ab[2ρ(∂t∂u)2−2T(∂x∂u)2+f(x,t)u]dx
有极值的必要条件(Euler-Lagrange 方程)
ρ
∂
2
u
∂
t
2
−
T
∂
2
u
∂
x
2
=
f
(
x
,
t
)
\rho\cfrac{∂^2u}{∂t^2}-T\cfrac{∂^2u}{∂x^2}=f(x,t)
ρ∂t2∂2u−T∂x2∂2u=f(x,t)
Rayleigh-Ritz 方法
对于微分方程,在多数实际情况下,往往只能求得近似解。在变分法的基础上,建立了实用的近似解法。
基本思路:讨论一元函数
y
(
x
)
y(x)
y(x) 泛函
J
[
y
]
=
∫
a
b
L
(
x
,
y
,
y
′
)
d
x
J[y]=\int_a^bL(x,y,y')\mathrm dx
J[y]=∫abL(x,y,y′)dx
的极值问题。可先选取一个合适的基函数序列
{
ϕ
i
(
x
)
}
\{\phi_i(x)\}
{ϕi(x)} 将函数
y
(
x
)
y(x)
y(x) 级数展开,设函数的
n
n
n 级近似解为
y
n
(
x
)
=
∑
i
=
1
n
c
i
ϕ
i
(
x
)
y_n(x)=\sum_{i=1}^nc_i\phi_i(x)
yn(x)=i=1∑nciϕi(x)
由此得
J
[
y
n
]
=
∫
a
b
L
(
x
,
y
n
,
y
n
′
)
d
x
J[y_n]=\int_a^bL(x,y_n,y_n')\mathrm dx
J[yn]=∫abL(x,yn,yn′)dx
它是线性组合系数
c
1
,
c
2
,
⋯
,
c
n
c_1,c_2,\cdots,c_n
c1,c2,⋯,cn 的函数,记为
J
(
c
1
,
c
2
,
⋯
,
c
n
)
J(c_1,c_2,\cdots,c_n)
J(c1,c2,⋯,cn) ,由多元函数取极值的必要条件知
∂
J
(
c
1
,
c
2
,
⋯
,
c
n
)
∂
c
i
=
0
(
i
=
1
,
2
,
⋯
,
n
)
\cfrac{∂J(c_1,c_2,\cdots,c_n)}{∂c_i}=0\quad (i=1,2,\cdots,n)
∂ci∂J(c1,c2,⋯,cn)=0(i=1,2,⋯,n)
这时关于
c
i
c_i
ci 的方程组,解出组合系数,从而可确定泛函极值的
n
n
n 级近似解。
非线性数学物理问题
孤立子:(soliton)
- KdV 方程
∂ u ∂ t + α u ∂ u ∂ x + ∂ 3 u ∂ x 3 = 0 \cfrac{∂u}{∂t}+\alpha u\cfrac{∂u}{∂x}+\cfrac{∂^3u}{∂x^3}=0 ∂t∂u+αu∂x∂u+∂x3∂3u=0
主要描述浅水中的表面波、含气泡的水中的声波、磁流体及等离子体中的声波等。
- Sin-Gordon 方程
∂ 2 u ∂ t 2 − ∂ 2 u ∂ x 2 + sin u = 0 \cfrac{∂^2u}{∂t^2}-\cfrac{∂^2u}{∂x^2}+\sin u=0 ∂t2∂2u−∂x2∂2u+sinu=0
主要用于描述晶体中的位错运动、约瑟夫森结中的磁通运动等。
- 非线性薛定谔方程
i ∂ u ∂ t + ∂ 2 u ∂ x 2 + β ∣ u ∣ 2 u = 0 \mathrm i\cfrac{∂u}{∂t}+\cfrac{∂^2u}{∂x^2}+\beta|u|^2u=0 i∂t∂u+∂x2∂2u+β∣u∣2u=0
主要用于描述二维平面电磁波的自聚焦、一维单色波的自调制,光纤中超短光脉冲的传播等。
- Toda 点阵方程
d q n d t = p n m d q n d t = exp ( q n − 1 − q n ) − exp ( q n − q n + 1 ) \begin{aligned} & \cfrac{\mathrm dq_n}{\mathrm dt}=\cfrac{p_n}{m} \\ & \cfrac{\mathrm dq_n}{\mathrm dt}= \exp(q_{n-1}-q_{n})-\exp(q_{n}-q_{n+1}) \end{aligned} dtdqn=mpndtdqn=exp(qn−1−qn)−exp(qn−qn+1)
主要用于描述晶格点阵中的声传播。