波动方程和行波法
第七章 行波法(一)

第七章 行波法
利用初值条件确定函数 F,G
u( x,0) ( x)
ut ( x,0) ( x)
F ( x) G ( x) ( x)
a[ F ( x) G( x)] ( x)
x
a[F ( x) G( x)] C ( )d
x0
其中
x
x1
x2
内,因此该三角区域称为
决定区域。
结论:达朗贝尔解表示沿x 轴正、反向传播的两列波速 为a波的叠加,故称为行波法。
第七章 行波法
影响区域、依赖区间、决定区域
波动是以一定的速度 a 向两个方向传播的。
如果在初始时刻 t=0,扰动仅仅在有限区间 [ x1 , x2 ] 上存在,则经过时间 t 后,扰动传到的范围为
x1 at x x2 at
第七章 行波法
无界弦振动的初值问题
2 2u 2 u x 2 a 2 x t u ( x, 0) ( x), u ( x, 0) ( x) t
第七章 行波法
2. 行波法的基本思想
这种方法是针对波动方程提出的。由于波动现象的普
1 过 x1 作斜率为 的直线 x x1 at a 1 过 x2 作斜率为 的直线 x x2 at a t 则 它们与区间 [ x1 , x2 ]
一起围成的三角形区域 中的任意一点 ( x, t ) 的 依赖区间都落在区间 [ x1 , x2 ]
x x1 at
x x2 at
遍性,对如何认识和解决波动问题,一直是物理学家和数 学家们长期探索的课题。 (1)波函数可写成位置和时间函数的分离形式,且波函数
是由无穷多个谐波分量叠加而成的,由此提出了分离变量
波动方程初值问题与行波法

1 x at 1 u d 2 2a x at 1
1 arctan( x at ) arctan( x at ) 2a
例4: 求二阶线性偏微分方程初值问题的解
uxx 2uxy 3u yy 0 2 u | 3 x , u y | y 0 0 y0
2 F 3 x G x 3 x F ' 3 x G ' x 0
1 F 3x G x C 3
9 2 F 3x x C ' 4 G x 3 x 2 C ' 4
P( x, t )
依赖区间
x at
x at
x
区间 [ x at , x at ] 为解的依赖区间。
2.决定区域 该区域中任一点(x, t )的依 赖区间都落在区间[c, d]内 部,因此解在此该区域中的 数值完全由区间[c, d]上的 初始条件决定。
t
x c at
x d at
例5 求二阶线性偏微分方程的通解
uxx 2sin xuxy cos xuyy 0.
2
解:特征方程为
dy
2
2sin xdxdy cos x dx 0
2 2
dy dy 1 sin x 1 sin x 0 dx dx
G(x-at)=G(x0+at-at)=G(x0)
u2 G ( x ) ( t 0)
O
at
u2 G ( x at ) ( t t0 )
x0
x x0 at
x
u1 F ( x at )
第七章 波动方程初值问题

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
波动方程和行波法剖析课件

目录 Contents
• 波动方程的基本概念 • 行波法的基本原理 • 波动方程的解析解法 • 波动方程的数值解法 • 行波法的应用实例
01
波动方程的基本概念
பைடு நூலகம்
波动方程的定义
波动方程
描述波动现象的基本数学模型,通常 用于描述物理场(如声场、电磁场、 水波等)随时间和空间的变化规律。
03
最后,通过迭代求解差分方程 ,得到波在每个网格点上的值 ,从而得到波的传播和演化过 程。
行波法的优缺点
优点
行波法简单易懂,易于编程实现,能够处理复杂的边界条件和初始条件,适用 于求解各种类型的波动方程。
缺点
行波法需要设定初始条件和边界条件,对于某些复杂的波动问题可能需要较高 的计算成本和精度要求。
水波传播的模拟
要点一
总结词
利用行波法模拟水波的传播,有助于研究水波的形成、演 化及对环境的影响。
要点二
详细描述
在水波传播的模拟中,行波法能够模拟水面的波动情况, 包括波浪的生成、传播和消散。通过调整参数,可以研究 不同条件下水波的传播规律,如风速、水深、地形等,对 于水文学、海洋学等领域具有重要意义。
03
波动方程的解析解法
分离变量法
将波动方程的解表示为若干个变量的 乘积或商的形式,以便分别求解。
VS
分离变量法是一种求解波动方程的常 用方法。通过假设波动方程的解可以 表示为若干个变量的乘积或商的形式 ,我们可以将一个复杂的偏微分方程 转化为若干个简单的常微分方程,从 而方便求解。
积分变换法
利用积分变换将波动方程化为易于求解的形式,再进行逆变换得到原方程的解。
地震学
用于模拟地震波的传播和反射,进行地震预 测和地球结构研究。
波动方程的特征线法

作变换 1 ( x, y ), 2 ( x, y ),
在区域Ω上作此变换下,可化简方程(1),甚至可求得其解. 此变换称为特征变换.
例1 一端固定的半无界弦的自由振动问题
2u 2u a2 0 ( t 0,0 x ), 2 2 x t u t 0 : u ( x ), ( x ) ( 0 x ), t 0 t t 0 x 0 : u 0.
举例
2u 2u a 2 2 , x R, t 0 t 2 x u ( x, 0) 1, xR 2 ut ( x, 0) x ,
例4:
例5:
2u 2u a2 2 , 2 x t u ( x,1) cos x, ut ( x,1) 0,
例2:
2u 2u 2 a2 2 0 t x u t 0 cos x, ut t 0 x
解:由达朗贝尔公式
( x at) ( x at) 1 x at u ( x, t ) ( )d 2 2a xat
cos(x at) cos(x at) 1 x at d 2 2a x at
此公式的意义在于把定解问 题的解表示为左、右行进波 相叠加,这种方法称为“行 波法”。
D’Alembert公式
1 1 x at u ( x, t ) ( x at ) ( x at ) ( )d 2 2a x at
注 : 当 ( x ) C 2 ( R ), ( x ) C 1 ( R )时, 初值问题( I )存在唯一的解 u( x , t ),由d ' Alembert 公式给出.
数理方程3.1 一维波动方程的初值问题

§3.1 一维波动方程的初值问题
21 5 x − y = C1 x + y = C , 2
特征线为
作变量代换
ξ = 5x − y η = x + y,
原方程化为
uξη = 0
其通解
u = F (5x − y ) + G(x + y )
利用初始条件可得
F (5x) + G(x) = 5x, −F ′ (5x) + G′ (x) = 0
utt = a2 uxx , u(x, 0) = 1 , ut (x, 0) = 0, 1 + 4 x2 − ∞ < x < +∞, t > 0, − ∞ < x < +∞
1 . 1+4x2
从静止开始
由达朗贝尔公式得
1 1 1 1 1 u(x, t) = [ϕ(x + at) + ϕ(x − at)] = · + · 2 2 2 1 + 4(x + at) 2 1 + 4(x − at)2
其中ϕ(x), ψ (x)分别表示初值位移和初始速度. 1. 泛定方程的通解 x 2 dx 2 采用第2.1节中的方法, 特征方程为( d (特征 dt ) − a = 0, 特征线方程为 dt = ±a, 其通解 线)为x + at = C1 , x − at = C2 , 作变换
ξ = x + at, η = x − at,
第三章
波动方程的初值问题与行波法
§3.1 一维波动方程的初值问题
本节思路: 无界弦的自由振动(utt = a2 uxx , −∞ < x < +∞): 经非奇异变换化为标准型后直接积分得通解, 代入初始条件得特解(达朗贝尔公式) 无界弦的受迫振动(utt = a2 uxx + f, −∞ < x < +∞): 由叠加原理分解为: 齐次问题+零初值的非齐次问题(由齐次化原理得解) 半无界弦的振动问题(utt = a2 uxx + f, 0 < x < +∞): 以某种方式延拓f 及初始函数, 转成无限长的弦的振动, 求出解后限制在半无界区域上.
数学物理方程:第3章 波动问题的行波法

第3章 波动问题的行波法§3.1 二阶线性方程的分类与化简本节讨论:①两个自变量方程的分类与化简,②多个自变量方程的分类与化简⒈ 两个自变量方程的分类与化简二阶方程的一般形式 二阶变系数方程可写为1112220(,)2(,,,,)(,)xx xy yy x y Lu x y a u a u a u x y u u u f x y =+++Φ= (3.1.1)式中:11a 、12a 、22a 为x 、y 的函数,0(,,,,)x y x y u u u Φ为低阶导数项。
公式关于二阶导数项为线性的,即称方程为准线性的。
若0(,,,,)x y x y u u u Φ关于u 及其x u 、y u 为线性的,则称方程为线性的。
方程的变换 为了简化上述方程,作可逆变换:(,)(,)x y x y ξξηη=⎧⎨=⎩, (,)0(,)J x y ξη∂=≠∂, (,)(,)x x y y ξηξη=⎧⎨=⎩(3.1.2) 代入方程中,不难得到:11122212(,,,,)(,)Lu A u A u A u u u u f ξξξηηηξηξηξη=+++Φ= (3.1.3)式中: 22111112222x x y yA a a a ξξξξ=++ (3.1.4) 12111222()x x x y y x y y A a a a ξηξηξηξη=+++ (3.1.5)22221112222x x y yA a a a ηηηη=++ (3.1.6) 我们化简的目的是使得二次项的项数尽量少,并且值尽量为简单(如0ij A =或1ij A =±)。
顾及ij A 的表达式,取关于z 的一阶非线性偏微分方程2211122220x x y y a z a z z a z ++= (3.1.7)若该方程有解),(1y x z ϕ=、),(2y x z ψ=,则110A =及220A =;公式大大简化了。
《数学物理方程》第四章§1

2/16
2u 2u a2 2 t 2 x
2u 2u a2 2 0 t 2 x
2 2 ( 2 a 2 2 )u 0 t x
0 1 0 a 2
dx 令 dt
2 a 2 0
a
x at x at
t t a a 1 1 x x
0 a a a 1 1 a 1 0 a 2 1 1
0 2a 2
《百科全书》不仅在于提供知识,而更重要的在 于改变读者的思想。 向前进,你就会产生信念 ————达朗贝尔
达朗贝尔脱下了微分学的神秘外衣 ————马克思
5/16
2u 2u a2 2 t 2 x
u( x , t ) = f1(x + at ) + f2(x – at )
u t 0 u ( x ), t
2a 2 0
3/16
0 1 0 a 2
0 2a 2
2a 2 0
( t , x ) ( , )
2u 2u a2 2 0 2 t x
2u 0
2u 4a 2 0
x , x [0,1 / 2] ( x ) 1 x , x [1 / 2,1] 0, 其它
随着时间的推移, u2 的图形以速度 a 向x 轴正方向 移动. 所以,u2表示一个以速度a 沿 x 轴正方向传播 的行波,称为右行波。
8/16
2u 2u a2 2 t 2 x u u t 0 ( x ), t
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
f ( x0 , y0 , z0 , t )
其中 f ( x0 , y0 , z0 , t ) 为已知函数。
35
第二类边界条件(Neuman 边界条件):
规定所研究物理量在边界外法线方向 n 上的
方向导数的数值.
u f n
u f ( x0 , y0 , z0 ) , n ( x0 , y0 , z0 )
36
第三类边界条件(混合边界条件 也叫 Robin边界条件 ):规定所研究物理量及其
外法向导数的线性组合在边界上的值
u Hun
( x0 , y0 , z0 )
f ( x0 , y0 , z0 , t )
u f H :常系数 u n
37
以上三类边界条件当 f 0 时,分别称为 第一、二、三类齐次边界条件。
22
应用微积分中值定理:
ux ( x dx, t ) T1ux ( x, t ) uxx dx
dy f ' ( x)dx
ux ( x dx, t ) T1ux ( x, t ) uxx dx
x Fdx dxutt
Tuxx dx Fdx dxutt
39
2 a u tt u xx 0 在这一点无意义.如果,将
l 分成 x x0 ,x x0 两段分别考虑,
在各段上,弦振动方程有意义,但它是一 根弦的两段,并不是各自振动的。从数学
上来讲,不可能在两端上分别列出定解问
题。两段可作为一个整体来研究,两段的 振动是相互关联的。
40
u
F(0,t)
15
即整根弦由相互牵连的质点组成,对每个
质点即每个小段可应用 F ma . 方法:将连续分布的介质离散化为多质点
系统,再取内部任一代表性的点进行研究。将
弦细分为许多极小的小段,取区间上 ( x, x dx) 小段为代表。无质量且柔软,故该段仅受到相
邻两段的拉力 T1 和 T2 .
16
④ 对弦的每一小段dx,沿x方向(纵向) 没有运动,沿 x方向所受合外力为零。任一
表征物理量的选择常常是建立一个新 方程的起点。 (一个或几个)。
5
3.寻找(猜测)物理过程所遵守的 物理定律或物理公理;
4.写出物理定律的表达式,即数学
模型。
6
1.1 弦振动方程
一、弦的横振动方程 二、定解条件的提出
三、三类定解问题
7
一、 弦的横振动方程(均匀弦的微小横振动) 演奏弦乐(二胡,提琴)的人用弓在弦上来回 拉动,弓所接触的是弦的很小的一段,似乎只能引 起这个小段的振动,实际上振动总是传播到整个弦,
29
2、初始条件 在求解含时间t变量的数理方程时,往往要追
溯到早些某个所谓“初始”时间的状况(“历
史” ),于是称物理过程初始状况的数学表达 式为初始条件。
30
如弦振动方程:
utt a uxx 0
2
其初始条件为: 同一时刻( t 0 )情况
u ut
t 0 t 0
( x) ( x)
α1
α2
xLeabharlann 41x x0 虽是折点,但它们连续,即
u( x0 0, t ) u( x0 0, t )
在 x0 ,力 F (t ) 应和张力平衡,即
F (t ) T sin 1 T sin 2 0
①
sin 1 tan 1 ux ( x0 0, t )
sin 2 tan 2 ux ( x0 0, t )
Tux ( x0 0, t ) Tux ( x0 0, t ) F (t )
②
①、②合称为衔接条件,这时振动问题适定。
42
再如,不同材料组成的杆的振动,在
衔接处的位移和能量相等,即:
u1
E1u1 x
x x0
u2
E1u2 x
x x0
x x0
x x0
u1 ( x, t ),u2 ( x, t ) :杆的两部分位移.
动。二维波动方程,如薄膜的横振动方程,管
道中小振动的传播,理想传输线的电报方程等
均可用上述波动方程描述。故称为一类方程,
即波动方程。(也是称其为泛定方程的远大) 可描述一类物理现象。流体力学与声学中推导 三维波动方程,这里不再一一推导。
27
二、定解条件的提出
1、必要性。导出方程后,就得对方程进行 求解。但是只有泛定方程不足以完全确定方程 的解,即不足以完全确定具体的物理过程,因 为具体的物理过程还与其初始状态及边界所受 的外界作用有关,因而必须找一些补充条件, 用以确定该物理过程。
E1 , E2 :两部分的杨氏模量.
43
静电场中,两种电介质的交界面
s
上电势应相等(连续),电位移矢量的法
向分量也应相等(连续),其衔接条件是:
u1 s u1 s u1 1 n u2 2 n
s
s
44
其中 u1 , u 2 代表两种电介质的电势,
1 , 2 代表两种电介质的介电常数,(设电
utt Tuxx F
23
即
utt a2uxx f
—— 弦的强迫横振动方程
T
其中: a
2
f ,
F
量纲分析:T : MLT 2 , : ML1
24
∴
2
T
MLT 2 2 2 : L T 1 ML
2
即
a :LT
2
a
:振动的传播速度 a
T
它与弦的张力的平方根成正比,与弦的
38
4、其它条件
⑴ 衔接条件 由于一些原因,在所研究的区域里出 现跃变点,泛定方程在该点失去意义。如 波动方程(弦),如果有横向力F (t )集中地 作用于 x0 点, 这就成了弦的折点。在点 x0 斜率 u x 的左极限 ux ( x0 0, t ) 不同于右极限
ux ( x0 0, t ),因而 u xx不存在,
13
3. 研究建立方程 ① 如图,选弦绷紧时(不振动)直线为 x 轴
u
F
1
2
T2
T1
s
0
A
x x x
B
x
14
② 弦离开平衡位置的位移记为 为表征物理量。
u ( x, t ),
③因弦的振动是机械振动,基本规律为:
F ma, 然而弦不是质点,故 F ma
对整根弦并不适用。但整根弦可以细分为许 多极小的小段,每个小段可以抽象为质点。
10
由于张力的作用,一个小段的振动必带动它
的邻段,邻段又带动它自己的邻段,这样一个 小段的振动必然传播到整个弦,这种振动的传 播现象叫作波。弦是轻质弦(其质量只有张力 的几万分之一)。跟张力相比,弦的质量完全
可以略去。
11
① 模型实际上就是:柔软轻质细弦(“没 有质量”的弦) ② 将无质量的弦紧绷,不振动时是一根直
∵
u ux tan x dx tan 1 ux x dx
20
小振动近似: x dx 与 x 两点间任一时刻横
向位移之差 u( x dx, t ) u( x, t ) 与 dx 相比是一 个小量,即
u 1 x
于是①、②化简为:
位移矢量分别为 D1 , D2 , 则
D1n s D2n
s
D E u
45
⑵ 自然边界条件
某些情况下,出于物理上的合理性等原因, 要求解为单值、有限,就提出自然边界条件, 这些条件通常都不是要研究的问题直接给出, 而是根据解的特性要求自然加上去,故称为自 然边界条件,如:
33
3、边界条件
求解方程时还需考虑边界状况(周边“环 境”)(边界状况将通过逐点影响所讨论的 整个区域),称物理过程边界状况的表达式 为边界条件,或称为边值条件。
边界条件在数学上分为三类:
34
第一类边界条件(Dirichlet边界条件):直 接规定所研究的物理量在边界上的数值
u f u ( x, y, z, t )
1
19
cos 2 1
22
2!
24
4!
1
sin 1 1
13
3!
15
5!
1 tan 1
sin 2 2
23
3!
25
5!
2 tan 2
ds (dx)2 (du )2 1 (ux ) 2 dx
线,取为 x 轴。 ③ 将弦上个点的横向位移记为 u u( x, t )
12
④ 已知:线密度 ( x, t ) (t ), 重量不计,
张力 T ( x, t ) 沿切线方向,不随x变化,弦中
各点的张力相等(小振动下T 与t 也无关). ⑤ 研究方法:连续介质,微积分思想, 任意性。
28
从物理角度看:泛定方程仅表示一般性(共
性),要为物体的运动个性化附加条件。 从数学角度看:微分方程解的任意性也需附 加条件。通解中含任意函数(解不能唯一确定 )。通过附加条件确定任意函数(常数),从 而确定解。这些附加条件就是前面所谈的问题 的“历史”与“环境”,即初始条件和边界条 件,统称为定解条件。
弦的各处都振动起来。振动如何传播呢?
8
1. 物理模型 实际问题:设有一根细长而柔软的弦,紧绷 于A,B两点之间,在平衡位置附近产生振幅极 为微小的横振动(以某种方式激发,在同一平
面内,弦上各点的振动方向相互平行,且与波
的传播方向(弦的长度方向)垂直),求弦上 各点的运动规律。