变分基础学习知识原理与变分法

1、第一章变分原理与变分法1.1关于变分原理与变分法（物质世界存在的基本守恒法则）一、大自然总是以可能最好的方式安排一切，似乎存在着各种安排原理：昼/夜，日/月，阴/阳，静止/运动 等矛盾/统一的协调体；对静止事物：平衡体的最小能量原理，对称 /相似原理；对运动事物：能量守恒，动量（矩）守恒，熵增原理等。变分原理是自然界静止（相对稳定状态）事物中的一个普遍适应的数学定律， 获称最小作用原理。Examples: 光线最短路径传播； 光线入射角等于反射角，光线在反射中也是光传播最短路径（Heron）； 光线折射遵循时间最短的途径（Fermat）；,E CAE EB AC CBSummary实际上光的传

2、播遵循最小能量原理；在静力学中的稳定平衡本质上是势能最小的原理。、变分法是自然界变分原理的数学规划方法 （求解约束方程系统极值的数学方 法），是计算泛函驻值的数学理论数学上的泛函定义定义：数学空间（集合）上的元素（定义域）与一个实数域间（值域）间 的（映射）关系特征描述法： J: X D R|J（x） r RExamples: 矩阵范数：线性算子（矩阵）空间 = 数域nII A II 1 = max aijj i 1max aiji j 1；IIAaij12 函数的积分：函数空间数域bJ fn(X)dXaNote:泛函的自变量是集合中的元素(定义域);值域是实数域。Discussi on： 判

3、定下列那些是泛函：f (x,y) ；3x+5y=2;f maxf(X)；(x x°)f(x)dx f(x°)试举另一泛函例子。物理问题中的泛函举例 弹性地基梁的系统势能i.梁的弯曲应变能:x'EJ (d-w)2dx0 dx2ii.弹性地基贮存的能量:iii.外力位能:I 2kw dx2 0io qwdxiv.系统总的势能:0得ej(器)20 dxqwdx; x 0也0dx泛函的提法：有一种梁的挠度函数(与载荷无关)，就会有一个对应的系 统势能。泛函驻值提法：在满足位移边界条件的所有挠度函数中，找一个 w(x),使 系统势能泛函取最小值。 最速降线问题问题：已知空间两

4、点A和B, A高于B，要求在两点间连接一条曲线，使 得有重物从A沿此曲线自由下滑时，从 A到B所需时间最短(忽略摩擦 力)。作法：i. 通过A和B作一垂直于水平面的平面，取坐标系如图。B点坐标(a, b),设曲线为 y = y(x),并已知：x = 0, y = 0 ； x = a, y = bii. 建立泛函：设P(x , y)是曲线上的点，P点的速度由能量守恒定律求得:2mv2 mgy v . 2gy命ds为曲线弧长的微分，有：ds 2dsds 1 y'.v ,.2gy dt dx重物从A点滑到B点的总时间:aT= 0 X 2 gyy'2dxdt2gy .,2gy泛函驻值提

5、法：在Owx< a的区间内找一个函数y(x)使其满足端点几何条 件并使T取最小值。 圆周冋题问题：在长度一定的闭曲线中，什么曲线所围成的面积最大。作法：i.假设所考虑的曲线用参数形式表示：x = x(s), y = y(s)s为参数。取si为曲线上的某一定点，则坐标表示xi=x(si), yi=y(si),因曲 线是圭寸闭的，必存在一个 s点使X2 = x(s2), y2 = y(s2)与点si(xi,yi)重合。ii. 该圭寸闭曲线的周长：L = s (彩2 (学)该曲线所围成的面积：R = dxdyiii. 转换R的表达式由Green公式：QP2()dxdy Pdx Qdyxys取P

6、 =-于，Q =扌，贝U: i22x y1 32i 32 R 2 xdy ydx i(xy' (s) yx'(s)ds2 s2 s泛函驻值的提法：等周问题即是在满足端点条件x(si) = x(s2),y(si) = y(s2)及周长一定勺，（ds）2 （ds）2l条件下，寻找一个曲线函数x（s）使泛sy（s）函R取驻值。 Discussion悬索线问题：已知空间中 A，B两点及一条长度L>AB的悬索，单位长的 质量为m。假设绳索的长度是不变的，并忽略绳索的弯曲刚度，把此绳索 的两端挂在A，B两点，求在平衡状态下绳索的形状。要求：列出悬索线应满足的泛函式及泛函驻值提法。提示

7、：绳索在平衡状态下，其势能应为最小值。1.2变分法（泛函驻值的计算方法）关于计算固体力学中的泛函、泛函极值的提法 这里所研究的泛函一般用积分显式表达，并不等于所有泛函都能用显 式积分表达。 所要研究的泛函都可表示成在一定区间或一定区域内的函数及其导数 （或偏导数）的积分形式，即：ba. i F（f（x）, f'（x）, f"（x）;x）dxab. 2 F（f （x,y）, fx（x,y）, fy（x,y）;x, y）dxdyc. 泛函中的可变化函数称为自变函数，或称宗量（argumen），x或y仅是积分变量，是被积函数的定义域。（被积函数是复合函数概念的推 广） 要说清楚一个

8、泛函的极值问题，应注意：a. 应把泛函本身讲清楚（即写出它的形式）；b. 还必须讲明白自变函数的性质，如：- 独立的自变函数的个数（导函数并不独立）；- 每个自变函数定义的区间/区域；- 这些自变函数应满足的条件（如：边界条件及其受约束的条件等）。c. 除了个别特殊情况外，一般情况下增加一个条件会使泛函极值及相应 的自变函数变化性质发生变化。如：极小值可能变大；极大值可能变 小;非极值的驻值可能成为极值。若干背景知识 泛函的驻值问题可以转化为等价的微分方程问题，变分法的理论计算就是完成这类工作。 本章内容沿袭此方法， 是要把问题的理论基础讲明确。 从近似解的角度出发，直接求解泛函的驻值，比解微

9、分方程更加方便， 也更为实用。特别计算机技术的发展，带来了大规模数值计算的可能性 (有限元的思想基础) 。 经 Euler, Lagrange, Dirichlet , Hilbert , Bernoulli 等数学先驱的卓越工作， 完成了的系统方法。 但把微分方程问题转换为泛函问题还很不成熟。在物理、力学中，即先 猜想一个泛函的驻值问题，再校对是否与原微分方程问题等价。 泛函驻值的计算(数值)先驱工作中以Ritz, Galerkin,Treft著名。关于变分法的一个预备定理若f(x)在a, b上连续，若对任意满足(a)= (b)=0的连续函数 x都有：bf(x) (x)dx 0a则f(x)在

10、a, b上处处为零。反证法：设xo为a, b中的点，在xo点f(xo)工0,可取f(xo)>O,v f(x)在区间上连续，必存在 xo的一个充分小邻域上f(x)>0, x 0- <x<xo+ 又T x为任意连续函数(满足边界条件)，可取x也在该邻域内大于零， 而在该邻域外恒等于零。所以有ba f (x) (x)dx o矛盾！即 f(x) 必须为零；同理可证小于零情况。 该定理可推广多元变量的函数问题。b1.2.1定积分 F(x,y, y)dx的驻值(变分)问题a目的：通过简单泛函的极值分析，获得建立变分法的基本概念、计算步骤(把变分解转化成微分方程)问题：在自变量x的区

11、间a，b 内决定一个函数y(x)，使它满足边界条件：y |x a ，y lx b 并使泛函：bV F(x, y, y )dx 取极值。a计算V方法1：先用变分观点解释GH曲线的增量设想已取得了一条曲线 GACH方程为：y= y (x)在GACH附近另取一条曲线GBDH，令该曲线无限接近GACH，其方程为:yi(x) y(x) y(x)y(x)是一个无穷小量，称为自变函数的变分(若 x不变，即为曲线纵坐标 的增量)(注意与函数微分的区别，这里函数的变分仍然是一个函数) 相应两条曲线，获得两个泛函值：bV a F(x,y, y )dxabV V F(x, y y, y y )dxa基本引理：(y)

12、 y证： y(x) y/x) y(x) ( y) y/x) y(x) y推广：(y) y另一条认识(y) y的思路：AC :y(xc)yg)y dxAB :W(Xb)ygyACD :力凶）y(xc)ycBD :y,XD)y,XB)y1dxy1y yy ydx ; y1%dxyiydxyi ybV F(x,y y,y y) F(x,y,y)dxa因为F（x, y,y）是x, y,y的连续可导函数（工程上一般如此），故y及 y很小时，V也很小，即 y, y 0 V 0取等式两端的一阶无穷小量，即:b FaFy dxy（可以从Tailor展开式去理解）V称为泛函V的一阶变分，简称变分，即泛函的一阶变

13、分是泛函增量中的 一阶小量部分（把自变函数的变分y作为一阶小量）所以，变分的运算服从 无穷小量的运算规则。计算V方法2：（把求泛函的极值转化成求普通函数的极值）记：y,x)yo(x)y(x) 01（yo及y固定）_bV （ ） a F（x,yoy,yo当V在yo上取极值，则相应于 极值条件：y）dxo的泛函值 V（）现在成为普通的函数V （ ）| o o （先不管该条件，现仅研究其导数计算）baF d(yo y dy)Ddxdy da上ya yydx上两式中出现,的项。y和y并不能独立变化，可设法把y项转换成只与y有关取分步积分:b Fbbydx uvdx uv L 取：u a yay dx（

14、上）ydx y|： a dx yy代入一阶变分式:V：上乎（上）ydx y|aa y dx yy要选定的函数满足边界条件，所以:b F d FV aG *ydx计算V 0若方括号内的函数在区间内不为 0,则可任选y使V大于零或小于零，即 使V不能获得极值，故需方括号的项为零。即： 2 （）0（ Euler 方程）y dx y此即与泛函驻值等价的微分方程。或：令V 0由变分基本定理：y任意连续函数，方括号中函数连续。-d（上）0 y dx yExample最速降线问题:（注不显含x）代入Euler方程，并乘以函数Q可得:FdFFFdFQ Q （） Q Q（Q ）0ydxyyydxy由于上 0（F中不显含x），上式中只要令Q y ,把上式配成全微分形式:dF、门£（f yy） 0这是因为:AfdxFdyy dxF dyy dx0）（代回原Euler方程，即得全微分）后两项由Q的假设 Q上Qy y由全微分方程F 十丫 C代入f的具体表达式:1.y(i y2)y(i y2) v令: y eta ntv1 cta n2tv