,化工应用数学 第五章,*,*,2024/11/15,1,5.1,引言,5.1.1,偏微分方程的定义,描述物理量在时、空域中变化规律的方程，,若含有未知函数的偏导数，则称之为偏微分方程。,第五章 偏微分方程与特殊函数,2023/9/2815.1 引言5.1.1 偏微分方程的定义,2024/11/15,2,5.1.2,偏微分方程的规定,（,1,）方程中出现的偏导数的,最高阶数,称为,方程的阶数,。,（,2,）若方程中没有未知函数及其偏导数的乘积或幂等非线性项称方程为线性的，反之统称成为非线性的。,在非线性方程中若仅对未知函数的所有最高阶偏导不是非线性的,拟线性的。,（,3,）不含未知函数及其偏导数的项称之为自由项。自由项为零的方程称为齐次方程，否则称为非齐次方程。,第五章 偏微分方程与特殊函数,2023/9/2825.1.2 偏微分方程的规定第五章 偏,2024/11/15,3,例,5.1.2,偏微分方程的规定,2023/9/283例5.1.2 偏微分方程的规定,2024/11/15,4,一般来说，求解,n,阶线性常微分方程的通解，必定包含有,n,个独立的任意函数，若存在,n,个边界条件，则可确定这,n,个常数，从而获得该方程满足边界条件的一个特定解。在偏微分方程中，通解具有特定形式的任意函数，且看以下例题。,(1),偏微分方程的定解问题,5.1.2,偏微分方程的规定,2023/9/284一般来说，求解n阶线性常微分方程的通解，,2024/11/15,5,(2),偏微分方程,是方程的通解，其中,f,是任意函数。,如,等都是该方程的特解。,这些不同的函数，都满足二维拉普拉斯方程,都可作为一维热传导方程,而函数,的通解。,2023/9/285 (2)偏微分方程,2024/11/15,6,由此可以得出两点结论：,偏微分方程的通解包含有任意函数，因此解偏微分方程，一般都不是先求通解，后由定解条件确定特解，而是直接求特解。,一个特定形式的偏微分方程可以描述许多物理现象的共性规律，可以有很多不同形式的特解。所以可称为“泛定方程”。,5.1.2,偏微分方程的规定,2023/9/286 5.1.2 偏微分方程的规定,2024/11/15,7,确定地描述某个系统的运动过程，除了反映运动一般规律的偏微分方程（泛定方程）外，还必须根据实际问题的模型提出定解条件。定解条件包括初始条件（当方程含有时间变量时）和边界条件（关于空间变量的约束条件）。,泛定方程加定解条件构成一个确定的物理过程的“定解问题”，此时问题才可能有确定的特解。,5.1.2,偏微分方程的规定,2023/9/287 确定地描述某个系统的运动过程，,2024/11/15,8,(3),偏微分方程的求解方法,1,、方程的建立,建模。,2,、求解,（,1,）解析解：,分离变量法,；,拉普拉斯变换法,等。,（,2,）数值解：尤拉法、龙格库塔法等。,5.1.2,偏微分方程的规定,2023/9/288(3)偏微分方程的求解方法5.1.2,2024/11/15,9,5.2,二阶偏微分方程分类,限两个自变量的二阶线性方程，未知函数,u,（,x,，,y,）,一般形式：,F,（,x,，,y,，,u,，,u,x,，,u,y,，,u,xx,，,u,yy,，,u,xy,）,=0,(5-4),线性形式：,Au,xx,+,2,Bu,xy,+Cu,yy,+Du,x,+Eu,y,+Gu+f=0,(5-5,）,A,，,B,，,C,，,D,，,E,，,f,，,G,是,x,，,y,的函数,当,A,，,B,，,C,，,D,，,E,，,G,是常数时，式（,5,5,）是二阶常系数线性偏微分方程，,f,f,（,x,，,y,）,为已知函数是自由项。,2023/9/2895.2 二阶偏微分方程分类限两个自变量的,2024/11/15,10,由参数,A,，,B,，,C,判断二阶线性方程的分类,设,M=M,（,x,，,y,）,为自变量域内的某一点，若在该点处有：,(1),B,2,AC,0,，则方程在该点处为双曲线型的，如：,u,xx,u,yy,0,（,5,6,）,(2),B,2,AC,0,，则方程在该点处为抛物线型，如：,u,y,u,xx,0,（,5,7,）,(3),B,2,AC,0,，则方程在该点处为椭圆型的，如：,u,xx,u,yy,0,（,5,8,）,2023/9/2810由参数A，B，C判断二阶线性方程的分类,2024/11/15,11,方程的类型在域内不一定是唯一的。如：,xu,xx,yu,yy,2yu,x,xu,y,0,B,2,-AC,0,xy,，,xy,0,，,M,在二，四象限双曲型,xy,0,，,M,在一，三象限椭圆型,xy,0(x,或,y,0),，,M,在,y,或者,x,轴上 抛物型,三类方程：（最典型的物理含义）,双曲线型：,u,tt,a,2,（,u,xx,u,yy,u,zz,）,波动方程,抛物线型：,u,t,a,2,（,u,xx,u,yy,u,zz,）,热传导方程,椭圆型：,u,xx,u,yy,u,zz,0,拉普拉斯方程,2023/9/2811方程的类型在域内不一定是唯一的。如：,2024/11/15,12,5.3,典型方程的建立,5.3.1.,拉格朗日方法,（一阶）（,将偏微分方程转化为常微分方程,）,2023/9/28125.3 典型方程的建立5.3.1.,2024/11/15,13,例,1,5.3,典型方程的建立,解：采用拉格朗日法,联立任意两个方程,2023/9/2813例 15.3 典型方程的建立解：采,2024/11/15,14,5.3,典型方程的建立,2023/9/28145.3 典型方程的建立,2024/11/15,15,5.3,典型方程的建立,其解为,即,方程的解为：,例,2,2023/9/28155.3 典型方程的建立其解为即方程的,2024/11/15,16,5.3.2,非线性（一阶）,5.3,典型方程的建立,2023/9/28165.3.2 非线性（一阶）5.3,2024/11/15,17,5.3.2,非线性（一阶）,解：,2023/9/28175.3.2 非线性（一阶）解：,2024/11/15,18,5.3.2,非线性（一阶）,2023/9/28185.3.2 非线性（一阶）,2024/11/15,19,5.3.3,拉普拉斯变换法,(1),将常微分方程转化为代数方程,(2),将偏微分方程转化为常微分方程,5.3,典型方程的建立,通过拉斯变换，可将两个自变量的偏微分方程转化为常微分方程。由于拉氏变换包括变量从零多无穷的积分，所以仅对从零到无穷有意义的自变量才有可能进行拉氏变换，一般来说只适用解初值问题。,2023/9/28195.3.3 拉普拉斯变换法(1),2024/11/15,20,5.3.3,拉普拉斯变换法,例,定解问题,用拉氏变化求解，需要首先确定对哪一个自变量施行拉氏变换，就本例而言，对,x,与,t,都可以。但对,x,偏导是二阶的，且没给出 ，所以选择对,t,作拉氏变换。,用,u(x,p),表示函数,u(x,t),关于,t,的拉氏变换，即,首先对两端作拉氏变换，另利用条件。可得新方程，过程如下：,2023/9/28205.3.3 拉普拉斯变换法例 定,2024/11/15,21,5.3.3,拉普拉斯变换法,同时对另一边界条件作拉氏变换,二阶常微分方程，其解为,由于 时。,u(x,t),应该有界，所以,u(x,p,),也应该有界，故,B=0,，,再有条件得：,2023/9/28215.3.3 拉普拉斯变换法同时对另一,2024/11/15,22,用拉氏变换解偏微分方程的要点是：,（,1,）首先确定对哪个自变量作拉氏变换。要求该自变量变化范围（,0,，），而且根据拉氏变换的微分性质,该变量必需具备上式有关的初值条件。如若有两个自变量都满足要求，那应取决于对哪个自变量求解过程最简单为准。,（,2,）除对方程作拉氏变换外，还要对凡在方程变化中没有用到的定解条件都要作拉氏变换，使其作为变换后新方程的定解条件。,5.3.3,拉普拉斯变换法,2023/9/2822 用拉氏变换解偏微分方程的要点是：,（,3,）最后得到定解问题的解的关键是对新方程之解作拉氏逆变换。当象函数较复杂时，运用查表和拉氏变换一章介绍的几种求逆变换方法也不得其解时，就只能运用拉氏变换的反演公式，通常用复变函数的围道积分法求解。,5.3.3,拉普拉斯变换法,（3）最后得到定解问题的解的关键是对新方程之解,2024/11/15,24,5.4,定解条件和定解问题,引言中指出,数学物理方程是具有某类共性的物理现象的泛定方程。在引言中我们也看到了，同一个泛定方程可以有多个不同函数的解。因此对实际的物理现象特性的讨论还需对其特定的,“,环境,”,和起始状态加以描述和限定,定解条件，结合泛定方程，便可确定定解问题的特解。,2023/9/28245.4 定解条件和定解问题,2024/11/15,25,5.4.1,初始条件（初值条件）,对于随着时间而发生变化的问题，必须考虑研究对象的初始时刻的状态，即初始条件。,凡泛定方程中只含,t,的一阶偏导数的只需要一个初始条件，,u,的初始分布。,泛定方程中含有,t,的二阶偏导数的则需要两个初始条件，初始分布和初始速度。,初始条件给出了整个系统的状态（,t=0,）。稳态过程因与,t,无关，则不存在初始条件,2023/9/2825 5.4.1 初始条件（初值,2024/11/15,26,5.4.2,边界条件,(1),第一边界条件,已知函数,直接给出在边界上的值（,s,上的动点）,如弦振动，长为 的弦两端固定，则边界条件为：,（,2,）第二类边界条件,已知导数,一维热传导（杆的导热），设杆的一端,x,a,绝热，则由外到内经过杆端的热量流速为零,2023/9/28265.4.2 边界条件(1)第一边界,2024/11/15,27,因,K,，,S,是常数，故,对于二维、三维应以边界的外法向导数表述,5.4.2,边界条件,（,2,）第二类边界条件,已知导数,（,3,）第三类边界条件,混合边界条件,给出边界上函数值与其法向导数构成的线性关系。如一维导热，杆端,x=a,处自由冷却，环境介质温度为,u,0,，则,2023/9/28275.4.2 边界条件（2）第二类边界,2024/11/15,28,（,3,）第三类边界条件,混合边界条件,杆端散发出的热流效率与端点温度与介质温度之差成正比,可改写,5.4.2,边界条件,对于长为 的杆两端自由冷却,第三类边界条件的一般形式,2023/9/2828（3）第三类边界条件混合边界条件5,2024/11/15,29,5.5,线性迭加原理,在讲如何用分离变量法求解偏微分方程的定解问题前，先介绍一下线性偏微分方程解的迭加原理。,定义线性偏微分算子,L,为,线性偏微分方程的一般形式,齐次线性偏微分方程,2023/9/28295.5 线性迭加原理 在,2024/11/15,30,设函数 是齐次线性偏微分方程 的特解，若级数 可逐项求偏微分，则该级数也是齐次线性偏微分方程,的解。,5.5,线性迭加原理,线性迭加原理,2023/9/2830 设函数,2024/11/15,31,5.6,分离变量法,对于多个自变量的偏微分方程定解问题的求解，,在可能的情况下，我们总设法使自变量的个数减少,。分离变量法就是基于这种想法产生的。,分离变量法也叫傅立叶方法，它利用变量分离形式的解法，将求解偏微分方程的定解问题化为求解常微分方程的固有值问题，,步骤是先找出一些满足边界条件的特解，然后利用迭加原理，作出这些解的线性组合，从而得到定解问题的解答。,分离变量法对定解条件尤其是边界条件的要求比较苛刻，一般只涉及较为规则的边界问题。下面通过各种例题来介绍分离变量法的具体应用。,2023/9/28315.6 分离变量法 对,2024/11/15,32,5.6,分离变量法,例,这是齐次方程，齐次边界条件的定解问题。,解,:,u(x,t),是其一个解函数。假设函数可以表示为各个自变量单元函数的乘积，代入方程