当前位置:文档之家› 波动方程求解法1

波动方程求解法1


=
ϕ
( x),
∂u ( x,1) ∂y

(
x)
(x ≥ 0)
(2.10)
解: 上述方程的特征线为 xy = c1,
做自变量变换 ξ = xy,η = x ,
x y
=
c2
y
标准型为
uξη

1


=
0.
(2.11)

ω = uη . (2.12)
25

ωξ

1

ω
=
0.
(2.13)
其通解为 ω = ξθ (η).
问题2: 在区间 [x1, x2 ] 上给定初始数据, 它能决定平
令 t′ = t −τ 并记 ω (x,t′;τ ) = ω(x,t′ + τ ;τ )
则定解问题(3)可化为以下的形式:
⎪⎧ ∂ 2ω
⎨ ∂t 2
=
a2
∂ 2ω
∂x 2
(t′ > 0),
⎪⎩ω (x,0;τ ) = 0,ωt (x,0;τ ) = f (x,τ ).
由达朗贝尔公式知其解为
∫ ω (x,t′;τ ) = 1 x+at′ f (ξ ,τ )dξ 2a x−at′
将此表达式代入(2.12)
可得
uη = ξθ (η).
积分可得
u(ξ ,η) = ξθ1(η) +θ1(ξ ).
进一步可得到通解为
Hale Waihona Puke u(x, y) =xyθ
2
(
x y
)
+
θ1
(
xy).
(2.14)
26
下面利用初始条件可得
θ1(x) + xθ2 (x) = ϕ(x), (2.15)
xθ1′(x) +
x 2
的解(其中 τ 为参数),
∫ 则 u(x,t) = t ω(x,t;τ )dτ 0
(4)
就是定解问题(2)的解.
利用含参变量积分的求导公式:
⎛ β (x)
⎞′ β (x)
∫ ∫ ⎜⎜⎝ α (x)
f (x, y)dy ⎟⎟⎠
=
α (x)
f ′(x, y)dy +
f (x, y)β ′(x) −
f (x, y)α ′(x).
θ
2
(
x)

x
3/


2
(
x)
=
ψ
(
x).
(2.16)
对(2.1θ5)1′(关x于) +x微2分1可x θ得2 (x) + xθ2′(x) = ϕ′(x).
用-x乘以上式再与(2.16)相加
由此可得
θ 2′ ( x)
=
1 2x
ϕ ′( x)

1 2

x
3

( x).
∫ ∫ θ2
(x)
=
1 2
x ϕ '(τ )dτ − 1
区间[x-at, x+at]称为点(x,t)的依赖区域。它的求法是过点 (x,t)作斜率为1/a,-1/a的两条直线与x轴交截而得的区间。
t
t
(x, t)
x‐at
x+at
x1
x
0
x2
x
0
x1
x2 x
a)点(x,t)的依赖区间
b)区间[x1,x2]的决定区域
图1.2 依赖区间、决定区域及影响区域
c)区间[x1,x2]的影响区域
⎪⎩u2t (x, 0) = 0,
(t > 0,
− ∞ < x < ∞),
定解问题解u为: u = u1 + u2.
齐次化原理
本小节先考虑如何求解非齐次的自由振动问题:
⎪⎪⎪⎨⎧u∂∂(t2xu2 ,0=)
a2
∂2u ∂x 2
= 0,
+
f
( x, t )
(t > 0,
− ∞ < x < ∞),
(2)
换回原变量, 则得
∫ ω(x,t;τ ) = 1 x+a(t−τ ) f (ξ ,τ )dξ
2a x−a(t −τ )
再代入(4)式就得到定解问题(2)的解为
∫ ∫ u(x,t) = 1 t x+a(t−τ ) f (ξ ,τ )dξdτ 2a 0 x−a(t −τ )
⎪⎧ ∂ 2u ⎨ ∂t 2
得到
c1+c2=0.
21
最后我们可得
u(x,t) = 1 (ϕ(x − at) + ϕ(x + at))
2
∫ + 1 x+atψ (ξ )dξ. 2a x−at
这个公式称为达朗贝尔公式。
22
举例,求解弦振动方程的柯西问题
⎧ ⎪⎪ ⎨ ⎪⎪⎩u
(
∂2u ∂t 2

∂2u ∂x2
=
0
x, 0) = x, ∂u(x, 0) ∂t
⎝0
⎠0
t
∫=
a
ω2 xx
(
x,
t;τ
)dτ
+
f
( x, t )
0
= a2uxx (x, t) + f (x, t)
0
u(x, 0) = ∫ω(x,t;τ )dτ = 0, 0
∫⎛ t

ut (x, 0) = ⎜ ωt (x,t;τ )dτ + ω(x,t;t) ⎟
= 0.
⎝0
⎠t =0
∫ 则 u(x,t) = t ω(x,t;τ )dτ 0 就是定解问题(2)的解.
xy 2
x/ y ϕ '(τ )dτ − xy τ
xy 2
x/ xy
y
ψ τ
(3τ/2)dτ
.
28
再将 θ2 (x)表达式代入(2.15) 可得
∫ ∫ θ1 (x) = ϕ(x) −
x 2
x ϕ '(τ )dτ + τ x0
x 2
x ψ (τ )dτ − c τ x0 3/2
x.
于是Cauchy问题(2.10)的解为
⎪⎩ut (x,0) = 0,
齐次化原理, 其基本思想是把非齐次方程的求解问题 转化为相应的齐次方程的情况来处理, 从而可以直接 利用前面有关齐次方程的结果来求解.
若 ω(x,t;τ ) 是定解问题
⎧ ∂ 2ω
⎪ ⎨
∂t
2
=
a2
∂2ω
∂x2
(t > τ ),
(3)
⎪⎩ω(x,τ ;τ ) = 0,ωt (x,τ ;τ ) = f (x,τ ),
所以由 ω (x,t;t) = 0,
⎛t
⎞′ t
ut (x,t) = ⎜ ∫ω(x,t;τ )dτ ⎟ = ∫ωt (x,t;τ )dτ + ω (x,t;t)
⎝0
⎠0
t
= ∫ωt (x,t;τ )dτ 0
⎛t
⎞′ t
∫ ∫ utt (x,t) = ⎜ ωt (x,t;τ )dτ ⎟ = ωtt (x,t;τ )dτ + ωt (x,t;t)
∫ ∫ u(x, y) = ϕ(xy) +
xy 2
x/ y ϕ '(τ )dτ − xy τ
xy 2
x/ xy
y
ψ τ
(3τ/2)dτ
.
29
依赖区域、决定区域和影响区域
问题1:Cauchy 问题(2.1)(2.2)的解u(x,t)在点(x,t)的值与x 轴上的哪些点的初始条件有关?
从达朗贝尔公式可看出,解u(x,t)在点(x,t)的值依赖于x轴 上区间[x-at, x+at]上的初始条件,而与其他点上的初始 条件无关。
x, 0)
− =
∂2u ∂x2
sin
= 2x x, ∂u(x,
∂t
(t > 0) = x
0,
−∞ < x < ∞) (−∞ < x < ∞)
由刚才的公式可得其解为:
∫ u(x,t) = 1 (sin(x − t) + sin(x + t)) + 1 x+t α dα
2
2 x−t
∫ ∫ + 1 t x+(t−τ ) 2ξ dξ dτ = sin x cos t + xt + xt2.
2a
G'(x) = 1 (aϕ'(x) +ψ (x)).
2a
(1) (2) (3)
20
再将以上两式关于x 积分一次就得到
∫ F (x) = 1 ϕ(x) − 1
2
2a

0

)dξ
+
c1,
∫ G(x) = 1 ϕ(x) + 1
2
2a

0

)dξ
+
c2.
其中c1与c2是常数。由 F (x) + G(x) = ϕ(x),
x∈R
x∈R
时,则与之相对应的Cauchy问题的解 u(x, t), u (x, t),满足
sup | u(x,t)-u (x,t) |< ε.
x∈R ,0≤t ≤T
24
例1,求解柯西问题
⎧ ⎪⎪
x2 ∂2u − y2 ∂2u = 0
∂x2
∂y 2
(x > 0, y > 1)
⎨ ⎪⎪⎩u ( x,1)
相关主题