数学物理方法分离变量法.ppt

《数学物理方法分离变量法.ppt》由会员分享，可在线阅读，更多相关《数学物理方法分离变量法.ppt（82页珍藏版）》请在三一办公上搜索。

1、本章中心内容,用分离变量法求解各种有界问题,第二章 分离变量法,我把数学看成是一件有意思的工作,而不是想为自己建立什么纪念碑。可以肯定地说,我对别人的工作比自己的更喜欢。我对自己的工作总是不满意。-拉格朗日,1,本章基本要求,掌握有界弦的自由振动解及其物理意义,着重掌握分离变量法的解题思路、解题步骤及其核心问题-本征值问题,2,分离变量法核心：,本章考虑问题（1）混合问题（2）边值问题,本章层次：,3,偏微分方程常微分方程,齐次方程+齐次边界条件,非齐次方程+齐次边界条件,非齐次方程+非齐次边界条件,分离变量法思路起源,物理上由乐器发出的声音可以分解为各种不同频率的单音，每种单音振动时形成正弦

2、曲线，可以表示成,4,2.1 齐次方程问题,特点：含两个变量的函数可以表示为两个分别只含一个变量的函数之积。,这个定解问题的特点是：偏微分方程是线性奇次的，边界条件也是奇次的。,研究两端固定的弦的自由振动,定解问题,解：,这是解的分离变量,5,研究两端固定的弦的自由振动,定解问题,（第一类齐次边界条件）,由前面思路，设,x,t 是相互独立的变量,（求非零解）,1、分离变量,代入方程中，,分离过程：,得出两个常微分方程：,代入边界条件：,6,高数中结论：,2、求解本征值问题,7,若有二阶常系数线性齐次方程,其中p、q为常数，则特征方程为,本方程特征方程r2+=0，由上面结论知，方程的解与的不同取

3、值有关，分情况讨论：,8,此时X（x）=0，只有零解，不合题意；,同样只有零解，不合题意；,C2是积分常数,9,则X（x）的一族非零解为,上解称为满足边界条件的固有解（特征解），称为固有值（特征值），sin函数称为固有函数（特征函数）。,固得到下面一族解：,A、B 是积分常数,3、解出时间函数，得到一族解,时间函数解,10,解方程,n=1,2,3,代入初始条件，有,一般情况下满足不了，怎么办？！,利用叠加原理！,11,4、通过初始条件，求出通解,此时要满足初始条件，则,12,则定解问题的最终解为,13,5、物理意义：,是驻波，（固有振动模式）,相邻节点之间距离等于半波长,波长=,14,15,波

4、腹,波节,（4）确定级数解中的待定常数（利用初始条件）,6、分离变量法概要：,（1）将偏微分方程化简为常微分方程（U=XT),（2）确定固有值和固有函数（利用边界条件）,（3）确定形式解（级数形式解）,16,17,例：求解,（第二类齐次边界条件）,解：,设,18,此时边界条件为：,相应的特征值问题为：,此时X（x）=0，只有零解，不合题意；,19,同样只有零解，不合题意；,20,则特征解为,将特征值代入T（t）的方程，解出,则u（x，t）的特解族为,同样很难满足初始条件，由叠加原理得,21,此时要满足初始条件，有,22,故定解问题的最终解为,2.2 有限长杆上的热传导,23,24,25,此特解

5、仍然很难满足初始条件，由叠加原理得级数解为,26,由初始条件有,2.3 二维拉普拉斯方程的定解问题,（1）圆域,因为边界形状是个圆周，圆域边界条件中x、y是不可直接分离的，故化为极坐标求解。,27,28,第一步：求满足齐次方程、周期边值条件和原点约束条件的变量分离形式的解,29,30,周期本征值问题,欧拉方程,第二步：求解周期本征值问题和欧拉方程,31,根据叠加原理，得到级数解,32,第三步：利用边界条件,利用傅立叶级数系数的求解公式,33,欧拉方程,常系数线性微分方程,附录：欧拉方程,34,欧拉方程的算子解法:,35,则由上述计算可知:,用归纳法可证,于是欧拉方程,转化为常系数线性方程:,3

6、6,（2）矩形域,37,38,39,叠加后的级数解为,40,k=1,2,3,k=0,1,2,3,41,k=0,1,2,3,k=0,1,2,3,42,2.4 非奇次方程的解法,研究一根弦在两端固定的情况下，受强迫力作用所产生的振动现象。,即考虑下列定解问题：,43,怎么办？！,很明显现在不能直接用前面的变量分离起手式进行分解，因为等式右边的非齐次尾巴没办法处理！,现在的情况下，弦的振动和两个原因有关，一是外力，二是初始状态。,有否经历过类似情景？是否有可借鉴的类似情况？,44,借用结论：,这里我们用一招移花接木！,全响应=零输入响应+零状态响应,零输入=初始状态引起振动，与外力无关；零状态=外力

7、引起振动，与初始状态无关,45,设解为：,初始状态原因（零输入）,外力原因（零状态）,46,（零输入响应）,（零状态响应）,47,对V（x，t），可直接用前面的变量分类法求出：,48,对W（x，t），如何求？,49,50,51,52,设,解法二,53,54,55,原方程的解为：,56,例 在环形域 内求解下列定解问题,解,由于求解区域是环形区域，所以改选用平面极坐标系，利用直角坐标与极坐标系之间的关系,57,将上述定解问题用极坐标表示出来：,利用上节求出的圆域拉普拉斯方程的本征函数，设解为,58,代入方程并整理得到：,比较两端的系数可得,59,再由边界条件得,通解为：,求解得,60,特解有,所

8、以有,代入边界条件有,原定解问题的解为,61,2.5 非齐次边界条件的处理,设有定解问题,边界条件非齐次，若用前面方法分离变量，由边界条件没有办法得到只与某个常微分方程有关的具体边界函数值。,怎么办？！,62,想办法把边界条件化为齐次！,设法作一代换将边界条件化为齐次的，令,所以要求,选取W(x,t)使V(x,t)的边界条件化为齐次的，即,63,一般这样的函数是很容易找到的，最简单的如选取关于x的线性函数：,代入w（x，t）要满足的边界条件，可求出：,64,此时关于V的定解问题为,因此只要做如下代换，V将满足齐次边界条件。,65,其中,关于V（x，t）的问题即前述非齐次方程、齐次边界条件问题。

9、,66,当边界条件不同时，方法一致（关键在与w（x，t）的选取），W（x，t）的形式不同。,常用的最简单的w（x，t）形式,67,通过上式可以求出W（x）的形式。,注：若f，u1，u2都与t无关，则可选取W（x）（与t无关），使V（x，t）同时满足齐次方程和齐次边界条件，此时 W（x）需满足：,68,此时u（x，t）=V（x，t）+W（x），则V（x，t）满足,69,例1：求,的形式解，其中A,B均为常数。,解：令,代入方程有,70,通过二次积分即边界条件求得：,则V的方程为：,71,利用分离变量法，带齐次边界的方程的解为,利用第二个初始条件,代入第一个初始条件有,即,72,由傅里叶系数公式可

10、得,因此，原定解问题的解为：,73,例2 求定解问题,其中b,u1均为常数。,解：令,代入方程有,74,分解为两个方程,（零输入响应）,（零状态响应）,75,对于问题（I），可以直接采用分离变量法求解。,代入有,由边界条件有,由此得到下面两个常微分方程,76,易求得特征值和特征函数为：,代入含T的方程有,77,它的通解为,从而问题（I）的解可表示为,其中Cn由初始条件确定为,78,故所求的解V（1）(x,t）为,对于问题（II），可以用特征函数法求解，将方程的自由项及解都按特征函数系来展开。,79,其中vn(t)满足,由此可解得,80,从而问题（II）的解为,原方程的解为,本章小结,分离变量法核心：偏微分方程常微分方程,（1）基础,齐次方程+齐次边界条件,（直接变量分离法求解）,（2）复杂化,非齐次方程+齐次边界条件,分解为,齐次方程+齐次边界条件+非零初始条件（零输入）,（直接变量分离法求解）,非齐次方程+齐次边界条件+零初始条件（零状态）,（固有函数法求解，设出级数解形式，求时间函数）,81,（3）近一步复杂化,非齐次方程+非齐次边界条件,1）做代换使边界条件齐次；,2）同（2）,综上（1）为（2）的基础，（2）为（3）的基础！,82,