变分法简介(简单_明了_易懂).doc

1§1 变分法简介作为数学的一个分支，变分法的诞生，是现实世界许多现象不断探索的结果，人们可以追寻到这样一个轨迹：约翰·伯努利（Johann Bernoulli，1667－1748）1696 年向全欧洲数学家挑战，提出一个难题：“设在垂直平面内有任意两点，一个质点受地心引力的作用，自较高点下滑至较低点，不计摩擦，问沿着什么曲线下滑，时间最短？”这就是著名的“最速降线”问题（The Brachistochrone Problem） 它的难处在于和普通的极大极小值求法不同，它是要求出一个未知函数（曲线） ，来满足所给的条件这问题的新颖和别出心裁引起了很大兴趣，罗比塔（Guillaume Francois Antonie de l'Hospital 1661-1704） 、雅可比·伯努利（Jacob Bernoulli 1654-1705） 、莱布尼茨（Gottfried Wilhelm Leibniz,1646-1716）和牛顿（Isaac Newton1642—1727）都得到了解答约翰的解法比较漂亮，而雅可布的解法虽然麻烦与费劲，却更为一般化后来欧拉（Euler Lonhard，1707～1783）和拉格朗日(Lagrange, Joseph Louis，1736-1813)发明了这一类问题的普遍解法，从而确立了数学的一个新分支——变分学。

有趣的是，在 1690 年约翰·伯努利的哥哥雅可比·伯努利曾提出著名的悬链线问题(The Hanging Chain Problem)向数学界征求答案，即，固定项链的两端，在重力场中让它自然垂下，问项链的曲线方程是什么在大自然中，除了悬垂的项链外，我們还可以观察到吊桥上方的悬垂钢索，挂着水珠的蜘蛛网，以及两根电线杆之间所架设的电线，这些都是悬链线（catenary） 伽利略（Galileo, 1564～1643）比贝努利更早注意到悬链线，他猜测悬链线是抛物线，从外表看的确象，但实际上不是惠更斯（Huygens, 1629 ～1695）在 1646 年（当时 17 岁），经由物理的论证，得知伽利略的猜测不对，但那时，他也求不出答案到 1691 年，也就是雅可比·伯努利提出悬链线问题的第二年，莱布尼兹、惠更斯（以 62 岁）与约翰·伯努利各自得到了正确答案，所用方法是诞生不久的微积分，具体说是把问题转化为求解一个二阶常微分方程解此方程并适当选取参数，得 （1）)(21axaxey即为悬链线悬链线问题本身和变分法并没有关系，然而这和最速降线问题一样都是贝努利兄弟间的相互争强好胜、不断争吵的导火索，虽然雅可比·贝努利在解决悬链线问题时略占下风，但他随后所证明的“悬挂于两个固定点之间的同一条项链，在所有可能的形状中，以悬链线的重心最低，具有最小势能” ，算是扳回了一局，俩兄弟扯平了！之所以提到悬链线问题，有两方面考虑，其一，这是有关数学史上著名的贝努利家族内的一个趣闻，而这是一个在变分法乃至整个数学物理领域有着巨大贡献的家族，其二，有关悬链线的得几个结论，可0)(122ydxadx2以用变分法来证明！ 现实中很多现象可以表达为泛函极小问题，我们称之为变分问题。

求解方法通常有两种：古典变分法和最优控制论我们这儿要介绍的基本属于古典变分法的范畴1.1 变分法的基本概念1.1.1泛函的概念设 为一函数集合，若对于每一个函数 有一个实数 与之对应，则称 是定S Stx)(JJ义在 上的泛函，记作 称为 的容许函数集)(txJSJ例如，在 上光滑曲线 y(x)的长度可定义为],[10（2）102xdxyJ考虑几个具体曲线，取 ，,10x若 ，则xy)( 102)((dxxJy若 y(x)为悬链线，则  10 10124))2( eedeeJxx对应 中不同的函数 y(x)，有不同曲线长度值 J，即 J 依赖于 y(x)，是定义在],[1C函数集合 上的一个泛函，此时我们可以写成0x)(xyJ我们称如下形式的泛函为最简泛函（3）ft dttFtJ0,)(&被积函数 包含自变量 ，未知函数 (t)及导数 (t)如，上述曲线长度泛函即为一最简Fx泛函1.1.2泛函极值问题考虑上述曲线长度泛函，我们可以提出下面问题：在所有连接定点 的平面曲线中，试求长度最小的曲线),(),(10yxBxA和即，求 ，使100)(,(][)( yxCyy 102)(xdJ取最小值。

此即为泛函极值问题的一个例子以极小值为例，一般的泛函极值问题可表述为，称泛函 在 取得极小值，如果对于任意一个与 接近的 ，)(txJSt0 )(0txStx)(都有 所谓接近，可以用距离 来度量，而距离可以 )(,0txd定义为 |})(|,)({|max)(,( 0000 txttttxdft &3泛函的极大值可以类似地定义其中 称为泛函的极值函数或极值曲线)(0tx1.1.3泛函的变分如同函数的微分是增量的线性主部一样，泛函的变分是泛函增量的线性主部作为泛函的自变量，函数 在 的增量记为)(tx0t )()(0txtx也称函数的变分由它引起的泛函的增量记作 00tJtJ如果 可以表为J )(,(),(00txrtxL其中 为 的线性项，而 是 的高阶项，则称 为泛函在 的变分，记作LxrL用变动的 代替 ，就有 )(0tJ)(t)(0tx)(txJ泛函变分的一个重要形式是它可以表为对参数 的导数：（4）0)(()( tttJ这是因为当变分存在时，增量),(),()()( xtrxtLtxtx 根据 和 的性质有Lr),(),(tt0limli00 xrx所以 )((li)(0JJxJ )(,(,,li0 xJLxrxL1.2 泛函极值的相关结论1.2.1泛函极值的变分表示利用变分的表达式（4） ，可以得到有关泛函极值的重要结论。

泛函极值的变分表示：若 在 达到极值（极大或极小） ，则)(txJ0t（5）证明：对任意给定的 ， 是变量 的函数，该函数在 处达到极值根x)(0xJ0据函数极值的必要条件知40)(0xJ再由（4）式，便可得到（5）式变分法的基本引理： ， ， ，有],[)(21C],[)(21xC0)(21x，210xdx则 ],[ ,0)( 2证明略1.2.2泛函极值的必要条件考虑最简泛函（3） ，其中 F 具有二阶连续偏导数，容许函数类 S 取为满足端点条件为固定端点（6）的二阶可微函数 （6）0)(xtffxt)(泛函极值的必要条件：设泛函（3）在 x(t)∈S 取得极值，则 x(t)满足欧拉方程（7）xxFdt&欧拉方程推导：首先计算（3）式的变分：0)((ttxJ ft dttxtxF0 0)(,, ft xx0 ])()([&对上式右端第二项做分部积分，并利用 ，有(0ftt，ff txtx dFdF00 ),),( &&所以 ftxxttJ0][利用泛函极值的变分表示，得0][0 ftxxdtFt&因为 的任意性，及 ，由基本引理，即得（7） 。

x)(0f（7）式也可写成（8）0xFFxtx &&通常这是关于 x(t)的二阶微分方程，通解中的任意常数由端点条件（6）确定1.2.3几种特殊形式最简泛函的欧拉方程(i) 不依赖于 ，即Fx),(xt5这时 ，欧拉方程为 ，这个方程以隐函数形式给出 ，但它一般0xF&0),(xtF)(tx不满足边界条件，因此，变分问题无解 (ii) 不依赖 ，即 ),(t&欧拉方程为0),(xtFd&将上式积分一次，便得首次积分 ，由此可求出 ，积分后得到可能1,cx ),(1ctx&的极值曲线族 dt1,(iii) 只依赖于 ，即Fx&)(xF这时 ，欧拉方程为0,,0xxtx&由此可设 或 ，如果 ，则得到含有两个参数的直线族 另外xxFx 21ctx若 有一个或几个实根时，则除了上面的直线族外，又得到含有一个参数 的直线族0x&，它包含于上面含有两个参数的直线族 中，于是，在 情ckt 21ctx)(xF&况下，极值曲线必然是直线族iv） 只依赖于 和 ，即Fx&),(F这时有 ，故欧拉方程为0xt&0xx&此方程具有首次积分为 1cFx事实上，注意到 不依赖于 ，于是有t。

0)()(  xxxxxx Fdtdtdt &&&&1. 3 几个经典的例子1.3.1 最速降线问题 最速降线问题 设 和 是铅直平面上不在同一铅直线上的两点，在所有连结 和AB A的平面曲线中，求一曲线，使质点仅受重力作用，初速度为零时，沿此曲线从 滑行至B的时间最短解 将 A 点取为坐标原点， B 点取为 B(x1,y1)，如图 1根据能量守恒定律，质点在曲6线 上任一点处的速度 满足（ 为弧长） A(0, 0) x)(xydtsmgyt21将 代入上式得 B(x1,y1)dxyds)('2y 图 1 最速降线问题xgyt'12于是质点滑行时间应表为 的泛函)(dxgyxyJ202'1端点条件为1)(,)(最速降线满足欧拉方程，因为yyF2'1)',(不含自变量 ，所以方程（8）可写作x0'''yFy等价于 )'('ydx作一次积分得12)'(c令 则方程化为,2'ctgy)cos1(2sin'112 cy又因 dcctgdydx )os1(2osi'1积分之，得 2)sin(2x7由边界条件 ，可知 ，故得0)(y02c).cos1(2inyx这是摆线（园滚线）的参数方程，其中常数 可利用另一边界条件 来确定。

1c1(yx）1.3.2 最小旋转面问题最小旋转面问题 对于 平面上过定点 和 的每一条光滑曲线xy),(1yxA),(2B，绕 轴旋转得一旋转体旋转体的侧面积是曲线 的泛函 ，易得)(xy xyJdxyxxyJ )(')(2)(12容许函数集可表示为 2121)()(][)( y, ,C|S 解 因 不包含 ，故有首次积分 "1yFx122' ''' cyy 化简得 21'c令 ，代入上式， shty' chttscy121由于 dttcdx1'积分之，得 21消去 ，就得到 t1cxhy这是悬链线方程，适当选择条件（令该悬链线过（0，1/a）点，且该点处的切线是水平的）就可得到（1） 本例说明，对于平面上过两个定点的所有光滑曲线，其中绕 轴旋转所得x旋转体的侧面积最小的是悬链线！1.3.3 悬链线势能最小1691 年，雅可比·伯努利证明：悬挂于两个固定点之间的同一条项链，在所有可能的形状中，以悬链线的重心最低，具有最小势能下面我们用变分法证明之考虑通过 A、B 两点的各种等长曲线令曲线 y＝f(x)的长度为 L，重心坐标为 ，),(yx8则 babadxydsL2)(1由重心公式有， Lxdyxba2)(1Lxdyyba2)(由于只需探讨曲线重心的高低，所以只对纵坐标的公式进行分析，注意到问题的表述，说明 L 是常数，不难看出重心的纵坐标是 y(x)的最简泛函，记作badxyxyJ2)(1)(此时对应的欧拉方程（8）可化为 0)(2令 解得 ，进而得dxyp0),1(22kp。

)]([1cxkhy此即为悬链线，它使重心最低，势能最小！大自然中的许多结构是符合。