变分法第八章

变分法第八章第1页，课件共42页，创作于2023年2月1泛函的概念最速落径问题,如图所示A、B两点不在同一铅垂线，也不在同一高度§8.1泛函与泛函的极值ABx(x,y,)我们知道，质点下落速率与下落高度间的关系为一质点在重力作用下无磨擦沿某曲线从A滑到B，求下滑的最短时间。或沿哪条曲线用时最短。所以第2页，课件共42页，创作于2023年2月T称为y(x)的泛函y(x)可取的函数种类，称泛函的定义域,泛函是函数的涵数（不指复合函数）一般地，C是函数的集合，B是实数（或复数）的集合，若对于C中的任一称元素y(x)，在B中均有一元素J与之对应，则称J为y(x)的泛函是函数。记为与通常函数的定义不同，泛函的值决定于函数的取形。即如上例中，T的变化决定于的变化，而非某一个自变量x的值进而某一个函数y的值。而是决定于函数集合C中的函数关系，即决定于函数的取形。第3页，课件共42页，创作于2023年2月通常，泛函多以积分形式出现，如称为泛函的核其中2泛函的极值与变分在泛函的概念下，最速落径问题归结为泛函的极值问题，所谓变分法,就是求泛函的极值问题。研究泛函极值问题的方法归为两类：直接法与间接法要讨论间接法，先讨论泛函的变分问题。第4页，课件共42页，创作于2023年2月设有连续函数即导数的变分等于变分的导数,变分微分运算可交换次序。将其微小变形为其中t是一个小参数，称为的变分，记为此时，函数相应变形为第5页，课件共42页，创作于2023年2月设对x,y,y’二阶可导，y’’连续中相对于y、y’作Tayler展开抵消t的0次项，保留t的1次项，略去t的高阶项有变分dy时，泛函J的变化为则函数可得第6页，课件共42页，创作于2023年2月上式称泛函J[y（x）]第一次变分，简称变分，记为3泛函极值的必要条件——欧拉方程设泛函J[y（x）]的极值问题有解，记为y=y(x)现在来推导此解y（x）满足的常微分方程设y=y（x）有变分，则可视为t的函数表示为

当t=0时

第7页，课件共42页，创作于2023年2月亦即,F(t)函数取极值。即取极值

这样，就把原来的泛函的极值问题转变成F(t)这种普通函数的极值问题。

令即将代入上式，得即第8页，课件共42页，创作于2023年2月

泛函取极值的必要条件是其变分为0，或者说，泛函J的极值函数y(x)必须是满足泛函的变分dJ=0的函数类所以泛函的极值问题称为变分问题在简单变分问题中，端点是固定的同乘t得即又因为（分步积分）第9页，课件共42页，创作于2023年2月欧拉(Euler)方程，泛函有极值的必要条件。所以，得ⅰ）单变量多函数的泛函以上为单变量单函数泛函极值问题的欧拉方程，较复杂的泛函欧拉方程可仿照上述方法导出。如与求多元函数的偏导数相似，分别对多函数泛函之某一函数取变分，其余函数保持不变。可得i=1,2,……n第10页，课件共42页，创作于2023年2月ⅱ）高阶导数的泛函取相应的欧拉方程为或写成ⅲ）多元函数的泛函取相应的欧拉方程为第11页，课件共42页，创作于2023年2月例1最速落径问题，即求解变分问题代入得解：由于欧拉方程变形为不显含x

第12页，课件共42页，创作于2023年2月求出偏导数，有通分并取平方取得令代入上式摆线的参数方程常数c1、c2由A、B位置决定第13页，课件共42页，创作于2023年2月4泛函的条件极值问题若变量函数y（x）受到附加条件的限制，则相应的极值问题，称为条件极值问题。典型的也是最重的限制是用积分形式表示的,如即所谓等周问题

均为常数，可仿照函数条件极值问题的Lagrange乘子法，即

其中

第14页，课件共42页，创作于2023年2月将附加条件乘以参数,确定特解l，求其变分，有这是通过a和b两点的y(x)在附加条件下,使泛函取极值的必要条件。则问题转化为一般的泛函变分问题，相应的欧拉方程为关于y(x)的二阶常微分方程，一般含三个参数，即l和两个积分常数，泛函取极值的必要条件。由来确定第15页，课件共42页，创作于2023年2月例2求的极值，其中y是归一化的，即得解：此泛函的条件极值问题，可转化为变分问题代入欧拉方程，有这里且已知第16页，课件共42页，创作于2023年2月的通解为代入归一化条件，得所以而泛函的极值为使泛函取极小值p2当n=1时，泛函满足条件第17页，课件共42页，创作于2023年2月5求泛函极值的直接方法(Ritz方法)从泛函自身出发，不经微分方程直接求出极值曲线，称为泛函极值问题的直接方法。Ritz方法—典型的直接方法：要点是不将其放在它全部定义域来考虑，而是在定义域的某一部分来考虑。使J转化为

设某种完备的函数系

试偿以其中的前几项来表示变分问题dJ=0的解

其中

为待定系数

的n元函数第18页，课件共42页，创作于2023年2月所以按多元函数求极值的方法，令不过这样得出的函数并非变分问题dJ=0的严格解

由于f的形式是我们预先选定的，比如即

由此解出

便确定出了函数y(x)而是近似解,记为yn(x)，严格解应为

Ritz法中函数系ji的选取至关重要，如何选取？

第19页，课件共42页，创作于2023年2月例3用Ritz方法求例2。即求采用试探解项的选取是为了满足解：以作为选取的函数系将其代入得下的变分问题。在约束条件且已知第20页，课件共42页，创作于2023年2月由即结果是第21页，课件共42页，创作于2023年2月把代入得显而易见，在c1=0时，J[y(x)]最小，最小值为10所以对比近似解，抛物线严格解，正弦曲线

且第22页，课件共42页，创作于2023年2月1）把偏微分方程的本征值问题或定解问题，与泛函的极值问题联系起来，使原来的方程是泛函的欧拉方程；2）用直接方法求出泛函的极值函数，由于此函数一定满足欧拉方程，所以，也一定满足原方程，即一定是原方程的解。用变分法求数理方程的基本原理本节以Helmhotz方程的本征值问题和Poisson方程的边值问题为例，讨论把上述问题转化为泛函极值问题或变分法的基本方法，然后来求解极值问题（用直接方法）。§8.2用分法求解数理方程第23页，课件共42页，创作于2023年2月（设u在区域t内有连续的二阶导数，l为参数,s为t的边界）取泛函令1本征问题与变分问题的关系Helmhotz本征值问题由第一格林公式则有或第24页，课件共42页，创作于2023年2月其中对应的欧拉方程为对于三元函数的泛函，其变分问题为所以第25页，课件共42页，创作于2023年2月即泛函中把代入欧拉方程，得欧拉方程变为极值问题的欧拉方程就是Helmhotz方程在边界条件下本征值问题而且，此泛函变分问题与泛函在附加条件就是说，Helmhotz方程的本征值问题，可归结为归一条件下J1[u]的极值问题。下的变分问题等价。第26页，课件共42页，创作于2023年2月所对应的泛函同样为若为第二类边界条件同样亦有即本征值问题若为第三类边界条件类似地有第27页，课件共42页，创作于2023年2月则本征值问题

记和边界条件下的极值问题可归结为在附加条件求泛函2泛函极值与本征值问题的关系仍以Helmhotz方程为例，先给出一重要结论：的最小值l0就是本征值问题泛函的最小本征值,而使泛函J1[u]在边界条件第28页，课件共42页，创作于2023年2月和附加条件u0就是该本征值问题对应本征值l0的本征函数。取得最小值的函数结论的证明：有最小值l0的极值函数，则有设u0是使泛函由边界条件知第29页，课件共42页，创作于2023年2月的欧拉方程为又，在附加条件下所以u0满足或代入J1[u0],有u0是本征函数。再证明：即l0是本征值，设l0是最小本征值。第30页，课件共42页，创作于2023年2月相应的本征函数为u1则有这与是的最小值相矛盾结论得证。有次小值l1的极值函数，类似地还可证明，若设u1是使泛函且满足边界条件和附加条件除此之外，还同时满足与u0正交的条件。即设第31页，课件共42页，创作于2023年2月相应的本征函数为u1满足依此类推，泛函取第i个极值的极值函数ui满足且满足边界条件和附加条件除此之外，还同时满足正交条件即由此得到的泛函的次极小值就是本征值问题的次极小值对于一系列本征值相应的本征函数为第32页，课件共42页，创作于2023年2月例用变分法求边界固定的圆膜横振动的本征振动。代入上式，得引入无量纲变量解：取平面极坐标，定解问题为令旋转对称

记得第33页，课件共42页，创作于2023年2月这是一个二阶常微分方程的本征值问题，用变分法对于任意的二阶常微分方程的本征值问题，形如在归一化条件能够证明，可转化归结为：及相应边界条件下求泛函的极值问题第34页，课件共42页，创作于2023年2月二方程对比在归一化条件有下，求泛函的极值问题所以方程的求解第35页，课件共42页，创作于2023年2月采用直接方法（Ritz方法）求解令代入归一化条件和泛函，得（如此取形使x=0处不出现尖点）算出各积分，得I,J

两个关于c1，c2的函数为第36页，课件共42页，创作于2023年2月由Lagrange乘子法，取极值的条件为其中:k=lb2c1、c2

非零解存在的条件是：解出k的两个解为，第37页，课件共42页，创作于2023年2月最小本征值为因l=b2/k将其代入c的方程和归一化条件解出本例的严格解可由分离变量法得出，结果为为最小本征值为相应的本征函数称为零阶Bessel函数，称为零阶Bessel函数的第一个零点。第38页，课件共42页，创作于2023年2月第一类边值问题3边值问题与变分问题的关系以Poisson问题为例来讨论s为t的边界取对取变分，有第39页，课件共42页，创作于2023年2月但由泛函取极值的条件为，所以相应的欧拉方程为（前式利