• 沒有找到結果。

解对初值和参数的连续依赖性

在文檔中 常微分方程 (頁 92-98)

第五章 高阶微分方程 73

5.3 解对初值和参数的连续依赖性

微分方程的解不仅决定于微分方程本身,而且也决定于解的初值。通常,微 分方程还包含某些参数。因此,我们需要考虑微分方程的解对初值和参数的依赖 性。

例如,线性单摆方程

d2x

dt2 + a2x = 0 满足初值条件

x(t0) = x0, x(t0) = v0

的解为

x = x0cos a(t− t0) +v0

a sin a(t− t0).

显然,它对于初值t0, x0, v0和参数a是连续的,而且也是连续可微的。注意,参数a =

g/l随重力加速度g和单摆长度l而定。由于初值x0, v0和常数g, l都是由测量得到的,

而任何测量都难免存在误差。所以上述解对x0, v0和常数a的连续性有明显的物理 意义:只要初值和参数的误差足够小,相应的单摆振动只有很小的偏差。

二体问题的解为

r = p

1 + e cos(θ− θ0), (5.29) 其中常数

e = C3 µ

C4+ ( µ

C3)2 ≥ 0, p = C32 µ > 0.

当0≤ e < 1时,它是椭圆;当e = 1时,它是抛物线;当e > 1时,它是双曲线。

我们将讨论一般n阶微分方程组的初值问题 dy

dx = f (x, y, λ), y(x0) = y0 (5.40)

的解y = φ(x; x0, y0; λ)关于初值(x0, y0)和λ的连续依赖性问题,这里λ = (λ1,· · · , λm) K ⊂ Rm。事实上,初值也可以合并为参数来看。为此,作变换

t = x− x0, u = y− y0,

其中t是新的自变量,而u = u(t)是未知函数。则初值问题(5.40)变成 du

dt = f (t + x0, u + y0, λ), u(0) = 0. (5.41)

§5.3 解对初值和参数的连续依赖性 93 注意,原来的初值(x0, y0)在(5.41)的微分方程中和λ一样以参数的形式出现,而

(5.41)的初值条件是固定不变的。因此,不失一般性,我们只讨论初值问题 (Eλ) : dy

dx = f (x, y, λ), y(0) = 0 (5.42) 的解y = φ(x, λ)对参量λ的依赖性,其中λm维的参数向量。

我们论证的思路是:先证明初值问题(Eλ)的Picard序列k(x, λ)}对参数λ的连 续性;再证明φk(x, λ)是一致收敛的,而且它的极限函数y = φ(x, λ)(Eλ)的解,从 而也就证明了有关的解对参数λ的连续性。令

∥y∥ = max{|y1|, |y2|, · · · , |yn|}

Theorem 5.3.1. 定理5.1.n维向量值函数f (x, y, λ)在区域 G :|x| ≤ a, ∥y∥ ≤ b, ∥λ − λ0∥ ≤ c 上是连续的,而且对y满足Lipschitz条件

∥f(x, y1, λ)− f(x, y2, λ)∥ ≤ L∥y1− y2∥,

其中常数L≥ 0。令正数M∥f(x, y, λ)∥在区域G的一个上界,而且令 h = min{a, b

M}.

则初值问题(Eλ)存在唯一解,并且解y = φ(x, λ)在区域 D :|x| ≤ h, ∥λ − λ0∥ ≤ c 上是连续的。

证明:由于这证明与第三章中Picard定理的证明非常类似,我们只列出要点:

(1)初值问题(Eλ)等价于积分方程 y =

x

0

f (x, y, λ)dx. (5.43)

ya(x) =

x

0

fa(s, y(s), λ)ds, a = 1,· · · , n.

(2)由此可以作Picard序列 yk+1(x, λ) =

x

0

f (x, yk(x, λ), λ)dx, k = 0, 1, 2,· · · (5.44)

其中y0(x, λ) = y0= 0和(x, λ)∈ D

§5.3 解对初值和参数的连续依赖性 95

|y3a(x, λ)− y2a(x, λ)| = |

x

0

[fa(s, y2(s, λ), λ)− fa(s, y1(s, λ), λ)]ds|

≤ |

x

0

L∥y2(s)− y1(s)∥ds|

≤ ML2|

x

0

|s|2 2 ds|

= M L2|x|3 6 .

因此,Picard序列ya= yak(x, λ)D上一致收敛。则极限函数 φa(x, λ) := lim

k→∞yak(x, λ), (x, λ)∈ D

D上连续。利用fa(x, y, λ)的连续性以及Picard序列yka(x, λ)的一致收敛性,在 yk+1(x, λ) =

x

0

f (x, yk(x, λ), λ)dx, k = 0, 1, 2,· · · (5.44) 中令k→ ∞就得到

φ(x, λ) =

x

0

f (s, φ(s, λ), λ)ds.

因此y = φ(x, λ)(Eλ)的解。并且同样可证明是唯一解。定理5.1由此得证。 2

Corollary 5.3.2. 推论:设n维向量值函数f (x, y)在区域 R :|x − x0| ≤ a, ∥y − y0∥ ≤ b 上连续,而且对y满足李氏条件。则微分方程初值问题

dy

dx = f (x, y), y(x0) = η 的解y = φ(x, η)在区域

Q :|x − x0| ≤ h

2, |η − y0| ≤ b 2 上是连续的,其中

h = min(a, b M), 而正数M∥f(x, y)∥在区域R上的一个上界。

我们知道,解的存在性可以从局部延伸到大范围。同样,解对初值(或参数)

的连续性也有类似的如下结论。

Theorem 5.3.3. 定理5.2.n维向量值函数f (x, y)(x, y)空间内的某个开区域G上 是连续的,而且对y满足局部李氏条件。假设y = ξ(x)是微分方程

dy

dx = f (x, y) (5.46)

的一个解,令它的存在区间为J。现在,在区间J内任取一个有界闭区间a≤ x ≤ b, 则存在常数δ > 0,使得对任何初值(x0, y0),

a≤ x0≤ b, ∥y0− ξ(x0)∥ ≤ δ, 柯西问题

(E) : dy

dx = f (x, y), y(x0) = y0

的解y = φ(x; x0, y0)也至少在区间a≤ x ≤ b上存在,并且它在闭区间 Dδ: a≤ x ≤ b, a ≤ x0≤ b, ∥y0− ξ(x0)∥ ≤ δ 上是连续的。

证明:我们仍利用Picard序列逼近来证明这个定理。仅指出其中不同于局部情 形的地方,省略细节推导。

注意到积分曲线段

Γ ={(x, y)|y = ξ(x), a ≤ x ≤ b}

G内的一个有界闭集。因此,利用有限覆盖定理可知,存在σ > 0,使得以Γ为中 心线的闭管状邻域

Σσ : a≤ x ≤ b, ∥y − ξ(x)∥ ≤ σ 坐落在开区域G内;并且f (x, y)在Σσ内有整体李氏常数L

(x0, y0)∈ Σσ。我们可以构造(E)的Picard序列:

φk+1(x; x0, y0) = y0+

x

x0

f (x, φk(x; x0, y0))dx, (5.47) 这里不同之处是我们选取(它是y = ξ(x)的平移)

φ0(x; x0, y0) = y0+ ξ(x)− ξ(x0). (5.48) 然后,要证明

∥φk(x; x0, y0)− ξ(x)∥ < σ (5.49)

§5.3 解对初值和参数的连续依赖性 97

在文檔中 常微分方程 (頁 92-98)