泛函和变分法PPT课件.ppt

计算物理,3/lesson/ComputationalPhysics,泛函和变分法,1,泛函和变分法,泛函和变分的基本概念最简泛函的极值问题其它类型泛函的极值问题泛函和变分用于微分方程边值问题,2,泛函和变分的基本概念(1/4),泛函的定义例(最短路径)：设C为定义在a,b上、满足条件y(a)=y1和y(b)=y2的、所有可微函数y(x)的集合。用L表示这样一段曲线的长(如右图所示)，L=Ly(x)问题：沿哪一条路径的路程最短函数的形式y(x)不同,例(捷线问题)：质点在重力作用下沿一条光滑的、从点A到点B的曲线运动，所需的时间T取决于曲线的形状(如右图所示)，T=Ty(x)问题：沿哪一条路径的下落时间最短函数的形式y(x)不同,3,泛函和变分的基本概念(2/4),定义：设C是函数(形式)的集合，B是实数集合；如果对C中的任一元素y(x)，在B中都有一个元素J与之对应，则称J为y(x)的泛函，记为Jy(x)泛函是函数的函数，以函数为自变量，而非普通变量最短路径：L=Ly(x),捷线问题：T=Ty(x),最简泛函：满足以下关系的泛函称为最简泛函,其中F(x,y,y)的称为核函数,4,泛函和变分的基本概念(3/4),函数的变分和泛函的变分定义：设y(x)是泛函Jy(x)的定义域内任意函数，如果y(x)变化为定义域内的另一新函数Y(x)，则Y(x)与y(x)之差dy=Y(x)-y(x)称为函数y(x)的变分函数变分和微分的比较变分和微分都是自变量x的函数微分是同一个函数y(x)，由于自变量x的取值不同而导致函数值y的变化；变分是由于函数形式的不同而导致函数值的变化函数求导和求变分可以交换次序,5,泛函和变分的基本概念(4/4),最简泛函的一阶和二阶变分,其中dJ称为泛函的一阶变分，d2J称为二阶变分泛函的极值条件就是一阶变分为零：dJ=0,6,最简泛函的极值问题(1/9),最简泛函的欧拉方程最简泛函的极值欧拉方程,欧拉方程的解仅仅对应极值函数，不关心泛函的大小通过变分运算等价于一定边界条件下的常微分方程例：如下泛函(不是最简泛函)的极值问题,等价于以下边界条件下的静电场中的泊松方程,7,最简泛函的极值问题(2/9),例：求以下最简泛函的极值问题,核函数和微分方程,满足边界条件的极值函数,例：求解最短路径问题,8,最简泛函的极值问题(3/9),例：求解捷线问题,9,最简泛函的极值问题(4/9),欧拉方程的其它算法如果F中不显含y,，不满足边界条件，则极值函数不存在,如果F中不显含y,如果F中不显含x,10,最简泛函的极值问题(5/9),例：再求解捷线问题,11,最简泛函的极值问题(6/9),例(最小旋转面)：光滑曲线以点A(x0,y0)和B(x1,y1)为端点(如右图)，求一条曲线使它绕Ox轴旋转时所得曲面的面积最小,以y(x)表示任意曲线，得旋转面面积,从欧拉方程的极值问题求曲线方程,12,最简泛函的极值问题(7/9),瑞利-里兹法的步骤选一组相对完备的基函数w0,w1,wn,，线性展开y,只取前面n项，作为y的近似，代入泛函，积分,Jy=I(a1,a2,an)按多元函数取极值方法,求解以上n个关于ai的方程，得到系数ai，代入展开式即可得到y的近似，再计算可得到Jy,取前面n+1项，重复以上2和3步，直至Jy收敛,13,最简泛函的极值问题(8/9),求解以下泛函的极值函数,取满足边界条件的基函数：wi=xi(1-x),只取前面n项，作为y的近似,14,最简泛函的极值问题(9/9),瑞利-里兹法的关键：选择合适的基函数幂函数：1,x,x2,=xi三角函数：1,cosx,sinx,cos2x,sin2x,其它：尽量同时满足边界条件,15,其它类型泛函的极值问题(1/4),依赖于多个函数的泛函泛函的一般形式,欧拉方程,例：求解以下泛函的极值问题,解：,16,其它类型泛函的极值问题(2/4),例：不均匀的介质中，折射率为n(x,y,z)，光的传播速度为c/n。求光从A(x0,y0,z0)到B(x1,y1,z1)的传播路径设G过A和B的某条光滑曲线：y=y(x),z=z(x)费马原理：光沿由A到B的所需时间最短的曲线行进,泛函的极值问题：要求T取极小值,17,其它类型泛函的极值问题(3/4),依赖于函数高阶导数的泛函泛函的一般形式,欧拉方程,例：求解以下泛函的极值问题,解：,18,其它类型泛函的极值问题(4/4),依赖于多元函数的泛函泛函的一般形式,欧拉方程,例：拉普拉斯方程的第三类边界问题,该定解问题转化为以下泛函的极值问题,19,泛函和变分用于(1/1),斯特姆-刘维型方程Ly=lr(x)y本征值：l1l2l3本征函数：y1(x),y2(x),y3(x),构成完备正交