理论教育 微分方程的变分问题改进:土建领域数值解法

微分方程的变分问题改进:土建领域数值解法

时间:2023-10-26 理论教育 版权反馈
【摘要】:本节最后我们将证明,如果L是对称正定线性算子,那么使得泛函式取得极小值的函数y就是微分方程式定解问题的解。泛函的变分计算可以通过微分方程求解来代替,而微分方程的求解也可通过泛函的变分计算来代替。

微分方程的变分问题改进:土建领域数值解法

在第1章介绍过微分算子和差分算子,本节对算子作进一步介绍。

1.内积

对于域Ω=D+Γ中的函数uv,称

uv的内积。内积是一个实数,容易验证内积具有以下性质:

uv>=<vu

C1u1+C2u2v>=C1u1v>+C2u2v

uu>≥0,且等式<uu>=0,当且仅当u≡0时成立。其中,uvu1u2Ω中的任意函数,C1C2为任意实数。

2.算子

若某一集合D与另一集合S建立了某种一一对应的关系,即D中任一元素u对应于S中的一个元素Lu,则称L为算子。集合D称为算子L定义域,集合S称为算子L的值域。对于拉普拉斯算子

定义域D是所有二次可微函数的集合,值域S是微分后生成的函数集合。

对于一般微分算子L,我们只讨论下列几种:

线性算子 若对于任意实数αβ,算子L具有性质

Lαu1+βu2)=αLu1)+βLu2

则称L为线性算子,其中,u1u2Dαu1+βu2D

对称算子 若对于线性算子L,定义域中的任意两个元素uv,都有

Lu),v>=<uLv)>

则称L为对称算子,其中,uvD

正定算子 若对于线性算子L定义域中的任意一个非零元素u,都有

Lu),u>>0

则称L为正定算子,其中,uD

算子方程 如果在算子L的值域S中给定一个元素f,则称等式

Lu)=f

为算子方程,其中,uD

考虑如下的二阶常微方程的边值问题

假定px)在区间[ab]上连续,且一阶可导,qx)和fx)在区间[ab]上连续,并且px)>0,qx)≥0。求函数yx),要求在区间[ab]上连续,且二阶可导。

现在我们尝试构造泛函

Jy)=<Ly),y>-2<fy> (2.3.16)

其中,

我们可以得到,算子L是对称正定的线性算子。先讨论其对称性,事实上,对于任意一个函数zx)∈D,应用分部积分公式以及式(2.3.15),有

可见,算子L是对称的。再看正定性,为此,令上式中的zx)=yx),有

可以看出,由于px)>0,qx)≥0,得

Ly),y>≥0

若<Ly),y>=0,即有

px)[y′x)]2+qx)[yx)]2=0

于是

y′x)=0

进而得

yx)=CC为常数)

边界条件yx=a=0,得

yx)=0

可见,算子L是正定的。(www.daowen.com)

本节最后我们将证明,如果L是对称正定线性算子,那么使得泛函式(2.3.16)取得极小值的函数yx)就是微分方程式(2.3.14)定解问题的解。

下面我们来证明泊松方程的算子L也是对称正定线性算子。对于二维情况

u|Γ=0 (2.3.18)

要求函数uxy)在求解域D内有连续二阶偏导数,在域D内和边界Γ上连续。设函数vxy)在求解域D内有连续二阶偏导数,根据内积的定义及格林公式,有

由于v满足定解问题的边界条件,故上式的边界积分项为零,于是

观察式(2.3.19)可见,等式右端关于uv是对称的,所以有

Lu),v>=<Lv),u

因此,泊松方程的算子是对称的。再看正定性,为此,令式(2.3.19)中的u=v,有

对于式(2.3.19),当<-Δuv>=0时,uxy)的一阶偏导数等于零,也就是uxy)为常数,根据该函数的边界条件得出uxyΓ=0,于是得出,泊松方程的算子为正定算子。

3.椭圆型方程边值问题和相应变分问题的等价性

定理L是对称正线性算子,若算子方程

Lu)=fuD) (2.3.21)

存在解u=u0,则u0所满足的充要条件是泛函

u=u0处取极小值。

证明 先证明必要性,即若u=u0是算子方程式(2.3.21),则有

Lu0)-f=0

D中任意函数

u=u0+εη

其中,η为任意连续函数;ε为一绝对值较小的参变量,且与变量xy等无关。于是

Ju)=Ju0+εη

=<Lu0+εη),u0+εη>-2<fu0+εη

L是对称正定算子,根据内积性质,上式可以展开

Ju)=<Lu0),u0>+<Lu0),εη>+<Lεη),u0>+<Lεη),εη>-2<fεη>-2<fu0

=Ju0)+<Lεη),εη>+<Lεη),u0>+<Lu0),εη>-2<fεη

=Ju0)+<Lεη),εη>+2<Lu0),εη>-2<fεη

=Ju0)+<Lεη),εη>+2<Lu0)-fεη

由于Lu0)-f=0,又从正定算子定义

Lεη),εη>0

于是

Ju)>Ju0

这说明泛函Ju),在u=u0时取极小值。

再证明充分性。因u=u0时,Ju)取极小值,于是

L是对称正定算子,有

Ju0+εη)=Ju0)+<Lεη),εη>+2<Lu0),εη>-2<fεη

=Ju0)+ε2Lη),η>+2εLu0)-fη

于是

由于η为任意连续函数,因此有

Lu0)-f=0

u0是方程式(2.3.21)的解。(定理证毕)

这样,就得到一个重要结论,对同一个物理过程,存在两种互为等价的数学提法,即在满足边界条件下的微分方程提法和变分提法。泛函的变分计算可以通过微分方程求解来代替,而微分方程的求解也可通过泛函的变分计算来代替。一般说来,若算子是正定对称,则采用里兹法,可以证明只要近似解构成的函数序列{un}是泛函的极小化序列,则un收敛于算子方程的精确解。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈