数学物理方法7数学物理方程的定解问题课件.ppt

1、2022-8-121参考书：参考书：R.Haberman著，郇中丹等译，著，郇中丹等译，实用偏微分方实用偏微分方程程 （原书第四版），机械工业出版社，（原书第四版），机械工业出版社，200720072022-8-1222022-8-123 物理现象物理现象 物理量物理量u 在空间和时间中的变在空间和时间中的变化规律，即物理量化规律，即物理量u在各个地点和各个时刻所取的值在各个地点和各个时刻所取的值之间的联系。之间的联系。数学语言描述 数学物理方程：从物理问题中导出的函数方程，数学物理方程：从物理问题中导出的函数方程，特别是偏微分方程和积分方程。特别是偏微分方程和积分方程。：二阶线性偏微分方程。

2、：二阶线性偏微分方程。例：牛顿第二定律反映的是力学现象的普遍规律，例：牛顿第二定律反映的是力学现象的普遍规律，跟具体条件无关。跟具体条件无关。2022-8-124三类典型的数学物理方程三类典型的数学物理方程双曲型方程双曲型方程波动方程为代表波动方程为代表抛物型方程抛物型方程扩散方程为代表扩散方程为代表椭圆型方程椭圆型方程泊松方程为代表泊松方程为代表退化为拉普拉斯方程退化为拉普拉斯方程2(,)ttxxua uf x t2(,)uauF x y z tt2(,)auF x y z0F0 u2022-8-1255 1 边界问题-边界条件体现边界状态的数学方程称为边界条件体现边界状态的数学方程称为边界

3、条件2 历史问题-初始条件体现历史状态的数学方程称为初始条件体现历史状态的数学方程称为初始条件例：一个物体做竖直上抛，一个物体斜抛。不同的初始条件例：一个物体做竖直上抛，一个物体斜抛。不同的初始条件 不同的运动状态，但都服从牛顿第二定律。不同的运动状态，但都服从牛顿第二定律。三、定解问题三、定解问题 在给定的边界条件和初始条件下，根据已知的物理规律，在在给定的边界条件和初始条件下，根据已知的物理规律，在给定的区域里解出某个物理量给定的区域里解出某个物理量u，即求即求u(x,y,z,t)。：边界条件和初始条件的总体。它反映了问题的边界条件和初始条件的总体。它反映了问题的 特殊性，即个性。特殊性，

4、即个性。：不带有边界和初始条件的方程称为泛定方程。不带有边界和初始条件的方程称为泛定方程。它反映了问题的共性。它反映了问题的共性。二、定解条件二、定解条件2022-8-1266具体问题求解的一般过程：具体问题求解的一般过程：1 1、根据系统的内在规律列出泛定方程、根据系统的内在规律列出泛定方程客观规律客观规律2 2、根据已知系统的边界状况和初始状况列出边界条件和、根据已知系统的边界状况和初始状况列出边界条件和 初始条件初始条件求解所必须的已知条件求解所必须的已知条件3 3、求解方法、求解方法 行波法、行波法、分离变量法分离变量法、积分变换法、格林、积分变换法、格林函数法和变分法函数法和变分法2

5、022-8-127建模步骤：建模步骤：（1 1）明确要研究的物理量是什么？）明确要研究的物理量是什么？从所研究的系统中从所研究的系统中，分析邻近部，分析邻近部分与它的相互作用。分与它的相互作用。（2 2）研究物理量遵循哪些物理规律？）研究物理量遵循哪些物理规律？（3 3）按物理定律写出数理方程（泛定方程）。）按物理定律写出数理方程（泛定方程）。2022-8-128 （如图）（如图）：沿：沿x轴绷紧的均匀柔软的细轴绷紧的均匀柔软的细弦，在平衡位置弦，在平衡位置(x轴轴)附近产生振幅极小的横向振动附近产生振幅极小的横向振动 建立与细弦上各点的振动规律相应的方程建立与细弦上各点的振动规律相应的方程

6、(1)弦弦不振动时不振动时静止于静止于x轴轴;(2)用用表示表示t时刻时刻弦上弦上x在在垂直于垂直于x轴轴方向上方向上的横向位移（的横向位移（偏离偏离）情况）情况 弦的横振动弦的横振动2022-8-129 选取不包括端点的一微元选取不包括端点的一微元x,x+dx弧弧B段段作为研究对象作为研究对象.研究对象研究对象:(4)(4)设单位长度上弦受力设单位长度上弦受力F(x,t)，线力密度为：，线力密度为：假设与近似：假设与近似：(1)(1)弦是柔软的弦是柔软的 (不抵抗弯曲不抵抗弯曲),),张力沿弦的切线方向张力沿弦的切线方向 (2)(2)振幅极小振幅极小,张力与水平方向的夹角张力与水平方向的夹角

7、 1 1和和 2 2 很小，很小，仅考虑仅考虑 1 1和和 2 2的一阶小量，略去二阶小量的一阶小量，略去二阶小量 (3)(3)弦的重量与张力相比很小，可以忽略弦的重量与张力相比很小，可以忽略(,)(,)/f x tF x t 质量线密度质量线密度，u(x)u+duu0 1 2T2T1xx+dxFB2022-8-1210B段段弦的原长近似为弦的原长近似为dx.振动拉伸后：振动拉伸后：222d(d)(d)d1(d/d)dsxuxuxx u(x)u+duu0 1 2T2T1xx+dxBFB段段的质量的质量：弦长弦长dx，质量线密度质量线密度，则，则B段段质量质量 m=dx物理规律：物理规律：用牛顿

8、运动定律分析用牛顿运动定律分析B段段弦的受力及运动状态弦的受力及运动状态：22ddttufmmut牛顿运动定律：牛顿运动定律：2022-8-1211沿沿x-方向方向：弦横向振动不出现弦横向振动不出现x方向平移，方向平移，得得力平衡方程力平衡方程2211coscos0TT沿垂直于沿垂直于x-轴方向轴方向：由牛顿运动定律得由牛顿运动定律得运动方程运动方程2211sinsin(,)d(d)ttTTF x txx u 1 21 20,cos1.，11sintanxxxuux 22dsintanxxxu 在微小振动近似下：在微小振动近似下：由由（1）式，弦中各点的张力相等式，弦中各点的张力相等u(x)u

9、+duu0 1 2T2T1xx+dxBF21TT（1 1）（2 2）2211dsinsinxxxxxTTT uu 2022-8-1212d(,)(,)dxxxxxxxttuuTF x tTuF x tux 2/aT 波动方程：波动方程：波速波速a(,)(,)/f x tF x t d(,)d(d)xxxxxttT uuF x txx u（）受迫振动方程受迫振动方程2(,)ttxxua uf x t单位质量弦所受单位质量弦所受外力，线力密度外力，线力密度令令一维波动方程一维波动方程2022-8-121322222uuafgtx一维波动方程一维波动方程-非齐次方程非齐次方程222220uuatx-

10、齐次方程齐次方程忽略重力和外力作用：忽略重力和外力作用：如考虑弦的重量：如考虑弦的重量：u(x)u+uu0 1 2T2T1xx+xBFdg x 沿沿x-方向，不出现平移方向，不出现平移2211coscosTT 沿垂直于沿垂直于x-轴方向轴方向2211sinsin(,)dd(d)ttTTF x txg xx u （1 1）（2 2）d02211ddsinsindddxxxxxxxuTTT uuT uTx 讨论：讨论：dd(,)dd(d)dttuTF x txg xx ux 2022-8-1214 扩散现象扩散现象：不均匀时，将出现物质从高浓度处不均匀时，将出现物质从高浓度处向低浓度处转移的现象，

11、称之为向低浓度处转移的现象，称之为扩散扩散。扩散定律即裴克定律：扩散定律即裴克定律：这是一条实验定律这是一条实验定律 数学建模：数学建模：建立空间各点建立空间各点的方程的方程 物理规律：物理规律：以扩散定律和粒子数守恒定律为研究基础以扩散定律和粒子数守恒定律为研究基础粒子数守恒粒子数守恒定律：定律：单位时间内流入某一体积的单位时间内流入某一体积的粒子数与流出这一体积的粒子数之差等于此体积粒子数与流出这一体积的粒子数之差等于此体积内的单位时间内粒子数的增加量内的单位时间内粒子数的增加量处理方法：处理方法：在在浓度不均匀的无源空间，划浓度不均匀的无源空间，划出出为研究对象，分析浓度变化规律为研究对

12、象，分析浓度变化规律。2022-8-1215浓度不均匀浓度不均匀:；u扩散流强弱（强度）扩散流强弱（强度）：用单位时间通用单位时间通过单位面积的物质的量过单位面积的物质的量 表示；表示；q扩散（裴克）实验定律：扩散（裴克）实验定律：qDuuuuDijkxyzxyz),(zyx(,)xdx ydy zdzdxdydzxqdxxq扩散系数扩散系数设定：设定：处理方法：处理方法：在在浓度不均匀的无源空间，划浓度不均匀的无源空间，划出出为研究对象，为研究对象，分析浓度变化规律分析浓度变化规律。扩散流强度与浓度梯度间关系：扩散流强度与浓度梯度间关系：采用裴克实验定采用裴克实验定律确定律确定体元体元内粒子

13、数：内粒子数：zyxtxyxuddd,2022-8-1216xyz),(zyx(,)xdx ydy zdzdxdydzxqdxxq单位时间沿单位时间沿x-方向净流入量方向净流入量()x dxxqqqdydzdxdydzx 同理沿同理沿y 和沿和沿z方向净流入量方向净流入量由粒子数守恒定律，有由粒子数守恒定律，有单位时间内向单位时间内向V的净流入量的净流入量d d d,qx y zyd d dqx y zzd d dd d dd d dqqqx y zx y zx y zxyz单位时间内单位时间内V内粒子数的增加量内粒子数的增加量d d dd d dd dd d ddqqqx y zx y zx

14、 y zxuxtyyzzd d dux y zt2022-8-1217如果如果扩散是均匀的扩散是均匀的,即即D是一常数，则可以令是一常数，则可以令D=a2 2，则有则有23(,)tuauF x y z t222222322220 xuuuuaauuatztyuut0()()()uuuuDDDtxxyyzz代入扩散定律代入扩散定律 如果所如果所研究的空间存在扩散源，源强度与研究的空间存在扩散源，源强度与u(x,y,z,t)无关，无关，且为且为F(x,y,z),这时扩散方程修改为这时扩散方程修改为 如果所如果所研究的空间存在源，源强度与研究的空间存在源，源强度与u(x,y,z,t)成正比，成正比，

15、即即F(x,y,z)=b2u(x,y,z)这时扩散方程修改为这时扩散方程修改为223(,)tuaub u x y z txqDuuqDx qqqudxdydzdxdydzdxdydzdxdydzxyzt2022-8-121823 auF密度场：密度在空间的分布构成一个标量场。密度场：密度在空间的分布构成一个标量场。23uauFt有扩散源时系统的密度场满足有扩散源时系统的密度场满足非齐次扩散方程非齐次扩散方程稳定状态：密度稳定状态：密度u 不随时间变化，则不随时间变化，则泊松方程泊松方程无扩散源：无扩散源：F=030u拉普拉斯方程拉普拉斯方程2022-8-1219所要研究的物理量：所要研究的物理

16、量：温度温度 (,)u x y z t物理规律物理规律：采用傅里叶实验定律：采用傅里叶实验定律热传导现象热传导现象：当导热介质中各点的温度分布不均匀时，有：当导热介质中各点的温度分布不均匀时，有热量从高温处流向低温处。热量从高温处流向低温处。数学建模数学建模：傅里叶定律：傅里叶定律：qku 温度不均匀温度不均匀:用用温度梯度温度梯度 表示表示；u传热的强弱即传热的强弱即热流强度热流强度：用单位时间内通过单用单位时间内通过单位面积的热量位面积的热量 表示；表示；q设定：设定：nuqknn 沿曲面法向流出热量：沿曲面法向流出热量：2022-8-1220nnqq 有限时间内有限时间内即时刻即时刻t1

17、 1到到t2 2通过闭曲面通过闭曲面S流入流入V的热量为的热量为 211ddttSQk uSt 高斯公式（矢量散度的体积分等于该矢量对包围该体积的面积分）高斯公式（矢量散度的体积分等于该矢量对包围该体积的面积分）2211213d dd dttttVVQku V tku V t ddd ddnnQuqSdtqSdtkS tn 热场MSdSV n处理方法：处理方法：在温在温度不均匀的无源空间，划出任一封闭曲面度不均匀的无源空间，划出任一封闭曲面S包包围的围的体积元体积元V（如图）（如图）。在在S 上选取任一足够小的微面元上选取任一足够小的微面元dS，在此在此面元范围内热流强度近似为常量。面元范围内

18、热流强度近似为常量。d dd duu nnku nS tk uS t ddSVk uSk uV 那么在那么在dt时间内从时间内从dS流入流入V的的为为(向为正向为正)：n nuqknn 2022-8-12212113d dttVQku V t 1(,)uy z tx2(,)uy z tx221(,)(,)dVQcu x y z tu x y z tV 12QQ流入的热量导致流入的热量导致V内的温度发生变化内的温度发生变化 22113d dd dttttVVukucV tV tt 3kcuut23uauft流入的热量：流入的热量：温度发生变化需要的热量温度发生变化需要的热量(c比热容，比热容，质

19、量密度质量密度)：21d dttVuct Vt 21d dttVucV tt 23auMSSV n如果物体内有热源，则温度满足如果物体内有热源，则温度满足非齐次热传导方程非齐次热传导方程2022-8-1222介质方程介质方程:DE ijk xyz其中其中:0)(2)ED （1高斯定理高斯定理:环路定理环路定理:物理规律物理规律：由电磁学可知，静电场满足静电学高斯定理、由电磁学可知，静电场满足静电学高斯定理、环路定理和介质方程。环路定理和介质方程。数学建模数学建模：建立电势：建立电势u(x,y,z)与电荷密度与电荷密度(x,y,z)的关系。的关系。由电场的高斯定理由电场的高斯定理/DE 物理问题

20、物理问题：在介电常数为在介电常数为的介质的介质空间空间，存在电荷分布，存在电荷分布(x,y,z)激发电场激发电场 形成电势分布形成电势分布u(x,y,z)。2022-8-12230若空间若空间无电荷，即电荷密度无电荷，即电荷密度，上式成为，上式成为 称这个方程为拉普拉斯方程拉普拉斯方程.由电场的环路定理，可知静电场是一个保守场由电场的环路定理，可知静电场是一个保守场.由保守由保守场的性质，场的性质，,且且，即，即：uE 3()Euuu3/u 22232220uuuuxyz进一步对电场取散度，有：进一步对电场取散度，有：泊松方程泊松方程 设电势为设电势为:u(x,y,z)。/E2022-8-12

21、24 7.1 3.4.2022-8-1225 数学物理方程的定解数学物理方程的定解 在给定的边界条件和初始条件下，根据已知的物理规律，在在给定的边界条件和初始条件下，根据已知的物理规律，在给定的区域里解出某个物理量给定的区域里解出某个物理量u，。：不带有边界和初始条件的方程称为泛定不带有边界和初始条件的方程称为泛定方程。方程。它反映了问题的共性。它反映了问题的共性。：边界条件和初始条件的总体。边界条件和初始条件的总体。它反映了问题的特殊性，即个性。它反映了问题的特殊性，即个性。2022-8-1226初始时刻的温度分布：初始时刻的温度分布：B B、热传导方程的初始条件、热传导方程的初始条件0(,

22、)|()tu x txC C、泊松方程和拉普拉斯方程的初始条件、泊松方程和拉普拉斯方程的初始条件A A、波动方程的初始条件波动方程的初始条件00|()()ttuxuxt描述系统的初始状态描述系统的初始状态系统各点的初位移系统各点的初位移系统各点的初速度系统各点的初速度波动方程含有时间的二阶导数，所以需二个初始条件波动方程含有时间的二阶导数，所以需二个初始条件 热传导方程含有时间的一阶导数，所以需一个初始条件热传导方程含有时间的一阶导数，所以需一个初始条件此类导方程不含时间的导数，所以不需要有初始条件此类导方程不含时间的导数，所以不需要有初始条件2022-8-1227和和 是空间坐标的函数是空间

23、坐标的函数(,)x y z(,)x y z 02 0222,(,)(),thlxxlu x thllxxll：初始条件给出系统在初始状态下物理量的分布，：初始条件给出系统在初始状态下物理量的分布，而不是某一位置处的情况。而不是某一位置处的情况。例例：一一根长为根长为l的弦，两端固定于的弦，两端固定于0和和l。在中点位置将弦沿着横向拉在中点位置将弦沿着横向拉开距离开距离h，如图所示，然后放手任其振动，试写出初始条件。如图所示，然后放手任其振动，试写出初始条件。l x l/2h解：解：初始时刻就是放手的那一瞬间，弦的形初始时刻就是放手的那一瞬间，弦的形状如图所示状如图所示,且弦处于静止状态，即有方

24、程且弦处于静止状态，即有方程00(,)ttu x t初始位移初始位移初始速度初始速度2022-8-1228 定义：定义：系统的物理量在边界上具有的情况。系统的物理量在边界上具有的情况。A.A.第一类（狄利克雷）边界条件第一类（狄利克雷）边界条件给出未知函数在给出未知函数在边界上的函数值边界上的函数值。例例2 2：两端固定的弦振动时的边界条件：两端固定的弦振动时的边界条件：00(,)xu x t0(,)x lu x t和和常见的线性边界条件分为三类：常见的线性边界条件分为三类：2022-8-1229例例3：细杆热传导：细杆热传导0 xlx 0(,)x lu x tu 细杆在细杆在x=l端的温度端

25、的温度随时间变化随时间变化,设温度变化规律为设温度变化规律为f(t),边边界的数理方程界的数理方程(,)()x lu x tf t细杆细杆x=l端的温度处于恒温状态端的温度处于恒温状态,边界的数理方程边界的数理方程000000,(,)(,)xyzu x y z tf xyzt边界2022-8-1230B.B.第二类（诺伊曼）边界条件第二类（诺伊曼）边界条件例例4 4：细杆热传导细杆热传导 我们用傅里叶我们用傅里叶(热传导热传导)定律来建立边界的数学物理方程定律来建立边界的数学物理方程.傅里叶实验定律傅里叶实验定律:单位时间内单位时间内,通过单位面积的热流为通过单位面积的热流为 0(,)xx a

26、ux t给出未知函数在边界上的给出未知函数在边界上的法线方向的导数之值法线方向的导数之值。000000,(,)(,)xyzu x y z tf xyztn边界细杆细杆x=a端点绝热的端点绝热的qku 设细杆沿设细杆沿x轴方向轴方向,则一维傅里叶实验定律改写为则一维傅里叶实验定律改写为 xxqku i其中其中u是所在位置处物体的是所在位置处物体的k是传热系数。是传热系数。(,)()xx aux tf t 细杆细杆x=a端点有热流端点有热流kf(t)流出的流出的2022-8-1231（3 3）.第三类（混合）边界条件第三类（混合）边界条件()qh un()xqh ui 一维傅里叶实验定律一维傅里叶

27、实验定律 xxqku i先引入两个基本物理定律：先引入两个基本物理定律：2022-8-1232例例5 5：写出写出导热细杆导热细杆l端端“自由自由”冷却的边界条件。冷却的边界条件。根据根据热传导定律热传导定律，在，在 x=l 处处：()xx lxx lkuqh u 0 xlx 负负x方向方向 n n正正x方向方向00()xxxkuh u)(/)(xx lHhuHuk 0()xxuHu 在在x=0 处：处：xxqku xxqku流出热流强度流出热流强度由由，此此流出热量与细杆和外界的温度差成流出热量与细杆和外界的温度差成正比，即正比，即即：即：讨论：如图情况讨论：如图情况xqu x2022-8-

28、1233例例6 6：细杆纵振动：细杆纵振动：端点与固定点弹性连接。应力为弹性力端点与固定点弹性连接。应力为弹性力胡克定律：胡克定律：xfESu弹性力：弹性力：fku 则在端点则在端点xkuESu0()xx lESuuk000000,()(,)边界xyzuuHf xyztn这些是最常见的线性边界条件，还有其它形式。这些是最常见的线性边界条件，还有其它形式。xlx k 系统中可能出现物理性质急剧变化的点系统中可能出现物理性质急剧变化的点(跃变点跃变点)。如两节。如两节具有不同的杨氏模量的细杆在具有不同的杨氏模量的细杆在 x=0 处连接，这一点就是跃变点。处连接，这一点就是跃变点。跃变点两边的物理过

29、程因此不同。但在跃变点，跃变点两边的物理过程因此不同。但在跃变点，某些物理量仍某些物理量仍然可以是连续的然可以是连续的，这就构成，这就构成衔接条件衔接条件。2022-8-1234120()sinsinF tTT这两个等式就是构成两段这两个等式就是构成两段衔接的衔接的是是。，横向力应与张力平衡：，横向力应与张力平衡：11022000sintan(,)sintan(,)xxuxtuxt 即即0000(,)(,)()xxTuxtTuxtF t ux0 x)(tF 1 20000()()u xu x位移极限值相同。位移极限值相同。0000()()xxuxux例例7 7横向力横向力F(t)集中作用于弦上

30、集中作用于弦上x0 0点点,使使x0点成为折点（如图）点成为折点（如图）。但二段是同一根弦，它们间相互关连。因此要建立此关系：但二段是同一根弦，它们间相互关连。因此要建立此关系：2022-8-1235例例 8 8 长为长为 l 的弦在的弦在 x=0 端固定，另一端端固定，另一端 x=l 自由，自由，且在初始时刻且在初始时刻 t=0=0 时处于水平状态，初始速度为时处于水平状态，初始速度为 x(l-x)，且已知弦作微小横振动，试写出此定解问题，且已知弦作微小横振动，试写出此定解问题.解解 （1 1）确定泛定方程：）确定泛定方程：x0 x 取弦的水平位置为取弦的水平位置为轴，轴，为原点，为原点，弦

31、作自由（无外力）横振动，所以泛定方程为齐次波动方程弦作自由（无外力）横振动，所以泛定方程为齐次波动方程 20ttxxua u(2)(2)确定边界条件确定边界条件 对于弦的固定端，显然有对于弦的固定端，显然有 u(x,t)|x=0=0,ux(x,t)|x=l=0 另一端自由，意味着弦的张力为零则另一端自由，意味着弦的张力为零则 0,xx lux t2022-8-1236(3)(3)确定初始条件确定初始条件 0t 根据题意，当根据题意，当时，弦处于水平状态，即初始位移为零时，弦处于水平状态，即初始位移为零 初始速度初始速度 综上讨论，故定解问题为综上讨论，故定解问题为00,tu x t0,tu x

32、 tx xlt200000000000(,),(),()ttxxxxx ltttua uxl tu x tux ttu x tux tx xlxl2022-8-1237例9 在均匀静电场中置入半径为R0 的导体球，若导体 球接有稳恒电池，使球与地保持电势差u0。试写出电势u满足的泛定方程与定解条件。解：选z轴沿均匀外电场的方向，见图1。00Rz0E0u1u(a)00EZdl(b)（图 1）2022-8-1238 设球内外电势分别用u0、u1 表示。（1）泛定方程。因为除球面上(R=R0)有自由电荷分布外，球内外的f=0，故（2）定解条件 球面上电势连续，即有边界条件2002100,0,uRRu

33、RR00Rz0E0u1u(a)00EZdl(b)（图 1）00100()()uRuRu2022-8-1239 现在计算上式从 到 的积分。由于在静电场中，上式的积分与积分的路线无关，故可取积分路线l为直线，如图(1)所示。将 Ecos作为常数提出积分号外，并将u(R0)=u0 代入，便有边界条件000000|cos|cos|cosRRRuuE RE RuuE R 无限远处：可以把导体表面有限的电荷分布产生的电势和电势u0看成点电荷和点电势源,由于点电荷在无限远处的贡献可以忽略不计，故可把目前问题简化为点电势在空间的分布问题。对于点电势，随着离开点势源的距离l的增加,电势是减少的,由图(1)可得

34、00ddcos dluEER 2022-8-1240 7.2 1.3.5.72022-8-1241(一一)线性二阶偏微分方程线性二阶偏微分方程(1)(1)偏微分方程偏微分方程 含有未知多元函数及其偏导数的方程，如含有未知多元函数及其偏导数的方程，如22222(,)0uuuuuF x yuxyxyx y 其中其中(,)u x y 是未知多元函数，是未知多元函数，而而 ,x y 是未知变量；是未知变量；,uuxy为为u的偏导数的偏导数.有时为了书有时为了书写方便，通常记写方便，通常记22,xyxxuuuuuuxyx7.3 7.3 数学物理方程的分类数学物理方程的分类*2022-8-1242(2)(

35、2)方程的阶方程的阶 偏微分方程中未知函数偏导数的最高阶数称偏微分方程中未知函数偏导数的最高阶数称为方程的阶为方程的阶(3)(3)方程的次数方程的次数 偏微分方程中最高阶偏导数的幂次数称为偏微分方程中最高阶偏导数的幂次数称为偏微分方程的次数偏微分方程的次数(4)(4)线性方程线性方程 一个偏微分方程对未知函数和未知函数的所一个偏微分方程对未知函数和未知函数的所有（组合）偏导数的有（组合）偏导数的幂次数幂次数都是一次的，就称为线性方程，都是一次的，就称为线性方程，高于一次以上的方程称为非线性方程高于一次以上的方程称为非线性方程(5)(5)准线性方程准线性方程 一个偏微分方程，如果仅对方程中所有最

36、一个偏微分方程，如果仅对方程中所有最高阶偏导数是线性的，则称方程为准线性方程高阶偏导数是线性的，则称方程为准线性方程(6)(6)自由项自由项 在偏微分方程中，不含有未知函数及其偏导数的在偏微分方程中，不含有未知函数及其偏导数的项称为自由项项称为自由项2022-8-1243（7 7）方程的通解）方程的通解如：二阶线性非齐次偏微分方程如：二阶线性非齐次偏微分方程 2xyuyx的的通解通解为为221(,)()()2u x yxyx yF xG y其中其中(),()F x G y是两个独立的任意函数是两个独立的任意函数若函数若函数 (),()F x G y的具体形式给定，则得到方程的特解。的具体形式给

37、定，则得到方程的特解。（8 8）方程的特解）方程的特解方程的解含有方程的解含有任意元素任意元素（即任意常数或任意函数）（即任意常数或任意函数）2022-8-12441112221220 xxxyyyxya ua ua ubub ucufa11,a12,a22,b1,b2,c,f 只是只是 x,y 的函数。的函数。f 0 方程为方程为齐次齐次的的;否则否则 ,为为非齐次非齐次的的.叠加原理叠加原理 定解问题的解可以看作几个部分的线性叠加，只要这些定解问题的解可以看作几个部分的线性叠加，只要这些部分各自所满足的泛定方程和定解条件的相应的线性叠部分各自所满足的泛定方程和定解条件的相应的线性叠加正好是

38、原来的泛定方程和定解条件。加正好是原来的泛定方程和定解条件。泛定方程、定解条件都是线性泛定方程、定解条件都是线性(二二)二阶线性偏微分方程的分类二阶线性偏微分方程的分类2022-8-1245作变换：作变换：为使变换非奇异，其雅克比行列式满足为使变换非奇异，其雅克比行列式满足,yyxx变换变换运算运算有有,u x y,u,x yx y即即 22222222xxyyyyxxxxxxxxxxxxxxxxxxxxxyxyxyyxxyxyxyyyyyyyyyyyuuuuuuuuuuuuuuuuuuuuuuuuuuuuuu()2022-8-1246采用新变量后的方程采用新变量后的方程1112221220A

39、 uA uA uB uB uCuF22111112221211122222221112221111222122111222122222xxyyxxxyyxyyxxyyxxxyyyxyxxxyyyxyAaaaAaaaAaaaBaaabbBaaabbCcFf()其中其中2022-8-1247注意注意A11和和A22形式相同，形式相同，和和用用z表示，如果表示，如果则有则有 （或（或 ）注意到注意到二阶线性偏微二阶线性偏微分方程的分方程的特征特征方程方程特征方程的根为：特征方程的根为：通过求解此微分方程可以得到变换函数（通过求解此微分方程可以得到变换函数（特征线特征线），从而线性偏微分方程得以简化。

40、，从而线性偏微分方程得以简化。2022-8-1248定义：定义：1.1.当判别式当判别式 以求得以求得两个实函数解两个实函数解 时，从特征方程可时，从特征方程可12(,)(,)x yCx yC及也就是说，偏微分方程也就是说，偏微分方程(1)(1)有有两条实的特征线两条实的特征线于是取于是取2(,)0uuuu 方程可化为：方程可化为：作为新的自变量，此时有：作为新的自变量，此时有：2022-8-1249或者进一步作变换或者进一步作变换,=于是有于是有,所以所以22222 uuu 又可以进一步将方程化为又可以进一步将方程化为22122(,)0uuuuu 形式形式我们前面建立的波动方程就属于此类型我

41、们前面建立的波动方程就属于此类型 这种类型的方程称为这种类型的方程称为双曲型方程双曲型方程，是双曲型方程的是双曲型方程的标准标准2022-8-12502 2当判别式当判别式 时：这时方程重根时：这时方程重根1211dyadxa0(,)x yC特征线为特征线为一条实特征线一条实特征线作变换作变换 (,)x y任意选取另一个变换，任意选取另一个变换，(,)x y只要它和只要它和 (,)x y彼此独立，即雅可比式彼此独立，即雅可比式(,)0(,)xy2022-8-1251方程可化为：方程可化为：22(,)0uuuu 此类方程称为此类方程称为抛物型方程抛物型方程热传导（扩散）方程就属于热传导（扩散）方

42、程就属于这种类型这种类型2022-8-12523.3.当判别式当判别式 面的讨论，只不过得到的面的讨论，只不过得到的 时：可以重复上时：可以重复上(,)x y和和 (,)x y是一是一对对共轭的共轭的复函数，或者说，复函数，或者说，两条特征线是一对两条特征线是一对共轭复函数共轭复函数族族:(,),(,)x yx y是一对共轭的复变量进一步引进两个新的实变量是一对共轭的复变量进一步引进两个新的实变量,=i()于是于是i,i2022-8-1253所以所以 22222 uuu 22222(,)0uuuuu 方程进一步化为方程进一步化为这种类型的方程称为这种类型的方程称为椭圆型方程椭圆型方程拉普拉斯拉

43、普拉斯(Laplace)(Laplace)方方程、泊松程、泊松(Poisson)(Poisson)方程和方程和Helmholtz Helmholtz 方程都属于这种类方程都属于这种类型型2022-8-1254小结：小结：1112221220 xxxyyyxya ua ua ubub ucuf =a212 a11a22判别式判别式 0 双曲型双曲型 =0 抛物线型抛物线型 x/a时，上式后两项无意义，必须将时，上式后两项无意义，必须将 u(x,t)延拓到这个范围延拓到这个范围(0,)0ut ，作奇延拓作奇延拓：()()xx ()()xx ()(0)()()(0)xxxxx ()(0)()()(0

44、)xxxxx 半无限长弦的自由振动半无限长弦的自由振动2022-8-126411(,)()()()22x atx atu x txatxatda x atx atx atatxxatxatdtx aau x txatatxdtx aa 11()()()(/)22(,)11()()()(/)22 2022-8-1265半波损失只有初始位移，没有初始速度只有初始位移，没有初始速度开始反射开始反射2022-8-1266一个端点自由一个端点自由22222()(,)0au x ttx (0)x 设初始条件为设初始条件为0()tux 和和0()ttux 边界条件边界条件00 xxu 应该是偶延拓应该是偶延

45、拓()(0)()()(0)xxxxx ()(0)()()(0)xxxxx 11(,)()()()22x atx atu x txatxatda 001122111222()()()(/)(,)()()()()(/)x atx atx atatxxatxatdtx aau x txatatxddtx aaa 2022-8-1267无半波损失无半波损失只有初始位移，没有初始速度只有初始位移，没有初始速度开始反射开始反射2022-8-1268（三）跃变点的反射（三）跃变点的反射 无限长杆，无限长杆，x0 两部分的杨氏模量和密度分别两部分的杨氏模量和密度分别为为 。x=0 是跃变点是跃变点。,IIII

46、IIYY 设有行波设有行波 从区域从区域 I 向向 x=0 点运动。到点运动。到x=0 产产生反射和透射生反射和透射。(,)(/)Iu x tf tx a取此波在取此波在 t=0 时刻抵达时刻抵达 x=0 .200,(0)()(0)IIttxxItIua uxxuf txa2000,(0)0,0(0)IIIIttxxIIIItttua uxuux衔接条件衔接条件0000,IIIxxIIIIIIxxxxuuY uY ux0 x,I II IY,IIY(/)If txa()IIxh ta()Ixg ta2022-8-1269区域区域 II II 中，只有透射波中，只有透射波(,)()(0)IIII

47、xux th txa0t()0,()0IIIIxxhhaa()0,()0(0)hh衔接条件衔接条件()()(),(0)11()()()IIIIIIf tg th ttYftgtYh taa()()(),(0)()()()IIIIIIf tg th tta Yf tg ta Yh t2()()(0)()()IIIIIIIIIIIIIIIIIIIIIa Yhfa Ya Ya Ya Ygfa Ya Y 区域区域 I I 中的行波：中的行波：(,)()(),(0)IIIxxux tf tg txaa0t()0Ixga()0,(0)g 2022-8-12700()()()()IIIIIIIIIIIIII

48、IIxtxag taa Ya Yxxf tta Ya Yaa0()()2()()IIIIIIIIIIIIIIIIIxtxah taa Yxxf tta Ya Yaa2221IIIIIIIIIIIIIIIIIIIIIa Ya Ya Ya Ya Ya Ya Y又又2IIIIIIIIIIIIa Ya Ya Ya Y22IIIIIIIIIa Ya Ya Y反射系数反射系数透射系数透射系数2022-8-1271 从从达朗贝尔公式达朗贝尔公式可以看出，波动方程的解，是初始可以看出，波动方程的解，是初始条件的演化。方程本身并不可能产生出超出初始条件的，条件的演化。方程本身并不可能产生出超出初始条件的，额外的

49、形式来。额外的形式来。而这种演化又受到边界条件的限制。而这种演化又受到边界条件的限制。这就说明了初始条件和边界条件在确定波动方程的这就说明了初始条件和边界条件在确定波动方程的解时的重要性。解时的重要性。11(,)()()()22x atx atu x txatxatda 达朗贝尔解表示沿达朗贝尔解表示沿x 轴正、反向传播的两列波的叠加轴正、反向传播的两列波的叠加2022-8-1272例例9 922222()(,)0au x ttx(0)x 设初始条件为设初始条件为0()tux 和0()ttux 边界条件边界条件0sinxxAutYS 定解问题定解问题20000,(0)(),()sinttxxt

50、ttxxua uxuxuxAutYS 20000,(0)(),()0ttxxtttxxua uxuxuxu 端点自由时的解端点自由时的解20000,(0)0,0sinttxxtttxxua uxuuAutYS 分解分解分析：分析：2022-8-1273解：通解为解：通解为12(,)()()u x tf xatfxat定解条件定解条件121200()()()()f xfxafxafx12()()sinAfatfattYS 求：求：12 0 (),(,)(),(,)ff 2022-8-1274求求2()?fxat121200()()()()f xfxafxafx21 02()(cos),aAcfY