变分原理与变分法.docx

第一章变分原理与变分法1.1关于变分原理与变分法（物质世界存在的基本守恒法则一、大自然总是以可能最好的方式安排一切，似乎存在着各种安排原理：昼/夜，日/月，阴/阳，静止/运动 等矛盾/统一的协调体；对静止事物：平衡体的最小能量原理，对称/相似原理；对运动事物：能量守恒，动量（矩）守恒，熵增原理等。变分原理是自然界静止（相对稳定状态）事物中的一 普遍适应的数学定律， 获称最小作用原理。Examples： 光线最短路径传播； 光线入射角等于反射角，光线在反射中也是光传播最短路径（Heron）； 光线折射遵循时间最短的途径（Fermat）；E CAE + EB > AC + CBSummary：实际上光的传播遵循最小能量原理；在静力学中的稳定平衡本质上是势能最小的原理。二、变分法是自然界变分原理的数学规划方法（求解约束方程系统极值的数学方法），是计算泛函驻值的数学理论数学上的泛函定义 定义：数学空间（集合）上的元素（定义域）与一个实数域间（值域）间的(映射)关系特征描述法： J: X u D R I J(x) = r e R Examples： 矩阵范数：线性算子（矩阵）空间 =数域IIAH1=寸max二ji=1aijAL=max £ a8 i j=1iaii =必j=1 i=1函数的积分：函数空间数域J = b f (x)dxf u Da nnNote:泛函的自变量是集合中的元素(定义域)；值域是实数域。Discussion：判定下列那些是泛函：|f | = max|f(x)| ;步)；3x+5y=2 ; f+w5 (x - x0) f (x )dx = f (x0)a< x<b°x-3x = 0,固支；x = 1,自由试举另一泛函例子。物理问题中的泛函举例弹性地基梁的系统势能i. 梁的弯曲应变能：n = -jlEJ(d-2w)2dxb 2 0 dx2ii. 弹性地基贮存的能量： n = - j lkw 2dxf 2 0iii. 外力位能：n i = - j1 qwdxiv. 系统总的势能：i d2wdwn = Jl 4-EJ( )2 + 4-kw2 - qwdx; x = 0 w = 0 = 0泛函的提法：有一种梁的挠度函数(与载荷无关)，就会有一个对应的系 统势能。泛函驻值提法：在满足位移边界条件的所有挠度函数中，找一个 w(x)，使 系统势能泛函取最小值。最速降线问题问题：已知空间两点A和B,A高于B，要求在两点间连接一条曲线，使得 有重物从A沿此曲线自由下滑时，从A到B所需时间最短(忽略摩擦力)。 作法：i. 通过A和B作一垂直于水平面的平面，取坐标系如图。B点坐标(& b), 设曲线为 y = y(x),并已知：x = 0,y = 0； x = a,y = bii. 建立泛函：设P(x , y)9曲线上的点，P点的速度由能量守恒定律求得：1 mv 2 = mgy n v =、：2 gy2 V命ds为曲线弧长的微分，有：ds 二 _ 土_ ds _ v 1 + y'2v , - 2 gy n dt - dxdtv2 gy v2 gy重物从A点滑到B点的总时间：T= j 八'i y '2 dx0 V2 gy泛函驻值提法：在0<x<a的区间内找一个函数y(x)使其满足端点几何条 件并使T取最小值。圆周问题问题：在长度一定的闭曲线中，什么曲线所围成的面积最大。作法：i. 假设所考虑的曲线用参数形式表示:x = x(s),y = y(s)s为参数。取s为曲线上的某一定点，则坐标表示x=x(s ),y =y(s ),因 '1111y(s )与点s (x y )重21曲线是封闭的，必存在一个s点使乂2= x(s), y =1, 1合。ii.该封闭曲线的周长:L = js2 档)2 +(半)2dss1 ds ds该曲线所围成的面积：R = jj dxdy Qiii.转换R的表达式由Green公式：jj (理-竺)dxdy j 七 Pdx + QdydxdysQi取 P=" y，号，则："=1-ydx1js2 (xy'(s) yx'(s)ds 2s1泛函驻值的提法：等周问题即是在满足端点条件x (s )y(s )2=x(s ),2y(s )1及周长一定% J(d)2 + 3)2 = L 条件下，寻找一个曲线函数J X(S)使泛七 ds ds y (s)函R取驻值。 Discussion悬索线问题：已知空间中A, B两点及一条长度L> ab的悬索，单位长的质 量为m。假设绳索的长度是不变的，并忽略绳索的弯曲刚度，把此绳索的 两端挂在A，B两点，求在平衡状态下绳索的形状。要求：列出悬索线应满足的泛函式及泛函驻值提法。提示：绳索在平衡状态下，其势能应为最小值。1.2变分法(泛函驻值的计算方法).关于计算固体力学中的泛函、泛函极值的提法 这里所研究的泛函一般用积分显式表达，并不等于所有泛函都能用显 式积分表达。 所要研究的泛函都可表示成在一定区间或一定区域内的函数及其导数 (或偏导数)的积分形式，即：a. n =F(f (x), f (x), f (x);x)dxb. n 2 = " F(f (x, y), fx (x, y), fy (x, y); x, y)dxdyQc. 泛函中的可变化函数称为自变函数，或称宗量(argument)，x 或 y 仅是积分变量，是被积函数的定义域。(被积函数是复合函数概念的 推广) 要说清楚一个泛函的极值问题，应注意：a. 应把泛函本身讲清楚(即写出它的形式)；b. 还必须讲明白自变函数的性质，如：- 独立的自变函数的个数(导函数并不独立)-每个自变函数定义的区间/区域- 这些自变函数应满足的条件(如：边界条件及其受约束的条件等)c. 除了个别特殊情况外，一般情况下增加一个条件会使泛函极值及相应 的自变函数变化性质发生变化。如：极小值可能变大；极大值可能变 小；非极值的驻值可能成为极值。.若干背景知识 泛函的驻值问题可以转化为等价的微分方程问题，变分法的理论计算就是完成这类工作。本章内容沿袭此方法，是要把问题的理论基础讲明确。 从近似解的角度出发，直接求解泛函的驻值，比解微分方程更加方便， 也更为实用。特别计算机技术的发展，带来了大规模数值计算的可能性(有限元的思想基础)。 经 Euler, Lagrange, Dirichlet Hilbert, Bernoulli 等数学先驱的卓越 工作，完成了的系统方法。 但把微分方程问题转换为泛函问题还很不成熟。在物理、力学中，即先 猜想一个泛函的驻值问题，再校对是否与原微分方程问题等价。 泛函驻值的计算(数值)先驱工作中以Ritz, Galerkin Treft著名。关于变分法的一个预备定理若f (x)在a, b上连续，若对任意满足(a)=(b)=0的连续函数x都有：j f (x)甲(x)dx = 0 a则 f (x)在a, b上处处为零。反证法：设x0为a,b中的点，在x0点f(x°)手0,可取f(x0)>0,.f (x)在区间上连续，必存在x的一个充分小邻域上f (x)>0, x - <x<x +x 为任意连续函数(满足边界条件),可取000又x也在该邻域内大于零，而在该邻域夕卜恒等于零o所以有j b f (x )甲(x) dx > 0a矛盾！即f (x)必须为零；同理可证小于零情况。该定理可推广多元变量的函数问题。1.2.1定积分bF3, y, y')dx的驻值(变分)问题 a目的：通过简单泛函的极值分析，获得建立变分法的基本概念、计算步骤(把 变分解转化成微分方程)问题：在自变量x的区间a，b 内决定一个函数y(x)，使它满足边界条件：y = a |, y = & I b 并使泛函：V = jb F(x, y, y ')dx 取极值。 a计算5 V方法1 ：先用变分观点解释G.H曲线的增量 设想已取得了一条曲线GACH方程为：y=y (x) 在GACH附近另取一条曲线GBDH，令该曲线无限接近GACH，其方程为：y1(x) = y (x) +5 y (x) 5y(x)是一个无穷小量，称为自变函数的变分(若 x 不变，即为曲线纵坐标 的增量)(注意与函数微分的区别，这里函数的变分仍然是一个函数) 相应两条曲线，获得两个泛函值：V = b F (x, y, y ')dxaV + AV = fb F(x, y + 5y, yr + 5y')dxa 基本引理：(5y) =5y'证： 5 y(x) = y (x) - y(x) n (5 y) = y'(x) - yf(x)三 5 y'11推广：(5y )" = 5y"另一条认识(5 y)r = d矿 的思路:A C :y (x ) = y (x ) + y 'dxA T B :y1(xB) = y(xA) +SyAC T D :y1(xD) = y(xC) + SyCB T D :y (x ) = y (x ) + y'dxy1=y+5 y；5y1 =矿产n (5y)'=5y -Sy1dx=y:-y'=5y'1 D 1 B 1AV = jb F (x, y + 5y, y' + 5y' ) F (x, y, y，)dxa因为F(x, y, yf)是x, y, y'的连续可导函数(工程上一般如此)，故Sy及Sy，很小时，AV也很小，即 Sy, Sy't 0 AV T 0取等式两端的一阶无穷小量，即：SV=j:尝"京y皿(可以从Tailor展开式去理解). S V称为泛函V的一阶变分，简称变分，即泛函的一阶变分是泛函增量中的 一阶小量部分(把自变函数的变分Sy作为一阶小量)所以，变分的运算服从 无穷小量的运算规则。计算SV方法2:(把求泛函的极值转化成求普通函数的极值)记：y (x) = y (x) + 宓 y(x)0 <£ < 1 ( y 及S y 固定)V 伍)=j :f (x, y0+eS y, y 0+eS y')dx当V 在y0上取极值，则相应于£ = 0的泛函值V(8)现在成为普通的函数 极值条件：V'(£ ) 1£ 0 = 0(先不管该条件，现仅研究其导数计算)SV =也 d£=jb竺土空 + 打(y0+£明dxd£=jb 竺 Sy +竺 Sy，dx d£a dyd£dy'd£a dydy'. 上两式中出现，Sy和Sy'并不能独立变化，可设法把Sy'项转换成只与S y有关的项。取分步积分：卜竺 5 ydx = -b a dy'u'vdx + uv |b取:aadFv =5 y8y'Jb 竺5y 'dx = -b± (当5ydx + 竺5y |b a Qy，a dx dy，dy， a代入一阶变分式:5V = J b 竺-d (dF )5ydx +dF 5y I a dy dx dy 'dy '要选定的函数满足边界条件，所以:5 y = 0 I ，5 y = 0 In 5 V = Jb 竺-(dF)5ydx a dydx dy'计算5 V = 0 若方括号内的函数在区间内不为 0，则可任选5y使5 V大于零或小于零， 使V不能获得极值，故需方括号的项为零。即：竺 一 （dF） = 0（Euler 方程）dy dx dy'此即与泛函驻值等价的微分方程。或：令5 V = 0由变分基本定理：5y任意连续函数，方括号中函数连续。n竺-d （竺）=0dy dx dy'1 + yf 2Example 最速降线问题：F =，一（注不显含x）.2 gy代入Euler方程，并乘以函数Q可得:dF d dF dF dF d dFQ - - Q (-) = Q - + Q - -(Q -) = 0dy dx dyf dy dyf dx dy!由于竺 =0 (F中不显含x)，上式中只要令Q = y,把上式酉己成全微分形式: dxd (F 卫 y') = 0dx8y，这是因为：d 日 _ °FdFdydFdy'8F_ 口dx dxdydxdy'dx( dx )(代回原Euler方程，即得全微分)一一 dFdF后两项由Q的假设= -fQ + = Q' dydy由全微分方程n F -y y，= C代入F的具体表达式:n 1= C n y(1 + y'2)=t'y (1 + y'2)令：y ' = c tan ty =1 + ctan21. 一， v=v sin21 = 2(1 - cos2t)、,dy2v sin t cos tn dx =dt = 2v sin2 tdt = v(1 - cos 2t)dtyfc tan tv上式积分得：x = (2t - sin 2t) + C 21v注意：y = (1 - cos2t)2引用初始条件：x=0,y=0, 只能有：C=0v令：目= ,2t =0 2x =四(0 - sin 0)八即为最速降线(圆滚线(渐开)方程。y =日(1 一 cos 0)Homework：在连接XY平面上两点的M0及M1的所有曲线中，要这样一条曲线使它 绕OX轴旋转成的曲面有最小表面积。