变分原理与变分法.pdf

《变分原理与变分法.pdf》由会员分享，可在线阅读，更多相关《变分原理与变分法.pdf（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

2、的途径（Fermat）；Summary:实际上光的传播遵循最小能量原理；在静力学中的稳定平衡本质上是势能最小的原理.二、变分法是自然界变分原理的数学规划方法（求解约束方程系统极值的数学方法），是计算泛函驻值的数学理论ABv1v2ECAE EB ACCB数学上的泛函定义定义：数学空间（集合）上的元素（定义域)与一个实数域间（值域）间的（映射)关系特征描述法:J：X X D D R R|J(x x)r R RExamples：矩阵范数：线性算子（矩阵）空间数域nnA1=maxaij；A maxji1iaij；A2(aij)21nn2j1j1 i1 函数的积分：函数空间数域J fn(x)dxabfn

3、 D DNote：泛函的自变量是集合中的元素（定义域）；值域是实数域。Discussion：判定下列那些是泛函：f(x,y)；3x+5y=2；(x x0)f(x)dx f(x0)f max f(x)；axbx 试举另一泛函例子.q(x)E、Jconstsx x物理问题中的泛函举例 弹性地基梁的系统势能x=0,固支；x=l,自由1ld2w2i.i.梁的弯曲应变能:bEJ(2)dx20dxii.ii.弹性地基贮存的能量：f1lkw2dx20l0iiiiii。外力位能：l qwdxiv iv。系统总的势能：d2w212 0 EJ(2)2kw qwdx;x 0w 0dxl12dw 0dx泛函的提法：有

4、一种梁的挠度函数（与载荷无关)，就会有一个对应的系统势能。泛函驻值提法：在满足位移边界条件满足位移边界条件的所有挠度函数中，找一个 w(x）,使系统势能泛函取最小值.最速降线问题问题：已知空间两点 A 和 B，A 高于 B，要求在两点间连接一条曲线，使得有重物从 A 沿此曲线自由下滑时，从 A 到 B 所需时间最短（忽略摩擦力）。作法：i i。通过 A 和 B 作一垂直于水平面的平面，取坐标系如图。B 点坐标(a,b),设曲线为 y=y(x),并已知:x=0，y=0；x=a，y=bii.ii.建立泛函：设 P(x，y）是曲线上的点,P 点的速度由能量守恒定律求得：21mv mgyv 2gy2命

5、 ds 为曲线弧长的微分,有：1 y2dsdsdx v 2gy dt dt2gy2gyAx重物从 A 点滑到 B 点的总时间:T=ay1 y22 gy0dxypB泛函驻值提法：在 0 xa 的区间内找一个函数 y（x)使其满足端点几何条件并使 T 取最小值。圆周问题问题:在长度一定的闭曲线中,什么曲线所围成的面积最大。作法:i.i.假设所考虑的曲线用参数形式表示：x=x（s),y=y(s）s 为参数。取 s1为曲线上的某一定点,则坐标表示 x1=x（s1)，y1=y（s1)，因曲线是封闭的，必存在一个 s2点使 x2=x（s2），y2=y（s2）与点 s1（x1,y1)重合。ii ii。该封闭

6、曲线的周长：L=s2s1(dx2dy2)()dsdsds该曲线所围成的面积：R=dxdyiiiiii。转换 R 的表达式由 Green 公式:s2QP()dxdy Pdx Qdys1xyyQPx1取 P=2，Q=2，则：xyR s122s1xdy ydx s12(xy(s)2s1yx(s)ds泛函驻值的提法：等周问题即是在满足端点条件 x(s1）=x(s2），y(s1）=y（s2)及周长一定s2s1dy 22(dx)()L条件下，寻找一个曲线函数dsdsx(s)y(s)使泛函R 取驻值。Discussion悬索线问题:已知空间中 A，B 两点及一条长度 LAB的悬索,单位长的质量为 m.假设绳

7、索的长度是不变的,并忽略绳索的弯曲刚度，把此绳索的两端挂在 A，B 两点，求在平衡状态下绳索的形状。要求要求：列出悬索线应满足的泛函式及泛函驻值提法.提示提示：绳索在平衡状态下，其势能应为最小值。1。2 变分法（泛函驻值的计算方法)关于计算固体力学中的泛函、泛函极值的提法 这里所研究的泛函一般用积分显式表达，并不等于所有泛函都能用显式积分表达。所要研究的泛函都可表示成在一定区间或一定区域内的函数及其导数（或偏导数）的积分形式，即：a a。1b.b.2aF(f(x),f(x),f(x);x)dxbF(f(x,y),fx(x,y),fy(x,y);x,y)dxdyc.c.泛函中的可变化函数称为自变

8、函数，或称宗量（argument），x 或 y 仅是积分变量，是被积函数的定义域。（被积函数是复合函数概念的推广）要说清楚一个泛函的极值问题,应注意:a.a.应把泛函本身讲清楚（即写出它的形式）；b.b.还必须讲明白自变函数的性质，如：独立的自变函数的个数(导函数并不独立）;-每个自变函数定义的区间/区域；-这些自变函数应满足的条件（如：边界条件及其受约束的条件等）。c.c.除了个别特殊情况外,一般情况下增加一个条件会使泛函极值及相应的自变函数变化性质发生变化。如：极小值可能变大;极大值可能变小；非极值的驻值可能成为极值。若干背景知识 泛函的驻值问题可以转化为等价的微分方程问题，变分法的理论计

9、算就是完成这类工作.本章内容沿袭此方法，是要把问题的理论基础讲明确。从近似解的角度出发，直接求解泛函的驻值，比解微分方程更加方便,也更为实用.特别计算机技术的发展，带来了大规模数值计算的可能性（有限元的思想基础）。经 Euler,Lagrange,Dirichlet，Hilbert，Bernoulli 等数学先驱的卓越工作,完成了的系统方法。但把微分方程问题转换为泛函问题还很不成熟。在物理、力学中,即先猜想一个泛函的驻值问题，再校对是否与原微分方程问题等价。泛函驻值的计算（数值)先驱工作中以 Ritz，Galerkin，Treft著名。关于变分法的一个预备定理都有：若 f（x）在a，b上连续，

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

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

12、x,y y,yy)dxab基本引理:(y)y(x)y(x)y证：y(x)y1(x)y(x)(y)y1推广：(y)y另一条认识(y)y的思路：yHA C：y(xC)y(xA)ydxA B：y1(xB)y(xA)yAC D：y1(xD)y(xC)yCdxB D:y1(xD)y1(xB)y1DBCAGadxy1 yyy ydx;y1 y1bdxx(y)by1ydx y y y1V F(x,y y,yy)F(x,y,y)dxa因为F(x,y,y)是x,y,y的连续可导函数(工程上一般如此),故y及y很小时,V也很小，即y,y 0V 0取等式两端的一阶无穷小量，即：V abFFy ydx（可以从 Tai

13、lor 展开式去理解）yyV称为泛函 V 的一阶变分，简称变分，即泛函的一阶变分是泛函增量中的一阶小量部分(把自变函数的变分y作为一阶小量）所以，变分的运算服从无穷小量的运算规则.计算计算V方法方法 2 2:（把求泛函的极值转化成求普通函数的极值）记：y1(x)y0(x)y(x)0 1（y0及y固定）y)dxV()F(x,y0y,y0ab当V在 y0上取极值，则相应于 0的泛函值V()现在成为普通的函数极值条件：V()|0 0（先不管该条件，现仅研究其导数计算）bF d(y y)bF y)dV()F d(y0F0V ddxdy ydxaadydydyy上两式中出现，y和y并不能独立变化，可设法

14、把y项转换成只与y有关的项。取分步积分:bababFFydx uvdx uv|bv ya取：u ayybdFFFy dx ()ydxy|baadx yyy代入一阶变分式：V abFd FF()ydxy|baydx yy要选定的函数满足边界条件，所以:y 0|xa，y 0|xbV abFdF()ydxydx y计算V 0若方括号内的函数在区间内不为 0，则可任选y使V大于零或小于零,即使V 不能获得极值，故需方括号的项为零。即：FdF()0（Euler 方程）ydx y此即与泛函驻值等价的微分方程。或：令V 0由变分基本定理：y任意连续函数,方括号中函数连续.FdF()0ydx y1 y2Exa

15、mple最速降线问题:F（注不显含 x）2gy代入 Euler 方程,并乘以函数 Q 可得：QFdFFFdFQ()QQ(Q)0ydx yyydxy由于F 0（F 中不显含 x）,上式中只要令Q y，把上式配成全微分形式:xdF(F y)0这是因为:dxydFF dyF dyFF（0）（代回原 Euler 方程，即得全微分）dxxy dxy dxxFF后两项由 Q的假设 Q QyyF F y C代入 F 的具体表达式：由全微分方程y1y(1 y2)C y(1 y2)v y v1 y2令：y ctanty dx vv2 vsin t(1cos2t)221 ctan tdy2vsintcostdt 2vsin2tdt v(1cos2t)dtyctant上式积分得：x 注意:y v(2t sin2t)C12v(1cos2t)2引用初始条件：x=0，y=0,只能有：C1=0v令：,2t 2x(sin)即为最速降线（圆滚线(渐开）)方程。y(1cos)Homework：在连接 XY 平面上两点的 M0及 M1的所有曲线中，要这样一条曲线使它绕 OX 轴旋转成的曲面有最小表面积.