第七章：偏微分方程第七章：偏微分方程一、 几个基本概念几个基本概念例：例：1、方程的阶数、方程的阶数方程中出现的最高阶导数的阶数即为方程的阶数方程中出现的最高阶导数的阶数即为方程的阶数一阶一阶二阶二阶 2、线性、非线性、拟线性、线性、非线性、拟线性 方程经过有理化并消去分式后，若方程中没有未知函方程经过有理化并消去分式后，若方程中没有未知函数及其偏导数的乘积或幂等非线性项，称该方程为线性数及其偏导数的乘积或幂等非线性项，称该方程为线性线性：线性：拟线性：拟线性： 在非线性方程中，如果未知函数的所有最高阶导数在非线性方程中，如果未知函数的所有最高阶导数不是非线性，则称此方程为拟线性不是非线性，则称此方程为拟线性 完全非线性：完全非线性：除拟线性之外的非线性方程除拟线性之外的非线性方程 二阶，线性二阶，线性二阶，拟线性二阶，拟线性二阶，完全非线性二阶，完全非线性 3、齐次、非齐次、齐次、非齐次 不含有未知函数及其偏导数的项不含有未知函数及其偏导数的项自由项：自由项：自由项为零的方程自由项为零的方程齐次方程：齐次方程： 自由项不为零的方程自由项不为零的方程非齐次方程：非齐次方程： 非齐次非齐次齐齐 次次 二、 二阶线性偏微分方程的分类二阶线性偏微分方程的分类，设设则二阶线性偏微分方程一般可表示为：则二阶线性偏微分方程一般可表示为： A,B,C,D,E,G,f 都是都是x,y的函数的函数 二阶线性方程的分类二阶线性方程的分类 ：：为方程中自变量域内任意一点，为方程中自变量域内任意一点，则分类判别条件如下：则分类判别条件如下：（（ⅰ）若在点）若在点M处有：处有： 则方程在该点处为双曲线型则方程在该点处为双曲线型 例如例如 （（ⅱ）若在点）若在点M处有：处有： 则方程在该点处为抛物线型则方程在该点处为抛物线型 例如例如 波动方程波动方程 热传导方程热传导方程 （（ⅲ）若在点）若在点M处有：处有： 则方程在该点处为椭圆型则方程在该点处为椭圆型 例如例如 例例1 1：判别下列方程的分类：判别下列方程的分类 当当 ，即，即 异号，异号，M点在二、四象限内，点在二、四象限内，该区域内方程是双曲型该区域内方程是双曲型当当 ，即，即 同号，同号，M点在一、三象限内，点在一、三象限内，该区域内方程是椭圆型该区域内方程是椭圆型当当x或或y为零，方程在为零，方程在x或或y轴上是抛物型轴上是抛物型 拉普拉斯方程拉普拉斯方程 三、三种典型方程的建立三、三种典型方程的建立建立理论模型的原则步骤建立理论模型的原则步骤 ①① 抽象出系统的物化模型并简化、假设抽象出系统的物化模型并简化、假设②② 确定输入、输出变量和模型参数，建立数学模型确定输入、输出变量和模型参数，建立数学模型③③ 模型模型 求解求解④ ④ 检验和修正所得的模型检验和修正所得的模型 1 1、均匀弦的微小横振动方程的建立、均匀弦的微小横振动方程的建立 设有一根均匀柔软的细弦，张紧后两端固定如图，设有一根均匀柔软的细弦，张紧后两端固定如图，给弦以扰动，使其产生振动。
给弦以扰动，使其产生振动确定弦上各点振动规律，即确定弦上各点振动规律，即 确定位移确定位移 满足的满足的方程 解：取一小微元段解：取一小微元段 ，分析受力情况，分析受力情况x 方向方向：u 方向方向：ABT1T20xx＋＋ΔΔxxuαα1 1αα2 2G(1)(2) 简化假设简化假设 （（1）弦是均匀的，所以质量均布，设单位长度弦的质量为）弦是均匀的，所以质量均布，设单位长度弦的质量为 kg/m （（5）弦是绝对柔软的，不能抗弯，因此弦上各点张力与该点切）弦是绝对柔软的，不能抗弯，因此弦上各点张力与该点切 线方向一致线方向一致2）弦的细小的，自重相比于张力显得很小，可以忽略弦的细小的，自重相比于张力显得很小，可以忽略3 3）振动方向与弦长方向相垂直，且振动保持在一固定平面内）振动方向与弦长方向相垂直，且振动保持在一固定平面内 （（4 4）振动是微小的，即弦上各点位移及弦的弯曲斜率很小）振动是微小的，即弦上各点位移及弦的弯曲斜率很小 微元段的质量：微元段的质量：G = 0x 方向无运动，即：方向无运动，即： 无外力的均匀弦微小横振动方程无外力的均匀弦微小横振动方程—— 齐次一维波动方程齐次一维波动方程x 方向方向：u 方向方向： ABT1T20xx＋＋ΔΔxxuαα1 1αα2 2F(x,t)x 方向方向：u 方向方向：假设振动过程中，假设振动过程中，除了张力外，除了张力外，还有其它外力作用，还有其它外力作用，设单位长度弦上的横向外力为设单位长度弦上的横向外力为弦的强迫振动方程弦的强迫振动方程 ——非齐次一维波动方程非齐次一维波动方程 2、热传导方程、热传导方程一根长为一根长为l的的均匀细杆均匀细杆侧面是绝热，横截面积足够小以至侧面是绝热，横截面积足够小以至在任何时刻都可以把断面上所有点的温度看作是相同。
在任何时刻都可以把断面上所有点的温度看作是相同 设杆的截面积为设杆的截面积为S，比热为，比热为C，导热系数为，导热系数为k，密度为，密度为ρ 试确定温度分布函数试确定温度分布函数 满足的方程满足的方程 x=0x=lxx+ΔΔx0热量衡算：输入－输出＋产生＝累计热量衡算：输入－输出＋产生＝累计 一维热传导方程一维热传导方程: : x处输入的热速率处输入的热速率: : x＋＋Δx处输出的热速率处输出的热速率: : 微元段累计的热速率微元段累计的热速率: : 一维热传导方程一维热传导方程一维波动方程一维波动方程（无外力）（无外力）（有外力）（有外力）（无内热源）（无内热源） 三维热传导方程三维热传导方程 若物体内部有一个热流，若物体内部有一个热流，T(x,y,z,t)其分布函数为其分布函数为M(x,y,z) 3、稳态方程（拉普拉斯方程）、稳态方程（拉普拉斯方程）热传导持续进行下去，如果达到稳定状态，温度的空间热传导持续进行下去，如果达到稳定状态，温度的空间分布不再变动，即分布不再变动，即 则方程则方程变为变为稳态浓度分布方程稳态浓度分布方程 稳态温度分布方程稳态温度分布方程 三种典型方程三种典型方程 是方程是方程其中其中f 是任意函数是任意函数 ,如如的通解的通解可以验证可以验证对于方程：对于方程：可以验证可以验证为其通解为其通解泛定方程泛定方程 暂停：休息暂停：休息 四、定解条件和定解问题四、定解条件和定解问题1、初始条件、初始条件（（I.C.)：： 初始时刻的状态初始时刻的状态 （一）定解条件（一）定解条件例例1 对于弦的微小横振动问题，假设初始速度为零，对于弦的微小横振动问题，假设初始速度为零， 初始位移符合正弦函数初始位移符合正弦函数 初始条件给定的是整个系统的状态，而不是某个局部初始条件给定的是整个系统的状态，而不是某个局部（如入口、出口等）的状态。
如入口、出口等）的状态特别注意：特别注意： 应是位置坐标的函数应是位置坐标的函数 （（1）第一类边界条件）第一类边界条件—已知函数已知函数直接给出未知函数直接给出未知函数 在边界在边界 上的值上的值 （狄里赫利（（狄里赫利（DirichletDirichlet）））） 例例2 一根弦长为一根弦长为l ，两端固定进行微小横振动，两端固定进行微小横振动 ，建立，建立其边界条件其边界条件2 2、边界条件（、边界条件（B.C.B.C.）） 例例3 细杆导热问题中，杆长细杆导热问题中，杆长l，两端分别保持温度，两端分别保持温度T1和和T2 ，建立边界条件建立边界条件 （（2 2）第二类边界条件）第二类边界条件—已知导数已知导数（牛曼（（牛曼（NoumannNoumann）条件））条件） 例例4 细杆导热问题中，杆长细杆导热问题中，杆长l，一端绝热，另一端有恒，一端绝热，另一端有恒定热流定热流q输入，试建立边界条件输入，试建立边界条件x=0x=l 热流输出怎样？？热流输出怎样？？ （（3 3）第三类边界条件）第三类边界条件——混合边界条件混合边界条件给出边界上函数值与其法向导数构成的线性关系给出边界上函数值与其法向导数构成的线性关系 （（RobinRobin条件）条件） 例例4 细杆导热问题中，杆长细杆导热问题中，杆长l，一端温度为，一端温度为T0，另一端，另一端与温度为与温度为T1的环境进行对流热交换，试建立边界条件。
的环境进行对流热交换，试建立边界条件x=0x=l （（4 4）积分）积分—微分边界条件微分边界条件 （（5）衔接条件）衔接条件内外层壁：内外层壁： （二）定解问题（二）定解问题 初值问题：初值问题：只有初始条件，没有边界条件的定解问题只有初始条件，没有边界条件的定解问题 边值问题：边值问题： 没有初始条件，只有边界条件的定解问题没有初始条件，只有边界条件的定解问题 混合问题：混合问题： 既有初始条件又有边界条件的定解问题既有初始条件又有边界条件的定解问题 定解问题定解问题 泛定方程（泛定方程（P.D.E) P.D.E) 定解条件定解条件 初始条件（初始条件（B.C.））边界条件边界条件（（I.C.）） 五、线性迭加原理五、线性迭加原理 所谓迭加：即几种不同因素综合作用于系统，产生所谓迭加：即几种不同因素综合作用于系统，产生的效果等于各因素独立作用产生的效果之总和的效果等于各因素独立作用产生的效果之总和 线性迭加原理：线性迭加原理： 则级数则级数 也是该方程之解也是该方程之解 假设函数假设函数是线性齐次微分方程是线性齐次微分方程的特解。
的特解i=1，，2……）） 六、分离变量法六、分离变量法 例例1 1、设两端固定的有界弦微小自由横振动过程，、设两端固定的有界弦微小自由横振动过程， 初始位移为初始位移为φ(x)，初始速度为，初始速度为ψ（（x））解：（解：（1 1）分离变量）分离变量 设设 （（1 1）） （（2 2）） （（3 3）） （（4 4）） 到到33页页 代入方程式（代入方程式（1 1）得）得 分离变量得分离变量得 = =设设 于是可得到两个常微分方程于是可得到两个常微分方程 到到35页页 = =+ 对于二阶线性常系数常微分方程：对于二阶线性常系数常微分方程： 先求特征方程的解：先求特征方程的解： 分三种情况分三种情况 （（ⅰⅰ）特征方程的根）特征方程的根r r1 1,r,r2 2为实数时为实数时 复复 习习 （（ⅱⅱ）特征方程的根）特征方程的根r r1 1=r=r2 2=r=r为重实数时为重实数时 （（ⅲⅲ）特征方程的根为复数，即）特征方程的根为复数，即 r r1 1=α+β=α+βi r r2 2= = α - β i 返回返回31页页 由边界条件式（由边界条件式（2 2）知）知 因为因为 所以只能所以只能 （（5 5）） （（6 6）） （（ⅰⅰ）设）设λλ>0>0，则方程式（，则方程式（5 5）通解为）通解为（（5 5）） 由由 由由 （（2 2）解定解问题）解定解问题 ——解固有值问题解固有值问题 到到37页页 解上述线性代数方程组得：解上述线性代数方程组得： 即即 所以所以 （（ⅱ）设）设λ =0，方程（，方程（5）通解为）通解为由由 由由 即即 所以所以 （（ⅲⅲ）设）设λ <0<0，，令令 ，，此时通解为此时通解为由由 由由 因为因为B不能为零，所以只能不能为零，所以只能所以所以－称为固有值（或本征值）－称为固有值（或本征值） ------称为固有函数（或本征函数）称为固有函数（或本征函数） （（7 7）） 返回返回36页页返回返回54页页 （（3 3）求解不构成本征问题的常微分方程的通解）求解不构成本征问题的常微分方程的通解（（6 6）） 将将 代入上式代入上式 （（8）） 将式（将式（7 7）与式）与式（（8））相乘，得到一组特解相乘，得到一组特解 其中，其中， 是任意常数是任意常数（（8））（（7）） （（4 4）应用线性叠加原理求）应用线性叠加原理求（（5 5）由傅里叶级数确定系数）由傅里叶级数确定系数（（9）） 暂停：休息暂停：休息 复习复习f(x)在〔－在〔－l，，l〕区间上的傅里叶展开式：〕区间上的傅里叶展开式： 在〔在〔0，，l〕区间上的只有余弦项的傅里叶展开式：〕区间上的只有余弦项的傅里叶展开式：在〔在〔0，，l〕区间上的只有正弦项的傅里叶展开式：〕区间上的只有正弦项的傅里叶展开式： 傅里叶展开是基于三角函数的正交性傅里叶展开是基于三角函数的正交性 根据正交性可直接确定系数：根据正交性可直接确定系数：0000（（10）） 至此，至此， （（9）～（）～（11）即构成了例）即构成了例1中定解问题的解中定解问题的解（（11）） 总结分离变量求解的关键点总结分离变量求解的关键点1、假设、假设代入代入齐次方程齐次方程进行分离变量；进行分离变量；2、分离、分离齐次边界条件齐次边界条件，组构本征值问题，并求解确定，组构本征值问题，并求解确定 本征值本征函数；本征值本征函数；3、求解另外一个常微分方程；、求解另外一个常微分方程；4、用线性叠加原理写出问题的解；、用线性叠加原理写出问题的解；5、代入初始条件，并用傅里叶展开法确定解中的常数。
代入初始条件，并用傅里叶展开法确定解中的常数 例例 2 2 一维导热问题（第三边值条件）一维导热问题（第三边值条件）长为长为l 的均匀细杆，其侧面（圆弧面）绝热，杆的一的均匀细杆，其侧面（圆弧面）绝热，杆的一端保持在端保持在0℃状态下，另一端则与温度为状态下，另一端则与温度为0℃环境介质环境介质进行自由热交换假设初始时刻温度分布为进行自由热交换假设初始时刻温度分布为φ(x)试确定杆上各点温度随时间变化规律确定杆上各点温度随时间变化规律 （（1））（（2））（（3））（（4））到第到第52页页 （（1 1）分离变量，设）分离变量，设 代入方程（代入方程（1 1）得：）得： 于是得：于是得： （（5 5）） （（6 6）） 比较波动问题得到两个常微分方程比较波动问题得到两个常微分方程 ：：返回上一页返回上一页 由边界条件式（由边界条件式（2 2）（）（3 3）知）知（（2 2）） 解本征值问题解本征值问题 （（7 7）） （（8 8）） 式（式（5）与式（）与式（7）、式（）、式（8）构成本征值问题）构成本征值问题 同前述同样方法讨论本征值同前述同样方法讨论本征值λλ的三种取值情况的三种取值情况 经讨论仅当经讨论仅当λλ<0时才有非零解时才有非零解 设设 于是方程式（于是方程式（5 5）通解为）通解为 由式（由式（7 7）得）得 由式（由式（8 8）得）得 则则 即即 得本征值得本征值(n =1，，2，，3… )得本征函数得本征函数 （(3(3）） 将本征值代入式（将本征值代入式（6 6）并求解）并求解（（4 4）由线性迭加原理得定解问题级数形式解）由线性迭加原理得定解问题级数形式解(Cn =BnAn) （（5 5）确定系数）确定系数由初始条件由初始条件由本证函数的正交性由本证函数的正交性 例例 3 3、求解稳态问题、求解稳态问题 解：（解：（1 1）分离变量）分离变量 设设 （（1 1）） （（2 2）） （（3 3）） 能否用分离变量法能否用分离变量法 求解？？？求解？？？ 能！！能！！ 代入方程式（代入方程式（1 1）得）得 分离变量得分离变量得 = =设设 于是可得到两个常微分方程于是可得到两个常微分方程 —由边界条件（由边界条件（3）） （（5 5）） （（6 6）） （（5 5）） （（2 2）解本征值问题）解本征值问题 凭经验当凭经验当λ <0<0，时有非零解，时有非零解令令 由由 由由 因为因为B不能为零能不能为零能 所以所以（（7 7）） （（3 3）求解不构成本征问题的常微分方程的通解）求解不构成本征问题的常微分方程的通解（（6 6）） 将将 代入上式代入上式 （（8）） （将式（将式（7 7）与式）与式（（8））相乘，得到一组特解相乘，得到一组特解 其中，其中， 是任意常数是任意常数 （（4 4）应用线性叠加原理求）应用线性叠加原理求（（5 5）由傅里叶级数确定系数）由傅里叶级数确定系数（（9）） （（10））（（9）） n 为偶数时为偶数时n 为奇数时为奇数时 1 1、只能用分离变量法直接求解齐次方程、、只能用分离变量法直接求解齐次方程、 齐次边界定解问题齐次边界定解问题2 2、本征值和本征函数由边界条件类型确定、本征值和本征函数由边界条件类型确定3 3、另外一个二阶常微分方程的形式及解、另外一个二阶常微分方程的形式及解 由泛定方程的类型确定由泛定方程的类型确定总总 结结 （（b b））（（c c））（（d d））（（a a）） 一维波动问题一维波动问题一维热传导问题一维热传导问题二维稳态问题二维稳态问题 暂停：休息暂停：休息 七、非齐次边界的处理七、非齐次边界的处理方法：方法： 将定解问题（将定解问题（A A）转化为齐次边界问题）转化为齐次边界问题 （（A））设设使使 满足相应的齐次边界条件，即满足相应的齐次边界条件，即（（1））（（2））（（3））（（4））（（5）） 则则设设由条件式（由条件式（6 6）确定：）确定： （（6））（（7））则则由条件式（由条件式（7 7）确定：）确定： 将以上结果代入将以上结果代入 （（B B）） （a） （（a a））（（b b））（（c c））（（d d））则则V（（x,t）满足的定解问题为）满足的定解问题为：： 设设 例例4 4 设一均匀细杆，初始时全杆有一均一温度设一均匀细杆，初始时全杆有一均一温度 ，然后，然后使其一端保持不变温度使其一端保持不变温度 ，另一端则有恒定的热流，另一端则有恒定的热流 输入，试求温度分布规律。
输入，试求温度分布规律 解：问题可由下述数学形式表述解：问题可由下述数学形式表述 （（1））（（2））（（3）） 设设 使使 则则 设设 则由式（则由式（6 6）求得）求得（（4））（（5））（（6））（（7）） 则则 代入原问题得代入原问题得 设设 W（（x,t）） 由经验，知由经验，知 唯唯 时有解时有解设设 则则 由边界，得由边界，得 所以所以 （=0，1，2……）（（n = 0，，1，，2……）） 于是固有函数为于是固有函数为 最后由初始条件定常数最后由初始条件定常数 暂停：休息暂停：休息 八、非齐次的泛定方程八、非齐次的泛定方程 设设 例例 5 5 两端固定的弦的强迫振动问题两端固定的弦的强迫振动问题 =0设设 解解: : 而而 （（1 1）） （（2 2）） （（3 3）） （（4 4）） （（5 5）） 转到转到82页页 （（b b））（（c c））（（d d））（（a a））返回上一页返回上一页 将式（将式（4），式（），式（5）代入式（）代入式（1）得）得 所以所以 将式（将式（3 3））代入所设：代入所设： 这里这里 （（5 5）） （（6 6）） （（7 7）） （（8 8）） 根据傅立叶展开根据傅立叶展开 于是有：于是有： 因已设因已设 设设 则则 （（6 6）） 对方程式（对方程式（6 6）两边拉氏变换）两边拉氏变换 由卷积定理知由卷积定理知= 暂停：休息暂停：休息 例例 6 6 用本征函数系展开法解非齐次热传导问题用本征函数系展开法解非齐次热传导问题 解：解： （（1 1）） （（2 2）） （（3 3）） 将上两式代入方程式（将上两式代入方程式（1 1）得）得＝＝0由初始条件式（由初始条件式（3 3）） 于是方程（于是方程（5）解为）解为所以：所以： 最终解为最终解为 （（4 4）） （（5 5）） 于是方程（于是方程（4）解为）解为 ★★ 重点掌握的内容：重点掌握的内容： 拉氏变换的定义和前拉氏变换的定义和前5 5条性质；三种常用求条性质；三种常用求解逆变换的方法；应用－求解常微分方程，积分解逆变换的方法；应用－求解常微分方程，积分- -微分方程的方法。
微分方程的方法★★ 本章主要内容：本章主要内容： 1、拉氏变换的概念和性质；、拉氏变换的概念和性质； 2、求拉氏逆变换的方法；、求拉氏逆变换的方法； 3、拉氏变化的应用，如何求解微分、积分方、拉氏变化的应用，如何求解微分、积分方 程、差分方程等程、差分方程等第五章：拉普拉斯变换第五章：拉普拉斯变换 第六章：场论初步第六章：场论初步★★ 本章重点内容：本章重点内容： 1、方向导数、通量、环量的概念及计算；、方向导数、通量、环量的概念及计算； 2、梯度、散度、旋度的概念及计算；、梯度、散度、旋度的概念及计算； 3、散度定理、斯托克斯定理；、散度定理、斯托克斯定理； 4、无源场、有势场（势函数）、调和场（调和、无源场、有势场（势函数）、调和场（调和函数）概念、性质及计算；函数）概念、性质及计算； 5、会用用场论方法简单的建模，重点是课上讲、会用用场论方法简单的建模，重点是课上讲的例题 第七章：偏微分方程第七章：偏微分方程★★ 本章重点内容：本章重点内容： 1、会用分离变量法求解二阶齐次方程、齐次边、会用分离变量法求解二阶齐次方程、齐次边界问题；界问题； 2、会建立简单的方程（类似三种典型方程）和、会建立简单的方程（类似三种典型方程）和定解条件；定解条件； 3、掌握方程的分类和线性性判别；、掌握方程的分类和线性性判别； 4、了解非齐次边界处理方法和非齐次方程的本、了解非齐次边界处理方法和非齐次方程的本征函数求解方法征函数求解方法 2012级化工数学期末考试信息级化工数学期末考试信息时间：时间：2014.6.30（第（第19周周一）周周一）？？:00—？？:00地点：地点：？楼？？楼？ （（ⅰⅰ）设）设λλ>0>0，则方程式（，则方程式（5 5）通解为）通解为由由 由由 解上述线性代数方程组得：解上述线性代数方程组得： 即即 所以所以 （（ⅱ）设）设λ =0，方程（，方程（5）通解为）通解为由由 由由 即即 （（ⅲⅲ）设）设λ <0<0，，令令 ，，此时通解为此时通解为由由 由由 所以综合有：所以综合有：------称为固有函数（或本征函数）称为固有函数（或本征函数） （（7 7）） 返回返回88页页 。