第四篇 分离变量(傅立叶级数)法1.pptVIP

* 第四章 分离变量法1 * u取非零解的条件为 因为ln0，一阶常微分方程的通解为 (A8) 因此得傅氏解 * 第四章 分离变量法1 * 由初始条件，得 傅里叶正弦级数的展开系数为 (A9) 第四章 分离变量法1 第四章 分离变量法1 * 第四章 分离变量法1 * 第四章 分离变量(傅立叶级数)法 §4.1 齐次方程的分离变量法 §4.2 非齐次方程和输运方程 §4.3 非齐次边界条件的处理 §4.4 泊松方程 （计划授课6学时） * 第四章 分离变量法1 * §4.1 齐次方程的分离变量法 0、分离变量法介绍 无限长(或半无限长)弦的解叫做行波解，由达朗贝尔公式给出。本章研究有限长弦的振动问题，用分离变量法(工具)，得到驻波解，它的叠加构成无穷级数形式的一般解。 分离变量法的实质就是将物理量u(x,y,z,t)按自变量写成若干个单自变量函数相乘的形式， 从而将偏微分方程(PDEs)的求解归结为若干个常微分方程(ODEs)的定解问题. * 第四章 分离变量法1 * (1) 应用分离变量法，令试探解为 (2) 将Eq.(2)代入Eq.(1)的泛定方程和边界条件，分别得到 (3a) (3b) 考虑有限长两端固定均匀弦的自由振动，则定解问题为 1、端点固定情况 * 第四章 分离变量法1 * 式(3a) 两边同除a2XT，得 第一个等号左边是时间t的函数，右边只与坐标x有关，左右相等的条件是二者相等且等于一个常数，记为-l. (5c) 另一方面，由式(3b) 得X(x)的边界条件 (4) 由式(4) 得到两个独立的二阶常微分方程： (5a) (5b) 式(5a)和(5b)又称为二阶常系数微分方程，下面根据常数l的符号，求出相应通解. * 第四章 分离变量法1 * 上式有两个根r＝?(-l)1/2，因此方程(5a)的通解为 其中系数c1和c2由边界条件X(0)=X(l)=0确定。 (6) (1) l=0的情况 方程(5a)为X’’=0，其通解为X(x)=c1x+c2，代入X(0)=X(l)=0中得 即X=0，u=XT=0为平凡解(trivial [’trivi?l]平庸的)，因此舍弃! 当l?0(即l0或l0)时，齐次ODE方程X’’+lX=0有e指数形式的通解X=erx，代回原方程中，得 * 第四章 分离变量法1 * (3) l0的情况 因为r＝?(-l)1/2= ?il1/2，由(6)式得通解为 代入X(0)=X(l)=0中，得 (7) (2) l0的情况 通解(6)式代入X(0)=X(l)=0中，得 此时X=0，即u=XT=0，这是一个平凡解，因此舍弃! * 第四章 分离变量法1 * 为了得到非平凡解，要求C2?0，因此必有 ，则 即 (8) 相应地 (9) 其中C2为积分常数。注意n?0，否则X(x)=0为平凡解。 常微分方程X’’+lX=0及边界条件X(0)=X(l)=0构成本征值问题，(8)式给出的特定的l值称为本征值，与之对应的非平凡解(9)式称为本征函数。 * 第四章 分离变量法1 * 下面由Eq.(5c)求解T(t). 对于每一个ln(或n)，齐次方程T’’+lna2T=0的通解T=ert, 其中待定系数r2=lna20, 因此得到 其中A, B是积分常数，由定解问题的初始条件确定. (10) 将Eqs. (9)和(10)代入Eq.(2)，得两端固定有限长弦上的驻波解 (又称傅里叶级数解) (11) 每一个正整数n对应一种驻波。不同的n给出两端固定弦的不同的本征振动模式：n=1称为基波，n1称为n次谐波. * 第四章 分离变量法1 * 驻波解(11)式又可写为(参看P95) 其中 ，本征振动模式n的圆频率 ，频率为 ，初始位相为 ，波数kn=np/l. (11’) 由频率得周期T(n)=1/fn=2p/wn，于是波长为T(n)a=2l/n. * 第四章 分离变量法1 * 2次谐波 基波 1次谐波 驻波n（第n个本征振动模式）有n个波腹，n+1个波节. 波腹 波腹 波腹 波腹 波腹 波腹 波节 波节 波节 波节 波节 波节 波节 波节 波节 * 第四章 分离变量法1 * 所有本征函数的线性叠加也是满足泛定方程和边界条件的解 (12) 其中An和Bn由初始条件(1b)确定，得 (13) 上式刚好是傅里叶正弦级数，系数为(P90) (14) * 第四章 分离变量法1 * 磁致伸缩换能器、鱼群探测器可归结为两端自由杆的纵振动，定解问题为 (15) 解：利用分离变数法，令试探解u(x,t)=T(t)X(x)，代入泛定方程，