32三维波动方程初值问题课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《32三维波动方程初值问题课件.ppt》由用户(ziliao2023)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 32 三维 波动 方程 初值问题 课件
- 资源描述:
-
1、3.2 三维波动方程初值问题三维波动方程初值问题三维齐次波动方程的球对称解三维齐次波动方程的球对称解三维齐次波动方程的泊松公式和三维齐次波动方程的泊松公式和球平均法球平均法泊松公式的物理意义泊松公式的物理意义三维非齐次波动方程的初值问题三维非齐次波动方程的初值问题和推迟势和推迟势2.2.三维波动方程初值问题三维波动方程初值问题基本思路:将三维问题转化为一维问题三维波动方程可描述声波、电磁波和光波等在空间中的传播,称为球面波。考虑初值问题其中 满足一定的光滑性条件。23300(),(,),0(2.1)(,),(,),(,)ttxxyyzztttua uuux y zR tux y z ux y
2、zx y zR2.1 2.1 三维齐次波动方程的球对称解三维齐次波动方程的球对称解,引入球坐标系 即sincos,sinsin,cosxryrzr(,),r 0,0,02,r 则方程(2.1)可化为2222222111 sin (2.2)sinsinttuuarrrruurr所谓球对称解球对称解,是指在球面上各点的值都相等的解(设球心为原点),即 与 和 无关。(,)(,)u x y z tu r t故当 u 是球对称函数时,方程(2.2)可化为22,0,0(2.3)ttrruuaurtrr或者等价地写成22()(2)(),ttttrrrrrruruaruuaru令 ru=v,则有 其通解可表
3、示为 2,ttrrva v()(),0,0,vF ratG ratrt其中F(r+at)是沿 r 负方向传播,为收敛波,G(r-at)是沿 r 正方向传播的行波,为发散波。从而,()()(,),0,0,F ratG ratu r trtr其中 F,G 是任意两个二阶连续可微函数。若考虑初始条件(,0)(),(,0)(),0,(2.4)tu rr u rrr则类似于半界弦的振动情况,可得初值问题(2.3)-(2.4)的解1()()()()21 (),;2(,)1()()()()21 (),00 .2at rr atr atr atatraratratratratrdaru r tratratrr
4、atarrtratd 2.2 2.2 三维齐次波动方程的泊松公式和球平均法三维齐次波动方程的泊松公式和球平均法(1)主要结果一维齐次波动方程的达朗贝尔解11(,)()()+()22x atx atu x txatxatda 可改写成1()21()2(,)+x atx atxxtataddatx ttttaut 其中 为初始位移 在 上的算术平均值,1()2x atx atdat,xat xat 为初始速度 在 上的算术均值,xat xat1()2x atx atdat 受此启发,在以M(x,y,z)为中心,以at为半径的球面上作初始函数 和 的平均值,分别为2 22 211(,),(,).44
5、MMatatSSdSdSa ta t 则问题(2.1)的解应该是(待证)2 22 211(,)(,)+(,)44MMatatSSu x y z ttdStdSta ta t 221(,)1(,)+,(2.5)44MMatatSSdSdSattat 其中 为M(x,y,z)为中心,以at为半径的球面。MatS为简化计算,将公式(2.5)在球坐标下化为累次积分,球面的方程为MatS2222()()()().xyzat设 为球面上的点,则(,)P sincos,sinsin,cos,xatyatzat三维齐次波动方程初值问题的三维齐次波动方程初值问题的Poisson公式公式2 2sindSa td
6、d 于是22 22 20022 22 2001(,)(,)sin41 +(,)sin4u x y z tta td dta tta td da t 200200(sin cos,sin sin,4cos)sin(sin cos,4 sin sin,cos)sin.(2.6)txatyatttzatd dxatyatzatd d (2)Poisson公式(5)的推导推导思路推导思路球平均法球平均法一般情况下,ru 未必满足一维波动方程。设法找一个与u有关的球对称函数 通过 把 u 求出来。,uu考虑 u 在球面 上的平均值,即MrS1211(,)(,),(2.7)44MMrSSu r tudSu
7、t dr 其中sincos,sinsin,0,0,02,cos,xryrrzr是球面 上的点的坐标,是单位球面上的面积元,且有 ,则MrS22sindSrd dr d 0(,)(,)lim(,)(0,)ru M tu x y z tu r tut球平均法下面证明 满足一维波动方程ru2(,)(,)(2.8)ttrrru r ta ru r t设 表示中心在 M 的半径为r的球域。对方程(2.1)的两边在 上积分,并利用高斯公式及(2.7),有MrBMrB2()()()MMrrttxxyyzzBBu dxdydzauuudxdydzd22(,)MMrrxyzSSuau u undSadSn122
8、224.MSuuar da rrr 另一方面,由于220MMrrttBSu dxdydzudSdt 1222222004,MrrSud dudtt 故有222220.ruuda rtr此式两端关于 r 求导,有22222().ur uartrr于是 满足22222,uaurtrrru即22222()().ruruatr所以(,)()()(2.9)ru r tF ratG rat即 满足一维波动方程。ru对(2.9)两边分别关于 r 和 t 求导,有()()(,)()(),ruuru r tF ratG ratrr1()()(),ruF ratG ratat将此二式相加,得()1()()1(,)
9、2(),ruruuuru r tF ratratrat令 有0,r(,)(0,)2().u x y z tutF at另一方面,在上式中取 t=0,有0()1()2()truruF rrat22011144MMrrSStrudSrudSrra tr01144MMrrtSStuudSdSrrar11.44MMrrSSdSdSrrar从而,用 at 取代 r,Poisson公式得证。定理定理1 1.若 则Poisson公式(2.5)表达的 u(x,y,z,t)在 内二阶连续可微,且为三维齐次波动方程初值问题的古典解。33(,),(,),x y zCx y zC3(0,)R 例例1 1.求解初值问题
10、233(),(,),0(,0),(,0)0,(,)ttxxyyzztua uuux y zR tu x y zxyz u x y zx y zR解解.由Poisson公式(2.6)得2001(,)4(sincossinsincos)sinu x y z ttxyztatd d 2001()sin4t xyzddt 22200(sincos)sinatdd 2200sincosatdd.xyz例例2 2.求解初值问题233(),(,),0(,0),(,0),(,)ttxxyyzztua uuux y zR tu x y zyz u x y zxzx y zR解解.法一法一.此处 由Poisson
11、公式(2.6)得200(,)1sin(sinsin)(cos)4u x y z ttyatzatd dt ,yzxz2001sin(sincos)(cos)4xatzatd d 20022 2002 21sin(sinsincos4sincos sin)sin(4cossincossincoscos).tyzzatyattta td dxzxatzata td d 由三角函数的周期性和正交性,有2200sincos0,dd 00cossincos0.dd 因此00(,)12sin2sin44u x y z tttyzdxzdt .yztxz法二法二.由于定解问题是线性的,故可由叠加原理,令12
12、3,uuuu其中 分别满足如下定解问题123,u u u1211010,00,tttxxttuaxR tuxuuxRz2220022,0,0,yytttttuuyzuayR tuyR3330032,0,0,0zztttttuuuazR tuzR由达朗贝尔公式可分别求得以上三个定解问题的解,为11,2x atx atuz dxzta 21()(),2uz yatz yatyz30,u 因此.uxztyz2.3 2.3 泊松公式的物理意义泊松公式的物理意义由泊松公式可见,定解问题(2.1)的解在M(x,y,z)点 t 时刻的值,由以 M 为中心,at 为半径的球面 上的初始值而确定。MatS这是由
13、于初值的影响是以速度 a 在时间 t 内从球面 上传播到 M 点的缘故。MatS设初始扰动限于空间某区域 内,(即在 外 ),记 d 和 D 分别为 M 点到区域 的最近和最远距离,则0,0(1)当 时,即 时,与 不相交,上的初始函数 为0,故u(M,t)=0。这说明扰动的前锋尚未达到 M 点。atddtaMatS,MatS(2)当 即 时,上的初始函数 不为0,故u一般不为0。这表明扰动正在经过M点。,datDdDtaa MatS,(3)当 时,即 时,与 也不相交,因而同样 u(M,t)=0,这说明扰动的阵尾已传过M点,M又恢复到静止状态。atDDtaMatS三维空间的初始局部扰动,在不
14、同的时间内对空间每一点发生影响,且波的传播有清晰的前锋和阵尾,这种现象物理上称为惠更斯(Huygens)原理或无后效现象。现实生活中声音的传播就是一例:从某处发出声音,经过一段时间后,才能听到,再经过一段时间之后恢复到静止状态。例例3 3.高空大气中有一半径为1的球形薄膜,薄膜内的压强超过大气的数值为 假定该薄膜突然消失,将会在大气中激起三维波,试求球外任意位置的附加压强P。0,P解解.设薄膜球心到球外任意一点的距离为d,则其定解问题为2303(),(,),0,1,(,0)(,0)0,(,)0,1,ttxxyyzztPaPPPx y zR tPdP x y zP x y zx y zRd当 时
15、,由泊松公式(2.5)有11datd 21(,)(,)4MatSP x y z tdSatt 22 202001sin4Pda tdatt 20212(1 cos)4Pa tat22 2202112142da tPa tatdat2021()14Padatatd 0()2Pdatd当atd+1 时,(,)0.P x y z t 2.4 2.4 三维非齐次波动方程的初值问题和推迟势三维非齐次波动方程的初值问题和推迟势考虑非齐次波动方程的初值问题23300()(,),(,),0(2.10)(,),(,),(,)ttxxyyzztttua uuuf x y z tx y zR tux y z ux
16、y zx y zR该问题可分解为下面两个问题23300(),(,),0(,(,),(,)2.11,),(,ttxxyyzztttva vvvx y zR tvx y z vx y zx y zR23300()(,),(,),(2.1200,(,)0,)ttxxyyzztttwawwwf x y z tx y zR twwx y zR问题(2.11)(2.12)的解分别为v(x,y,z,t),w(x,y,z,t),由叠加原理,问题(2.10)的解 u=v+w.对于问题(2.12),齐次化原理同样成立,即()201(,)(,)4Ma ttSfw x y z tdSdat 作代换 上式为 20(,)
展开阅读全文