这篇文章是为了回答一位知友的提问而作的向未了解过变分法这个概念的读者介绍它的文章,因此本文只涉及一些变分法的入门知识。在从极值条件导出欧拉方程时只涉及泛函变分的Lagrange定义,这样一来就不必引入关于分离泛函线性主部的讨论。行文时并不过分追求数学的严密性并只假定读者具有基础的微积分尤其是复合函数求取偏导数的知识。

1.泛函与泛函极值

提起变分法这个名词就不得不先来谈谈泛函的概念,撇开数学上的定义不谈让我们先来看几个例子:

首先就是题主问到的最速降线问题,它是Galileo于1630年首先提出来的。如章首图所示在平面直角坐标系里给定点

,在区间

上有一条曲线

连结这两个点。现在让我们在

点放上一个小球,这样它就由于受到重力作用沿着曲线运动起来,运动到

点的时间很容易计算出来:

(如果读者觉得这不容易的话阅读这篇文章并不是个明智之举)。这里我们看到所用的时间

是依赖于函数

的具体形式的,现在的问题是如何选择

的形式能够使

取得最小值?

再来举一个吹肥皂泡的例子:想要玩吹泡泡这个游戏看来我们需要准备一些工具——一个由铁丝围起来的封闭形状。我们的生活经验会毫不犹豫地告诉我们肥皂液在铁丝上形成一个平面而不是凹凸不平的曲面(略去重力的影响)。在物理上这件事的解释是液体的膜具有一种叫做表面张力的性质,某种程度上读者可以把肥皂膜理解为橡皮膜,它有使自己绷紧达到表面积最小的趋势,显然平面比凹凸不平的面具有更小的面积。不过如果我们的铁丝圈本来不在一个平面上会怎么样呢?这样一来无论如何肥皂膜再不能是平面了。我们假设那个铁丝圈由曲线方程

决定,肥皂膜的方程由定义在

上(闭曲线

在坐标平面上的投影)的函数

表示,它的面积就可以表示成二重积分

,现在的问题是

取何种形式能使得这积分式最小?

最后来举一个力学上的例子。对于某些确定的力学系统我们可以用一些变量来决定它的位形。比如说一个弹簧振子,决定它的位形只需要用其偏离平衡位置的坐标就好了(设为

),如果我们想知道这振子的运动规律只需要知道

依赖于时间的形式即

。那这件事如何办到呢?这可以从一个物理原理——最小作用量原理出发。实验表明存在这样的一个函数

(动能减势能),真实的运动方程使得积分

在给定边界条件(初末态的位置)时取得最小值。现在我们该怎样求得

呢?

我们看到上面的这些例子里我们所关注的变量都依赖于一个函数的表达形式,这种以函数为自变量的函数就被称作是泛函。而要解决这些问题就是要定出当泛函取得极值时自变量函数的形式。

2.依赖于一个一元函数的固定边界最简泛函极值问题

下面就让我们来考虑一种特定形式的泛函并且来想一想如何求出它的极值。设

是依赖于三个独立变量

的函数,设待求解函数

应满足在端点处取固定值,记作

,泛函则形如

这种特定形式的泛函极值问题就被称为固定边界的最简泛函极值问题。类似于求取普通的函数极值我们有微分法,求取泛函极值的一种特定的方法被称作变分法。我们马上就将看到这方法具体是如何实施的。

对于给定的函数

,我们可以在其上加上一个改变量

使它变为另一个函数

,这被加上的改变量称作是函数

处的变分。为了使

仍满足固定边界条件我们取

。值得指出的是变分与微分是很不同的,后者是由于自变量的改变而引起的函数值的改变,而前者指自变量取相同值时函数的变化量。显然这时导数的变分

,后面我们将把

直接记作

。现在在变分

前乘上一个数

并构造以此数为自变量的函数:

.

如果在这个式子中

已经使得泛函取极值,那在

的情况下

任取为不等于0的数都会使得函数

的值变大(或变小),换句话说

是函数

的极值点,此时应有

。于是:

上面的推导中用到了求导与积分交换顺序和复合函数求偏导数的数学定理。对积分的第二项使用分部积分法:

.

由于变分

在边界处的取值为0,所以积分号外的项值为0。而

是能任意选取的函数所以若要积分式恒为0必有其前面的系数

恒成立,这个微分方程就是泛函取极值时应满足的微分方程,它是瑞士数学家Euler于1736年得到的。

Euler方程是求解泛函极值问题的极好工具,因为它将定义在整个区间上的问题转化成了只需考虑函数临近部分关系的微分方程问题。在笔者看来能做到这点的原因是与最简泛函的形式密切相关的,如果泛函在整体上具有某种取得最小值的性质,由于积分是各个分段区间上的求和,整体的最值条件要求被积函数在分段的小区间上也必须满足取最值的条件,因此本来是在有限区间上的问题就能被我们转化成微分方程了。

Euler方程可以被推广到依赖于多元函数的情形,对于文章第一部分提到的二元函数情形,极值曲面应满足的微分方程是

3.应用举例

【例1】求解最速降线问题。

【解】在这个问题里被积表达式不显含

,Euler方程等价于

(请读者试着推导) ,也就是说

为一常数。

取参数

使得

,于是

。而

,积分求得

的表达式为

。带入边界条件可以定出任意常数的值,稍微化简一下得到:

这是一条摆线的方程。

【例2】(弦振动方程)考虑一根长为

的两端固定的弦,设其线密度为

,弦上具有恒定张力

(由于微振动故假设张力在振动过程中不改变)。求解弦线做微小振动时位移满足的微分方程。弦振动问题

【解】设时刻

时位于

处质元的纵向位移为

,遵照本文之前提到的最小作用量原理,先来表示弦线具有的动能和势能。

在弦线形变不大时将它的伸长量精确到一阶,即

。势能化作

现在泛函

应取极值, 利用上文提到的推广到二元函数的Euler方程就得到弦振动方程

.

熟悉这方程的读者想必知道

即是弦线上波动的传播速度。

结合泛函极值_泛函极值与变分法相关推荐

  1. 结合泛函极值_泛函极值及变分法教程.doc

    PAGE PAGE 53 第二章 泛函极值及变分法(补充内容)2.1 变分的基本概念2.1.1 泛函和变分泛函是一种广义的函数,是指对于某一类函数{y(x)}中的每一个函数y(x),变量J有一值与之对 ...

  2. 结合泛函极值_泛函极值及变分法讲义.doc

    第二章 泛函极值及变分法(补充内容) 2.1 变分的基本概念 2.1.1 泛函和变分 泛函是一种广义的函数,是指对于某一类函数{y(x)}中的每一个函数y(x),变量J有一值与之对应,或者说数J对应于 ...

  3. 结合泛函极值_泛函的极值

    相关文献 0 引 言在许多工程中,特别在可靠性最优化问题上,被广为实践应用的是被抽象为表示为分式的线性泛函数G(f)=M(f)N(f)(1)式中M.N分别表示在某一分布函数f(x)(x∈R(n),n& ...

  4. python广义极值_广义极值(GEV)极大似然拟合的奇异pdf

    首先,我认为您可能希望将位置参数固定在0.在 其次,你的数据中有0,结果拟合可能在x=0处有{}pdf,例如对于GEV拟合或Weibull拟合. 因此,拟合实际上是正确的,但是当绘制pdf(包括x=0 ...

  5. 结合泛函极值_第2章泛函的极值.doc

    第2章泛函的极值 第2章 泛函的极值 , , 如果, 当 (或者说)时, 有 那么, 我们称在处是连续的, 记为. 2.1.2 函数的可微性 更进一步, 如果存在, 使得 那么我们称在处是可微的, 或 ...

  6. 05 ,函数的极值,最值,极值定理,极值与导 0 的关系 :

    1 ,极值,极值点 : 定义 : 去心邻域内的最大值或者最小值 2 ,最值 : 定义域内的最大值或者最小值 3 ,极值 vs 最值 : 极值 : 局部性概念 最值 :全局性概念 比较 : 4 ,极值定 ...

  7. 结合泛函极值_(二) 泛函的极值

    极值的概念 函数 \(f(x)\) 在 \(x_0\) 处取得极小值,是指当 \(x\) 在 \(x_0\) 点及其附近 \(|x - x_0| < \varepsilon\) 时,恒有 \(f ...

  8. python迭代法求极值_用Python实现最速下降法求极值的方法

    对于一个多元函数 ,用最速下降法(又称梯度下降法)求其极小值的迭代格式为 其中 为负梯度方向,即最速下降方向,αkαk为搜索步长. 一般情况下,最优步长αkαk的确定要用到线性搜索技术,比如精确线性搜 ...

  9. 判断极值_拐点的三个充要条件以及一个必要

    极值 一阶可导点是极值点的必要条件 设 f(x) 在 x=x0 处可导,且在点 x0 处取得极值,则必有 又是费马定理 判断极值的第一充分条件 判断极值的第二充分条件 判断极值的第三充分条件 设 f( ...

最新文章

  1. linux硬盘满了问题排查
  2. c语言中count的头文件,求助C语言大佬 , 只会写到一个.c文件里 ,不会用.h头文件...
  3. SqlServer 行转列
  4. 放置奇兵 算法 月度活动 破碎时空记录 1-3关
  5. php如何求同列元素之和_求得这个数组中各个元素之和
  6. Tensorflow深度学习应用(进阶篇)
  7. 机器学习与计算机视觉(数据集的选择)
  8. 导航猫(NaviCat for MySql)建立表的方法
  9. select top 1 oracle,在ORACLE中SELECT TOP N的实现方法
  10. VS2015安装使用番茄助手Visual Assist
  11. 水经注地图发布服务中间件服务配置功能说明
  12. 前端特效——简单下雪(纯css)
  13. 请启封,您的双态IT大会邀请函
  14. 微信小程序自定义tabbar
  15. aics6圆角插件怎么安装_AI圆角插件!助你更快设计字体 Xtream Path1.4的安装使用教程...
  16. 2019辽宁公务员考试行测常识大全:公务员常识40000问(四十八)(2)
  17. Windows锁屏图片文件存放位置
  18. Linux服务器开发,手把手设计实现epoll
  19. Ackerman阿克曼(单舵轮)前轮驱动并转向移动机器人的建模及控制算法
  20. 如何使用ghost制作系统镜像

热门文章

  1. PayPal轮询收单系统升级之PayPalme亲友支付
  2. 怎么免费把pdf转换成jpg图片
  3. sklearn - Dimensionality reduction
  4. c/c++ -nan(ind) NAN
  5. python 知乎 合并 pdf_实例4:用Python提取不同PDF文件中的页面合并进新的PDF文件...
  6. oracle 下bdump占用大,【学习笔记】Oracle oradebug 利用oradebug释放被删除文件空间的案例...
  7. C# Aspose.Words 用法
  8. eNSP之简单WLAN拓扑配置
  9. 【NLP】BERT语言模型
  10. 查看已结束的中国大学MOOC课程