第七章 一维波动方程的解题方法及习题答案
复变函数第7章
(7.2) (7.3) (7.4)
分离变量法 求解步骤:
1、分离变量: 2、求解特征值问题 3、求解定解问题
1、分离变量
i)分离变量形式解:设形式解为u(x,t)=T(t)X(x).
ii) 分离方程:将形式解代入泛定方程(7.2)得 TX = a2TX T X (为一常数) 即 2 a T X
r1 r2 r1 r2 r r i
二阶常系数微分方程的通解:
y C1e r1x C2e r2 x rx rx y C1e C2 xe y e x (C cos x C sin x) 1 2
X X 0 2、求解特征值问题 X (0) X (l ) 0
3、求解定解问题
由叠加原理得,方程(7.2)满足边界条件(7.3)的解为
n at n at n x u ( x, t ) (Cn cos Dn sin )sin l l l n 1
n x ( x) u ( x, 0) Cn sin l n 1 ( x) u ( x, 0) D n a sin n x t n l l n 1
X X 0 故得特征值问题 X (0) X (l ) 0
注:① 的值为该常微分方程边值问题的特征值(或本征值或固有值) ② 相应的非平凡解称为特征函数(或本证函数或固有函数) ③ 求特征值和特征函数的问题称为特征值问题(或本证函数问题或固有函数问题)
补充:
特征值问题是二阶常系数微分方程的求解问题, 所以考虑二阶常微分方程:y"+py'+qy=0的通解. 特征方程:r2+pr+q=0 的根分三种情况:
第七章 波动方程初值问题
x1 x0 at
即, f1(x - at) 表示波速为 a 的右行波
同理可知, f2(x + at) 表示波速为 a 的左行波. 因此,行波解为左行波与右行波的叠加. 三. 半无界弦的自由振动
utt a 2 uxx 0 u x0 0 u t 0 ( x ), ut
二. 行波解的物理意义 行波法的通解为:
u( x, t ) f1 ( x at ) f 2 ( x at )
对 f1(x - at),在 t0 时刻,x0 位置的波动位移为:
f1 ( x0 at0 )
若在t0+Δt 时刻, x1位置的波动位移也为 f1 ( x0 at0 ) 则:
t 0
a f1 ( x at ) x
f 2 ( x at ) t 0 a x
t 0
a f1 '( x ) a f 2 '( x ) y ( x )
对上式积分:
1 x x0 y ( )d [ f1 ( x ) f1 ( x0 )] [ f2 ( x ) f2 ( x0 )] (2) a
(1)
t 0
y ( x ) a f1 '( x ) a f 2 '( x )
1 x x0 y ( )d f1 ( x ) f 2 ( x ) c a
(2)
1 1 x c f1 ( x ) 2 [ ( x ) a x0 y ( )d ] 2 由 (1) (2) (x > 0) 解得: x f ( x ) 1 [ ( x ) 1 y ( )d ] c 2 2 a x0 2
数学物理方法课件第七章-----行波法
变量代换
x at
x at
2 u( , ) 0
a a u ( x, t ) 0 x t x t
u f1 ( ) f 2 ( )
行波法解题要领
• 行波法的提法来自于研究行进波。其解题要领为: • (1)引入特征变换,把方程化为变量可积的形式,从 而得到方程的通解; • (2)使用定解条件确定通解中的任意函数(对于常微 分方程为常数),从而得到其特解。 • 注意:由于偏微分方程求解较难,大部分偏微分方程 的通解均不易获得,使用定解条件确定其任意函数或 常数也绝非易事,故行波法也有其较大的局限性。但 是对于研究波动问题,行波法自有其独特的优点(实际 上我们主要只使用它研究波动问题)。因此行波法是求 解数学物理方程的基本的和主要的方法之一。
utt a u xx , ( Ⅰ )u |t 0 ( x) u | ( x) t t 0
2
- x
① ② ③
其中 ( x)和 ( x)为已知函数。
§7.1 行波法—一维波动方程的达朗贝尔解
解: 1 )做特征变换,求定解问题Ⅰ中方程①的通 () 一、达朗贝尔公式 dx 2 ①的特征方程为: ( ) a2 0 算符分解 dt ①式 dx dx a a u 0 x0 x 即( a )(t a) t dt dt 从而得到两簇特征线 (积分后得到 )如下: x a( ) t 坐标变换: x at c1 , x at c2 做特征变换 x at x at ④
§7.1 行波法—一维波动方程的达朗贝尔解
一、达朗贝尔公式 利用复合函数求导法则,有 u u u u u x x x
第七章一维波动方程的解题方法与习题答案
第七章一维波动方程的傅里叶解小结及习题答案第二篇数学物理方程——物理问题中的二阶线性偏微分方程及其解法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轴。
第七章 一维波动方程的解题方法及习题答案
第二篇 数学物理方程——物理问题中的二阶线性偏微分方程及其解法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)立假设:抓主要因素,舍弃次要因素,将问题“理想化”---“无理取闹”(物理趣乐)。
复变函数第二部分课后答案
⎧ utt = a 2u xx (1 < x < 2, t > 0) ⎪ ⎪ u (0, t ) = u (l , t ) = 0(t ≥ 0) ⎪ (0 ≤ x ≤ 1) ⎧ hx ⎨ ⎪ u ( x, 0) = ⎨ h(2 − x) (1 ≤ x ≤ 2) ⎩ ⎪ ⎪ ⎩ut ( x, 0) = 0
1
2
解:其付氏解为:
∞ u (r ,θ ) = A0 + ∑ ( An cos nθ + B n sin nθ )r n 2 n =1
,
α sin ϕ An = 1 n ∫02π f (ϕ )cos nϕdϕ = 1 2π A cos nϕ dϕ = nA π −α π ∫0 πl 其中:
= 2 A sin nα nπ
u rr + r u r + r uθθ = 0 。
⎧ + 1u + 1 u =0 ⎪u rr r r r 2 θθ ⎪ ⎨ ⎧ A, θ < α , (− π ≤ θ ≤ π ) ⎪u (1,θ ) = ⎪ ⎨ ⎪ 0, θ ≥ α ⎪ ⎩ ⎩ 2、 求解狄利克雷问题 , 其中 A,α 为
已知常数。
∞
0
2 ∞ − a 2 µ 2t e π ∫0
sin x π dx = x 2。 sin µ cos( µ x)d µ µ
u ( x, t ) = u (0, 0) =
2 sin µ e0 cos(0) d µ = 1 ∫ π µ ,
即:
2 ∞ sin µ dµ =1 π ∫0 µ
2 ∞ sin x ∫0 x dx = 1 令 x = µ ,则有: π ∞ sin x π dx = ∫ 0 x 2 得证。 即:
数学物理方法答案(完整版)
高等数学 第四册(第三版) 数学物理方法 答案(完整版)第七章 一维波动方程的傅氏解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 αβ<<<试求其傅氏解。
物理 波动例题(供学习参考)
(2) 质点振动的最大速度。
解 (1)a. 比较法(与标准形式比较) 标准形式 波函 ) 0 ] T 50 0.10 y 0.04 cos 2π ( t x) 2 2 2 T 0.04 s A 0.04 m 50 2 20 m u 500 m/s 0.10 T
3
解:设两质点的间距 x xb xa ,并设波的 表式为: 2 y A cos(2 t x 0 ) (1)
b 2 10 1.0
2
0.2 0
3
(3)
21
相位差为
5 0.1 b a ( xb xa ) 3 2 6
x 2 10 cos 100 t 200 2
2
将 x 4m 代入波的表达式,该处质点的振动位移
y 2 10 cos 100t 2 该质点在 t 2 s 时的振动速度为
2
dy 2 v 2 10 100 sin 200 6.28 m / s dt 2 16
2
所以波源首次回到平衡位置的时间为: y
2 / 3 3 / 2
t
3 2 1 2 3 t t 0 ( s) 100 120
t=0
12
例题5: 一平面简谐波沿Ox轴的负方向传播,波长 为 ,P处质点的振动规律如图所示。(1)求P处质 点的振动方程;( 2 )此波的波动方程;( 3 )若图 中d = /2,求坐标原点O处质点的振动方程。 y m
yo A cos t 2
14
《电磁场与电磁波》课后习题解答(第七章)
第7章习题解答【7.1】 解:设第一个分子的球心位置为原点,即0d (d 为分子直径)处 依题意任意时刻都要满足%5)10()0(0≤-E d d E E (1)其中E 是空间变化的电场,其形式为)exp(0ikx E -=E ,ck ω=,则(1)式变为%5)210exp(1≤--cfdi π (2) 可以求出 15151019.11056.1215⨯≈⨯≤f 所以频率上限的数量级为1510【7.2】解p V k ω=p pg p g p kdV dV d V V V dk dk V d ωωω===+ 1pg pp V V V d ωω=-22()1p i o rcc V n n ωωαω==-+0i n → p V c ∴= g p V V c ==即 2g p V V c ⋅=【7.3】解(1)波数681221501022310k f πππ===⨯⨯⨯⨯=⨯(rad/m ) 相速81.510p v ===⨯ (m/s )波长 21kπλ==(m )波阻抗60ηπ==(Ω) (2)均匀平面波的平均坡印廷矢量26z m S 0.26510z e e -==⨯平均 (W/m 2)得 31010m E -=⨯(V/m )当t = 0,z = 0时33sin 10100.8668.66103m E E π--⎛⎫==⨯⨯=⨯ ⎪⎝⎭(V/m )(3) t = 0.1s μ后210sin 23E ft kz ππ-⎛⎫=-+ ⎪⎝⎭267310sin 21501011028.66103z πππ---⎛⎫=⨯⨯⨯⨯-+=⨯ ⎪⎝⎭得 1sin 3028.66103z πππ-⎛⎫+-=⨯ ⎪⎝⎭15z =(m )【7.4】 解:电磁波的频率为8820310********v f λ-⨯===⨯⨯(Hz ) 在无损耗媒质中的波长为 12810vfλ-==⨯ (m ) 故波速为12888102510210v f λ-==⨯⨯⨯=⨯=(m/s )而无损耗媒质的本征阻抗为505000.1E H η==== (Ω) 联解以下两式:8210=⨯500= 得 1.99, 1.13r r με==【7.5】 解: 803100.2c f fλ⨯===故 883101510()0.2f Hz ⨯==⨯ 而 0.09vfλ== 故 880.090.091510 1.3510(/)v f m s =⨯=⨯⨯=⨯ 又v ===故 2882(/)(310/1.3510) 4.94r c v ε==⨯⨯=【7.6】 解:由题意知 7610ωπ=⨯0.8k π==106016E Hηππ====联解6100.8ππ⨯= 和60π= 得 8,2r r εμ==【7.7】 解:因4101σωε=<<,为低损耗媒质。
A7.行波法求解一维波动方程的初值问题—半无界问题
行波法求解一维波动方程的初值问题—半无界问题行波法求解一维波动方程的两个基本公式:1.达朗贝尔(d'Alembert )公式:⎰+-+-++=at x at x d aat x at x t x u ξξψφφ)(21))()((21),(; 2.Kirchhoff 公式:⎰⎰⎰----+-++-++=t t a x t a x at x at x d d f a d a at x at x t x u 0)()(),(21)(21))()((21),(τττξτξξξψφφ半无界弦的振动问题对于半无界域上波动方程初值问题的讨论,需要根据端点所处的物理状态不同分别加以讨论。
1. 端点固定(1)齐次端点条件 考虑定解问题.0,0,0,00),0(),()0,(),()0,(),,(2≥+∞<≤>+∞<<⎪⎩⎪⎨⎧===+=t x t x t u x x u x x u t x f u a u t xx tt ψφ求解上述问题的基本思路是以某种方式延拓函数,,,ψφf 使其在0<<∞-x 也有定义,这样把半无界区域+∞<≤x 0上的问题转变为+∞<<∞-x 上的初值问题。
然后利用达朗贝尔公式,求出在+∞<<∞-x 上的解),(t x u 。
同时使此解),(t x u 满足0),0(=t u 。
这样当x 限制在+∞<≤x 0上就是我们所要求的半无界区域+∞<≤x 0上的解。
由微积分知识可知,如果一个连续可微函数)(x g 在),(+∞-∞上是奇函数,则必有0)0(=g 。
因此,要使解),(t x u u =满足0),0(=t u ,只要),(t x u 是x 的奇函数便可。
因此对函数ψφ和,f 关于x 作奇延拓。
我们定义)()(),,(x x t x F ψΦ和如下:⎩⎨⎧≥<--≥≥⎩⎨⎧<≥--=ψ<≥⎩⎨⎧--=Φ.0,0),,(,0,0),,(),(.0,0),(),()(.0,0),(),()(t x t x f t x t x f t x F x x x x x x x x x x ψψφφ 显然函数在和)()(),,(x x t x F ψΦ+∞<<∞-x 上是奇函数。
大学物理2-1第七章(波动光学)习题答案
习 题 七7-1 如图所示,O S O S 21=。
若在O S 1中放入一折射率为n ,厚度为e 的透明介质片,求O S 1与O S 2之间的光程差。
如果1S 和2S 是两个波长为λ的同相位的相干光源,求两光在O 点的相位差。
[解] O S 1与O S 2的几何路程相等 光程差为()e n 1-=δ 位相差为()e n 122-==∆λπδλπϕ7-2 一束绿光照射到两相距 0.6mm 的双缝上,在距双缝2.5m 处的屏上出现干涉条纹。
测得两相邻明条纹中心间的距离为2.27mm ,试求入射光的波长。
[解] 由杨氏双缝干涉知,dD x λ=∆ 所以5448m 10448.55.21060.01027.2733=⨯=⨯⨯⨯=∆=---D xd λÅ7-3 如图所示,在双缝干涉实验中,21SS SS =,用波长为λ的单色光照S ,通过空气后在屏幕E 上形成干涉条纹。
已知点P 处为第3级干涉明条纹,求1S 和2S 到点P 的光程差。
若整个装置放于某种透明液体中,点P 为第4级干涉明条纹,求该液体的折射率。
[解] 1S 和2S 到P 点的光程差满足λλδ312==-=k r r 整个装置放置于液体中,1S 和2S 到P 点的光程差满足()λδ412=-=r r nλλ43=n 所以得到 33.134==n7-4 如习题7-1图所示,1S 和2S 是两个同相位的相干光源,它们发出波长λ=5000Å的光波,设O 是它们中垂线上的一点,在点1S 与点O 之间的插入一折射率n =1.50的薄玻璃,点O 恰为第4级明条纹的中心,求它的厚度e 。
[解] 在O 点是第4级明条纹的中心 光程差 λδ4=-=e ne所以 410414⨯=-=n e λÅ7-5 初位相相同的两相干光源产生的波长为6000Å的光波在空间某点P 相遇产生干涉,其几何路径之差为6102.1-⨯m 。
如果光线通过的介质分别为空气(11=n )、水(=2n 1.33)或松节油(=3n 1.50)时,点P 的干涉是加强还是减弱。
7-1 一维波动方程的达朗贝尔公式 chen
t
( x, t )
依赖区间
O x − at
x + at x
退出
数学物理方程与特殊函数 主页 上一页 下一页
t
x=
t
t
x1
x=
O
x1
x2
x
O
x1
x2
决定区间
影响区间
在区间[ x1 , x2 ]上给定初始条件,就可以在其决定 区间域中决定初值问题的解.
数学物理方程与特殊函数 主页 上一页 下一页 退出
退出
因为在特征线 x − at = C 2 上,右行波 u2 = f 2 ( x − at ) 的振幅取常数值 f 2 (C 2 ) ,在特征线 x + at = C1 上,左行波 u 1 = f1 ( x + at )的振幅取常数值 f1 (C1 ) ,且这两个数值随 特征线的移动(即常数 C i ( i = 1, 2) 的改变)而改变,所以波 动实际上是沿特征线传播的.
⎧ uxx + 2uxy − 3u yy = 0 ⎪ ⎨ 2 ⎪u | y=0 = 3 x , uy | y=0 = 0 ⎩
y > 0, −∞ < x < +∞
−∞ < x < +∞
先确定所给方程的特征线.为此写出它的特征方程: 它的两族积分线为:
(dy )2 − 2dxdy − 3(dx )2 = 0
主页上一页下一页退出数学物理方程与特殊函数由此可以看出在xt平面上斜率为的两族直线常数对一维波动方程的研究起着重要的作用我们称其为一维波动方程ttxxttxx主页上一页下一页退出数学物理方程与特殊函数因为在特征线上右行波的振幅取常数值在特征线上左行波的振幅取常数值且这两个数值随特征线的移动即常数的改变而改变所以波动实际上是沿特征线传播的
第七章数理方程教材
一般说来,任何一个本征解都不能单独满足初 始条件,因此本征解并不是定解问题的解。
为了获得满足初始条件的解,通常要将本征解 进行线性叠加,从而形成如下的通解式:
可以证明,通解式既满足微分方程,又满足边 值条件。若要使其满足初始条件,那么
(x)和(x)的傅氏展开
根据以上初始条件,可以进一步确定通解式中 待定常数
再假设初始条件为 那么完整的定解问题为:
小结:
1. 定解问题: 描述物理现象的偏微分方程+定解条件; 2. 微元法建立偏微分方程: 在系统中任选一微元,将有
关的物理定律用于这一微元,建立它的运动方程.然 后取趋向于无穷小的极限,保留最低阶小量,略去高 阶小量,就可得到所需的偏微分方程; 3. 定解条件: 边界条件+初始条件(+附加条件);
则w(x)必须满足条件:
求解以上定解问题很容易求出:
根据 v(x,t)定解问题中的初始条件,就可以 确定待定系数
§7.5 有阻尼的波动问题 例10 两端固定弦的小阻尼振动问题
f (x,t)x
弦在振动过程中所受阻力一般正比于速率。 ( , 为常数)
类似于本章例1的推导可以得到:
(阻尼因子)
解:采用分离变量法,设
代入边界条件后得: ,若要使 ,那么
相应的本征函数为: 因此该问题的本征解为:
管乐器中空气的本征振动角频率为: 当n=0时,对应于最低频率ν 0(基频)。 当n>1时,相应的本征振动频率是n次谐频。
管乐器声音中只有奇次谐频,没有偶次谐频。
分离变量法解题的四步:
1. 设具有分离变量法的试探解,并代入偏微分方程和边界条 件,从而化为几个常微分方程(必需有一个方程构成本征 值问题)和相应的边界条件;
第七章一维波动方程的傅氏解
第七章 一维波动方程的傅氏解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就得到这个系统的波动方程:
在这个例子中,波速。
数理方法 第七章 一维波动方程的付氏解
第二节 齐次方程混合问题的分离变数解法 §7.2.1分离变数法 例:两端固定的弦的自由振动
utt
a2uxx ,(0
x
l,t
0)
u(0,t) 0,u(l,t) 0,(t 0)
(7.1) (7.2)
u(
x
,
0)
(
x
),
ut
(
x
,
0)
(
x
),
(0
x l)
(7.3)
代表的是驻波,称为弦的本征振动,其波矢量为kn ,
圆频率为n ,1 称为基频.其余称为谐频(泛音)
波长为:n
2
a
n
2
kn
2l n
n次谐波(泛音)波节的位置在:
sin n x 0 x m l m n ,(m 0,1, 2,L , n)
l
n2
相邻波节间距 l / n n / 2
初始条件通常通过未知函数及其导数在自变量的 一个特定点的值给出。
如:
y f ( x, y),
(1).
y
|x0
2
y py qy f ( x)
(2).
y
|x0
3
y
|x0
4
2.弦振动方程的定解条件
2.1 初始条件 类似于常微分方程定解过程的初值。
nl
l
0
(
)sin kn
d
2
n a
l
0
2021年数理方法课件 精美PPT 07第7章 一维波动方程的傅里叶解
两边除以 ρΔx, 然后取极限Δx→0:
utt ( x, t) =
T0
uxx(x, t) +
f
( x, t )
➢ 弦振动的泛定方程 ut t = a2ux x + f
u(x,t) ➔ x 处的质元在 t 时刻相对平衡位置的位移 f(x,t) ➔ t 时刻 x 处单位质量所受的横向外力
a = T0 / : 弦中横波的波速 T0 ➔ 初始张力,ρ➔质量线密度
• 整个系统初始状况的表达式称为初始条件 • 对弦振动,泛定方程为 ut t = a2 2u + f
需给出弦在初始时刻 t=0 的位移和速度: u( x,0) = ( x), ut ( x,0) = ( x)
• 泛定方程出现时间的 n 阶偏导数时需要 n 个初始条件
• 对物理量的稳态分布,无初始条件
lxntalnbtalnatxunnn??????sinsincos1??????????????lndxlxnxula0sin02??lxnblanxunnt????sin001????????????????202sincos14ldxlxnxnlh??????????????????????lxlxlhhlxxlhxu22220200ll2xux0h分离变量法得出解的一般形式bn0lxntalnnnhtxun????????sincossin????????12228222sincos14lnnnlhan??????????????????2020cos1coslldxlxnlnxlnlxni????????????20sinldxlxnxi??对奇数n计算积分202sin2cos2llnlxnnlnl??????????????22sinlnn??????2sin822????nnh???回顾
一维波动方程的解法
一维波动方程的解法波动现象是自然界和人类生活中广泛存在的一种现象,它具有许多重要的物理意义,例如声波、光波等。
一维波动方程是描述波动的重要方程之一,本文将介绍一维波动方程的解法。
一、一维波动方程的基本形式和意义一维波动方程的基本形式为:$$\frac{\partial^2 u}{\partial t^2}-c^2\frac{\partial^2 u}{\partialx^2}=0$$其中,$u(x,t)$表示波动的幅度,$c$表示波速。
这个方程描述了介质中的一种波动现象:波动传播速度为$c$,波动在媒质中沿$x$轴方向的传播,波动的幅度随时间$t$的变化而变化。
在声波和电磁波中,$u$分别是空气压力和电场强度,$c$分别是声速和光速。
二、1. 分离变量法分离变量法是一种基本的解法,其思想是将波动方程中的未知函数$u(x,t)$表示成仅包含$x$的函数和仅包含$t$的函数的乘积形式:$$u(x,t)=X(x)T(t)$$将$u(x,t)$代入一维波动方程中,得到:$$\frac{T''(t)}{c^2T(t)}=\frac{X''(x)}{X(x)}=-\lambda^2$$其中,$\lambda$是一个常数。
由此可得到两个关于未知函数的简单微分方程:$$T''(t)+\lambda^2c^2T(t)=0$$和$$X''(x)+\lambda^2X(x)=0$$其中,第一个微分方程的解为:$$T(t)=A\cos(\lambda ct)+B\sin(\lambda ct)$$其中,A、B是常数。
第二个微分方程的解为:$$X(x)=C\cos(\lambda x)+D\sin(\lambda x)$$其中,C、D是常数。
因此,一维波动方程的通解为:$$u(x,t)=\sum_{n=1}^\infty a_n\cos(\lambda_n x)\cos(\lambda_n ct)+b_n\sin(\lambda_n x)\sin(\lambda_n ct)$$其中,$\lambda_n=n\pi/L$,$L$为介质的长度,$a_n$和$b_n$是待定常数。
一维波动方程推导
①下行波中质点的运动方向与所受力的方向始终一致;
②上行波中质点的运动方向与所受力的方程始终相反。
第1第3十三页页,/共共444页。3页
(23)
(24)
3 通过计算可以分离出桩身各截面
上行波和下行波的值,具体如下:
F F F
(25)
V V V
将式(20)和式(24)带入式(25),即得
1 F ( F ZV )
,即得著名的一维波动方程
c2
2 x
x 2
2u t 2
第第6六页页,/共共444页3。页
(1) (2)
(3)
1.2.2一维波动方程的解
• 求解一维波动方程有多种方法,常用的有行波法、分离 变量法、特征线法,这里主要介绍基桩检测常用的行波 法。
作变量代换: x ct
x ct
(4)
第第7七页页,/共共444页3。页
2003。
E. 土阻力的发挥在一定范围内与桩土之间的相对位移成正比,因此测试时,桩必须获得永久性 的贯入度,如果桩没有被打动,或者贯入度极小,则得到的承载力仅仅是一种激发值,就好 像不做到破坏的静载试验一样。
第第2二4十页四页/,共共444页3。页
1.2.7 去除土阻尼的影响
CASE法分析时,假设土的动阻力集中在桩尖,桩周不存在动阻力。动阻力与桩土相对运动 速度和阻尼有关,CASE法中一般假设动阻力与速度和阻尼成正比,即
的性质是否变化。
(2)应力波反射和透射能力大小取决于两种介质波阻抗的差异情况。
两种介质波阻抗相差愈大,反射能力愈大,透射能力愈小;波阻抗相 等,只有透射没有反射。
(3)波的性质:压力波,拉力波的变换
(4)当桩身缩径、夹泥、离析、断桩等缺陷时,Z2<Z1,入射波与反射 波同号;桩身存在扩径时,Z2>Z1,入射波与反射波反号;
物理学7章习题解答
[物理学7章习题解答]7-2 一个运动质点的位移与时间的关系为m ,其中x的单位是m,t的单位是s。
试求:(1)周期、角频率、频率、振幅和初相位;(2) t = 2 s时质点的位移、速度和加速度。
解(1)将位移与时间的关系与简谐振动的一般形式相比较,可以得到角频率s 1, 频率, 周期, 振幅,初相位.(2) t = 2 s时质点的位移.t = 2 s时质点的速度.t = 2 s时质点的加速度.7-3 一个质量为2.5 kg的物体系于水平放置的轻弹簧的一端,弹簧的另一端被固定。
若弹簧受10 n的拉力,其伸长量为5.0 cm,求物体的振动周期。
解根据已知条件可以求得弹簧的劲度系数,于是,振动系统的角频率为.所以,物体的振动周期为.7-4求图7-5所示振动装置的振动频率,已知物体的质量为m,两个轻弹簧的劲度系数分别为k1 和k2。
解 以平衡位置o 为坐标原点,建立如图7-5所示的坐标系。
若物体向右移动了x ,则它所受的力为.根据牛顿第二定律,应有,改写为.所以,.7-5 求图7-6所示振动装置的振动频率,已知物体的质量为m ,两个轻弹簧的劲度系数分别为k 1 和k 2。
解 以平衡位置o 为坐标原点,建立如图7-6所示的坐标系。
当物体由原点o 向右移动x 时,弹簧1伸长了x 1 ,弹簧2伸长了x 2 ,并有.物体所受的力为,式中k 是两个弹簧串联后的劲度系数。
由上式可得 , .于是,物体所受的力可另写为,由上式可得 ,所以. 图7-5 图7-6装置的振动角频率为,装置的振动频率为.7-6仿照式(7-15)的推导过程,导出在单摆系统中物体的速度与角位移的关系式。
解由教材中的例题7-3,单摆的角位移θ与时间t的关系可以写为θ = θ0 cos (ω t+ϕ) ,单摆系统的机械能包括两部分, 一部分是小物体运动的动能,另一部分是系统的势能,即单摆与地球所组成的系统的重力势能.单摆系统的总能量等于其动能和势能之和,即,因为, 所以上式可以化为.于是就得到,由此可以求得单摆系统中物体的速度为.这就是题目所要求推导的单摆系统中物体的速度与角位移的关系式。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第二篇 数学物理方程——物理问题中的二阶线性偏微分方程及其解法Abstracts:1、根据物理问题导出数理方程—偏微分方程;2、给定数理方程的附加条件:初始条件、边界条件、物理条件(自然条件,连接条件),从而与数理方程一起构成定解问题;3、方程齐次化;4、数理方程的线性导致解的叠加。
一、数理方程的来源和分类(状态描述、变化规律)1、来源I .质点力学:牛顿第二定律F mr =r r && 连续体力学2222()(,)(,)0(()0;v 1()0(Euler eq.).u r t a u r t t v t v v p f t ρρρ⎧⎧∂⎪⎪-∇=⎨⎪∂⎪⎪⎩⎪∂⎪+∇⋅=⎨∂⎪∂-⎪+⋅∇=+=⎪∂⎪⎪⎩r r r r r r r r &弹性定律弦弹性体力学杆 振动:波动方程);膜流体力学:质量守恒律:热力学物态方程: 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 σρτρσ⎧⋅=⇒∇⋅=⋅=⋅⇒∇⨯=⎪⎪⎪⋅=⇒∇⋅=⋅=+⋅⇒∇⨯=+⎨⎪=-∇=∇⨯⎪⇒⇒⎪⎩⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰r r r r r r r r r &&r r r r r r r r r r r &&r r r r 已已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∂=-∇+∂h h2. 分类物理过程方 程 数学分类 振动与波波动方程222210u u a t ∂∇-=∂ 双曲线 输运方程20u k u t ⎧∂-∇=⎨∂⎩能量:热传导质量:扩 散 抛物线 稳态方程 Laplace equation 20u ∇=椭圆型 二、数理方程的导出推导泛定方程的原则性步骤:(1)定变量:找出表征物理过程的物理量作为未知数(特征量),并确定影响未知函数的自变量。
(2)立假设:抓主要因素,舍弃次要因素,将问题“理想化”---“无理取闹”(物理趣乐)。
(3)取局部:从对象中找出微小的局部(微元),相对于此局部一切高阶无穷小均可忽略---线性化。
(4)找作用:根据已知物理规律或定律,找出局部和邻近部分的作用关系。
(5)列方程:根据物理规律在局部上的表现,联系局部作用列出微分方程。
Chapter 7 一维波动方程的傅里叶解第一节 一维波动方程-弦振动方程的建立弦横振动方程的建立(一根张紧的柔软弦的微小振动问题)(1)定变量:取弦的平衡位置为x 轴。
表征振动的物理量为各点的横向位移),(t x u ,从而速度为t u ,加速度为tt u .(2)立假设:①弦振动是微小的,1<<α,因此,sin tan ααα≈≈,1cos ≈α,又tan u x αα∂=≈∂Q ,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 T u 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 =∇+(非齐次方程)定解问题第一节从物理问题和相应的物理定律导出了其所满足的偏微分方程,但总是选择物体内部,不含端点或边界,对一小部分来讨论其运动状况,仅反映了物体内部各部分之间的相互联系,且在区域内部相邻之间、相继时刻之间的这种联系(规律)通常与周围环境(边界上)和初始时刻对象(体系)所处的状态无关。
仅有方程还不足以确定物体的运动,因为外界的作用通常是通过物体边界“传”到内部的;一个方程可能有多个解,通解中含若干任意常数(函数),初始条件和边界条件就是确定它们的条件。
求一个微分方程的解满足一定初始条件和边界条件的问题称为定解问题:泛定方程& ⎧⎪⎪⎨⎪⎪⎩初始条件边界条件定解条件衔接条件自然条件。
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 :求解数理方程定解问题的方法有分离变量法、行波法、积分变换法、变分法、复变函数论等,这些方法各有千秋。
分离变量法普遍适用,在其使用条件下,自然导致了问题的核心—本征值问题。
求解常微分方程:一般先求通解,再用初始/边界条件定其参数;求解偏微分方程,即使求得通解,亦难于由定解条件来定解(含任意函数)—本征值问题可解决此类问题。
利用分离变量法求解齐次弦振动方程的混合问题分离变量法:把二元函数(,)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 λ ()0)()(=+''x X x X λ ()第二步,将(,)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] ()这样就完成了分离变量法求解偏微分方程定解(亦定界)问题的前两步:分离变量。