变分法第八章

1、变分法第八章第1页，共42页，2022年，5月20日，4点54分，星期二1 泛函的概念最速落径问题,如图所示A、B两点不在同一铅垂线，也不在同一高度8.1泛函与泛函的极值ABx(x,y,)我们知道，质点下落速率与下落高度间的关系为一质点在重力作用下无磨擦沿某曲线从A滑到B，求下滑的最短时间。或沿哪条曲线用时最短。所以第2页，共42页，2022年，5月20日，4点54分，星期二T称为y(x)的泛函y(x)可取的函数种类，称泛函的定义域,泛函是函数的涵数（不指复合函数）一般地， C是函数的集合， B是实数（或复数）的集合，若对于C中的任一称元素y(x) ，在B中均有一元素J与之对应，则称J为y(x

2、) 的泛函是函数。记为与通常函数的定义不同，泛函的值决定于函数的取形。即如上例中，T的变化决定于 的变化，而非某一个自变量x的值进而某一个函数y的值。而是决定于函数集合C中的函数关系，即决定于函数的取形。第3页，共42页，2022年，5月20日，4点54分，星期二通常，泛函多以积分形式出现，如称为泛函的核其中2 泛函的极值与变分在泛函的概念下，最速落径问题归结为泛函的极值问题，所谓变分法,就是求泛函的极值问题。研究泛函极值问题的方法归为两类：直接法与间接法要讨论间接法，先讨论泛函的变分问题。第4页，共42页，2022年，5月20日，4点54分，星期二设有连续函数即导数的变分等于变分的导数,变分

3、微分运算可交换次序。将其微小变形为其中t是一个小参数，称为 的变分，记为此时，函数相应变形为第5页，共42页，2022年，5月20日，4点54分，星期二设对x, y, y二阶可导，y连续中相对于y、y作Tayler展开抵消t的0次项，保留t的1次项，略去t的高阶项有变分dy 时，泛函J的变化为则函数可得第6页，共42页，2022年，5月20日，4点54分，星期二上式称泛函 J y（x）第一次变分，简称变分，记为3 泛函极值的必要条件欧拉方程设泛函 J y（x）的极值问题有解，记为y = y(x)现在来推导此解y（x）满足的常微分方程设y=y（x）有变分 ， 则可视为t 的函数 表示为 当t=0

4、时 第7页，共42页，2022年，5月20日，4点54分，星期二亦即, F(t)函数取极值。即 取极值 这样，就把原来的泛函的极值问题转变成F(t)这种普通函数的极值问题。 令即将代入上式，得即第8页，共42页，2022年，5月20日，4点54分，星期二 泛函取极值的必要条件是其变分为0，或者说，泛函J的极值函数y(x)必须是满足泛函的变分dJ=0的函数类所以泛函的极值问题称为变分问题在简单变分问题中，端点是固定的同乘t 得即又因为（分步积分）第9页，共42页，2022年，5月20日，4点54分，星期二欧拉(Euler)方程，泛函有极值的必要条件。所以，得）单变量多函数的泛函以上为单变量单函数

5、泛函极值问题的欧拉方程，较复杂的泛函欧拉方程可仿照上述方法导出。如与求多元函数的偏导数相似，分别对多函数泛函之某一函数取变分，其余函数保持不变。可得i=1,2, n第10页，共42页，2022年，5月20日，4点54分，星期二）高阶导数的泛函取相应的欧拉方程为或写成）多元函数的泛函取相应的欧拉方程为第11页，共42页，2022年，5月20日，4点54分，星期二例1 最速落径问题，即求解变分问题 代入得解：由于 欧拉方程变形为 不显含x 第12页，共42页，2022年，5月20日，4点54分，星期二求出偏导数，有 通分并取平方取得 令代入上式摆线的参数方程常数c1 、c2由A、B位置决定第13页

6、，共42页，2022年，5月20日，4点54分，星期二4 泛函的条件极值问题若变量函数 y（x）受到附加条件的限制，则相应的极值问题，称为条件极值问题。典型的也是最重的限制是用积分形式表示的,如即所谓等周问题 均为常数，可仿照函数条件极值问题的Lagrange乘子法，即 其中 第14页，共42页，2022年，5月20日，4点54分，星期二将附加条件乘以参数,确定特解l，求其变分，有这是通过a和b两点的y(x)在附加条件下,使泛函取极值的必要条件。则问题转化为一般的泛函变分问题，相应的欧拉方程为关于y(x)的二阶常微分方程，一般含三个参数，即l和两个积分常数，泛函取极值的必要条件。由来确定第15

7、页，共42页，2022年，5月20日，4点54分，星期二例2 求 的极值，其中y是归一化的，即 得解：此泛函的条件极值问题，可转化为变分问题 代入欧拉方程，有 这里且已知第16页，共42页，2022年，5月20日，4点54分，星期二的通解为代入归一化条件，得所以而泛函的极值为使泛函取极小值 p2当n=1时，泛函满足条件第17页，共42页，2022年，5月20日，4点54分，星期二5 求泛函极值的直接方法(Ritz 方法)从泛函自身出发，不经微分方程直接求出极值曲线，称为泛函极值问题的直接方法。Ritz 方法典型的直接方法：要点是不将其放在它全部定义域来考虑，而是在定义域的某一部分来考虑。使J转

8、化为 设某种完备的函数系 试偿以其中的前几项来表示变分问题 dJ = 0 的解 其中 为待定系数 的n 元函数第18页，共42页，2022年，5月20日，4点54分，星期二所以按多元函数求极值的方法，令不过这样得出的函数并非变分问题dJ = 0的严格解 由于f的形式是我们预先选定的，比如即 由此解出 便确定出了函数y(x)而是近似解,记为yn(x)，严格解应为 Ritz法中函数系ji的选取至关重要，如何选取？ 第19页，共42页，2022年，5月20日，4点54分，星期二例3 用Ritz方法求例2。即求采用试探解项的选取是为了满足解：以 作为选取的函数系 将其代入得下的变分问题。在约束条件且已

9、知第20页，共42页，2022年，5月20日，4点54分，星期二由即结果是第21页，共42页，2022年，5月20日，4点54分，星期二把代入得显而易见，在c1=0时，Jy(x)最小，最小值为10所以对比近似解，抛物线严格解，正弦曲线 且第22页，共42页，2022年，5月20日，4点54分，星期二1）把偏微分方程的本征值问题或定解问题，与泛函的极值问题联系起来，使原来的方程是泛函的欧拉方程；2）用直接方法求出泛函的极值函数，由于此函数一定满足欧拉方程，所以，也一定满足原方程，即一定是原方程的解。用变分法求数理方程的基本原理本节以Helmhotz方程的本征值问题和Poisson方程的边值问题为

10、例，讨论把上述问题转化为泛函极值问题或变分法的基本方法，然后来求解极值问题（用直接方法）。8.2用分法求解数理方程第23页，共42页，2022年，5月20日，4点54分，星期二（设u在区域t内有连续的二阶导数，l为参数, s为t的边界）取泛函令1 本征问题与变分问题的关系Helmhotz本征值问题由第一格林公式则有或第24页，共42页，2022年，5月20日，4点54分，星期二其中对应的欧拉方程为对于三元函数的泛函，其变分问题为所以第25页，共42页，2022年，5月20日，4点54分，星期二即泛函中把代入欧拉方程，得欧拉方程变为极值问题的欧拉方程就是Helmhotz方程在 边界条件下本征值问

11、题而且，此泛函变分问题与泛函 在附加条件就是说，Helmhotz方程的本征值问题，可归结为归一条件下J1u的极值问题。 下的变分问题等价。第26页，共42页，2022年，5月20日，4点54分，星期二所对应的泛函同样为若为第二类边界条件 同样亦有即本征值问题若为第三类边界条件 类似地有第27页，共42页，2022年，5月20日，4点54分，星期二则本征值问题 记和边界条件 下 的极值问题可归结为在附加条件 求泛函2 泛函极值与本征值问题的关系仍以Helmhotz方程为例，先给出一重要结论：的最小值l0就是本征值问题泛函的最小本征值,而使泛函J1u在边界条件第28页，共42页，2022年，5月2

12、0日，4点54分，星期二和附加条件u0就是该本征值问题对应本征值l0的本征函数。取得最小值的函数结论的证明：有最小值l0 的极值函数，则有设u0是使泛函由边界条件知第29页，共42页，2022年，5月20日，4点54分，星期二的欧拉方程为 又，在附加条件下所以u0满足 或代入J1u0,有u0是本征函数。再证明：即l0 是本征值，设l0是最小本征值。第30页，共42页，2022年，5月20日，4点54分，星期二相应的本征函数为u1 则有这与是的最小值相矛盾结论得证。有次小值l1 的极值函数，类似地还可证明，若设u1是使泛函且满足边界条件和附加条件除此之外，还同时满足与u0正交的条件。即设第31页

13、，共42页，2022年，5月20日，4点54分，星期二相应的本征函数为u1 满足依此类推，泛函取第i个极值的极值函数ui 满足且满足边界条件和附加条件除此之外，还同时满足正交条件即由此得到的泛函的次极小值就是本征值问题的次极小值对于一系列本征值相应的本征函数为第32页，共42页，2022年，5月20日，4点54分，星期二例 用变分法求边界固定的圆膜横振动的本征振动。 代入上式，得引入无量纲变量解：取平面极坐标，定解问题为 令 旋转对称 记得第33页，共42页，2022年，5月20日，4点54分，星期二这是一个二阶常微分方程的本征值问题，用变分法对于任意的二阶常微分方程的本征值问题，形如在归一化

14、条件能够证明，可转化归结为：及相应边界条件下求泛函的极值问题第34页，共42页，2022年，5月20日，4点54分，星期二二方程对比在归一化条件有下，求泛函的极值问题所以方程的求解第35页，共42页，2022年，5月20日，4点54分，星期二采用直接方法（Ritz方法）求解 令 代入 归一化条件和泛函，得（如此取形使x=0处不出现尖点）算出各积分，得I,J 两个关于c1，c2的函数为第36页，共42页，2022年，5月20日，4点54分，星期二由Lagrange乘子法，取极值的条件为 其中:k=lb2c1 、c2 非零解存在的条件是：解出k的两个解为， 第37页，共42页，2022年，5月20日，4点54分，星期二最小本征值为 因l=b2/k将其代入c的方程和归一化条件解出本例的严格解可由分离变量法得出，结果为为最小本征值为相应的本征函数称为零阶Bessel函数，称为零阶Bessel函数的第一个零点。第38页，共42页，2022年，5月20日，4点54分，星期二第一类边值问题3 边值问题与变分问题的关系以Poisson问题为例来讨论s为t的边界取 对取变分，有第39页，共42页，2022年，5月20日，4点54分，星期二 但由泛函取极值的条件为，所以相应的欧拉方程为（前式利用格林