数学物理方法变分法

《数学物理方法变分法》由会员分享，可在线阅读，更多相关《数学物理方法变分法（42页珍藏版）》请在装配图网上搜索。

1、数学物理方法变分法数学物理方法变分法 从前面的定解问题的解法中，我们容易想到由于边界形从前面的定解问题的解法中，我们容易想到由于边界形状较为复杂，或由于泛定方程较为复杂，或由于其它各种条状较为复杂，或由于泛定方程较为复杂，或由于其它各种条件发生变化，将使得定解问题难以严格解出，因此又发展了件发生变化，将使得定解问题难以严格解出，因此又发展了一些切实可用的近似方法，通过本章的学习我们会看到近似一些切实可用的近似方法，通过本章的学习我们会看到近似解的价值一点也不低于严格解的价值事实上，我们应该已解的价值一点也不低于严格解的价值事实上，我们应该已经注意到，从推导数学物理方程时难免要作一些简化假定，经

2、注意到，从推导数学物理方程时难免要作一些简化假定，定解条件本身也带有或多或少的近似性，前面所谓的严格解定解条件本身也带有或多或少的近似性，前面所谓的严格解其实也是某种程度的近似其实也是某种程度的近似数学物理方法变分法数学物理方法变分法 如果某个定解问题不能严格解出，但另一个与它差别如果某个定解问题不能严格解出，但另一个与它差别甚微的定解问题能严格解出，那么就可以运用微扰法求近甚微的定解问题能严格解出，那么就可以运用微扰法求近似解量子力学教科书中一般都要介绍微扰法，限于课时，似解量子力学教科书中一般都要介绍微扰法，限于课时，这里就不再重复介绍这里就不再重复介绍近似解法涉及近似解法涉及：变分法，有

3、限差分法和模拟法等变分法，有限差分法和模拟法等 变分法变分法是研究求解泛函极值（极大或极小）的方法，变是研究求解泛函极值（极大或极小）的方法，变分问题即是分问题即是求泛函的极值问题求泛函的极值问题把定解问题转化为变分把定解问题转化为变分问题，再求变分问题的解问题，再求变分问题的解数学物理方法变分法数学物理方法变分法变分法的优点变分法的优点:(2)(2)变分法易于变分法易于实现数学的统一化实现数学的统一化因为一般而言，数因为一般而言，数学物理方程的定解问题都可以转化为变分问题尤其是前学物理方程的定解问题都可以转化为变分问题尤其是前面介绍的斯特姆刘维尔本征值问题可转化为变分问题，面介绍的斯特姆刘维

4、尔本征值问题可转化为变分问题，变分法提供了施刘型本征值问题的本征函数系的完备性变分法提供了施刘型本征值问题的本征函数系的完备性等结论的证明；等结论的证明；(1)(1)变分法在物理上可以归纳定律因为几乎所有的变分法在物理上可以归纳定律因为几乎所有的自然定律都能用变分原理的形式予以表达；自然定律都能用变分原理的形式予以表达；数学物理方法变分法数学物理方法变分法(3)(3)变分法是解数学物理定解问题常用的近似方法，变分法是解数学物理定解问题常用的近似方法，其基本思想是把数学物理定解问题转化为变分问题由其基本思想是把数学物理定解问题转化为变分问题由直接解变分问题发展了一些近似解法，其中最有用的直接解变

5、分问题发展了一些近似解法，其中最有用的是里茨是里茨 （RitzRitz）法）法 由于里茨法中的试探函数的选由于里茨法中的试探函数的选取较为麻烦，计算系数矩阵也十分困难，随着计算机取较为麻烦，计算系数矩阵也十分困难，随着计算机的展，又迅速发展了一种有限元法；的展，又迅速发展了一种有限元法；(4)(4)变分法的应用不仅在经典物理和工程技术域，而变分法的应用不仅在经典物理和工程技术域，而且在现代量子场论，现代控制理论和现代信息理论等且在现代量子场论，现代控制理论和现代信息理论等高技术领域都有十分广泛的应用高技术领域都有十分广泛的应用数学物理方法变分法数学物理方法变分法有限差分法有限差分法：有限差分法

6、把定解问题转化为代数方程，：有限差分法把定解问题转化为代数方程，然后通过电子计算机求定解问题的数值解然后通过电子计算机求定解问题的数值解模拟法模拟法：即用一定的物理模型来模拟所研究的定解问题，：即用一定的物理模型来模拟所研究的定解问题，而在模型上实测解的数值而在模型上实测解的数值 变分法变分法是这些方法中最为重要和切实有效的方法，是这些方法中最为重要和切实有效的方法，已经广泛应用于科学研究和工程计算之中，限于篇幅故已经广泛应用于科学研究和工程计算之中，限于篇幅故本书主要详细介绍经典变分法的基本概念和理论本书主要详细介绍经典变分法的基本概念和理论数学物理方法变分法数学物理方法变分法定义：定义：变

7、分法变分法 变分问题变分问题 变分法就是求泛函极值的方法变分问题即是求变分法就是求泛函极值的方法变分问题即是求泛函的极值问题泛函的极值问题一、泛函一、泛函 变分法研究的对象是泛函，泛函是函数概念的推广变分法研究的对象是泛函，泛函是函数概念的推广为了说明泛函概念先看一个例题：为了说明泛函概念先看一个例题：数学物理方法变分法数学物理方法变分法 考虑著名的考虑著名的最速降线落径问题最速降线落径问题。如图。如图13.1 13.1 所示，所示，已知已知A A和和B B为不在同一铅垂线和不同高度的两点，要求为不在同一铅垂线和不同高度的两点，要求找出找出A A、B B间的这样一条曲线，当一质点在重力作用下沿

8、间的这样一条曲线，当一质点在重力作用下沿这条曲线无摩擦地从这条曲线无摩擦地从A A滑到滑到B B时，所需的时间时，所需的时间T T最小最小 y x A B(x,y)图 19.1 图图13.1数学物理方法变分法数学物理方法变分法我们知道，此时质点的速度是我们知道，此时质点的速度是 d2dsgyt因此从因此从 A A滑到滑到B B所需的时间为所需的时间为21+ddd22BAtBBtAAysTtxgygy即为即为21+()d2BAyT y xxgy（13.1.1）数学物理方法变分法数学物理方法变分法yxT()y x()y x()T y x式中式中 代表对代表对求一阶导数求一阶导数 我们称上述的我们称

9、上述的为为的的泛函泛函，而称，而称为可取的函数类，为泛函为可取的函数类，为泛函的定义域。简单地说，泛函就是函数的函数（不是复合函数的定义域。简单地说，泛函就是函数的函数（不是复合函数的那种含义）的那种含义）一般来说，设一般来说，设C C是是函数的集合函数的集合，B B是是实数或复数的集合实数或复数的集合，如果对于如果对于C C的任一元素的任一元素 ()y x在在B B中都有一个元素中都有一个元素J与之对应，与之对应，J()y x()JJ y x则称则称为为的泛函，记为的泛函，记为必须注意，泛函不同于通常讲的函数决定通常函数值的必须注意，泛函不同于通常讲的函数决定通常函数值的数学物理方法变分法数

10、学物理方法变分法因素是自变量的取值，而决定泛函的值的因素则是函数的取因素是自变量的取值，而决定泛函的值的因素则是函数的取形如上面例子中的泛函形如上面例子中的泛函T T的变化是由函数的变化是由函数 ()y xx（即从（即从A A到到B B的不同曲线）的不同曲线）值，也不取决值，也不取决所引起的它的值既不取决于某一个所引起的它的值既不取决于某一个本身的变化本身的变化于某一个于某一个 yyx值，而是取决于整个集合值，而是取决于整个集合C C中中与与的函数关系的函数关系定义：泛函定义：泛函 泛函的核泛函的核 泛函通常以积分形式出现，比如上面描述的最速降线泛函通常以积分形式出现，比如上面描述的最速降线落

11、径问题的式（落径问题的式（13.1.113.1.1）更为一般而又典型的泛函定义为）更为一般而又典型的泛函定义为()(,)dbaJ y xF x y yx(13.1.2)其中其中 (,)F x y y称为称为泛函的核泛函的核数学物理方法变分法数学物理方法变分法 二、泛函的极值二、泛函的极值变分法变分法对于不同的自变量函数对于不同的自变量函数 ()y x，与此相应的泛函，与此相应的泛函 ()J y x也有不同的数值找出一个确定的自变量函数也有不同的数值找出一个确定的自变量函数 ()y x，使泛函，使泛函 ()J y x 具有极值（极小或极大），这种泛函的极小值与极大具有极值（极小或极大），这种泛函

12、的极小值与极大值统称为泛函的极值值统称为泛函的极值引入泛函的概念后，对于上述的最速降线落径问题变引入泛函的概念后，对于上述的最速降线落径问题变为泛函为泛函()J y x的极小值问题物理学中常见的有光学的极小值问题物理学中常见的有光学中的中的费马费马(Fermat)(Fermat)原理，分析力学中的哈密顿原理，分析力学中的哈密顿(Hamiton)(Hamiton)原理等，都是泛函的极值原理等，都是泛函的极值问题问题数学物理方法变分法数学物理方法变分法即直接分析所提出的问题；另一类叫间接法，即把问题转化为即直接分析所提出的问题；另一类叫间接法，即把问题转化为求解微分方程为讨论间接方法，先介绍变分和

13、泛函的变分求解微分方程为讨论间接方法，先介绍变分和泛函的变分三、三、变分变分 定义：定义：变分变分 如果我们将泛函取极值时的函数（或函数曲线）定义为如果我们将泛函取极值时的函数（或函数曲线）定义为 ();y x并定义与函数曲线并定义与函数曲线 ()y x邻近的曲线（或略为变形的邻近的曲线（或略为变形的定义：定义：变分法：所谓的变分法变分法：所谓的变分法就是求泛函极值的方法就是求泛函极值的方法研究泛函极值问题的方法可以归为两类：一类叫直接法，研究泛函极值问题的方法可以归为两类：一类叫直接法，数学物理方法变分法数学物理方法变分法曲线）作为比较曲线，记为曲线）作为比较曲线，记为(,)()()y xy

14、 xx其中其中 是一个小参数；是一个小参数；()x是一个具有二阶导数的任意是一个具有二阶导数的任意选定函数，规定选定函数，规定 它在一个小范围内变化，这限制主要保证泛它在一个小范围内变化，这限制主要保证泛函在极值处连续在研究泛函极值时，通常将函在极值处连续在研究泛函极值时，通常将 ()x固定，固定，而令而令变化，这样规定的好处在于：建立了由参数变化，这样规定的好处在于：建立了由参数 到泛函到泛函 ()J y x值之间的对应关系，因此泛函值之间的对应关系，因此泛函 ()J y x就成为了参数就成为了参数 的普通函数原来泛函的极值问题就成为的普通函数原来泛函的极值问题就成为数学物理方法变分法数学物

15、理方法变分法普通函数对普通函数对 的求极值的问题同时，函数曲线的求极值的问题同时，函数曲线 ()y x的变分定义为的变分定义为0(,)|()dyy xx(13.1.3)(13.1.3)因此可得因此可得 ()dyx(13.1.4)(13.1.4)这里这里 ,y代表对代表对x求一阶导数求一阶导数 所以所以 ddyyx(13.1.5)(13.1.5)即变分和微分可以交换次序即变分和微分可以交换次序 数学物理方法变分法数学物理方法变分法()dbaFFJyyxyy(13.1.6)(13.1.6)在极值曲线在极值曲线 ()y x附近，泛函附近，泛函 ()J y x的增量，定义为的增量，定义为(,)()JJ

16、 y xJ y x(13.1.7)(13.1.7)依照上述约定，当依照上述约定，当 0时，泛函增量时，泛函增量 J的线性的线性主要部分定义为泛函的变分，记为主要部分定义为泛函的变分，记为 四、四、泛函的变分泛函的变分定义：定义：泛函的变分泛函的变分 泛函的增量泛函的增量 变分问题变分问题泛函的变分定义为泛函的变分定义为0|dJJ(13.1.8)数学物理方法变分法数学物理方法变分法 在求一元或多元函数的极值时，微分起了很大的作用；同样在求一元或多元函数的极值时，微分起了很大的作用；同样在研究泛函极值问题时，变分起着类似微分的作用因此，通常在研究泛函极值问题时，变分起着类似微分的作用因此，通常称泛

17、函极值问题为变分问题；称求泛函极值的方法为变分法称泛函极值问题为变分问题；称求泛函极值的方法为变分法解解 1721()()dJ y xy exyx1172711111771111111()()d(2)dd 2dd2d|d 2dJ y xyexyxxy y eyxxy y xey xxy y x eyxxy y x注意：最后一步利用了一般在边界上函数变分为零的事实，即注意：最后一步利用了一般在边界上函数变分为零的事实，即 711|0ey例例 1 计算泛函的变分计算泛函的变分数学物理方法变分法数学物理方法变分法13132 2 泛函的极值泛函的极值 泛函的极值问题，一般来说是比较复杂的因为它与泛泛函

18、的极值问题，一般来说是比较复杂的因为它与泛函包含的自变量个数，未知函数的个数以及函数导数的阶数函包含的自变量个数，未知函数的个数以及函数导数的阶数等相关另外，在求泛函极值时，有的还要加约束条件，且等相关另外，在求泛函极值时，有的还要加约束条件，且约束条件的类型也有不同，等等下面我们首先讨论泛函的约束条件的类型也有不同，等等下面我们首先讨论泛函的极值的必要条件极值的必要条件 数学物理方法变分法数学物理方法变分法一、一、泛函的极值的必要条件泛函的极值的必要条件欧拉拉格朗日方程欧拉拉格朗日方程 设设 ()J y x的极值问题有解的极值问题有解()yy x（13.2.1）现在推导这个解所满足的常微分方

19、程，这是用间接法现在推导这个解所满足的常微分方程，这是用间接法研究泛函极值问题的重要一环设想这个解有变分研究泛函极值问题的重要一环设想这个解有变分 ()x则则 ()()J y xx可视为参数可视为参数 的函数的函数 ()()().J y xx而当而当 0时，时，数学物理方法变分法数学物理方法变分法()()()y xxy x对应于式对应于式(13.2.1),(13.2.1),即为即为 ()()J y xx取极值于是原来的泛函极值取极值于是原来的泛函极值问题，就化为一个求普通函数问题，就化为一个求普通函数 ()的极值问题由函数的极值问题由函数取极值的必要条件，有取极值的必要条件，有0d|0d即有即

20、有 0|0J（13.2.213.2.2）数学物理方法变分法数学物理方法变分法 1.泛函表示为一个自变量，一个函数及其一阶导数泛函表示为一个自变量，一个函数及其一阶导数的积分形式的积分形式泛函表示为一个自变量，一个函数及其一阶导数的积分形式，泛函表示为一个自变量，一个函数及其一阶导数的积分形式，（13.1.213.1.2）()(,)baJ y xF x y y dx若考虑两端固定边界的泛函问题若考虑两端固定边界的泛函问题：积分是在区域内通过两点积分是在区域内通过两点 1122(,),(,)x yxy的任意曲线进行的，其中的任意曲线进行的，其中 12,xa xb数学物理方法变分法数学物理方法变分法

21、泛函中泛函中 y为为(,)()()y xy xx由于两端固定，所以要求由于两端固定，所以要求 ()0,()0ab，即，即|0,|0 x ax byy由由(13.1.8)(13.1.8)，有，有 0()()|d ()d()d d dbabaJ y xxJFFxxxyyFFyyxyy（13.2.313.2.3）数学物理方法变分法数学物理方法变分法式式(13.2.3)(13.2.3)的积分号下既有的积分号下既有 y，又有，又有 y，对第二项，对第二项应用分部积分法可使积分号下出现应用分部积分法可使积分号下出现yd|()ddbbaaFFFJyy xyyxy(13.2.4)(13.2.4)根据（根据（1

22、7.2.217.2.2）,所以所以 0|0JJd,再根据再根据（13.2.413.2.4）故有）故有d|()d0dbbaaFFFJyy xyyxy（13.2.513.2.5）数学物理方法变分法数学物理方法变分法 因为因为|0,|0 x ax byy并且并且 y是任意的，所以是任意的，所以 d()0dFFyxy (13.2.6)(13.2.6)上式上式(13.2.6)(13.2.6)称为欧拉（称为欧拉（EulerEuler）拉格朗日（）拉格朗日（LagrangeLagrange）方程，简称为方程，简称为E-LE-L方程方程 此即泛函取极值的必要条件即泛函此即泛函取极值的必要条件即泛函 J的极值函

23、数的极值函数 ()y x必须是满足泛函的变分必须是满足泛函的变分 0J的函数类的函数类 ()y x因此，因此，数学物理方法变分法数学物理方法变分法把泛函的极值问题称为变分问题把泛函的极值问题称为变分问题 注明：注明：E-LE-L方程是泛函取极值的必要条件，而不是充分方程是泛函取极值的必要条件，而不是充分条件如果讨论充分条件，则要计算二阶变分，并考虑其正、条件如果讨论充分条件，则要计算二阶变分，并考虑其正、负值负值,但对于实际问题中，当泛函具有明确的物理涵义，极值的但对于实际问题中，当泛函具有明确的物理涵义，极值的存在性往往间接地在问题的提法中就可以肯定，所以极值的存存在性往往间接地在问题的提法

24、中就可以肯定，所以极值的存在性是不成问题的，只要解出在性是不成问题的，只要解出E-LE-L方程，就可以得到泛函的极方程，就可以得到泛函的极值值 E-L E-L方程除了上面给出的形式方程除了上面给出的形式(13.2.6)(13.2.6)之外，之外，另外还有四种特殊情况：另外还有四种特殊情况：数学物理方法变分法数学物理方法变分法(1)(1)F不显含不显含 x(,),FF y y且且 0Fx因为因为ddd()()()dddFFFFFFFyyyxyxyxyyxy若若 0,y E-LE-L方程等价于方程等价于 FFycy （13.2.713.2.7）数学物理方法变分法数学物理方法变分法(2)(2)F不依

25、赖于不依赖于 y(,),FF x y且且 0Fy则则E-LE-L方程化为方程化为d()0,dFFcxyy （13.2.813.2.8）(3)(3)F不依赖于不依赖于 y(,),FF x y且且 0Fy则则E-LE-L方程化为方程化为0Fy（13.2.913.2.9）数学物理方法变分法数学物理方法变分法由此可见由此可见 F仅为仅为 x的函数的函数 (4)(4)F关于关于 y是线性的：是线性的：(,)(,)(,)F x y yf x yg x y y则则E-LE-L方程化为方程化为0fgyx （13.2.1013.2.10）对于含有一个自变量，多个变量函数，以及有较高阶变量对于含有一个自变量，多个

26、变量函数，以及有较高阶变量函数导数的泛函，类似上面的推导可得如下结论：函数导数的泛函，类似上面的推导可得如下结论：数学物理方法变分法数学物理方法变分法2.泛函表示为多个函数的积分形式泛函表示为多个函数的积分形式1122()(,)dbnnaJ y xF x y y yyyyx|0,|=0 (1,2,)ix aix byyin则与此泛函极值问题相应的则与此泛函极值问题相应的E-LE-L方程为方程为d()0 (1,2,)diiFFinyxy（13.2.11）数学物理方法变分法数学物理方法变分法3.泛函的积分形式中含有高阶导数泛函的积分形式中含有高阶导数()()(,)dbnaJ y xF x y y

27、yyx(1)()()()0ny ay ay a(1)()()()0ny by by b与此泛函极值问题相应的与此泛函极值问题相应的E-LE-L方程为方程为22()ddd()()(1)()0dddnnnnFFFFyxyxyxy（13.2.12）数学物理方法变分法数学物理方法变分法4.泛函的积分形式中含有多元函数泛函的积分形式中含有多元函数(,)u x y,x y设设为为的二元函数，则的二元函数，则22111212(,)d d(,)(,)(,)(,)0 xyxyxyJF x y u u ux yu x yu xyu x yu x y 与此泛函极值问题相应的与此泛函极值问题相应的E-LE-L方程为方

28、程为()()0 xyFFFuxuyu（13.2.1313.2.13）数学物理方法变分法数学物理方法变分法212yFgy不显含不显含 x，故其，故其E-LE-L方程为（方程为（13.2.713.2.7）式）式0221122yyFFyygycyygy令令 02cgc故有故有 221(1)yyc例例2 2 试求解最速降线落径问题，即变分问题试求解最速降线落径问题，即变分问题21d02BAyxgy解目前，我们只能用间接方法来求解，由于解目前，我们只能用间接方法来求解，由于数学物理方法变分法数学物理方法变分法令令 121cc分离变量得到分离变量得到1dd ycyxy再令再令 12sin2cy代入上式得到

29、代入上式得到112dsind(1 cos)d22cxc即得到即得到121(sin)2(1c o s2)cccxy数学物理方法变分法数学物理方法变分法此即为摆线的参数方程，积分常数可由初始位置此即为摆线的参数方程，积分常数可由初始位置 （图（图13.113.1的的A,BA,B两点）决定两点）决定13.2.213.2.2泛函的条件极值问题泛函的条件极值问题 在许多泛函的极值问题中，变量函数还受到一些附加条件在许多泛函的极值问题中，变量函数还受到一些附加条件的限制，其中最常见和重要的一种是以积分形式表示的限制的限制，其中最常见和重要的一种是以积分形式表示的限制条件条件(,)dbaG x y yxl（

30、13.2.14）数学物理方法变分法数学物理方法变分法即所谓的等周问题：即所谓的等周问题：01()(,)d,(),()(,)d babaJ y xF x y yxy ayy byG x y yxl(13.2.15)(13.2.15)（注：这种问题之所以称为等周问题，是因为在历史上起源注：这种问题之所以称为等周问题，是因为在历史上起源于求一条通过两点，长度固定为于求一条通过两点，长度固定为l l的曲线的曲线 (),y x使面积使面积 ()dbaSy xx取极大值）取极大值）数学物理方法变分法数学物理方法变分法其中其中 01,l yy为常数此类问题可以仿照普通函数的为常数此类问题可以仿照普通函数的条

31、件极值问题的拉格朗日乘子法即将附加条件条件极值问题的拉格朗日乘子法即将附加条件(13.2.14)(13.2.14)乘以乘以参数，求其变分后，加到泛函取极值的必要条件中得到参数，求其变分后，加到泛函取极值的必要条件中得到(;,)(;,)d0baF x y yG x y yx于是问题转化为不带条件的由上式所表示的变分问题于是问题转化为不带条件的由上式所表示的变分问题 其对应的其对应的E-LE-L方程为方程为数学物理方法变分法数学物理方法变分法d()0dFGFGyyxyy这是通过这是通过 a和和 b两点的两点的 ()y x之下使泛函取极值的必要条件它实际上是一个关于之下使泛函取极值的必要条件它实际上

32、是一个关于 在附加条件（在附加条件（13.2.1413.2.14）()y x的二阶常微分方程其通解中含有三个参数，即的二阶常微分方程其通解中含有三个参数，即和两个积分和两个积分常数它们可由条件常数它们可由条件 01(),()y ayy by（13.2.1413.2.14）来确定）来确定 .和附加条件和附加条件数学物理方法变分法数学物理方法变分法 例例3 3 求求 120()()dJ y xyx的极值，其中的极值，其中 y是归一化的，即是归一化的，即 120d1yx，且已知，且已知 (0)0,(1)0.yy 解解 本题是求泛函的条件极值问题，可化为变分问题本题是求泛函的条件极值问题，可化为变分问

33、题1220()d0yyx 对应的对应的E-LE-L方程为方程为0yy其通解为其通解为cossin (0)yAxBx数学物理方法变分法数学物理方法变分法代入附加条件代入附加条件 (0)0,(1)0.yy得到得到()sin()(1,2,)nnyxcn xn代入归一化条件得到代入归一化条件得到1220sin()d1ncn x x 于是得到于是得到 2nc ，故原极值问题的解为，故原极值问题的解为2sin()nyn x 而题中要求的泛函而题中要求的泛函 120()dyx的极值为的极值为数学物理方法变分法数学物理方法变分法12222202 cos()dnn xxn当当 1n 时，极值函数时，极值函数 1

34、()2siny xx 使得泛函数取得最小值使得泛函数取得最小值 2例例4 4 求泛函求泛函 20(2 cos)dJ yyyxx在条件在条件 (0)0,()0yy下的极值曲线下的极值曲线.解解 此时此时 xyyyyxFcos2),(2 则偏导数则偏导数 数学物理方法变分法数学物理方法变分法yFxFyy2,cos2.对应的对应的EulerEuler方程为方程为0cos xy其通解为其通解为 21cosCxCxy代入边界条件可得代入边界条件可得 12C 12C所以极值曲线为所以极值曲线为 2cos1yxx数学物理方法变分法数学物理方法变分法泛函极值问题的求解泛函极值问题的求解,通常有两种结果：通常有

35、两种结果：（i i）解析解）解析解 由变分法得到的由变分法得到的E-LE-L方程求解，一般来说，是很困难的方程求解，一般来说，是很困难的但在分析力学中往往还是采用这一办法来求解因为历史悠但在分析力学中往往还是采用这一办法来求解因为历史悠久，它自有一套办法久，它自有一套办法（iiii）近似解）近似解 所谓近似解即由泛函本身出发，而不需求解所谓近似解即由泛函本身出发，而不需求解E-LE-L方程，方程，直接求得所需要的解直接求得所需要的解极值曲线极值曲线 因此，常常称它为研究泛函极值问题的直接法因此，常常称它为研究泛函极值问题的直接法数学物理方法变分法数学物理方法变分法例例5 假设大气的光折射率假设

36、大气的光折射率 n只依赖于高度只依赖于高度 y利用费马（利用费马（FermatFermat）原理导出在大气中光线轨迹的微分方程；）原理导出在大气中光线轨迹的微分方程；解解（1 1）根据费马原理：光线的实际路径上，光程的变分为零）根据费马原理：光线的实际路径上，光程的变分为零21d0ttn l（13.3.1）其中其中 n为介质中的光折射率，为介质中的光折射率，dl为沿光线进行方向的路程为沿光线进行方向的路程元上述问题也可表示为如下泛函极值问题：元上述问题也可表示为如下泛函极值问题：2()()1dbaI y xn yyx（13.3.2）由于由于 (,)F x y y不显含不显含 x，根据公式，根据公式(13.2.7)(13.2.7)，可得首次积分，可得首次积分02()1n ycy （13.3.3）