第一章 变分原理与变分法1.1 关于变分原理与变分法（物质世界存在的基本守恒法则）一、 大自然总是以可能最好的方式安排一切，似乎存在着各种安排原理：昼/夜，日/月，阴/阳，静止/运动 等矛盾/统一的协调体； 对静止事物：平衡体的最小能量原理，对称/相似原理；对运动事物：能量守恒，动量（矩）守恒，熵增原理等。

变分原理是自然界静止(相对稳定状态)事物中的一个普遍适应的数学定律，获称最小作用原理。

Examples :① 光线最短路径传播；② 光线入射角等于反射角，光线在反射中也是光传播最短路径（Heron ）； ③CB AC EB AE +>+Summary : 实际上光的传播遵循最小能量原理；在静力学中的稳定平衡本质上是势能最小的原理。

二、变分法是自然界变分原理的数学规划方法（求解约束方程系统极值的数学方法），是计算泛函驻值的数学理论数学上的泛函定义定义：数学空间（集合）上的元素（定义域）与一个实数域间（值域）间的（映射）关系特征描述法：{ J :R x R D X ∈=→⊂r J )(|}Examples :① 矩阵范数：线性算子（矩阵）空间数域‖A ‖1 = ∑=ni ij ja 1max ；∑=∞=nj ij ia A 1max；21)(1122∑∑===nj ni ij a A② 函数的积分： 函数空间数域D ⊂=⎰n ban f dxx f J )(Note : 泛函的自变量是集合中的元素（定义域）；值域是实数域。

Discussion ：① 判定下列那些是泛函：)(max x f f b x a <<=；x y x f ∂∂),(； 3x+5y=2； ⎰+∞∞-=-)()()(00x f dx x f x x δ ② 试举另一泛函例子。

物理问题中的泛函举例① 弹性地基梁的系统势能i. 梁的弯曲应变能： ⎰=∏l b dx dxw d EJ 0222)(21ii. 弹性地基贮存的能量： dx kw lf⎰=∏0221 iii. 外力位能： ⎰-=∏l l qwdx 0iv. 系统总的势能：000;})({221222021===-+=∏⎰dxdww x dx qw kw dxw d EJ l 泛函的提法：有一种梁的挠度函数（与载荷无关），就会有一个对应的系统势能。

泛函驻值提法：在满足位移边界条件的所有挠度函数中，找一个w (x ),使系统势能泛函取最小值。

② 最速降线问题问题：已知空间两点A 和B ,A 高于B ，要求在两点间连接一条曲线，使得有重物从A 沿此曲线自由下滑时，从A 到B 所需时间最短（忽略摩擦力）。

作法：i. 通过A 和B 作一垂直于水平面的平面，取坐标系如图。

B 点坐标(a , b ),设曲线为y = y (x ),并已知：x = 0,y = 0；x = a ,y = b ii. 建立泛函：x设P (x , y )是曲线上的点，P 点的速度由能量守恒定律求得：gy v mgy mv 2221=⇒= 命ds 为曲线弧长的微分，有：dx gy y gy dsdt gy v dt ds 2'1222+==⇒== 重物从A 点滑到B 点的总时间：T =dx gyy a⎰+022'1泛函驻值提法：在0≤x ≤a 的区间内找一个函数y (x )使其满足端点几何条件并使T 取最小值。

③ 圆周问题问题：在长度一定的闭曲线中，什么曲线所围成的面积最大。

作法：i. 假设所考虑的曲线用参数形式表示： x = x (s ), y = y (s )s 为参数。

取s 1为曲线上的某一定点，则坐标表示x 1=x (s 1),y 1=y (s 1),因曲线是封闭的，必存在一个s 2点使x 2 = x (s 2),y 2 = y (s 2)与点s 1(x 1,y 1)重合。

ii. 该封闭曲线的周长： L =ds dsdy ds dx s s⎰+2122)()(该曲线所围成的面积：R = ⎰⎰Ωdxdyiii. 转换R 的表达式 由Green 公式：⎰⎰⎰+=∂∂-∂∂21)(s s Qdy Pdx dxdy yP x Q Ω取P =-2y ,Q =2x, 则：1=∂∂-∂∂yP x Q ∴⎰⎰-=-=2121))(')('(2121sss sds s yx s xy ydx xdy R泛函驻值的提法：等周问题即是在满足端点条件x (s 1) = x (s 2), y (s 1) = y (s 2)及周长一定L s s dsdy ds dx =+⎰2122)()( 条件下，寻找一个曲线函数⎩⎨⎧)()(s y s x 使泛函R 取驻值。

④ Discussion悬索线问题：已知空间中A ，B 两点及一条长度L >AB 的悬索，单位长的质量为m 。

假设绳索的长度是不变的，并忽略绳索的弯曲刚度，把此绳索的两端挂在A ，B 两点，求在平衡状态下绳索的形状。

要求：列出悬索线应满足的泛函式及泛函驻值提法。

提示：绳索在平衡状态下，其势能应为最小值。

1.2 变分法（泛函驻值的计算方法）● 关于计算固体力学中的泛函、泛函极值的提法① 这里所研究的泛函一般用积分显式表达，并不等于所有泛函都能用显式积分表达。

② 所要研究的泛函都可表示成在一定区间或一定区域内的函数及其导数（或偏导数）的积分形式，即： a. ⎰=ba dx x x f x f x f F ));("),('),((1∏b. ⎰⎰=Ω∏dxdy y x y x f y x f y x f F y x ),);,(),,(),,((2c. 泛函中的可变化函数称为自变函数，或称宗量（argument ），x 或y 仅是积分变量，是被积函数的定义域。

（被积函数是复合函数概念的推广）③ 要说清楚一个泛函的极值问题，应注意： a. 应把泛函本身讲清楚（即写出它的形式）； b. 还必须讲明白自变函数的性质，如：- 独立的自变函数的个数（导函数并不独立）； - 每个自变函数定义的区间/区域；- 这些自变函数应满足的条件（如：边界条件及其受约束的条件等）。

c. 除了个别特殊情况外，一般情况下增加一个条件会使泛函极值及相应的自变函数变化性质发生变化。

如：极小值可能变大；极大值可能变小；非极值的驻值可能成为极值。

●若干背景知识① 泛函的驻值问题可以转化为等价的微分方程问题，变分法的理论计算就是完成这类工作。

本章内容沿袭此方法，是要把问题的理论基础讲明确。

② 从近似解的角度出发，直接求解泛函的驻值，比解微分方程更加方便，也更为实用。

特别计算机技术的发展，带来了大规模数值计算的可能性（有限元的思想基础）。

③ 经Euler ,Lagrange ,Dirichlet ,Hilbert ,Bernoulli 等数学先驱的卓越工作，完成了①的系统方法。

④ 但把微分方程问题转换为泛函问题还很不成熟。

在物理、力学中，即先猜想一个泛函的驻值问题，再校对是否与原微分方程问题等价。

⑤ 泛函驻值的计算（数值）先驱工作中以Ritz ,Galerkin ,Treft 著名。

关于变分法的一个预备定理 若f (x )在[a ,b ]上连续，若对任意满足 ϕ(a )= ϕ(b )=0 的连续函数ϕ(x ),都有：⎰=ba dx x x f 0)()(ϕ则 f (x )在[a ,b ]上处处为零。

反证法：设x 0为[a ,b ]中的点，在x 0点f (x 0)≠0，可取f (x 0)>0,∵ f (x )在区间上连续，必存在x 0的一个充分小邻域上f (x )>0, x 0-ε<x<x 0+ε 又∵ ϕ(x )为任意连续函数（满足边界条件），可取ϕ(x )也在该邻域内大于零，而在该邻域外恒等于零。

所以有⎰>badx x x f 0)()(ϕ矛盾！即)(x f 必须为零；同理可证小于零情况。

该定理可推广多元变量的函数问题。

1.2.1 定积分⎰'ba dx y y x F ),,(的驻值（变分）问题目的：通过简单泛函的极值分析，获得建立变分法的基本概念、 计算步骤（把变分解转化成微分方程）问题：在自变量x 的区间[ a ，b ]内决定一个函数y (x )，使它满足边界条件：a x y ==|α ，b x y ==|β 并使泛函：⎰'=badx y y x F V ),,( 取极值。

计算V δ方法1：先用变分观点解释G .H 曲线的增量x● 设想已取得了一条曲线GACH 方程为：y= y (x )●在GACH 附近另取一条曲线GBDH ，令该曲线无限接近GACH ，其方程为：)()()(1x y x y x y δ+=●)(x y δ是一个无穷小量，称为自变函数的变分（若x 不变，即为曲线纵坐标的增量）（注意与函数微分的区别，这里函数的变分仍然是一个函数） ●相应两条曲线，获得两个泛函值：dx y y x F V ba ⎰'=),,(⎰'+'+=∆+badx y y y y x F V V ),,(δδ●基本引理： y y '='δδ)(证： y x y x y y x y x y x y '≡'-'='⇒-=δδδ)()()()()()(11 推广： y y ''=''δδ)(另一条认识y y '='δδ)(的思路：C A →： dx y x y x y A C '+=)()( B A →： A A B y x y x y δ+=)()(1D C →： C C D y x y x y δ+=)()(1D B →： dx y x y x y B D 111)()('+= y y y δ+=1 dx y y dx y y 11;'='=⇔δδ )('⇒y δ＝y y y dxyy '='-'=-δδδ11 ⎰'-'+'+=∆badx y y x F y y y y x F V )},,(),,({δδ因为),,(y y x F '是y y x ',,的连续可导函数（工程上一般如此），故y y 'δδ及很小时，V ∆也很小，即 0,→'y y δδ 0→∆V 取等式两端的一阶无穷小量，即：dx y yFy y F V ba][''∂∂+∂∂=⎰δδδ （可以从Tailor 展开式去理解）●V δ称为泛函V 的一阶变分，简称变分，即泛函的一阶变分是泛函增量中的一阶小量部分（把自变函数的变分y δ作为一阶小量）所以，变分的运算服从无穷小量的运算规则。