变分法的应用
- 我们将寻找如下形式泛函的极小化解:
J ( y ) = ∫ a b L ( x , y ( x ) , y ′ ( x ) ) d x , J(y) = \int_a^b L(x, y(x), y'(x)) \, dx, J(y)=∫abL(x,y(x),y′(x))dx,
该泛函定义在满足边界条件的光滑函数 y y y的集合 D D D上:
y ( a ) = α 和 y ( b ) = β . y(a) = \alpha \quad \text{和} \quad y(b) = \beta. y(a)=α和y(b)=β.
对于这个问题,可行变分类 A A A由在区间端点处为0的光滑函数 h h h组成,即
h ( a ) = h ( b ) = 0. h(a) = h(b) = 0. h(a)=h(b)=0.
理解起来还好吧
为了专注于变分法的主线,不去细究一般的细节定义,我们将假设 L , y L, y L,y和 h h h足够光滑。为了计算泛函极值,需要采取如下步骤:
a. 对于 y ∈ D y \in D y∈D和 h ∈ A h \in A h∈A,计算 Gâteaux 变分 δ J ( y , h ) \delta J(y, h) δJ(y,h)。
b. 求解极值方程 δ J ( y , h ) = 0 \delta J(y, h) = 0 δJ(y,h)=0。
c. 在极值中识别极小值。
在微积分中,(c)通常通过二阶导计算完成。变分法中有二阶导数测试,但我们将忽略它,转而利用几何或物理直觉来确定极小值。
- 现在我们准备好进行步骤 (a) 和 (b)。对于给定 y ∈ D y \in D y∈D和任意 h ∈ A h \in A h∈A,
δ J ( y , h ) = d d ϵ J ( y + ϵ h ) ∣ ϵ = 0 = d d ϵ ∫ a b L ( x , y + ϵ h , y ′ + ϵ h ′ ) d x ∣ ϵ = 0 = ∫ a b [ L y ( x , y , y ′ ) h + L y ′ ( x , y , y ′ ) h ′ ] d x . \delta J(y, h) \\ = \left. \frac{d}{d\epsilon} J(y + \epsilon h) \right|_{\epsilon=0}\\ = \left. \frac{d}{d\epsilon} \int_a^b L(x, y + \epsilon h, y' + \epsilon h') \, dx \right|_{\epsilon=0} \\ = \int_a^b [L_y(x, y, y')h + L_{y'}(x, y, y')h'] \, dx. δJ(y,h)=dϵdJ(y+ϵh) ϵ=0=dϵd∫abL(x,y+ϵh,y′+ϵh′)dx ϵ=0=∫ab[Ly(x,y,y′)h+Ly′(x,y,y′)h′]dx.
通过分部积分去除 h ′ h' h′上的导数,
δ J ( y , h ) = ∫ a b [ L y ( x , y , y ′ ) − d d x L y ′ ( x , y , y ′ ) ] h d x + L y ′ ( x , y , y ′ ) h ∣ a b . \delta J(y, h) = \int_a^b \left[ L_y(x, y, y') - \frac{d}{dx} L_{y'}(x, y, y') \right] h \, dx + L_{y'}(x, y, y')h \bigg|_a^b. δJ(y,h)=∫ab[Ly(x,y,y′)−dxdLy′(x,y,y′)]hdx+Ly′(x,y,y′)h ab.
由于 (2) 中的边界项消失,剩下
δ J ( y , h ) = ∫ a b [ L y ( x , y , y ′ ) − d d x L y ′ ( x , y , y ′ ) ] h d x . \delta J(y, h) = \int_a^b \left[ L_y(x, y, y') - \frac{d}{dx} L_{y'}(x, y, y') \right] h \, dx. δJ(y,h)=∫ab[Ly(x,y,y′)−dxdLy′(x,y,y′)]hdx.
因此,如果 y y y满足,
L y ( x , y ( x ) , y ′ ( x ) ) − d d x L y ′ ( x , y ( x ) , y ′ ( x ) ) = 0 , L_y(x, y(x), y'(x)) - \frac{d}{dx} L_{y'}(x, y(x), y'(x)) = 0, Ly(x,y(x),y′(x))−dxdLy′(x,y(x),y′(x))=0,
则对所有可行变分
h
h
h有
δ
J
(
y
,
h
)
=
0
\delta J(y, h) = 0
δJ(y,h)=0。
如果
y
y
y满足
L
y
(
x
,
y
(
x
)
,
y
′
(
x
)
)
−
d
d
x
L
y
′
(
x
,
y
(
x
)
,
y
′
(
x
)
)
=
0
,
L_y(x, y(x), y'(x)) - \frac{d}{dx} L_{y'}(x, y(x), y'(x)) = 0,
Ly(x,y(x),y′(x))−dxdLy′(x,y(x),y′(x))=0,
称
y
y
y满足欧拉方程.
我们将求解两点边值问题:
{
L
y
(
x
,
y
,
y
′
)
−
d
d
x
L
y
′
(
x
,
y
,
y
′
)
=
0
,
对于
a
<
x
<
b
,
y
(
a
)
=
α
,
y
(
b
)
=
β
.
\begin{cases} L_y(x, y, y') - \frac{d}{dx} L_{y'}(x, y, y') = 0, & \text{对于 } a < x < b, \\ y(a) = \alpha, \\ y(b) = \beta. \end{cases}
⎩
⎨
⎧Ly(x,y,y′)−dxdLy′(x,y,y′)=0,y(a)=α,y(b)=β.对于 a<x<b,
- 例子两点之间直线最短:考虑上次笔记中定义的弧长泛函 J J J。对于定义域
D = { y ∈ C 1 [ 0 , 1 ] ∣ y ( 0 ) = 0 , y ( 1 ) = 1 } , D = \{ y \in C^1[0, 1] \mid y(0) = 0, y(1) = 1 \}, D={y∈C1[0,1]∣y(0)=0,y(1)=1},
我们设
J ( y ) = ∫ 0 1 1 + y ′ ( x ) 2 d x . J(y) = \int_0^1 \sqrt{1 + y'(x)^2} \, dx. J(y)=∫011+y′(x)2dx.
J ( y ) J(y) J(y)是从 (0, 0) 到 (1, 1) 的光滑曲线 y = y ( x ) y = y(x) y=y(x)的长度。这两点间的最短路径是 J J J的极小化解。这个极小化解应在 J J J的极值中寻找。正如我们刚刚看到的,这些极值点是边值问题的解:
{ L y − d d x L y ′ = 0 , y ( 0 ) = 0 , y ( 1 ) = 1. \begin{cases} L_y - \frac{d}{dx} L_{y'} = 0, \\ y(0) = 0, \\ y(1) = 1. \end{cases} ⎩ ⎨ ⎧Ly−dxdLy′=0,y(0)=0,y(1)=1.
拉格朗日函数是
L ( y ′ ) = 1 + y ′ 2 , L(y') = \sqrt{1 + y'^2}, L(y′)=1+y′2,
因此欧拉方程为
d d x L y ′ = d d x ( y ′ ( x ) 1 + y ′ ( x ) 2 ) = 0. \frac{d}{dx} L_{y'} = \frac{d}{dx} \left( \frac{y'(x)}{\sqrt{1 + y'(x)^2}} \right) = 0. dxdLy′=dxd(1+y′(x)2y′(x))=0.
由此可得
y ′ ( x ) 1 + y ′ ( x ) 2 = C , \frac{y'(x)}{\sqrt{1 + y'(x)^2}} = C, 1+y′(x)2y′(x)=C,
对于某个常数 C C C。由 (11) 可知 y ′ ( x ) = A y'(x) = A y′(x)=A,其中 A 2 = C 2 1 − C 2 A^2 = \frac{C^2}{1 - C^2} A2=1−C2C2。因此
y ( x ) = A x + B , y(x) = Ax + B, y(x)=Ax+B,
对于常数 A A A和 B B B。根据边界条件, A = 1 A = 1 A=1且 B = 0 B = 0 B=0。因此 J J J的唯一极值是 y ( x ) = x y(x) = x y(x)=x。显然,这也是所寻求的极小化解 y ∗ ( x ) y^*(x) y∗(x)。
- 与拉格朗日函数 L L L对应的哈密顿量为
H = − L ( x , y , y ′ ) + y ′ L y ′ ( x , y , y ′ ) . H = -L(x, y, y') + y' L_{y'}(x, y, y'). H=−L(x,y,y′)+y′Ly′(x,y,y′).
- 二阶微分方程 F ( x , y , y ′ , y ′ ′ ) = 0 F(x, y, y', y'') = 0 F(x,y,y′,y′′)=0的一阶积分 g ( x , y , y ′ ) g(x, y, y') g(x,y,y′)是一个守恒量,即如果 y y y满足微分方程,则
g ( x , y ( x ) , y ′ ( x ) ) = 常数。 g(x, y(x), y'(x)) = \text{常数}。 g(x,y(x),y′(x))=常数。
我们通常可以为欧拉方程写下一阶积分。
a. 如果 L = L ( y , y ′ ) L = L(y, y') L=L(y,y′),则哈密顿量
H = − L ( y , y ′ ) + y ′ L y ′ ( y , y ′ ) , H = -L(y, y') + y' L_{y'}(y, y'), H=−L(y,y′)+y′Ly′(y,y′),
是欧拉方程的一阶积分。
b. 如果 L = L ( x , y ′ ) L = L(x, y') L=L(x,y′),则
L y ′ ( x , y ′ ) , L_{y'}(x, y'), Ly′(x,y′),
是一阶积分。
c. 如果 L = L ( x , y ) L = L(x, y) L=L(x,y),则欧拉方程简化为代数方程
L y ( x , y ) = 0 。 L_y(x, y) = 0。 Ly(x,y)=0。
-
示例:(最速降线问题)一个质量为 m m m的小球,最初处于静止状态,沿着曲线 y = y ( x ) y = y(x) y=y(x)从 ( a , α ) (a, \alpha) (a,α)滑动到 ( b , β ) (b, \beta) (b,β),其中 a < b a < b a<b且 α > β \alpha > \beta α>β。作用在小球上的唯一力是重力。我们将使用变分法来找到使小球的传输时间 T T T最小化的曲线 y ∗ y^* y∗。设 s s s是沿着曲线的弧长, S S S是曲线的总长度。则
d s = d x 2 + d y 2 = 1 + y ′ ( x ) 2 d x ds = \sqrt{dx^2 + dy^2} = \sqrt{1 + y'(x)^2} dx ds=dx2+dy2=1+y′(x)2dx
因此, d s d x = 1 + y ′ ( x ) 2 \frac{ds}{dx} = \sqrt{1 + y'(x)^2} dxds=1+y′(x)2。(12)
设 v v v为小球的速度。小球最初处于静止状态。因此,在 t = 0 t = 0 t=0时没有动能,总能量等于势能 m g α mg\alpha mgα。由于能量守恒, 1 2 m v 2 + m g y = m g α \frac{1}{2}mv^2 + mgy = mg\alpha 21mv2+mgy=mgα。
因此,
v = d s d t = 2 g ( α − y ) 1 + y ′ ( x ) 2 . ( 13 ) v = \frac{ds}{dt} = \sqrt{\frac{2g(\alpha - y)}{1 + y'(x)^2}}. (13) v=dtds=1+y′(x)22g(α−y).(13)
根据(12)和(13),
T = ∫ 0 T d t = ∫ 0 S d t d s d s = ∫ a b d t d s d s d x d x = ∫ a b 1 1 + y ′ ( x ) 2 2 g ( α − y ( x ) ) 1 + y ′ ( x ) 2 d x ( 14 ) T = \int_0^T dt = \int_0^S \frac{dt}{ds} ds = \int_a^b \frac{dt}{ds} \frac{ds}{dx} dx = \int_a^b \frac{1}{\sqrt{1 + y'(x)^2}} \sqrt{\frac{2g(\alpha - y(x))}{1 + y'(x)^2}} dx (14) T=∫0Tdt=∫0Sdsdtds=∫abdsdtdxdsdx=∫ab1+y′(x)211+y′(x)22g(α−y(x))dx(14)
积分(14)给出了 T T T作为 y y y的泛函。 T T T的极值函数是满足边界条件(2)的欧拉方程的解 y y y。拉格朗日量是
L ( y , y ′ ) = 1 + y ′ 2 2 g ( α − y ) L(y, y') = \frac{\sqrt{1 + y'^2}}{2g(\alpha - y)} L(y,y′)=2g(α−y)1+y′2
哈密顿量是一个第一积分。因此,
H ( y , y ′ ) = − 1 + y ′ 2 2 g ( α − y ) + y ′ 2 2 g ( α − y ) = C , ( 15 ) H(y, y') = -\frac{\sqrt{1 + y'^2}}{2g(\alpha - y)} + \frac{y'^2}{2g(\alpha - y)} = C, (15) H(y,y′)=−2g(α−y)1+y′2+2g(α−y)y′2=C,(15)
对于常数 C C C。将(15)乘以 1 2 g ( α − y ) \frac{1}{\sqrt{2g(\alpha - y)}} 2g(α−y)1,将 2 g 2g 2g吸收到一个新的常数 A A A中,并对得到的方程两边平方,得到 1 = A 2 ( 1 + y ′ 2 ) ( α − y ) 1 = A^2(1 + y'^2)(\alpha - y) 1=A2(1+y′2)(α−y)。(16)
因此, y y y满足一阶微分方程
y ′ = ± 1 − A 2 ( α − y ) A 2 ( α − y ) . ( 17 ) y' = \pm \sqrt{\frac{1 - A^2(\alpha - y)}{A^2(\alpha - y)}}. (17) y′=±A2(α−y)1−A2(α−y).(17)
由于小球应该沿着曲线下滑,我们在(17)中取负平方根。接下来,设 A 2 ( α − y ) = sin 2 θ A^2(\alpha - y) = \sin^2\theta A2(α−y)=sin2θ,(18)
其中 0 ≤ θ ≤ π / 2 0 \leq \theta \leq \pi/2 0≤θ≤π/2。从这里,再做一些微积分运算,我们得到
x = 1 A 2 ( θ − sin 2 θ ) + B , ( 19 ) x = \frac{1}{A^2}(\theta - \sin 2\theta) + B, (19) x=A21(θ−sin2θ)+B,(19)
其中 B B B是一个常数。(18)和(19)是摆线的参数方程。常数 A A A和 B B B由边界条件(2)确定。
-
显然,定义在(19)中的函数 y y y是 T T T的极值函数。我们最后需要证明,就是 y y y是一个最小值函数。此略