第十一章变分法
§11.1泛函和泛函的极值
11.1.1泛函
1696年,Basel大学Bernoulli提出,最速落径:
v=dsdt =2gy − − − √ ∴总下降时间
T=∫ t 2 (B) t 1 (A) dt=∫ B A dsv ,由于ds=dx 2 +dy 2 − − − − − − − − √ =1+(y ′ ) 2 − − − − − − − √ dx
∴T=∫ B A 1+(y ′ ) 2 − − − − − − − √ dx2gy − − − √
即T=T[y(x)]=∫ B A 1+(y ′ ) 2 − − − − − − − √ dx2gy − − − √ ,−是函数的函数
1.定义
泛函是函数的函数。即J=J[y(x)]。其中,J∈B:复(实);数集y(x)∈C:复数集.
注意:(1)不同于普通函数
(2)不同于复合函数
(3)定义域:y(x)可取类
2.典型表达方式:J[y(x)]=∫ b a F(x,y,y ′ )dx→(1)
F(x,y,y ′ )泛函的核
3.最速落径问题:求泛函T[y(x)]的极小值问题。
11.1.2泛函的极值
类似于上述的求极值问题,如光学中费马原理,力学中最小作用力问题。在物理学中很多,我们将会看到求泛函极值可归结为两种方法.
1.求泛函极值方法
(1)直接方法:从泛函直接求极值曲线
(2)间接法:化为解微分方程−与变分问题有联系
为此,建立有关变分概念。
2.变分
(1):函数的变分:若y(x)−→ 微变 y(x)+tη(x),t为小参数,则记δ(y)=tη(x)(2)–称tη(x)为y(x)的变分.注意:δy不同于dy,dy有一取极值过程,δy不取极限,略去了高阶小量.此时,
y ′ (x)=lim Δx→0 ΔyΔx →lim Δx→0 Δ(y+tη)Δx =y ′ (x)+tη ′ (x)
(2)变分与微分可交换次序:
δ(y ′ )=tη ′ (x)=ddx [tη(x)]=ddx δ(y)→y ′ (x)的变分
即:δ(y ′ )≡ddx δ(y)
(3):函数的变分
若(1)中F∈C 2 ,y∈C 2 ,则当y→y+tη
ΔJ=J[y(x)+tη]−J[y]
=∫ b a [F(x,y+tη,y ′ +tη ′ )−F(x,y,y ′ )]dx
=∫ b a [∂F∂y tη+∂F∂y ′ tη ′ +t的高阶小量]dx
即δJ=∫ b a (∂F∂y δy+∂F∂y ′ δy ′ )dx→(4)−称为泛函J[y(x)]的第一次变分
11.1.3泛函极值的必要条件
设J[y(x)]的极值函数为y(x),且有变分tη(x),则J[y(x)+tη(x)]=ϕ(t)(因为y(x)已经设定),由一元函数取极值条件,则J极值条件
→ϕ(t)极值条件→ϕdt | t=0 =0
∂J[y(x)+tη(x)]∂t | t=0 =0,即,∫ b a ∂F(x,y+tη,y ′ +tη ′ ]∂t | t=0 dx=0
∫ b a [∂F∂(y+tη) η+∂F∂(y ′ +tη ′ ) η ′ ]| t=0 dx=0,
即:∫ b a [∂F∂y η+∂F∂y ′ η ′ ]dx=0→∫ b a [∂F∂y δy+∂F∂y ′ δy ′ ]dx=0(5)
1.泛函极值的问题是变分问题
对照(4):极值条件δJ=0有
∵∫ b a ∂F∂y ′ δy ′ dx=∫ b a ∂F∂y ′ ddx (δy)dx
=∂F∂y ′ δy| b a −∫ b a ddx (∂F∂y ′ )δydx=−∫ b a ddx (∂F∂y ′ )δydx
代入(5):∫ b a [∂F∂y −ddx (∂F∂y ′ )]δydx=0
∴∂F∂y −ddx (∂F∂y ′ )=0→Euler方程
2.泛函取极值的条件−y(x)满足Euler方程
(1)对于J[y(x)]=∫ b a F(x,y,y ′ )dx
有∂F∂y −ddx (∂F∂y ′ )=0→二阶常微分方程⇒y ′ ∂F∂y ′ −F=C
(2)对于J[y 1 (x),y 2 (x),⋯,y n (x)]=∫ b a F(x;y 1 ,y 2 ,⋯,y n ;y ′ 1 ,y ′ 2 ,⋯,y ′ n )dx
有∂F∂y i −ddx (∂F∂y ′ i )=0,(i=0,1,2,⋯,n)
(3)对于J[y ′ (x),y ′ (x)⋯y (n) (x)]=∫ b a F(x;y;y ′ ,y ′′ ⋯y (n) )dx
有∂F∂y −ddx (∂F∂y ′ )+d 2 dx 2 (∂F∂y ′′ )−d 3 dx 3 (∂F∂y ′′′ )=0
(4)多元函数J[u(x,y,z)]=∫ b a F(x,y,z;u;u x ,u y ,u z )dxdydz
∂F∂u −∂∂x (∂F∂u x )−∂∂y (∂F∂u y )−∂∂z (∂F∂u z )=0
于是求泛函数极值问题→解Euler方程问题
3.若F不显含x则
∂F∂y −ddx (∂F∂y ′ )=0→y ′ ∂F∂y ′ −F=C
考虑ddx [y ′ ∂F∂y ′ −F]=y ′′ ∂F∂y ′ +y ′ ddx (∂F∂y ′ )−∂F∂x −∂F∂y y ′ −∂F∂y ′ y ′′
=−y ′ [∂F∂y −ddx (∂F∂y ′ )]−∂F∂x
∵F不显示含x,即∂F∂x =0⇒ddx [y ′ ∂F∂y ′ −F]=0,
∴积分可得y ′ ∂F∂y ′ −F=C(B)
例1.求最速落径
T[y(x)]=∫ x B x A 1+(y ′ ) 2 − − − − − − − √ dx2gy − − − √
∴δ∫ x B x A 1+(y ′ ) 2 − − − − − − − √ dx2gy − − − √ =0
F=1+(y ′ ) 2 − − − − − − − √ dx2gy − − − √ =12g − − √ 1+(y ′ ) 2 − − − − − − − √ y √ 不显含x
于是Eulerf方程
y ′ ⋅∂∂y ′ 1+(y ′ ) 2 − − − − − − − √ y √ −1+(y ′ ) 2 − − − − − − − √ y √ =C
(y ′ ) 2 [1+(y ′ ) 2 ]y − − − − − − − − − √ −1+(y ′ ) 2 − − − − − − − √ y √ =C
⇒(y ′ ) 4 [1+(y ′ ) 2 ]y +1+(y ′ ) 2 y −2y ′2 y =C 2
1[1+(y ′ ) 2 ]y =C 2 ,令1C 2 =C 1
得y ′ =C 1 −y − − − − − √ y √ ,→x−C 2 =∫y √ C 1 −y − − − − − √ dy
令y=C 1 sin 2 θ2 ,则⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ x=C 1 2 (θ−sinθ)+C 2 y=C 1 2 (1−cosθ)
是由半径为C 1 2 的圆周上一固定点运动产生的.
在途中x轴下方滚动.存在一条且仅一条通过原点及点(x b ,y b )的摆线.适当选择C 1 和C 2 可以给出这条摆线.
11.1.4泛函的条件极值
⎧ ⎩ ⎨ J[y(x)]=∫ b a F(x,y,y ′ )dx∫ b a G(x,y,y ′ )dx=l ,y(a)=y 0 ,y(b)=y 1
拉格朗日(Lagrange)乘子法:
考虑:δ∫ b a [F(x,y,y ′ )+λG(x,y,y ′ )]dx=0
则→∂F∂y +λ∂G∂y −ddx [(∂F∂y ′ )+λddx (∂G∂y ′ )]=0
积分常数C 1 ,C 2 和λ可由附加条件定出
例2.⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ J[y(x)]=∫ 1 0 y ′2 dx∫ 1 0 y 2 dx=1y(0)=0,y(1)=1
考虑δ∫ 1 0 (y ′2 +λy 2 )dx=0
不显含x,也可推出一阶Euler方程,此处直接用二阶Euler也不困难:2λy−ddx (2y ′ )=0
即:y ′′ −λy=0
y=C 1 e λ √ x +C 2 e −λ √ x
由{y(0)=0y(1)=1 得:y n =C n sinnπx(n=1,2,⋯)
再由∫ 1 0 y 2 dx=1得:C n =±2 √ ,∴y n =±2 √ sinnπx
J[y(x)]=∫ 1 0 [ddx ±2 √ sinnπx] 2 dx=n 2 π 2
极小值为J[y 1 (x)]=π 2
11.1.5里兹方法
1.对于J[f(x)],
令y(x)=f(φ 1 (x),φ 2 (x)⋯φ n (x);C 1 ,C 2 ⋯C n )
则J[f(x)]=φ(C 1 ,C 2 ,⋯,C n )
于是当∂φ∂C i =0时,J[f(x)]取极值,i=1,2,⋯,n
2.注意
(1)f是近似解f(x)=lim n→∞ f(φ 1 ,φ 2 ,⋯,φ n ;C 1 ,C 2 ,⋯,C n )
(2)适当选φ,f
注:一般为多项式三角式为它们的线性组合,满足边界条件
3.例:求J[f(x)]=∫ 1 0 y ′2 dx(的极小值)(1)
其中{∫ 1 0 y 2 dx=1y(0)=0,y(1)=0 (2)
解:选φ n (x)=C n x n
令y(x)=x(x−1)(C 0 +C 1 x)=C 1 x 3 +(C 0 −C 1 )x 2 −C 0 x(3)
代入(1):ϕ=J[y(x)]=∫ 1 0 [3C 1 x 2 +2(C 0 −C 1 )x−C 0 ] 2 dx
=13 (C 2 0 +C 0 C 1 +dfrac25C 2 1 )(4),
即ϕ=ϕ(C 0 ,C 1 )
代入(2):∫ 1 0 y 2 dx=130 (C 2 0 +C 0 C 1 +27 C 2 1 )=1(5)
即ψ=130 (C 2 0 +C 0 C 1 +27 C 2 1 )−1=0
由拉格朗日乘子法:若要求y(x)=f(x 1 ,x 2 ,⋯,x n )在m个约束条件:
g k (x 1 ,x 2 ,⋯,x n )=0,(k=1,2,⋯,m)
下的极值只需考虑F=y+∑ k=1 m λ k g k
使⎧ ⎩ ⎨ ∂F∂x =0g k (x i )=0 i=1,2,⋯,n即可
∴应由⎧ ⎩ ⎨ ⎪ ⎪ ∂F∂C i =∂ϕ∂C i +λ∂ψ∂C i =0φ=0 联立消去λ来求y
但在此不必,由(5):C 2 0 +C 0 C 1 =30−27 C 2 1 (6)
代入(4):ϕ=J[y(x)]=13 (30−27 C 2 1 +25 C 2 1 )
=23 (15+235 C 2 1 ),C 1 =0有极小值,此时J[f(x)]=10
C 1 =0代入(6):C 0 =±30 − − √ ,代入(4)
∴y(x)=±30 − − √ x(x−1)=±30 − − √ x(x−1)
=±30 − − √ [(x−12 ) 2 −14 ]
(x−12 ) 2 =1±30 − − √ y+14 =±21230 − − √ (y±30 − − √ 4 )
令h=12 ,k=−(±30 − − √ 4 ),p=1230 − − √
⇒(x−h) 2 =±2p(y−k)
顶点:(h,k);焦点:(h,k±p2 );准线:y=−k±p2
§11.2用变分法解数理方程
11.2.1变分法
1.思路
泛函极值←泛函的Euler方程的解
↑
直接方法求解
2.步骤
(1)写出定解问题所对应的泛函
(2)用直接法(里兹法)求该泛函极值,此即写方程的解.