我的数模讲义.doc

1、最优控制建模现实生活中很多现象可以表达为求泛函的极大与极小值问题，我们称之为泛函极值问题。其求解方法通常有三种：变分法(古典变分法,近代变分法)，动态规划方法，最大值原理- 它们都属于最优控制的范畴。关于最优控制的近似计算方法用得较多的有：无约束最优控制问题的梯度方法（或者最速下降方法） ，Newton 方法,以及有约束最优控制问题的罚函数方法.读者可参考有关文献.另外，需要指出来的是最优控制的近似计算方法远较有限维非线性规划问题的计算方法复杂。迄今为止，最优控制的计算方法的研究，在深度和广度方面都远不如非线性规划的最优化计算方法的研究，一些最优控制的计算方法，如最优控制的罚函数方法，其收敛速

2、度等问题都还未解决. 我在这主要要介绍的是古典变分法方法（不准备介绍近代变分方法）.古典变分法简介作为数学的一个分支，变分法的诞生，是现实世界许多现象不断探索的结果，人们可以追寻到这样一个轨迹：约翰伯努利（Johann Bernoulli，16671748）1696 年向全欧洲数学家挑战，提出一个难题：“设在垂直平面内有任意两点，一个质点受地心引力的作用，自较高点下滑至较低点，不计摩擦，问沿着什么曲线下滑，时间最短？”这就是著名的“最速降线”问题（The Brachistochrone Problem） 。它的难处在于和普通的极大极小值求法不同，它是要求出一个未知函数（曲线） ，来满足所给的条

3、件。这问题的新颖和别出心裁引起了很大兴趣，罗比塔（Guillaume Francois Antonie de lHospital 1661-1704） 、雅可比伯努利（Jacob Bernoulli 1654-1705） 、莱布尼茨（Gottfried Wilhelm Leibniz,1646-1716）和牛顿（Isaac Newton16421727）都得到了解答。约翰的解法比较漂亮，而雅可布的解法虽然麻烦与费劲，却更一般化。后来欧拉（Euler Lonhard，17071783）和拉格朗日(Lagrange, Joseph Louis，1736-1813)发明了这一类问题的普遍解法，从而确

4、立了数学的新分支-变分学。有趣的是，在 1690 年约翰伯努利的哥哥雅可比伯努利曾提出著名的悬链线问题(The Hanging Chain Problem)向数学界征求答案，即，固定项链的两端，在重力场中让它自然垂下，问项链的曲线方程是什么。在大自然中，除了悬垂的项链外，我们还可以观察到吊桥上方的悬垂钢索，挂着水珠的蜘蛛网，以及两根电线杆之间所架设的电线，这些都是悬链线（catenary） 。 伽利略（Galileo, 15641643）比贝努利更早注意到悬链线，他猜测悬链线是抛物线，从外表看的确象，但实际上不是。惠更斯（Huygens, 16291695）在 1646 年（当时 17 岁）

5、，经由物理的论证，得知伽利略的猜测不对，但那时，他也求不出答案。到 1691 年，也就是雅可比伯努利提出悬链线问题的第二年，莱布尼兹、惠更斯（以 62 岁）与约翰伯努利各自得到了正确答案，所用方法是诞生不久的微积分，具体说是把问题转化为求解一个二阶常微分方程解此方程并适当选取参数，得 （1）)(21axaxey即为悬链线。悬链线问题本身和变分法并没有关系，然而这和最速降线问题一样都是贝努利兄弟间的相互争强好胜、不断争吵的导火索，虽然雅可比贝努利在解决悬链线问题时略占下风，但他随后所证明的“悬挂于两个固定点之间的同一条项链，在所有可能的形状中，以悬链线的重心最低，具有最小势能” ，算是扳回了一局

6、，俩兄弟扯平了！之所以提到悬链线问题，有两方面考虑，其一，这是有关数学史上著名的贝努利家族内的一个趣闻，而这是一个在变分法乃至整个数学物理领域有着巨大贡献的家族，其二，有关悬链线的得几个结论，可以用变分法来证明！ 1.1 变分法的基本概念1.1.1 泛函的概念设 为一函数集合，若对于每一个函数 有一个实数 与之S Stx)(J对应，则称 是定义在 上的泛函，记作 。 称为 的容许函JSJ数集。0)(122ydxadx例如，在 上光滑曲线 y(x)的长度可定义为,10x（2）02xdxyJ考虑几个具体曲线，取 ，1,若 ，则xy)( 02)(dxxJy若 y(x)为悬链线，则 10 10124)

7、2( eedeeJ xxx对应 中不同的函数 y(x)，有不同曲线长度值 J，即 J 依,1C赖于 y(x)，是定义在函数集合 上的一个泛函，我们可以写成,10xC)(yJ我们称如下形式的泛函为最简泛函（3）ft dtxtFtxJ0)(,()(被积函数 包含自变量 ，未知函数 (t)及导数 (t)。上述曲线长度x泛函即为一最简泛函。1.1.2 泛函极值问题考虑上述曲线长度泛函，我们可以提出下面问题：在所有连接定点 的平面曲线中，试求长度最小的曲),(),(10yxBxA和线。即，求 ，使100)(,)(,) yxCy 102(xdyJ取最小值。此即为泛函极值问题的一个例子。以极小值为例，一般的

8、泛函极值问题可表述为，称泛函 在 取得极小值，如果对于任意一个与)(txJSt)(0接近的 ，都有 。所谓接近，可以用距离)(0tx )(0txJt来度量，而距离可以定义为,d|)(|,)(|max)(,( 0000 txttttxdft 泛函的极大值可以类似地定义。其中 称为泛函的极值函数或极)0值曲线。1.1.3 泛函的变分如同函数的微分是增量的线性主部一样，泛函的变分是泛函增量的线性主部。作为泛函的自变量，函数 在 的增量记为)(tx0t)(0tx也称函数的变分。由它引起的泛函的增量记作 )()(00txJtJ如果 可以表为J)(,(),(00trtxL其中 为 的线性项，而 是 的高阶

9、项，则称 为泛函在 的变LxrL)(0tx分，记作 。用变动的 代替 ，就有 。)(0tJ)(tx)(0t)(tJ泛函变分的一个重要形式是它可以表为对参数 的导数： （4）0)()( txtJtxJ这是因为当变分存在时，增量),(),()()( xtrxtLtxJtxJ 根据 和 的性质有Lr),(),(xtxt0limli00 rr所以)(lim)(0xJxJxJ )(,(),(),(li0 xJLxrxL 1.2 泛函极值的相关结论1.2.1 泛函极值的变分表示利用变分的表达式（4） ，可以得到有关泛函极值的重要结论。泛函极值的变分表示：若 在 达到极值（极大或极小） ，)(txJ)(0t

10、则 （5）0)(txJ证明：对任意给定的 ， 是变量 的函数，该函数在x)(0xJ处达到极值。根据函数极值的必要条件知0)(0xJ再由（4）式，便可得到（5）式。变分法的基本引理：若 ，对于 ，,)(21xC,)(21xC，有 ，则 。0)(21x210xd, ,0)( 21x证明略。1.2.2 泛函极值的必要条件考虑最简泛函（3） ，其中 F 具有二阶连续偏导数，容许函数类S 取为满足端点条件为固定端点（6）的二阶可微函数。， （6）0)(xtffxt)(泛函极值的必要条件：设泛函（3）在 x(t)S 取得极值，则 x(t)满足欧拉方程（7）0xxFdt欧拉方程推导：首先计算（3）式的变分：

11、0)(txtJ ft dttxtF0 0)(, ft xx0 )()(对上式右端第二项做分布积分，并利用 ，有)(0ftxt，ff txtx dFdF00 ,(),(所以 ftxxttJ0利用泛函极值的变分表示，得00 ftxxdtFt因为 的任意性，及 ，由基本引理，即得（7） 。x)(0f（7）式也可写成（8）xFFxtx 通常这是关于 x(t)的二阶微分方程，通解中的任意常数由端点条件（6）确定。1.2.3 几种特殊形式最简泛函的欧拉方程(i) 不依赖于 ，即Fx),(xtF这时 ，欧拉方程为 ，这个方程以隐函数形式给出0x 0x，但它一般不满足边界条件，因此，变分问题无解。 )(tx(

12、ii) 不依赖 ，即Fx),(xtF欧拉方程为0),(xtFd将上式积分一次，便得首次积分 ，由此可求出 ，1),(cxt ),(1ctx积分后得到可能的极值曲线族 dtcx1,(iii) 只依赖于 ，即F)(xF这时 ，欧拉方程为0,0xxtx由此可设 或 ，如果 ，则得到含有两个参数的直线族xxFx。另外若 有一个或几个实根时，则除了上面的直线族21ct0外，又得到含有一个参数 的直线族 ，它包含于上面含有两ccktx个参数的直线族 中，于是，在 情况下，极值曲线21tx)(xF必然是直线族。（iv） 只依赖于 和 ，即Fx),(x这时有 ，故欧拉方程为0xt 0xxF此方程具有首次积分为

13、 1cx事实上，注意到 不依赖于 ，于是有Ft。0)()( xxxxxx FdtdtFdt 1. 3 几个经典的例子1.3.1 最速降线问题 最速降线问题 设 和 是铅直平面上不在同一铅直线上的两AB点，在所有连结 和 的平面曲线中，求一曲线，使质点仅受重力AB作用，初速度为零时，沿此曲线从 滑行至 的时间最短。AB解 将 A 点取为坐标原点，B 点取为 B(x1,y1)，如图 1。根据能量守恒定律，质点在曲线 上任一点处的速度 满足（ 为弧长）)(xydtssA(0, 0) xy B(x1,y1)图 1 最速降线问题mgydts21将 代入上式得 于是质点滑行时间应表xs)(2dxgydt2

14、1为 的泛函 )(xyyyJx20)(端点条件为 最速降线满足欧拉方程，1,由于 不含自变量 ，所以方程（ 8）可写作yyF21),(x0yF等价于 )(ydx作一次积分得12)(cy令 则方程化为,2ctgy)cos1(2sin112 cy又因 dcctgdydx )os1(2os2in1积分之，得 2)sin(2cx由边界条件 ，可知 ，故得0)(y0).cos1(2inyx这是摆线（园滚线）的参数方程，其中常数 可利用另一边界条件1c来确定。1(yx）1.3.2 最小旋转面问题最小旋转面问题 对于 平面上过定点 和 的每xy),(1yxA),(2yxB一条光滑曲线 ，绕 轴旋转得一旋转体。旋转体的侧面积是曲线)(xy的泛函 ，易得)(xyJ dxyxxyJ )(1)(2)(12容许函数集可表示为 2121)()()( yx, ,yxC|yxS 解 因 不包含 ，故有首次积分 “1F