数理方程第讲市公开课金奖市赛课一等奖课件.pptx

1第1页第1页2第二章 分离变量法2.1 有界弦自由振动第2页第2页3在高等数学中我们知道一个普通函数f(x)经常能够展开成级数.比如,幂级数形式就是:f(x)=a0+a1x+a2x2+a3x3+其中无穷多个函数v0(x)=1,v1(x)=x,v2(x)=x2,等等,构成了级数展开一个函数系.而三角级数形式就是f(x)=a0+a1sinx+b1cosx+a2sin2x+b2cos2x+其中无穷多个函数1,sinx,cosx,sin2x,cos2x,也构成了级数展开一个函数系.第3页第3页4因此,普通而言,一个函数f(x)能够在一个函数系v0(x),v1(x),v2(x),下展开成级数形式为f(x)=a0v0(x)+a1v1(x)+a2v2(x)+那么,一个二元函数u(x,t),将t固定住视为常数,看作x函数,则也能够在函数系v0,v1,v2,下展开成级数形式 u(x,t)=a0(t)v0(x)+a1(t)v1(x)+a2(t)v2(x)+其中每一项都是两个一元函数乘积ai(t)vi(x),这样构成二元函数我们称之为可分离变量.而假如级数中每一项都是线性偏微分方程解,则此级数也就是线性偏微分方程解.第4页第4页5讨论两端固定弦自由振动定解问题:设u(x,t)=X(x)T(t)则第5页第5页6代入方程(2.1)得X(x)T(t)=a2X(x)T(t)或此式左端仅是x函数,右端仅是t函数,普通情况不也许相等,除非它们均为常数,令此常数为-l,则有这样能够得到两个常微分方程:第6页第6页7再利用边界条件(2.2),由于u(x,t)=X(x)T(t),X(0)T(t)=0,X(l)T(t)=0.但T(t)0,假如T(t)=0,这种解称为平凡解,因此X(0)=X(l)=0(2.6)因此,要求方程(2.1)满足条件(2.2)变量分离形式解,就先要求解下列常微分方程边值问题第7页第7页8要拟定l取何值时(2.5)才有满足条件(2.6)非零解,又要求出这个非零解X(x).这样问题称为常微分方程(2.5)在条件(2.6)下特性值问题,使问题(2.5),(2.6)有非零解l称为该问题特性值,相应非零解X(x)称为它特性函数.下面分l0三种情况来讨论,将得出结论l0和l=0不能成立.第8页第8页9第9页第9页10而方程X(x)+lX(x)=0特性方程为r2+l=0当l0时,特性根为方程通解为第10页第10页111 设l0,并令l=b2,b为非零常数.此时方程(2.5)通解为 X(x)=A cos bx+B sin bx,由条件(2.6)得A=0B sin bl=0由于B不能为零,因此sin bl=0,即从而第13页第13页14(2.5),(2.6)一系列特性值及相应特性函数为:将上式中特性值代入到(2.4)得第14页第14页15其通解为:因此可分离变量方程特解为其中 是任意常数.第15页第15页16为满足初始条件(2.3),求出原问题解,将(2.10)中所有函数un(x,t)叠加起来:第16页第16页17将初始条件(2.3)代入上式得:第17页第17页18复习高等数学中周期为2l傅立叶级数:假如周期为2l周期函数f(x)为奇函数,则有其中系数bn为:第18页第18页19第19页第19页20第20页第20页21 解：令u(x,t)=X(x)T(t)是齐次方程和齐次边界条件非零解,则有第21页第21页22方程特解为第22页第22页23这时l=10,并给定a2=10000.这个问题傅里叶级数形式解可由(2.11)给出.其系数按(2.12)式为Dn=0,第23页第23页24因此,所求解为第24页第24页25解题中惯用到积分表内容:第25页第25页26分析一下级数形式解(2.11)物理意义.先固定t,看看任意指定期刻波是什么形状;再固定x,看该点振动规律.(2.11)中一项:其中第26页第26页27第27页第27页28第28页第28页29第29页第29页30某一时刻n=1,2,3驻波形状xOulxOulxOuln=1n=2n=3第30页第30页31综合上述,可知u1(x,t),u2(x,t),un(x,t),是一系列驻波,它们频率,位相与振幅都随n不同而不同.因此一维波动方程用分离变量法解出结果u(x,t)是由一系列驻波叠加而成,而每一个驻波波形由特性函数确定,它频率由特性值确定.这完全符合实际情况.因为人们在考察弦振动时,就发觉许多驻波,它们叠加又能够组成各种各样波形,因此很自然地会想到用驻波叠加表示弦振动方程解.这就是分离变量法物理背景,因此分离变量法也称为驻波法.第31页第31页322.2 有限长杆上热传导第32页第32页33设有一均匀细杆,长为l,两端点坐标为x=0与x=l,杆侧面是绝热,且在端点x=0处温度是零摄氏度,而在另一端x=l处杆热量自由发散到周围温度地零度介质中去,已知初始温度分布为j(x).求杆上温度改变规律,也就是要考虑下列定解问题:第33页第33页34用分离变量法来解此问题,设u(x,t)=X(x)T(t),代入方程(2.13)得上式左端不含有x,右端不含有t,只有当两端均为常数时才也许相等.令此常数为-b2,则有第34页第34页35从而得到两个线性常微分方程解方程(2.16)得X(x)=A cos bx+B sin bx,由边界条件(2.14)可知X(0)=0,X(l)+hX(l)=0.(2.17)从X(0)=0得A=0,从X(l)+hX(l)=0得b cos bl+h sin bl=0(2.17)a第35页第35页36第36页第36页37ygy=agy=tan gg1g2g3-g1-g2第37页第37页38于是得到无穷多个特性值及相应特性函数再由(2.16)解得得到一组满足边界条件特解为其中Cn=AnBn第38页第38页39由于方程(2.13)与边界条件(2.14)都是齐次,因此仍满足方程与边界条件.最后考虑u(x,t)能否满足初始条件(2.15),从(2.22)式得现在希望它等于已知函数j(x),首先要问在0,l上定义函数j(x)能否展开成上式形式,另一方面要问系数Cn如何拟定.前者答案是必定(不证).主要讨论后者.第39页第39页40不难证实令于是在两端乘上sin bkx,然后在0,l上积分得即将(2.24)代入(2.32)式即得原定解问题解.第40页第40页41分离变量法主要环节为:一一,首先将偏微分方程定解问题通过度离变量转化为常微分方程定解问题,这对线性齐次偏微分方程是能够做到.二二,求特性值问题，即拟定特性值与特性函数。当边界条件是齐次时,求特性函数就是求一个常微分方程满足零边界条件非零解.第41页第41页42三三,定出特性值、特性函数后,再解其它常微分方程,把得到解与特性函数乘起来成为un(x,t),这时un(x,t)中还包括着任意常数.四四,最后为了使解满足其余定解条件,需要把所有un(x,t)叠加起来成为级数形式,这时级数中一系列任意常数就由其余定解条件拟定.在这最后一步工作中,需要把已知函数展开为特性函数项级数,这种展开合理性将在2.6中叙述.第42页第42页432.3 圆域内二维拉普拉斯方程定解问题第43页第43页44 一个半径为r0薄圆盘,上下两面绝热,圆周围缘温度分布为已知,求达到稳恒状态时圆盘内温度分布.这时温度分布应满足拉普拉斯方程 2u=0由于边界形状是个圆周,它在极坐标下方程为r=r0,因此在极坐标系下边界条件可表为既然边界条件用极坐标形式表示出来很简朴,因此就在极坐标系下求解这个定解问题.第44页第44页45因r,q取值范围分别是0,r0与0,2p,而圆内包括中心温度有限,且(r,q)与(r,q+2p)事实上表示同一点,温度应当相同,即应当有|u(0,q)|+(2.27)u(r,q)=u(r,q+2p)(2.28)现在来求满足方程(2.25)及条件(2.26),(2.27),(2.28)解.先令 u(r,q)=R(r)F(q),第45页第45页46代入方程(2.25)得即令比值为常数l即得两个常微分方程F+lF=0,r2R+rR-lR=0.再由条件(2.27)及(2.28)可得|R(0)|+,F(q+2p)=F(q).(2.29)第46页第46页47因此得到两个常微分方程定解问题先解哪一个要看哪一个能够定出特性值.由于条件(2.29)满足可加性(即所有满足(2.29)函数加起来依旧满足(2.29),因此只能先解问题(2.30).第47页第47页48采用与2.1中同样办法能够得到当l0时,取l=b 2,这时(2.30)解为Fb(q)=abcosbq+bbsin bq,且为使F(q)以2p为周期,b必须是整数n,n=1,2,3,则可将上面得到解表示成Fn(q)=ancos nq+bnsin nq.第48页第48页49高等数学复习:求解欧拉方程x2y+xy-n2y=0(a)作变换x=et或t=ln x,则有代入(a)得通解为y=Cent+De-nt=Cxn+Dx-n(n0)和y=Ct+B=Cln x+D(n=0)第49页第49页50对于非齐次欧拉方程x2y+xy-n2y=axm特解形式应当是y*=Cxm,将之代入上式可拟定常数C.第50页第50页51至此,已经定出了特性值问题(2.30)特性值b2n=n2,特性函数Fn(q).接下去是解(2.31).其中方程是欧拉(Enler)方程,它通解为R0=c0+d0lnr,当l=0;Rn=cnrn+dnr-n,当l=n2(n=1,2,3,)为了确保|R(0)|+,只有dn=0(n=0,1,2,),即Rn=cnrn(n=0,1,2,).因此利用叠加原理,方程(2.25)满足条件(2.27),(2.28)解能够表示为级数第51页第51页52第52页第52页53将这些系数代入(2.32)式即得所求解.为了以后应用起来以便,还能够将解(2.32)写成另一个形式.为此,将(2.34)式系数代入(2.32)式通过简化后可得第53页第53页54利用下面已知恒等式第54页第54页55可将(2.35)中解u(r,q)表示为公式(2.36)称为圆域内泊松公式圆域内泊松公式.它作用在于把解写成了积分形式,这样便于作理论上研究.第55页第55页56例例 解下列定解问题A为常数.解解 利用公式(2.34)并注意三角函数正交性代入(2.32)即得所求解为第56页第56页572.4 非齐次方程解法第57页第57页58研究一根弦在两端固定情况下,受逼迫力作用所产生振动现象.即要考虑下列定解问题:由于非齐次方程解经叠加后普通不再是原方程解，因此不能用分离变量法直接求解非齐次方程定解问题。但是依据齐次方程(2.11)解，且它边值条件和方程(2.37)同样，故先假设定解问题(2.37)(2.39)解u(x,t)能够展开成下列傅立叶级数形式第58页第58页59并且把定解数据f(x,t)，和 都按固有函数系 展开第59页第59页60其中显然，u(x,t)满足边界条件(2.38)，因此只需u(x,t)再满足方程(2.37)和初值条件(2.39)就得到原问题解，把上面展开式分别代入方程(2.37)和初始条件(2.39)，可得第60页第60页61比较上面三个展开式系数可得由于(2.41)相应齐次方程解为第61页第61页62由高等数学可知：假如Cy1(x)是齐次线性方程解，那么能够利用变换y=uy1(x)(这变换是把齐次方程解中任意常数C换成未知函数u(x)而得到)去解非齐次线性方程。这一办法也适合用于高阶线性方程。下面就二阶情形来做讨论。假如已知齐次方程那么，我们能够用下列常数变易法去求非齐次方程通解为第62页第62页63令要拟定未知函数v1(x)及v2(x)使(3)式所表示函数y满足非齐次方程(2)。为此，对(3)式求导，得由于两个未知函数v1、v2只需满足一个关系式，因此可要求它们再满足一个关系式。从 上述表示式可看出，为了使 表示式中不含 和 ，可设第63页第63页64从而再求导，得把 、代入方程(2)，得整理得注意到y1及y2是齐次方程(1)解，故上式即为第64页第64页65联立方程(4)和(5)，在系数行列式时，可解得对上两式积分(假定f(x)连续)，得于是得非齐次方程(2)通解为第65页第65页66则利用常微分方程中参数变易法，即设为方程(2.41)解，则待定函数an(t),bn(t)由下面方程组拟定求得第66页第66页67把它们代入(2.43)，就得到方程(2.37)通解再由初始条件(2.39)，可拟定故问题(2.41)，(2.42)解为(2.44)把(2.44)代入(2.40)，可知(2.40)是非齐次问题(2.37)(2.39)解，完整形式为第67页第67页68(2.45)从解形式上看，可分为两部分：等号右端第一个级数项表示初始位移和初始速度对弦振动影响；第二个级数项表示外力f(x,t)对弦振动影响。若弦所受外力为零，则(2.45)式就是齐次问题(2.1)(2.3)解(2.11)。第68页第68页69上述这种解法是把方程非齐次项以及解按对应齐次方程一族固有函数展开，伴随方程与边界条件不同，固有函数族也就不同，但总是把非齐次方程解按对应固有函数展开，这种方法又称固有函数法。此方法对其它类型方程也是合用。第69页第69页702.5 非齐次边界条件处理第70页第70页71 前面所讨论定解问题解法,无论方程是齐次还是非齐次,边界条件都是齐次.假如碰到非齐次边界条件情况,应当如何处理?总原则是设法将边界边界条件化成齐次条件化成齐次.现在以下列定解问题为例,阐明选取代换办法.以弦振动为例，设含有非齐次边界条件弦振动定解问题为令第71页第71页72其中v(x,t)满足和u(x,t)相同边界条件(2.47)则当x=0或x=l时这样关于V(x,t)定解问题边值条件就是齐次，相应辅助函数v(x,t)也容易找到，对于第一边值问题，普通可设v(x,t)=a(t)x+b(t)，代入(2.47)可得第72页第72页73即把v(x,t)代入(2.49)，可得再把上式代入(2.46)(2.48)，则定解问题转化为第73页第73页74重复2.4节做法，就可得到V(x,t)，进而求出u(x,t).对方程和边界条件都是非齐次，且f,u1,u2都与t无关,则可取适当v(x)(也与t无关),使V(x,t)方程与边界条件同时都化为齐次,这样做就能够省掉下面对V(x,t)要进行解非齐次方程繁重工作.这种v(x)终归怎么找,将在下面例题中阐明.第74页第74页75例例1 求下列定解问题:形式解,其中A,B均为常数.第75页第75页76解解 这个定解问题特点是:方程及边界条件都是非齐次.依据上述原则,首先应将边界条件化成齐次,由于方程(2.63)自由项及边界条件都与t无关,因此我们有也许通过一次代换将方程与边界条件都变成齐次,详细做法下列:令 V(x,t)=u(x,t)-W(x),代入方程(2.63)得第76页第76页77为了使这个方程及边界条件同时化成齐次,选W(x)满足(2.66)是一个二阶常系数线性非齐次方程边值问题,它解能够通过两次积分得第77页第77页78再由(2.65)可知函数V(x,t)为下列定解问题:解.采用分离变量法,可得(2.67)满足齐次边界条件(2.68)解为第78页第78页79利用(2.69)中第二个条件可得Dn=0.于是定解问题(2.67),(2.68),(2.69)解表示为代入(2.69)中第一个条件得即第79页第79页80由傅里叶级数系数公式得第80页第80页81因此,原定解问题解为其中Cn由(2.71)拟定.第81页第81页82对边界条件不全是第一类,本节方法依然合用,不同只是辅助函数v(x,t)形式。第82页第82页831v(x,t)=u2(t)x+u1(t)2v(x,t)=u1(t)(x-l)+u2(t)3v(x,t)=(1/2l)(u2(t)-u1(t)x2+u1(t)x注意以上v(x,t)选取不是唯一。第83页第83页84以上各节阐明了如何用分离变量法来解定解问题,其主要环节小结下列:一,依据边界形状选取适当坐标系,选取原则是使在此坐标系中边界条件表示式最为简朴.圆、圆环、扇形等域用极坐标系较以便.二,若边界条件是非齐次,则无论方程是否为齐次,必须先作函数代换使其化为含有齐次边界条件问题,然后再求解.三,非齐次方程,齐次边界条件定解问题(无论初始条件如何)能够用特性函数法求解.第84页第84页852.6 积分变换法第85页第85页86积分变换通过特定积分运算，把一个函数变成另一个函数变换。积分变换法是通过积分变换，将数学模型转化，从而简化定解问题一个求解办法。如通过积分变换将偏微分方程化为常微分方程，于是求解问题得以简化。尤其对于无限或半无限区域上定解问题，采用积分变换有固定程序求解，更为以便。第86页第86页87这里函数f(x)经过上述积分运算变成另一函数F(s)就称为一个积分变换，其中k(x,s)称为积分变换核，当选取不同积分变换核和积分域时，就得到不同积分变换。下面介绍积分变换(傅里叶变换和拉普拉斯变换)在求解偏微分方程定解问题中应用。普通，含参变量s积分第87页第87页88一、傅里叶变换(Fourier)1.傅里叶变换定义若函数f(x)在 上满足：逐段光滑；绝对可积，即 收敛则称 为函数f(x)傅里叶变换，简称傅氏变换，记为 能够推出 称为傅氏逆变换，记为第88页第88页892.傅里叶变换性质 线性性质 设 ，a,b为任意常数，则 同样逆变换也成立，即微分性质 类似地，逆变换有第89页第89页90积分性质卷积定理 卷积定义：若已知函数f1(x),f2(x)，则积分称为函数f1(x)与f2(x)卷积,记为f1(x)*f2(x),即第90页第90页91显然卷积定理：假定则或第91页第91页92二、拉普拉斯变换 傅氏变换要求进行变换函数一定要满足绝对可积，这样条件是比较强，许多简朴函数如1，xn,ex和sinx等都不满足在内绝对可积；另外，在工程应用中许多以时间t为自变量函数仅在 上有定义。因此，傅氏变换应用范围受到很大限制。为了克服傅氏变换缺点，就需要适当地把傅氏变换加以改造，从而导出拉普拉斯变换，简称拉氏变换。第92页第92页931.拉氏变换定义 设函数f(t)，当 时有定义，且积分在s某一区域内收敛，其中 是复参量，则由此拟定函数称为f(t)拉普拉斯变换式。记作第93页第93页94并称函数F(s)为f(t)拉氏变换，称f(t)为F(s)拉氏逆变换，记作显然若F(s)是f(t)拉氏变换，则可推出2.拉氏变换性质 线性性质 若a,b为常数，则第94页第94页95且逆变换也有微分性质 设f(t)在 上连续，则积分性质第95页第95页96推论:卷积定义：若f1(t)，f2(t)满足拉氏变换存在条件，则积分称为f1(t)，f2(t)卷积，记为f1(t)*f2(t)，即实质上，拉氏变换卷积与傅氏变换卷积定义是一致，当 时，若f1(t)，f2(t)第96页第96页97满足条件则从傅氏变换卷积可得到拉氏变换卷积，即卷积定理：设f1(t)，f2(t)满足拉氏变换存在定理中条件，且第97页第97页98则或三、用傅氏变换法求解定解问题 例：求解一维齐次热传导方程定解问题第98页第98页99解：首先对于未知函数u(x,t)及初始条件中函数 关于x作傅氏变换，记然后，对方程两边关于x作傅氏变换，并利用微分性质得即第99页第99页100这是一个含参数 一阶常微分方程，对初始条件也作同样变换得解常微分方程初值问题，其解两端关于 作傅氏逆变换，左端为而右端依据卷积定理知第100页第100页101由于 ，并查傅氏变换表得故定解问题解第101页第101页102四、用拉氏变换求解定解问题解：这个问题显然不能用傅里叶变换来求解了，由于x,t变换范围都是 ，下面用拉氏变换来解，从x,t变换范围来看，对x,t都能取拉氏变换，但由于方程中含有 ，而在x=0处未给出 值，故不能第102页第102页103对x取拉氏变换，而对t来说，由于方程中只出现关于t一阶偏导数，只要知道当t=0时u值就够了，这个值已由给出，故我们采用关于t拉氏变换。用U(x,p),F(p)分别表示函数u(x,t),f(t)关于t拉氏变换，即首先，对方程两端取拉氏变换，并利用条件，可得第103页第103页104再对条件取同样变换，可得方程是关于U(x,p)线性二阶常系数常微分方程，它通解为由于当 时，u(x,t)应当有界，因此第104页第104页105U(x,p)也应当有界，故C2=0。再由条件可得C1=F(p)，从而可得由拉氏变换表可查到其中余误差函数LL第105页第105页106再依据拉氏变换微分性质可得因此第106页第106页107最后由拉氏变换卷积性质可得原问题解第107页第107页108选择变换原则:1.自变量改变范围，由于傅氏变换要求自变量在 内改变，拉氏变换要求自变量在 内改变，因此要依据自变量改变范围，作为选择变换条件之一。2.定解条件形式，由于拉氏变换微分性质：能够看出，对函数关于某自变量取拉氏变换时，必须在定解条件中给出当自变量等于0时函数值及相关导数值。第108页第108页109用积分变换法求解定解问题过程:一,依据自变量改变范围以及定界条件详细情况,选取适当积分变换,然后对方程两端取变换,把一个含两个自变量偏微分方程化为含一个参量常微分方程.二,对定解条件取相应变换,导出新方程定解条件.三,解所得常微分方程,求得原定解问题解变换式(即象函数).四,对所得变换式取逆变换,得到原定解问题解.第109页第109页110作业习题二41页开始第1(1),2(2),3(1),6(1),7,10题第110页第110页