行波法与积分变换法课件

第三章第三章 行波法与积分变换法行波法与积分变换法行波法(求解无界区域内波动方程定解问题)积分变换法(无界或有界区域)3.1 一维波动方程的达朗贝尔公式3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式考虑代换利用复合函数求导法则得3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式同理有:代入方程，得到 在上式中对 积分,得(是 的任意可微函数)3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式再将此式对 积分,其中 都是任意二次连续可微函数.利用初始条件，确定两个函数的具体形式。由第二式得.其中3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式由由,解得代入通解表达式，得达朗贝尔达朗贝尔(DAlembert)(DAlembert)公式公式.3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式图 3-1 u2xt=0u2xu2xt=1/2u2xt=1t=2考虑 的物理意义随着时间t 的推移u2的图形以速度a 向x轴正向移动.3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式物理意义物理意义:随着时间 t 的推移,的图形以速度 a 向 x 轴正方向移动,也就是说,它表示一个以速度a 向x 轴正方向行进的波,称为右行波右行波右行波右行波.同样道理,以速度a 向x 轴负方向传播的行波,称为左行波左行波左行波左行波.3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式在 平面上斜率为 的两族直线 ,对一维波动方程的研究起到重要作用,称这两族直线为一维波动方程的特征线特征线特征线特征线,变换称为特征变换特征变换特征变换特征变换,行波法也叫特征线法特征线法特征线法特征线法.3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式的积分曲线,这个常微分方程称为它的特征方程特征方程特征方程特征方程.一维波动方程的两族特征线恰好是常微分方程3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式一般的二阶线性偏微分方程它的特征方程为(*)这个常微分方程的积分曲线称为偏微分方程(*)的特征曲线特征曲线特征曲线特征曲线.记称其为二阶线性偏微分方程的判别式判别式双曲型方程椭圆型方程抛物型方程3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式可以证明，当 时，有两条相异的实特征线因此特征线法对双曲型方程都是有效的，沿着特征线做自变量替换 总可以把双曲型方程化为 从而得到方程的通解3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式例 求下面问题的解:(3.1)解:特征方程 两族积分曲线为 做特征变换 3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式代入方程化简得:3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式它的通解为其中 ,是两个二次连续可微函数.于是原方程的通解为代入初始条件 ,，得 第二式的两端得关于 积分得解得所求问题的解为 3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式解 特征方程为特征曲线为 例 求方程的一般解.3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式所以，做变换则原方程可以变为 其中 ,是任意的二次连续可微函数.于是，方程的通解为3.1 3.1 一维波动方程的达朗贝尔公式一维波动方程的达朗贝尔公式3.2 三维波动方程的泊松公式研究波在空间传播问题.三维波动方程的初值问题 3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式一、球对称情形球坐标系 3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式 若 仅是 r 的函数,则是r 和 t 的函数,此时称定解问题是球对称球对称的。球对称波动方程进一步有3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式对对球对称球对称问题球对称情形下，三维波动方程边值问题可化为 3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式这个问题我熟悉！由达朗贝尔公式3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式二二.一般情况一般情况令 表示 在球面 上的平均值。其中M=M(x,y,z),是球面 上的点,3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式二二.一般情况一般情况令3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式表示以 M 为中心的单位球面，表示 上的面积元素，表示单位球面上的面积元素，即而3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式以下推导 所满足方程及初始条件。3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式进一步有：两边关于 r 求导，得 得由3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式即可得：由3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式由初值条件和 的表达式，有：其中 分别是函数 在 上的球平均值。满足如下定解问题：3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式方程的通解为利用初始条件有其中是两个二次连续可微的任意函数3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式所以解方程组得3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式将 延拓到r0的范围内。并且同理 也是偶函数利用3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式所以 3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式由于 ，只考虑 的情形利用洛必达法则3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式即简记成3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式三维波动方程的泊松公式三、泊松公式的物理意义 从泊松公式出发，解释波在三维空间的传播现象.设 且，1.在任一固定点 的振动情况 设 ，由 沿以 M 为中心，at 为半径的球面的曲面 积分所决定。3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式M 点处于静止状态，说明 T 的振动尚未达到 M 点。当 时，为空集，所以 当 时，不为空集，所以M点处于振动状态,表明 T 的振动已传到 M 点。当 时，为空集，说明振动已 传过 M 点，M 点仍回复到静止状态。3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式2.在某固定时刻 ，初始时刻的振动所传播的范围 设 ，T 是半径为 R 的球体。由Poisson公式，只有与 M 相距为 的点上的初始扰动能够影响 的值，故 P 点的初始扰动，在时刻 只影响到以 P 为球心，以 为半径的球面 当 P 在 T 内移动时，球面族的包络面所围成的区域即为 T 内各点的振动在 时刻所传播的区域，称为 T 在时刻 的影响区域影响区域。3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式 总之，三维空间中有限区域 T上的初始振动，有着清晰的前阵面和后阵面，对空间的任一点，振动传过后，仍回复到平衡状态，这种只在有限时间内引起振动的现象称为 Huygens Huygens 原原 理理。在 足够大时，包络面以T 的心o(T)为心，分别以 和 为半径的球面所夹部分。故 时刻的影响区域为 的球壳，球面 是振动到来的前峰，称为波的 前前 阵阵 面面,球面 是振动传过后的后沿，称为波的后阵后阵面面。3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式R3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式解 例.设已知三维波动问题中的初位移，初速度分别为：,求解相应的Cauchy问题。3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式三.降维法及二维波动方程考虑二维波动方程的初值问题 设 解 为 ,令 ，则3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式由泊松公式 球面 在平面 上投影 为 设其上面积微元为 ，则由投影关系有：其中 v 表示 dS 的单位法向量与 之夹角，3.2 3.2 三维波动方程的泊松公式三维波动方程的泊松公式又上、下两球面的投影有对称关系，故柱面波3.3 积分变换法 3.3 3.3 积分变换法积分变换法 常见的两种积分变换常见的两种积分变换-傅立叶变换傅立叶变换-拉普拉斯变换拉普拉斯变换.如果 满足上面的条件，我们可以定义傅立叶逆变换为:如果函数 在 上绝对可积，它的傅立叶变换定义如下有时把 记为 。一.傅立叶变换反演公式反演公式3.3 3.3 积分变换法积分变换法傅立叶变换的性质:1）线性性质线性性质 设 f,g 是绝对可积函数，是任 意复常数，则 2）微分性质微分性质 设 f,绝对可积函数，则 3）乘多项式乘多项式 设 f,x f 绝对可积，则 3.3 3.3 积分变换法积分变换法4）伸缩性质伸缩性质 设 f(x)绝对可积，则 6）卷积性质卷积性质 设f,g 是绝对可积函数，令 则5）平移性质平移性质 设 f(x)绝对可积，则 3.3 3.3 积分变换法积分变换法例例 用积分变换法解方程：解：作关于 x 的傅立叶变换,方程可变为 设3.3 3.3 积分变换法积分变换法可解得 由于即则3.3 3.3 积分变换法积分变换法从而方程的解 3.3 3.3 积分变换法积分变换法 例例 用积分变换法解方程解：作关于 的傅立叶变换。设方程变为 3.3 3.3 积分变换法积分变换法用常数变易法可解得 而则3.3 3.3 积分变换法积分变换法利用反演公式有3.3 3.3 积分变换法积分变换法傅立叶逆变换是一种把分析运算化为代数运算的傅立叶逆变换是一种把分析运算化为代数运算的有效方法有效方法,但但1.傅立叶变换要求原象函数在傅立叶变换要求原象函数在R上绝对可积上绝对可积,大部大部分函数不能作傅立叶变换分函数不能作傅立叶变换2.傅立叶变换要求函数在整个数轴上有定义傅立叶变换要求函数在整个数轴上有定义,研究研究混合问题时失效混合问题时失效.3.3 3.3 积分变换法积分变换法二二.拉普拉斯变换拉普拉斯变换定义定义:f(t)定义在 上，若其满足下列条件1.f(t)分段光滑；2.存在常数 M 和 使得 3.则称f(t)为初始函数,称为f(t)的增长指数.反例反例3.3 3.3 积分变换法积分变换法定理定理:设f(t)是一以 为增长指数的初始函数,则经变换得到的函数F(p)是 上的解析函数.上述变换称为拉普拉斯变换3.3 3.3 积分变换法积分变换法例例 3.3 3.3 积分变换法积分变换法反演公式反演公式：在 f(t)的每一个连续点均有 其中，3.3 3.3 积分变换法积分变换法基本性质基本性质:1）线性性质 设 f,g 的拉普拉斯变换分别 为L(f),L(g ),是任意复常数，则 2）微分性质 假设 ,则3.3 3.3 积分变换法积分变换法6）卷积性质定义4）延迟性质5）伸缩性质则3）积分性质3.3 3.3 积分变换法积分变换法例 设 求解常微分方程的初值问题 解 对 进行拉普拉斯变换,设 ,则3.3 3.3 积分变换法积分变换法于是原方程变为由上式得:对 进行拉普拉斯逆变换,得3.3 3.3 积分变换法积分变换法解 问题归结为求解下列定解问题:例 一条半无限长的杆，端点温度变化已知，杆的初始温度为0，求杆上温度分布规律。3.3 3.3 积分变换法积分变换法对 t 进行拉普拉斯变换怎么变换？为什么？知道 的值了方程通解为 表示温度，当 时，一定有界，所以 亦有界，从而 .对 t 进行拉普拉斯变换，设于是原问题变为3.3 3.3 积分变换法积分变换法由边值条件可知 ，即 对 p 进行拉普拉斯逆变换，有3.3 3.3 积分变换法积分变换法于是查表得而易证3.3 3.3 积分变换法积分变换法所以于是3.3 3.3 积分变换法积分变换法 例 设 求解下面定解问题解 对 进行拉普拉斯变换,则原方程 变为即3.3 3.3 积分变换法积分变换法由条件 得 解得对 取拉普拉斯逆变换，得 3.3 3.3 积分变换法积分变换法数学物理方程数学物理方程+定解条件定解条件解解常微分方程常微分方程+定解条件定解条件解解积分变换积分变换逆变换逆变换3.3 3.3 积分变换法积分变换法如何使用积分变换法求解定解问题：如何使用积分变换法求解定解问题：1)1)选选取取恰恰当当的的积积分分变变换换，对对某某个个（某某些些）自自变变量量作作积积分分变变换换，得得到到象象函函数数的的含含参参变变量量的的常常微微分分方方程；程；2 2）对对部部分分定定解解条条件件取取相相应应的的积积分分变变换换,导导出出象象函函数方程的定解条件；数方程的定解条件；3 3）解关于象函数的定解问题）解关于象函数的定解问题,求出象函数；求出象函数；4 4）将象函数取积分逆变换，即得原定解问题的解）将象函数取积分逆变换，即得原定解问题的解.3.3 3.3 积分变换法积分变换法1）傅立叶变换的取值范围是 ，拉普拉斯变换的取值范围是 。2）注意定解条件的形式。假如对 进行拉普拉斯变换，而原方程为 阶方程，则定解条件中应出现 需要注意需要注意3.3 3.3 积分变换法积分变换法