第七章 一维波动方程的解题方法及习题答案
- 格式:doc
- 大小:1008.50 KB
- 文档页数:23
一维波动方程是描述一维波动现象的偏微分方程,其一般形式为:
∂²u/∂t²= c²∂²u/∂x²
其中,u(x,t)是波动的位移或振幅,c是波速。
解一维波动方程的一般步骤是将其转化为一个简单的常微分方程或特殊的偏微分方程,然后通过求解这个方程得到波动的解析表达式。
这里介绍两种求解方法:
分离变量法
假设u(x,t) 可以表示为两个函数的乘积形式:u(x,t) = X(x)T(t),代入原方程得到:
X''(x)/X(x) = (1/c²) T''(t)/T(t)
由于左边的方程只涉及x,右边的方程只涉及t,因此两边必须都等于一个常数k²,即:X''(x)/X(x) = k²
T''(t)/T(t) = k²c²
分别解上面两个方程,得到:
X(x) = A sin(kx) + B cos(kx)
T(t) = C sin(ckt) + D cos(ckt)
其中,A、B、C、D 是待定常数,可以根据边界条件和初值条件确定。
将上述两个函数代回原方程,得到波动的解析表达式:
u(x,t) = Σ[An sin(nπx/L) + Bn cos(nπx/L)] sin(ncπt/L)
其中,An、Bn 是待定常数,L 是波动区间的长度,n 为正整数。
第七章一维波动方程的傅里叶解小结及习题答案第二篇数学物理方程——物理问题中的二阶线性偏微分方程及其解法Abstracts:1、根据物理问题导出数理方程—偏微分方程;2、给定数理方程的附加条件:初始条件、边界条件、物理条件(自然条件,连接条件),从而与数理方程一起构成定解问题;3、方程齐次化;4、数理方程的线性导致解的叠加。
一、数理方程的来源和分类(状态描述、变化规律)1、来源I.质点力学:牛顿第二定律Fmr连续体力学弦2u(r,t)弹性体力学杆振动:22波动方程);au(r,t)0(2t(弹性定律)膜流体力学:质量守恒律:(v)0;t热力学物态方程:v1(v)vpf0(Eulereq.).tII.麦克斯韦方程DddD;EdlBdsEB;Bd0B0;Hdl(jD)dsHjD.Eu,BA,u,A满足波动方程。
Lorenz力公式力学方程;Maxwelleqs.+电导定律电报方程。
III.热力学统计物理热传导方程:扩散方程:Ttt2kT2D0;0.特别:稳态(0t):20(Laplaceequation).IV.量子力学的薛定谔方程:2u2.iuVut2m2.分类物理过程方程数学分类振动与波波动方程2u 12u22at双曲线输运方程能量:热传导质量:扩散ut20ku抛物线1稳态方程Laplaceequation 2u0椭圆型二、数理方程的导出推导泛定方程的原则性步骤:(1)定变量:找出表征物理过程的物理量作为未知数(特征量),并确定影响未知函数的自变量。
(2)立假设:抓主要因素,舍弃次要因素,将问题“理想化”---“无理取闹”(物理趣乐)。
(3)取局部:从对象中找出微小的局部(微元),相对于此局部一切高阶无穷小均可忽略---线性化。
(4)找作用:根据已知物理规律或定律,找出局部和邻近部分的作用关系。
(5)列方程:根据物理规律在局部上的表现,联系局部作用列出微分方程。
Chapter7一维波动方程的傅里叶解第一节一维波动方程-弦振动方程的建立1.弦横振动方程的建立(一根张紧的柔软弦的微小振动问题)(1)定变量:取弦的平衡位置为x轴。
高等数学 第四册(第三版) 数学物理方法 答案(完整版)第七章 一维波动方程的傅氏解1. 今有一弦,其两端被钉子钉紧,作自由,它的初位移为: 2.(01)()(2)(12)hx x x h x x ϕ≤<⎧=⎨-≤≤⎩,初速度为0,试求其付氏解,其中h 为已知常数。
解:所求问题是一维波动方程的混合问题:2(12,0)(0,)(,)0(0)(01)(,0)(2)(12)(,0)0tt xx t u a u x t u t u l t t hx x u x h x x u x ⎧=<<>⎪==≥⎪⎪≤≤⎧⎨=⎨⎪-≤≤⎩⎪⎪=⎩,根据前面分离变量解法得其傅氏解为:1(,)(cossin )sin n n n n at n at n xu x t C D l l l πππ∞==+∑。
其中,122201228()sin [sin (2)sin ]222l n n n n hC d h d h d l l n πξπξπξϕξξξξξξπ==+-=⎰⎰⎰,0n D =,于是所求傅氏解为:2218(,)cos sin n h n at n xu x t n l l πππ∞==∑2.将前题之初始条件改为:(1)(10)()(1)(01)h x x x h x x ϕ+-≤≤⎧=⎨-≤≤⎩,试求其傅氏解。
解:所求问题为一维波动方程的混合问题:211((1)sin (1)sin n n l l l h d h d πξπξξξξξ--=++-⎰⎰n c 012222211(sinsinsin )n n n h d d d πξπξπξξξξξ--=++⎰⎰⎰2282sin h n n ππ=22821(,)sin cossinh n n at n x lln n u x t ππππ∞=∴=∑。
3今有一弦,其两端0x =和x l =为钉所固定,作自由摇动,它的初位移为0。
初速度为[](2()0(2,c x x x βϕβ≤≤⎧=⎨∉⎩,其中c 为常数,0,l αβ<<<试求其傅氏解。
波动方程习题答案波动方程是物理学中一种重要的方程形式,描述了波动现象的传播和演变规律。
在学习波动方程时,经常会遇到一些习题,需要我们进行解答。
本文将针对波动方程的一些习题进行解答,帮助读者更好地理解和应用波动方程。
1. 一维弦上的波动问题考虑一根固定在两端的弦,当弦被扰动后,波动会沿着弦传播。
对于一维弦上的波动问题,可以通过波动方程进行描述。
波动方程的一般形式为:∂²u/∂t² = v²∂²u/∂x²其中,u表示弦上的位移,t表示时间,x表示弦上的位置,v表示波速。
对于给定的初值条件和边界条件,可以通过求解波动方程得到弦上的位移分布。
具体的求解方法包括分离变量法、傅里叶级数法等。
2. 反射和折射问题当波动遇到界面时,会发生反射和折射现象。
这些现象可以通过波动方程的解来进行分析。
对于反射问题,我们可以考虑一个波从一个介质传播到另一个介质的情况。
根据波动方程和边界条件,可以求解出反射波和入射波的幅度和相位关系。
对于折射问题,我们可以考虑一个波从一个介质传播到另一个介质时发生的折射现象。
根据波动方程和边界条件,可以求解出折射波的传播方向和折射角。
3. 球面波问题当波动在空间中传播时,会形成球面波。
球面波是一种特殊的波动形式,可以通过波动方程进行描述。
对于球面波问题,我们可以考虑一个点源在空间中产生的波动。
根据波动方程和边界条件,可以求解出球面波的传播速度和幅度衰减规律。
4. 波动方程的应用波动方程在物理学和工程学中有着广泛的应用。
例如,在声学中,可以通过波动方程来描述声波的传播和衰减规律。
在光学中,可以通过波动方程来描述光波的传播和干涉现象。
此外,波动方程还可以应用于地震学、电磁学、量子力学等领域。
在地震学中,可以通过波动方程来研究地震波的传播和地壳的结构。
在电磁学中,可以通过波动方程来描述电磁波的传播和电磁场的分布。
在量子力学中,可以通过波动方程来描述粒子的波动性质和量子态的演化。
第七章 一维波动方程的傅氏解1.今有一弦,其两端被钉子钉紧,作自由振动,它的初始位移为()()⎩⎨⎧≤<-≤≤=2x 1 ,21x 0 ,x h hx x ϕ,初速度为0,试求其傅氏解,其中h 为已知常数。
解法1,(1)此问题归结为定解问题:()()()()()()⎪⎩⎪⎨⎧≤≤==≥==<<=(3))(0 0,,0,(2) )0( 0,,0,0(1))(0 2L x x x u x x u t t L u t u L x u a u t xx tt ψϕ其中()()⎩⎨⎧≤<-≤≤=2x 1 ,21x 0 ,x h hx x ϕ(2)求其定解问题的傅氏解,应用分离变量法a 变量分离 设方程的解为()()x X x T u =代入(1)得X X aT T ''"=令 λ==X X aT T ''"于是得()()⎪⎩⎪⎨⎧===+0000"L X X x X λ (()00=X ,()0=L X 由(2)得到)(5) 02"= +T a T λB 解特征值问题()()(4) 0000"⎪⎩⎪⎨⎧===+L X X x X λ(a )0<λ 时,方程(4)的通解为()x x Be Ae x X λλ---+=其中AB 为任意常数. 当0=x 时 ()B A X +==00当L x =时, ()x x Be Ae L X λλ---+==0所以 ⎪⎩⎪⎨⎧=+=+---00x x Be Ae B A λλ因为011≠=∆---x x e e λλ 所以A=B=0 ,所以当0<λ时特征问题(4)只有平凡解.(b)当0=λ时特征问题也有平凡解. (c) 当0>λ时,方程0"=+x X λ的通解为()x B x A x X λλsin cos += 代入边界条件(2)得()00100=⋅+⋅==B A x ,所以A=0()0sin 0sin 1cos =⇒=+=L B L B A L x λλλ 为使()0≠x X 应有0≠B ,所以0sin =L λ 所以 ,3,2,1 ==n n L πλ满足上面等式的λ值称为特征值,记为n λ既 ,2,1 222==n Ln n πλ 相应与的函数()x L n B x X n n πsin =称为特征函数有时记为()x Ln x X n πsin = C 解不构成特征值问题的常微分方程02"=+T a T n λ 其通解为()t La n D t L a n C t T n n n ππsin cos += 于是我们得到方程(1)满足边界条件(2)的可分离变量的一系列特解:()()()x L n t L a n D t L a n C x X t T t x u n n n n n πππsin sin cos ,⎪⎭⎫ ⎝⎛+== d.由叠加原理得形式解:()∑∞=⎪⎭⎫ ⎝⎛+=1sin sin cos ,n n n x L n t L a n D t L a n C t x u πππ e.由Fourier 级数确定n n D C ,()()∑∞===0sin0,n n x Ln C x x u πϕ ()()∑∞===1sin 0,n nt x Ln L a n D x x u ππψ 所以()ξξπξϕd Ln L C L n ⎰=0sin 2()()⎰⎰==L n d Ln a n d L n a n L L D 0sin 2sin 2ξξπξψπξξπξψπ ()2sin 82sin 22sin 22221021ππππn n h xdx n x h dx n hx C n =⎥⎦⎤⎢⎣⎡-+=⎰⎰ 因为()0=x ψ 所以0=n D 所以()∑∞==1222sin 2cos 2sin 8,n x n t a n n n h t x u ππππ 解法2 设傅氏解为 ()∑=⎪⎭⎫ ⎝⎛+=1sin sin cos ,n n n n x L n t L a n D t L a n C t x u πππ 其中()() ,2,1sin 2sin 200=⎪⎪⎩⎪⎪⎨⎧==⎰⎰n d L n a n D d L n L C L n L n ξπξξψπξπξξϕ 由此计算出()2sin 82sin 22sin 22221021ππππn n h xdx n x h dx n hx C n =⎥⎦⎤⎢⎣⎡-+=⎰⎰ 0=n D 故()∑∞==1222sin 2cos 2sin 8,n x n t a n n n h t x u ππππ 2.将前题之初始条件该为()()()⎩⎨⎧≤≤+<<-+=11101 1x x h x x h x ϕ 试求其傅氏解。
一维波动方程可用如下的方式推导:一列质量为m的小质点,相邻质点间用长度h的弹簧连接。
弹簧的劲度系数(又称“倔强系数”)为k:
其中u(x) 表示位于x的质点偏离平衡位置的距离。
施加在位于x+h处的质点m上的力为:
其中代表根据牛顿第二定律计算的质点惯性力,代表根据
胡克定律计算的弹簧作用力。
所以根据分析力学中的达朗贝尔原理,位于x+h处质点的运动方程为:
式中已注明u(x) 是时间t的显函数。
若N个质点间隔均匀地固定在长度L = N h的弹簧链上,总质量M = N m,链的总体劲度系数为K = k/N,我们可以将上面的方程写为:
取极限N, h就得到这个系统的波动方程:
在这个例子中,波速。
第二篇 数学物理方程——物理问题中的二阶线性偏微分方程及其解法Abstracts:1、根据物理问题导出数理方程—偏微分方程;2、给定数理方程的附加条件:初始条件、边界条件、物理条件(自然条件,连接条件),从而与数理方程一起构成定解问题;3、方程齐次化;4、数理方程的线性导致解的叠加。
一、数理方程的来源和分类(状态描述、变化规律)1、来源I .质点力学:牛顿第二定律F mr = 连续体力学2222()(,)(,)0(()0;v 1()0(Euler eq.).u r t a u r t t v t v v p f t ρρρ⎧⎧∂⎪⎪-∇=⎨⎪∂⎪⎪⎩⎪∂⎪+∇⋅=⎨∂⎪∂-⎪+⋅∇=+=⎪∂⎪⎪⎩弹性定律弦弹性体力学杆 振动:波动方程);膜流体力学:质量守恒律:热力学物态方程: II.麦克斯韦方程;;00;().,,,D D E l B s E B B B H l j D s H j D E u B A u A σρτρσ⎧⋅=⇒∇⋅=⋅=⋅⇒∇⨯=⎪⎪⎪⋅=⇒∇⋅=⋅=+⋅⇒∇⨯=+⎨⎪=-∇=∇⨯⎪⇒⇒⎪⎩⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰d d d d d d d 满足波动方程。
Lorenz 力公式力学方程;Maxwell eqs.+电导定律电报方程。
III. 热力学统计物理220;0.T k T t D t ρρ∂⎧-∇=⎪⎪∂⎨∂⎪-∇=⎪∂⎩热传导方程:扩 散方程:特别: 稳态(0t ρ∂=∂):20ρ∇= (Laplace equation). IV . 量子力学的薛定谔方程:22.2u i u Vu t m∂=-∇+∂稳态方程 Laplace equation 20u ∇= 椭圆型二、数理方程的导出推导泛定方程的原则性步骤:(1)定变量:找出表征物理过程的物理量作为未知数(特征量),并确定影响未知函数的自变量。
(2)立假设:抓主要因素,舍弃次要因素,将问题“理想化”---“无理取闹”(物理趣乐)。
(3)取局部:从对象中找出微小的局部(微元),相对于此局部一切高阶无穷小均可忽略---线性化。
(4)找作用:根据已知物理规律或定律,找出局部和邻近部分的作用关系。
(5)列方程:根据物理规律在局部上的表现,联系局部作用列出微分方程。
Chapter 7 一维波动方程的傅里叶解第一节 一维波动方程-弦振动方程的建立7.1.1 弦横振动方程的建立(一根张紧的柔软弦的微小振动问题)(1)定变量:取弦的平衡位置为x 轴。
表征振动的物理量为各点的横向位移),(t x u ,从而速度为t u ,加速度为tt u .(2)立假设:①弦振动是微小的,1<<α,因此,sin tan ααα≈≈,1cos ≈α,又tan u x αα∂=≈∂,1<<∂∂∴xu ;②弦是柔软的,即在它的横截面内不产生应力,则在拉紧的情况下弦上相互间的拉力即张力),(t x T 始终是沿弦的切向(等价于弦上相互间有小的弹簧相连);③所有外力都垂直于x 轴,外力线密度为),(t x F ;④设弦的线密度(细长)为),(t x ρ,重力不计。
(3)取局部:在点x 处取弦段d x ,d x 是如此之小,以至可以把它看成质点(微元)。
质量微元:x t x d ),(ρ;微弧长:x x x u u x s d d 1d d d 222≈⎪⎭⎫ ⎝⎛∂∂+=+=(即这一小段的长度在振动过程中可以认为是不变的,因此它的密度()t x ,ρ不随时间变化,另外根据Hooke 定律F k x δδ=-可知,张力),(t x T 也不随时间变化,我们把它们分别记为()x ρ和)(x T .(4)找作用:找出弦段所受的力。
外力:x t x F d ),(,垂直于x 轴方向;张力变化:()()d cos |cos |(d )()x x x T T T x x T x αα+-=+-,x 方向紧绷,()()()()()d d sin |sin |||d x x x x x x x x x x T T Tu Tu Tu x αα++-=-=,垂直于x 轴方向。
(5)列方程:根据牛顿第二定律0)()d (=-+x T x x T ,因x 方向无位移,故T x T x x T ==+)()d (.()x Tu x t x F x Tu x t x F xu x xx x x tt d d ),(d d ),(d )(+=+=ρ 即,),(t x f u Tu xx tt =-ρ,其中ρ),(),(t x F t x f =是单位质量所受外力。
如果弦是均匀的,即ρ为常数,则可写ρT a =为弦振动的传播速度,则自由振动(0f ≡): 20tt xx u a u -=(齐次方程)。
小结1:对于弦的横振动、杆的纵振动方程(一根弹性均匀细杆的微小振动问题)、薄膜的横振动方程(张紧的柔软膜的微小振动问题),在不受外力情况下,其振动的微分方程为:22tt u a u =∇(齐次方程)其中a 为振动的传播的速度。
当单位质量所受外力为f 时,其振动微分方程为:22tt u a u f =∇+(非齐次方程)7.1.2 定解问题第一节从物理问题和相应的物理定律导出了其所满足的偏微分方程,但总是选择物体内部,不含端点或边界,对一小部分来讨论其运动状况,仅反映了物体内部各部分之间的相互联系,且在区域内部相邻之间、相继时刻之间的这种联系(规律)通常与周围环境(边界上)和初始时刻对象(体系)所处的状态无关。
仅有方程还不足以确定物体的运动,因为外界的作用通常是通过物体边界“传”到内部的;一个方程可能有多个解,通解中含若干任意常数(函数),初始条件和边界条件就是确定它们的条件。
求一个微分方程的解满足一定初始条件和边界条件的问题称为定解问题:泛定方程& ⎧⎪⎪⎨⎪⎪⎩初始条件边界条件定解条件衔接条件自然条件。
1. 初始条件00(,)()(,)().t t t u x t x u x t x ϕψ==⎧=⎪⎨=⎪⎩,即已知初位移)(x ϕ和初速度)(x ψ 2. 边界条件i. 第一类边界条件-狄利克雷条件(Dirichlet 边界条件):直接给出了未知函数在边界上的值。
ii. 第二类边界条件-诺依曼条件(Neumann 边界条件):给出未知函数在边界上法向导数的值。
自由端点边界(端点不受外力,自由振动,意味着弦张力在振动方向无分量)属于此类,边界条件为(0,)0(,)0或x x u t u l t ==iii. 第三类边界条件-罗宾条件:给出未知函数和其边界法向导数在边界上的线性关系。
弹性支撑边界(端点受到弹簧的约束而无外力)属于此类,边界条件为:(,)(,)000x u t hu t -=Note :初始条件和边界条件是场运动规律的极限。
例1.对弦的横振动问题导出下列情况的定解条件:弦的两端点0=x 和l x =固定,用手将弦上的点(0)x c c l =<<拉开使之与平衡位置的偏离为h (l h <<),然后放手。
解:两端固定,所以边界条件为:(0,)0,(,)0u t u l t ==由点c x =的初始位移求出其他点的初始位移,它们是两段直线方程,容易求得:(0)(,0)()() ()h x x c c u x x h l x c x l l c ϕ⎧≤≤⎪⎪==⎨⎪-≤≤⎪-⎩, , 显然,初速度为零:(,0)0t u x =第二节 齐次方程混合问题的傅里叶解——分离变量法 本征值问题Abstract :求解数理方程定解问题的方法有分离变量法、行波法、积分变换法、变分法、复变函数论等,这些方法各有千秋。
分离变量法普遍适用,在其使用条件下,自然导致了问题的核心—本征值问题。
求解常微分方程:一般先求通解,再用初始/边界条件定其参数;求解偏微分方程,即使求得通解,亦难于由定解条件来定解(含任意函数)—本征值问题可解决此类问题。
7.2.1 利用分离变量法求解齐次弦振动方程的混合问题分离变量法:把二元函数(,)u x t 表示为两个一元函数相乘(,)()()u x t X x T t =⋅;然后带入函数的二阶偏微分齐次方程20tt xx u a u -=,把偏微分方程化为两个常微分方程;把偏微分方程的边界条件转化为常微分方程的边界条件。
题型I :方程和边界条件都是齐次的,而初始条件是非齐次的。
例题1:下面以两端固定弦的自由振动为例(第一类齐次边界条件):()20000 0,0; 0,(); ().tt xx x x l t t t u a u x l u u u x u x ϕψ====⎧-=<<⎪⎪==⎨⎪==⎪⎩ 注意这里的边界条件。
第一步, 分离变量,将二阶偏微分方程转化为两个常微分方程。
设)()(),(t T x X t x u =[取此特解形式,可得驻波解:()T t 是振荡函数,而与x 无关,()X x 是幅度函数,与t 无关],将此)()(),(t T x X t x u =代入泛定方程,即得2()()()().X x T t a X x T t ''''=等式两端除以)()(2t T x X a ,就有)()()()(2x X x X t T a t T ''=''. 注意在这个等式中,左端只是t 的函数,与x 无关,而右端只是x 的函数,与t 无关。
因此,左端和右端相等,就必须共同等于一个既与x 无关、又与t 无关的常数。
令这个常数为λ-(参数),即,λ-=''='')()()()(2x X x X t T a t T . 由此得到两个常微分方程:0)()(2=+''t T a t T λ (7.1)0)()(=+''x X x X λ (7.2)第二步,将(,)u x t 原来的边界条件转化为()X x 的边界条件。
将此(,)()()u x t X x T t =代入边界条件,得0)()0(=t T X ,0)()(=t T l X ,转化为()X x 的边界条件:0)0(=X ,0)(=l X [因为)(t T 不可能恒为0,否则),(t x u 恒为0] (7.3)这样就完成了分离变量法求解偏微分方程定解(亦定界)问题的前两步:分离变量。
在这两步中,假设所要求的是变量分离形式的非零解)()(),(t T x X t x u =,导出了函数)(x X 应该满足的常微分方程和边界条件,以及)(t T 所满足的常微分方程。
分离变量之所以能够实现,是因为原来的偏微分方程和边界条件都是齐次的(可分离变量)。
第三步,求解本征值问题上面得到的函数)(x X 的常微分方程定解问题,称为本征值问题。