行波法和达朗贝尔公式

行波法与达朗贝尔公式

我们已经熟悉常微分方程的常规解法:先不考虑任何附加条件,从方程本身求出通解,通解中含有任意常数(积分常数),然后利用附加条件确定这些常数。偏微分方程能否仿照这种办法求解呢? (一)达朗贝尔公式

试研究均匀弦的横振动方程(7-1-6)、均匀杆的纵振动方程(7-1-9)、理想传输线方程(7-1-14),它们具有同一形式

∂⎫⎛∂2

-a ⎪ u =0, 22

∂ x ⎭⎝∂ t

2

2

∂⎫⎛∂∂⎫⎛∂

+a -a ⎪ ⎪ u =0.

∂ x ⎭⎝∂ t ∂ x ⎭⎝∂ t

(7-4-1)

(1)通解

方程(7-4-1)的形式提示我们作代换

x =a (ξ+η), t =ξ-η,

(7-4-2)

因为在这个代换下,

∂∂ ξ∂

=

∂ t ∂∂ ξ∂ t

+

∂ x ∂∂ ξ∂ x

=

∂∂ t

+a

∂∂ x

,

∂⎫⎛∂

= + =- -a ⎪, ∂ η∂ η∂ t ∂ η∂ x ∂ x ⎭ ⎝∂ t

∂ t ∂∂ x ∂

方程(7-4-1)就成为 修改为

(∂/∂ξ ∂η) u =0

2

。但为了以后的书写便利,把代换(7-4-2)

1⎧

x =(ξ+η), ⎪⎪2⎨

⎪t =1(ξ-η), ⎪2a ⎩

⎧ξ=x +at ,

⎩η=x -at .

在此代换下,方程(7-4-1)化为

∂u ∂ξ ∂η

2

=0,

(7-4-3)

就很容易求解了。 先对

η

积分,得

∂ u ∂ξ

=f (ξ)

(7-4-4)

其中

f

是任意函数。再对

u =

ξ

积分,就得到通解

f (ξ) d ξ+f 2(η) =f 1(ξ) +f 2(η)

=f 1(x +at ) +f 2(x -at ),

(7-4-5)

其中

f 1

f 2

都是任意函数。

式(7-4-5)就是偏微分方程(7-4-1)的通解。不同于常微分方程的情况,式中出现任意函数而不是任意常数。 通解(7-4-5)具有鲜明物理意义。以

f 2(x -at )

而论,改用以速度 a 沿

x 正方向移动的坐标轴 X ,则新旧坐标和时间之间的关系为

⎧X =x -at , ⎨

⎩T =t ,

f 2(x -at ) =f 2(X ),

与时间 T 无关。这就是说,函数的图象在动坐标系中保持不变,亦即是随着动

f (x +at )

坐标系以速度 a 沿 x 正方向移动的行波。同理,1是以速度 a 沿

x 负方向移动的行波。

这样,偏微分方程(7-4-1)描写以速度 a 向两方向传播的行波。

(2)函数

f 1

f 2

的确定

f 1

通解(7-4-5)中的函数 与

f 2

可用定解条件确定。

我们假定所研究的弦、杆、传输线是“无限长的” (此词的真正含义见§7.2(二)末),这就不存在边界条件,设初始条件是

u

t =0

=ϕ(x ), u t

t =0

=ψ(x ). (-∞

(7-4-6)

以一般解(7-4-5)代入初始条件,得

⎧f 1(x ) +f 2(x ) =ϕ(x ),

⎨' ' af (x ) -af (x ) =ψ(x ); ⎩12

⎧f 1(x ) +f 2(x ) =ϕ(x ), ⎪⎨1 x

⎪f 1(x ) -f 2(x ) =⎰ x ψ(ξ) d ξ+f 1(x 0) -f 2(x 0) .

a ⎩

由此解得

f 1(x ) =⎪⎪⎨

⎪f (x ) =2⎪⎩

1212

ϕ(x ) +ϕ(x ) -

ψ(ξ) d ξ⎰2a

x 0

1

x

+-

1212

[f [f

1

(x 0) -f 2(x 0) ], (x 0) -f 2(x 0) ].

ψ(ξ) d ξ⎰2a

x 0

1

x

1

x +at

以此代回(7-4-5)即得满足初始条件(7-4-6)的特解

u (x , t ) =

12

[ϕ(x +at ) +ϕ(x -at ) ]+

12a

x -at

ψ(ξ) d ξ.

(7-4-7)

这叫作达朗贝尔公式

作为第一个例子,设初速为零即

(x 1, x 2)

ψ(x ) =0

,而初始位移

u 0

ϕ(x )

只在区间

上不为零,于

x =(x 1+x 2) /2

处达到最大值 ,如图7-14a 所示。

x -x 1⎧

2u ⎪0x -x , x 1

21

x 2-x x 1⎪2u , ϕ(x ) =⎨0

x 2-x 1

⎪⎪⎪

⎩0. x

≤x ≤+x 22

x 1+x 2

2

≤x ≤x 2

x 1 或 x >x 2

达朗贝尔公式(7-4-7)给

u (x , t ) =12

12

ϕ(x +at ) +

ϕ(x -at )

,即初始位移

(图7-14b 最下一图的粗线所描画)分为两半(该图细

线),分别向左右两方以速

度 a 移动(图7-14b 由下而上各图的细线所描画),这两个行波的和(图7-14b 由下而上的粗线所描画)给出各个时刻的波形。 作为第二个例子,设初始位移为零即 区间

(x 1, x 2)

ϕ(x ) =0

,而且初速

ψ(x )

也只在

上不为零,

ψ(x ) =⎨

⎧ 常数 ψ0 ( x 在 (x 1, x 2) 上), ⎩ 0 ( x 不在 (x 1, x 2) 上).

达朗贝尔公式(7-4-7)给出

u (x , t ) =

⎰2a

1

x +at

-∞

ψ(ξ) d ξ-

⎰2a

1

x -at

-∞

ψ(ξ) d ξ

=ψ(x +at ) -ψ(x -at ),

这里 ψ 指的是(图7-15)

ψ(x ) =

⎰2a

1

x +at

x -at

ψ(ξ) d ξ

⎪0 (x ≤x 1), ⎪⎪1

=⎨(x -x 1) ψ0 (x 1≤x ≤x 2),

⎪2a ⎪1

(x 2-x 1) ψ0 (x 2≤x ). ⎪⎩2a

于是,作出

+ψ(x )

-

ψ(x )

两个图形,

让它们以速度 a 分别向左右两方移动(图

7-16由下而上各图的细线所描画),两者的和(图7-16由下而上各图的粗线)就描画出各个时刻的波形。

在图7-14b 中,波已“通过”的地区,振动消失而弦静止在原平衡位置;在图7-16中,波已“通过”的地区,虽然振动已消失,但偏离了原平衡位置。 (二)端点的反射

研究半无限长弦的自由振动。半无限长的弦具有一个端点。 先考察端点固定的情况,即定解问题。

u tt -a u xx =0, (0

2

(7-4-8)

⎧⎪u t =0=ϕ(x ),

⎨⎪⎩u t t =0=ψ(x ),

u

x =0

(0≤x

(7-4-9)

=0.

(7-4-10)

注意初始条件(7-4-9)里的

ϕ(x )

ψ(x )

必须宗量 x ≥0 才有意义,这是

因为在 x

(t >x /a )

,达朗贝尔公式里的

ϕ(x -at )

x -at

ϕ(ξ) d ξ

失去意义,

公式也就不能应用。

参照§5.2习题4,不妨把这根半无限长弦当作某根无限长弦的 x ≥0 的部分。按照(7-4-10),这无限长弦的振动过程中,点 x =0 必须保持不动。这是说,无限长弦的位移 和初始速度

ψ(x )

u (x , t )

应当是奇函数,因而无限长弦的初始位移

Φ(x )

都应当是奇函数,即

⎧ϕ(x ) (x ≥0), ⎧ψ(x ) (x ≥0),

Φ(x ) =⎨ψ(x ) =⎨

⎩-ϕ(-x ) (x

通常采用“延拓”一词把(7-4-11)说成“把 奇延拓到整个无界区间,分别成为

ϕ(x ) 和

ψ(x )

从半无界区间 x ≥0

Φ(x ) ψ(x )

”。现在完全可以应用达郎

贝尔公式(7-4-7)求解无限长弦的自由振动,它的 x ≥0 的部分正是我们所考

察的半无限长弦。根据(7-4-7),

u (x , t ) =

12

[Φ(x +at ) +Φ(x -at ) ]+

⎰2a

1

x +at

x -at

ψ(ξ) d ξ.

把(7-4-11)代入上式,

⎧1

[ϕ(x +at ) +ϕ(x -at ) ]+⎪⎪2

u (x , t ) =⎨

⎪1[ϕ(x +at ) -ϕ(at -x ) ]+⎪⎩2

⎰2a ⎰2a

1

1

x +at

x -at

ψ(ξ) d ξ, (t ≤ ψ(ξ) d ξ. (t ≥

x a x a

) )

x +at

x -at

(7-4-12)

为了阐明(7-4-12)的物理意义,图7-17描画了只有初始位移而没有初始速度的情况。最下一图右半边用实线描画出分别向左右两方移动的波,左半边用细线描画出其奇延拓,奇延拓的波也分别向左右两方移动。在这图中,端点还没有引起什么影响。由下而上各图按着时间顺序描画了波的传播情况,粗线为合成的波形,端点 x =0 确实保持不动。由图可见,端点的影响表现为反射波。这反射波的相位跟入射波相反,这就是所谓半波损失。

再考察半无限长杆的自由振动,杆的端点自由,这个定解问题是

u tt -a u xx =0, (0

2

(7-4-13)

⎧⎪u t =0=ϕ(x ),

⎨⎪⎩u t t =0=ψ(x ), (0≤

x

u x

x =0

=0.

(7-4-15)

同样,不妨把这根半无限长杆当作某根无限长杆的 x ≥0 的部分。按照

u

(7-4-15),这无限长杆的振动过程中,在点 x =0 的相对伸长 x 必须保持为

零。这是说,无限长杆的位移 移

Φ (x )

u (x , t )

应当是偶函数,因而无限长杆的初始位

和初始速度

ψ (x )

都应当是偶函数,即

⎧ψ(x ) (x ≥0),

ψ(x ) =⎨

⎩ψ(-x ) (x

⎧ϕ(x ) (x ≥0),

Φ(x ) =⎨

⎩ϕ(-x ) (x

(7-4-16)

这就是“把 成为

Φ (x )

ϕ (x ) 和

ψ (x )

从半无界区间 x ≥0 偶延拓到整个无界区间分别

ψ (x )

”。现在,应用达朗贝尔公式(7-4-7)求解无限长杆的自

由振动,

u (x , t ) =

12

[Φ(x +at ) +Φ(x -at ) ]+

12a

x +at

x -at

Ψ(ξ) d ξ.

把(7-4-16)代入上式,

1 x +at x ⎧1

[ϕ(x +at ) +ϕ(x -at ) ]+ ψ(ξ) d ξ, (t ≤) ⎰ x -at ⎪22a a ⎪⎪1

u (x , t ) =⎨[ϕ(x +at ) +ϕ(at -x ) ]

⎪2

1 x +at 1 x -at x ⎪

⎪ +2a ⎰ 0 ψ(ξ) d ξ+2a ⎰ 0 ψ(ξ) d ξ. (t >a ) ⎩

(7-4-17)

自由端点的影响可以仿照图7-17加以阐明,这也是一种反射波。不同的是反射波的相位跟入射波相同,没有半波损失。 (三)定解问题是一个整体

从偏微分方程(7-4-1)解出达朗贝尔公式(7-4-7)的过程,与读者所熟悉的常微分方程的求解过程是完全类似的。

但是,很可惜,绝大多数偏微分方程很难求出通解;即使已求得通解,用定

解条件确定其中待定函数往往更加困难。

在本章的开头已指出,从物理的角度来说,问题的完整提法是在给定的定解条件下求解数学物理方程。现在我们要指出,除了达朗贝尔公式一类极少的例外,从数学的角度来说,不可能先求偏微分方程的通解然后再考虑定解条件,必须同时考虑偏微分方程和定解条件以进行求解。

这样,不管从物理上说还是从数学上说,定解问题是一个整体。 (四)定解问题的适定性

定解问题来自实际,它的解答也应回到实际中去。为此,应当要求定解问题(1)有解,(2)其解是唯一的,(3)解是稳定的,解的存在性和唯一性这两个要求明白易懂。至于第三个要求即稳定性说的是:如果定解条件的数值有细微的改变,解的数值也只作细微的改变。

为什么要求稳定呢?由于测量不可能绝对精密,来自实际的定解条件不免带有细微的误差,如果解不是稳定的,那么它就很可能与实际情况相去甚远,没有价值。

定解问题如果满足以上三个条件,就称为适定的。非适定的定解问题应当修改其提法,使其成为适定的。 以达朗贝尔解(7-4-7)为例。如果 属于具有二阶连续导数的函数类),

ϕ (x ) ∈C

1

2

(这些记号的意思是

ϕ (x )

ψ (x ) ∈C

,不难直接验证它确实满足方程

(7-4-1)和条件(7-4-6)。这是说,解是存在的。

在推得达朗贝尔公式(7-4-7)的过程中,没有对所求解的 u 作过任何假定和限制,凡满足方程(7-4-1)和条件(7-4-6)的解必可表为(7-4-7)。这是说,解是唯一的。

最后,证明达朗贝尔解(7-4-7)的稳定性。设有相差很细微的两组初始条件

⎧⎪u t =0=ϕ 1(x ),

⎨⎪⎩u t t =0=ψ1(x );

⎧⎪u t =0=ϕ 2(x ),

⎨⎪⎩u t t =0=ψ2(x );

(7-4-18) (7-4-19)

ϕ 1-ϕ 2

u 2

则相应的两解 和 相差

u 1-u 2 ≤

12

ϕ 1(x +at ) -ϕ 2(x +at ) +

12

ϕ 1(x -at ) -ϕ 2(x -at )

+

12

⎰2a

δ+

1

x +at

x -at

ψ1(ξ) -ψ2(ξ) d ξ

12a

2at δ=(1+t ) δ.

12

δ+

两解的差确是细微的。

行波法与达朗贝尔公式

我们已经熟悉常微分方程的常规解法:先不考虑任何附加条件,从方程本身求出通解,通解中含有任意常数(积分常数),然后利用附加条件确定这些常数。偏微分方程能否仿照这种办法求解呢? (一)达朗贝尔公式

试研究均匀弦的横振动方程(7-1-6)、均匀杆的纵振动方程(7-1-9)、理想传输线方程(7-1-14),它们具有同一形式

∂⎫⎛∂2

-a ⎪ u =0, 22

∂ x ⎭⎝∂ t

2

2

∂⎫⎛∂∂⎫⎛∂

+a -a ⎪ ⎪ u =0.

∂ x ⎭⎝∂ t ∂ x ⎭⎝∂ t

(7-4-1)

(1)通解

方程(7-4-1)的形式提示我们作代换

x =a (ξ+η), t =ξ-η,

(7-4-2)

因为在这个代换下,

∂∂ ξ∂

=

∂ t ∂∂ ξ∂ t

+

∂ x ∂∂ ξ∂ x

=

∂∂ t

+a

∂∂ x

,

∂⎫⎛∂

= + =- -a ⎪, ∂ η∂ η∂ t ∂ η∂ x ∂ x ⎭ ⎝∂ t

∂ t ∂∂ x ∂

方程(7-4-1)就成为 修改为

(∂/∂ξ ∂η) u =0

2

。但为了以后的书写便利,把代换(7-4-2)

1⎧

x =(ξ+η), ⎪⎪2⎨

⎪t =1(ξ-η), ⎪2a ⎩

⎧ξ=x +at ,

⎩η=x -at .

在此代换下,方程(7-4-1)化为

∂u ∂ξ ∂η

2

=0,

(7-4-3)

就很容易求解了。 先对

η

积分,得

∂ u ∂ξ

=f (ξ)

(7-4-4)

其中

f

是任意函数。再对

u =

ξ

积分,就得到通解

f (ξ) d ξ+f 2(η) =f 1(ξ) +f 2(η)

=f 1(x +at ) +f 2(x -at ),

(7-4-5)

其中

f 1

f 2

都是任意函数。

式(7-4-5)就是偏微分方程(7-4-1)的通解。不同于常微分方程的情况,式中出现任意函数而不是任意常数。 通解(7-4-5)具有鲜明物理意义。以

f 2(x -at )

而论,改用以速度 a 沿

x 正方向移动的坐标轴 X ,则新旧坐标和时间之间的关系为

⎧X =x -at , ⎨

⎩T =t ,

f 2(x -at ) =f 2(X ),

与时间 T 无关。这就是说,函数的图象在动坐标系中保持不变,亦即是随着动

f (x +at )

坐标系以速度 a 沿 x 正方向移动的行波。同理,1是以速度 a 沿

x 负方向移动的行波。

这样,偏微分方程(7-4-1)描写以速度 a 向两方向传播的行波。

(2)函数

f 1

f 2

的确定

f 1

通解(7-4-5)中的函数 与

f 2

可用定解条件确定。

我们假定所研究的弦、杆、传输线是“无限长的” (此词的真正含义见§7.2(二)末),这就不存在边界条件,设初始条件是

u

t =0

=ϕ(x ), u t

t =0

=ψ(x ). (-∞

(7-4-6)

以一般解(7-4-5)代入初始条件,得

⎧f 1(x ) +f 2(x ) =ϕ(x ),

⎨' ' af (x ) -af (x ) =ψ(x ); ⎩12

⎧f 1(x ) +f 2(x ) =ϕ(x ), ⎪⎨1 x

⎪f 1(x ) -f 2(x ) =⎰ x ψ(ξ) d ξ+f 1(x 0) -f 2(x 0) .

a ⎩

由此解得

f 1(x ) =⎪⎪⎨

⎪f (x ) =2⎪⎩

1212

ϕ(x ) +ϕ(x ) -

ψ(ξ) d ξ⎰2a

x 0

1

x

+-

1212

[f [f

1

(x 0) -f 2(x 0) ], (x 0) -f 2(x 0) ].

ψ(ξ) d ξ⎰2a

x 0

1

x

1

x +at

以此代回(7-4-5)即得满足初始条件(7-4-6)的特解

u (x , t ) =

12

[ϕ(x +at ) +ϕ(x -at ) ]+

12a

x -at

ψ(ξ) d ξ.

(7-4-7)

这叫作达朗贝尔公式

作为第一个例子,设初速为零即

(x 1, x 2)

ψ(x ) =0

,而初始位移

u 0

ϕ(x )

只在区间

上不为零,于

x =(x 1+x 2) /2

处达到最大值 ,如图7-14a 所示。

x -x 1⎧

2u ⎪0x -x , x 1

21

x 2-x x 1⎪2u , ϕ(x ) =⎨0

x 2-x 1

⎪⎪⎪

⎩0. x

≤x ≤+x 22

x 1+x 2

2

≤x ≤x 2

x 1 或 x >x 2

达朗贝尔公式(7-4-7)给

u (x , t ) =12

12

ϕ(x +at ) +

ϕ(x -at )

,即初始位移

(图7-14b 最下一图的粗线所描画)分为两半(该图细

线),分别向左右两方以速

度 a 移动(图7-14b 由下而上各图的细线所描画),这两个行波的和(图7-14b 由下而上的粗线所描画)给出各个时刻的波形。 作为第二个例子,设初始位移为零即 区间

(x 1, x 2)

ϕ(x ) =0

,而且初速

ψ(x )

也只在

上不为零,

ψ(x ) =⎨

⎧ 常数 ψ0 ( x 在 (x 1, x 2) 上), ⎩ 0 ( x 不在 (x 1, x 2) 上).

达朗贝尔公式(7-4-7)给出

u (x , t ) =

⎰2a

1

x +at

-∞

ψ(ξ) d ξ-

⎰2a

1

x -at

-∞

ψ(ξ) d ξ

=ψ(x +at ) -ψ(x -at ),

这里 ψ 指的是(图7-15)

ψ(x ) =

⎰2a

1

x +at

x -at

ψ(ξ) d ξ

⎪0 (x ≤x 1), ⎪⎪1

=⎨(x -x 1) ψ0 (x 1≤x ≤x 2),

⎪2a ⎪1

(x 2-x 1) ψ0 (x 2≤x ). ⎪⎩2a

于是,作出

+ψ(x )

-

ψ(x )

两个图形,

让它们以速度 a 分别向左右两方移动(图

7-16由下而上各图的细线所描画),两者的和(图7-16由下而上各图的粗线)就描画出各个时刻的波形。

在图7-14b 中,波已“通过”的地区,振动消失而弦静止在原平衡位置;在图7-16中,波已“通过”的地区,虽然振动已消失,但偏离了原平衡位置。 (二)端点的反射

研究半无限长弦的自由振动。半无限长的弦具有一个端点。 先考察端点固定的情况,即定解问题。

u tt -a u xx =0, (0

2

(7-4-8)

⎧⎪u t =0=ϕ(x ),

⎨⎪⎩u t t =0=ψ(x ),

u

x =0

(0≤x

(7-4-9)

=0.

(7-4-10)

注意初始条件(7-4-9)里的

ϕ(x )

ψ(x )

必须宗量 x ≥0 才有意义,这是

因为在 x

(t >x /a )

,达朗贝尔公式里的

ϕ(x -at )

x -at

ϕ(ξ) d ξ

失去意义,

公式也就不能应用。

参照§5.2习题4,不妨把这根半无限长弦当作某根无限长弦的 x ≥0 的部分。按照(7-4-10),这无限长弦的振动过程中,点 x =0 必须保持不动。这是说,无限长弦的位移 和初始速度

ψ(x )

u (x , t )

应当是奇函数,因而无限长弦的初始位移

Φ(x )

都应当是奇函数,即

⎧ϕ(x ) (x ≥0), ⎧ψ(x ) (x ≥0),

Φ(x ) =⎨ψ(x ) =⎨

⎩-ϕ(-x ) (x

通常采用“延拓”一词把(7-4-11)说成“把 奇延拓到整个无界区间,分别成为

ϕ(x ) 和

ψ(x )

从半无界区间 x ≥0

Φ(x ) ψ(x )

”。现在完全可以应用达郎

贝尔公式(7-4-7)求解无限长弦的自由振动,它的 x ≥0 的部分正是我们所考

察的半无限长弦。根据(7-4-7),

u (x , t ) =

12

[Φ(x +at ) +Φ(x -at ) ]+

⎰2a

1

x +at

x -at

ψ(ξ) d ξ.

把(7-4-11)代入上式,

⎧1

[ϕ(x +at ) +ϕ(x -at ) ]+⎪⎪2

u (x , t ) =⎨

⎪1[ϕ(x +at ) -ϕ(at -x ) ]+⎪⎩2

⎰2a ⎰2a

1

1

x +at

x -at

ψ(ξ) d ξ, (t ≤ ψ(ξ) d ξ. (t ≥

x a x a

) )

x +at

x -at

(7-4-12)

为了阐明(7-4-12)的物理意义,图7-17描画了只有初始位移而没有初始速度的情况。最下一图右半边用实线描画出分别向左右两方移动的波,左半边用细线描画出其奇延拓,奇延拓的波也分别向左右两方移动。在这图中,端点还没有引起什么影响。由下而上各图按着时间顺序描画了波的传播情况,粗线为合成的波形,端点 x =0 确实保持不动。由图可见,端点的影响表现为反射波。这反射波的相位跟入射波相反,这就是所谓半波损失。

再考察半无限长杆的自由振动,杆的端点自由,这个定解问题是

u tt -a u xx =0, (0

2

(7-4-13)

⎧⎪u t =0=ϕ(x ),

⎨⎪⎩u t t =0=ψ(x ), (0≤

x

u x

x =0

=0.

(7-4-15)

同样,不妨把这根半无限长杆当作某根无限长杆的 x ≥0 的部分。按照

u

(7-4-15),这无限长杆的振动过程中,在点 x =0 的相对伸长 x 必须保持为

零。这是说,无限长杆的位移 移

Φ (x )

u (x , t )

应当是偶函数,因而无限长杆的初始位

和初始速度

ψ (x )

都应当是偶函数,即

⎧ψ(x ) (x ≥0),

ψ(x ) =⎨

⎩ψ(-x ) (x

⎧ϕ(x ) (x ≥0),

Φ(x ) =⎨

⎩ϕ(-x ) (x

(7-4-16)

这就是“把 成为

Φ (x )

ϕ (x ) 和

ψ (x )

从半无界区间 x ≥0 偶延拓到整个无界区间分别

ψ (x )

”。现在,应用达朗贝尔公式(7-4-7)求解无限长杆的自

由振动,

u (x , t ) =

12

[Φ(x +at ) +Φ(x -at ) ]+

12a

x +at

x -at

Ψ(ξ) d ξ.

把(7-4-16)代入上式,

1 x +at x ⎧1

[ϕ(x +at ) +ϕ(x -at ) ]+ ψ(ξ) d ξ, (t ≤) ⎰ x -at ⎪22a a ⎪⎪1

u (x , t ) =⎨[ϕ(x +at ) +ϕ(at -x ) ]

⎪2

1 x +at 1 x -at x ⎪

⎪ +2a ⎰ 0 ψ(ξ) d ξ+2a ⎰ 0 ψ(ξ) d ξ. (t >a ) ⎩

(7-4-17)

自由端点的影响可以仿照图7-17加以阐明,这也是一种反射波。不同的是反射波的相位跟入射波相同,没有半波损失。 (三)定解问题是一个整体

从偏微分方程(7-4-1)解出达朗贝尔公式(7-4-7)的过程,与读者所熟悉的常微分方程的求解过程是完全类似的。

但是,很可惜,绝大多数偏微分方程很难求出通解;即使已求得通解,用定

解条件确定其中待定函数往往更加困难。

在本章的开头已指出,从物理的角度来说,问题的完整提法是在给定的定解条件下求解数学物理方程。现在我们要指出,除了达朗贝尔公式一类极少的例外,从数学的角度来说,不可能先求偏微分方程的通解然后再考虑定解条件,必须同时考虑偏微分方程和定解条件以进行求解。

这样,不管从物理上说还是从数学上说,定解问题是一个整体。 (四)定解问题的适定性

定解问题来自实际,它的解答也应回到实际中去。为此,应当要求定解问题(1)有解,(2)其解是唯一的,(3)解是稳定的,解的存在性和唯一性这两个要求明白易懂。至于第三个要求即稳定性说的是:如果定解条件的数值有细微的改变,解的数值也只作细微的改变。

为什么要求稳定呢?由于测量不可能绝对精密,来自实际的定解条件不免带有细微的误差,如果解不是稳定的,那么它就很可能与实际情况相去甚远,没有价值。

定解问题如果满足以上三个条件,就称为适定的。非适定的定解问题应当修改其提法,使其成为适定的。 以达朗贝尔解(7-4-7)为例。如果 属于具有二阶连续导数的函数类),

ϕ (x ) ∈C

1

2

(这些记号的意思是

ϕ (x )

ψ (x ) ∈C

,不难直接验证它确实满足方程

(7-4-1)和条件(7-4-6)。这是说,解是存在的。

在推得达朗贝尔公式(7-4-7)的过程中,没有对所求解的 u 作过任何假定和限制,凡满足方程(7-4-1)和条件(7-4-6)的解必可表为(7-4-7)。这是说,解是唯一的。

最后,证明达朗贝尔解(7-4-7)的稳定性。设有相差很细微的两组初始条件

⎧⎪u t =0=ϕ 1(x ),

⎨⎪⎩u t t =0=ψ1(x );

⎧⎪u t =0=ϕ 2(x ),

⎨⎪⎩u t t =0=ψ2(x );

(7-4-18) (7-4-19)

ϕ 1-ϕ 2

u 2

则相应的两解 和 相差

u 1-u 2 ≤

12

ϕ 1(x +at ) -ϕ 2(x +at ) +

12

ϕ 1(x -at ) -ϕ 2(x -at )

+

12

⎰2a

δ+

1

x +at

x -at

ψ1(ξ) -ψ2(ξ) d ξ

12a

2at δ=(1+t ) δ.

12

δ+

两解的差确是细微的。


相关文章

  • 从阿贝尔变换看定积分分部积分公式
  • 从阿贝尔变换看定积分分部积分公式 刘鹏飞 数学与应用数学专业 05级基地班 指导老师 尹小玲 2006年9月 摘要:通过深入了解阿贝尔变换的几何意义, 分析它与定积分存在某种联系; 经过进一步探讨, 得到由阿贝尔变换可以推导出定积分分部积分 ...查看


  • 高等数学大纲(物理类)
  • <高等数学>教学大纲 课程名称:高等数学 适用层次.专业:理科.工科各专业 学 时:320学时 学 分:20学分 课程类型:通识教育平台课 课 程 性 质:必修课 一.课程的教学目标与任务 高等数学是理.工.管等相关专业的第一基 ...查看


  • 617 数学分析
  • 617 数学分析 三.考试形式一)试卷满分及考试时间 本试卷满分为150分,考试时间为180分钟. (二)答题方式 答题方式为闭卷.笔试.试卷由试题和答题纸组成,所有题目的答案必须写在答题纸相应的位置上.考生不得携带具有存储功能的计算器. ...查看


  • 东华大学2012研究生入学考试 数学物理方法 考试大纲
  • 东华大学2012研究生入学考试 数学物理方法 考试大纲 一.考试性质和范围 <数学物理方法>介绍的是研究古典物理问题的数学方法,目的在于为后继专业课提供必要的数学基础和工具,巩固和深化在大学数学课程中所学到的数学知识,对学生应用 ...查看


  • 马尔可夫决策过程实例讲解
  • Machine Learning 16-Reinforcement Learning 之前我们学过3个部分的内容:监督学习.学习理论.半监督学习.现在我们来学习第四部分:自增强学习. 在监督学习中,给定了训练集以及对应的标签y ,算法要做的 ...查看


  • 注册会计师考试:财管失败(经验总结)
  • 注会财管的考试,对于很多人来说是非常难通过的,有人连续考5.6年都通不过,这其中的原因除了教材内容比较多,并且难度大以外,归纳起来,还有其他的一些原因. 第一个原因是学习时分不清主次,眉毛胡子一把抓. 注会财管教材(2012年)总共20章, ...查看


  • LTE站点功率提升计算方法
  • 一.华为 1. 功率计算公式: TDL 的RS 功率计算见下面公式,PsinglesAntena 是单通道可用功率,即单通道最大功率-单通道TDS 功率.Nrb=100,PB=1. 单通道TDS 功率= MAX(TS0上各信道功率之和, 小 ...查看


  • 中国产业内贸易指数放大效应的测度研究
  • 理论新探 中国产业内贸易指数放大效应的测度 靖 飞1,2 (1. 渤海大学管理学院,辽宁锦州121013:2. 中国社会科学院农村发展研究所,北京100005) 摘要:文章通过对中国产业内贸易水平的实证分析,发现在测度农产品等初级产品的产业 ...查看


  • 高等数学教学大纲模板
  • [高职教学大纲模板] 江西工商职业技术学院 系 < 高等数学 >课程教学大纲 一.课程的性质与任务 (一)本课程的性质 (高等数学是高等职业院校经济类专科科各专业学生的一门必修的重要基础理论课,它是为培养我国社会主义现代化建设所 ...查看


热门内容