由2.1节可知,时标上的线性矩阵方程(1.4.1)的解具有多项式表达,且其系数完全由与 ω 相伴的动态解决定。本节将与 ω 相伴的动态解记为 f ω ( t ),并讨论 f ω ( t )的显式表达。
本节假设 ω ( t )具有如下形式
其中
a
j
为互不相同的复数,
m
j
为非负整数且满足
n
+1。
由于
,我们首先考虑微分算子(
D
-
al
)
k
。
如(1.4.6)所示,通过简单演算易得
满足
我们知道,函数集合
时是线性无关的
[92]
。事实上,对于一般的时标
,也是线性无关的。
引理2.2.1
设
为一时标,
上的基本广义指数多项式为
则函数集合
线性无关。
证明:
首先证明对于固定的
,集合
线性无关。
任取 n 个不同的函数
不妨设 k 1 < k 2 <…< k n 。
若
用
作用(2.2.4)得
故 λ n =0。以此类推可得
即
线性无关。
其次证明函数集合
线性无关。
任取两两不同的
n
(
n
>1)个数
,设
用算子
作用(2.2.7)有
故有
λ
n
=0。同理可以证明
λ
1
=
λ
2
=…
λ
n
-1
=0,所以函数族
线性无关。
下面用数学归纳法证明函数集合
线性无关。
任取两个函数
。
设
若
a
=
b
,不妨设
k
<
m
.用
作用(2.2.10)式两端有
,故
λ
2
=0。再由(2.2.10)式得
λ
1
=0。
若
a≠b
,用
作用(2.2.10)式两端有
故 λ 2 =0。再由(2.2.10)式得 λ 1 =0。
所以 g a,k , g b,m 线性无关。
假设集合 G 中任意 n 个函数线性无关。
任取
n
+1个两两互不相同的函数
。
设
若
,用算子
作用(2.2.11)式两端有
从而 λ n +1 =0,由归纳假设得 λ 1 = λ 2 =…= λ n =0。
若
中至少有两个数相等,不妨设
a
1
=
a
2
=…
a
m
,且
k
1
<
k
2
,…<
k
m
,其中1<
m
≤
n
+1。
用算子
作用(2.2.11)式两端有
从而
λ
m
=0,再由归纳假设有
。
综上,函数集合
线性无关。
我们把由集合
生成的复向量空间记为
ε
,称
ε
中的元素为广义指数多项式。
易见,函数组
的任意线性组合都是(
D
-
aI
)
j
核空间中的元素。因此若
ω
如(2.2.1)所示,则集合
生成的向量空间包含在
ω
(
D
)的核空间内。
在集合 S 上定义可交换的卷积运算*。
若 a ≠ b ,定义
其中
若 a = b ,定义
定义线性运算为
对
。
则卷积运算可线性扩张到空间 ε 上。
定义 ε 上的线性泛函Φ为
事实上,即为
处的值。
引理2.2.2
设Φ为
ε
上的线性泛函(2.2.17)。那么对任意
f
,有
证明: 只需验证
再由线性条件即得(2.2.18)。
(2.2.19)式的左端为
(2.2.19)式的右端为
两端 g b,j ,0≤ j ≤ m 的系数均相等,故(2.2.19)式成立。
与 ω 相伴的动态解 f ω 可以写成指数多项式的卷积形式,即
下面给出定理1.4.2的证明。
定理1.4.2的证明 :由卷积的定义易知由 S 生成的 ε 的子空间关于卷积运算*是封闭的,因此 ω ( D ) f ω =0。只需验证 f ω 满足(1.4.3)的初值条件即可。由(1.2.14)及(2.2.15),我们有
由(2.2.16)知
因此,任意两个指数多项式的卷积在 t =0点的值均为零。从而有
下面我们来确定 f ω 的初值。
若
ω
只有一个
n
+1重根
a
,则
。令0≤
k
≤
n
,由于
,有
从而有
若 ω 至少有两个不同的根,设 n >1,0≤ k ≤ n 。
f ω 以写成
其中
由引理2.2.2有
,故
按(2.2.28)依次迭代得
往证
对 n 用数学归纳法证明(2.2.30)。
当
n
=1时,
。若
k
=0,则
若 k =1,则
由(2.2.31)及(2.2.32)得
。
假设(2.2.30)在 n 时成立。对 n +1,设
其中
由(2.2.2)得
。故
对0≤ k ≤ n +1,由(2.2.29)式及(2.2.30)式有
综上,(2.2.30)式成立,故
是与ω相伴的动态解。
我们还可以考虑线性常矩阵▽动力学方程 [7]
其中
A
为
m
阶常矩阵,
为矩阵值函数。
称方程(2.2.34)的解为▽矩阵指数函数,记为
。
与Δ动力学方程(1.4.1)的讨论类似,我们可以对应得到▽动力学方程(2.2.34)的求解方法。这里只列出一些主要的内容。
设ω同(1.4.2)。称初值问题
的解为与ω相伴的▽动态解,其中 D 代表▽导算子。
定义时标
上的▽多项式序列为
其中
是方程
的解。
定义时标
上的基本▽指数多项式为
定义卷积运算
为
其中
我们可以利用与ω相伴的▽动态解将矩阵方程(2.2.34)的解表示成 A 的多项式。
定理2.2.1
设
为ω的Horner多项式;
是与ω相伴的▽动态解。
则▽动力学方程(2.2.34)的解(即▽矩阵指数函数)为
证明: 与定理1.4.1类似,略。
由(2.2.41)和(2.2.42)可见,▽矩阵指数函数的系数完全由与ω相伴的▽动态解决定。
我们同样可以得到与ω相伴▽动态解的显式表达式。
定理2.2.2 设ω同(1.4.2),令
则
是与ω相伴的▽动态解。
证明: 证明与定理1.4.2类似,略。
由定理2.2.1和定理2.2.2即可得到▽动态方程(2.2.34)的解的显式表达。