向前增长的方向不适合从遥远过去累积外力,但反向传播时它会衰减。指数二分把两类方向分别处理:稳定部分从过去积分,不稳定部分从未来积分。这样即使零解不吸引,也能找到全实线上唯一有界的受迫响应。
形式陈述
设 A : R → R n × n 连续,U ( t , s ) 是线性系统 理路 线性常微分方程组 Linear system of ordinary differential equations 形如 x′=A(t)x+b(t) 的向量值一阶线性方程组。 x ′ = A ( t ) x 的演化矩阵,即 U ( s , s ) = I 。称系统在全实线具有一致指数二分,若存在连续矩阵 P ( t ) 、常数 K ≥ 1 , α > 0 ,使
(1) P ( t ) 2 = P ( t ) , U ( t , s ) P ( s ) = P ( t ) U ( t , s ) , 以及
(2) ‖ U ( t , s ) P ( s ) ‖ ≤ K e − α ( t − s ) ( t ≥ s ) , (3) ‖ U ( t , s ) ( I − P ( s ) ) ‖ ≤ K e − α ( s − t ) ( t ≤ s ) . P ( t ) 是投影,但不要求正交。式(1)表示稳定方向随动力学一致搬运;式(2)控制它向前衰减;式(3)控制其补方向向后衰减。允许 P = I 或 P = 0 ,并非强制两部分都非零。
有界响应定理。 若上述条件成立,g : R → R n 连续有界,则 x ′ = A ( t ) x + g ( t ) 恰有一个全线有界解:
(4) x ∗ ( t ) = ∫ − ∞ t U ( t , s ) P ( s ) g ( s ) d s − ∫ t ∞ U ( t , s ) ( I − P ( s ) ) g ( s ) d s . 两项绝对收敛,且
(5) ‖ x ∗ ‖ ∞ ≤ 2 K α ‖ g ‖ ∞ . 若某部分为零,相应积分和对应的一个 K / α 可直接省去。这里的有界是同时控制过去和未来,不只是某个初值的向前有界。
直觉
两个积分的符号从哪里来
稳定部分在时间 s 注入后,可向后来的时刻 t 传播并指数减弱,因此累积整个过去。补方向若这样做会放大过去误差;改用终端条件,沿未来时刻 s 向回传播到 t 。常数变易中将终点移回当前时刻,正好带来式(4)的负号。
用式(2)–(3)积分 K e − α | t − s | ‖ g ‖ ∞ ,立即得到绝对收敛及式(5)。在任意紧的 t 区间上,A ( t ) 有界,因而同一个指数控制也允许对两端积分求导。上限为 t 的积分贡献 P ( t ) g ( t ) ,带负号的下限为 t 的积分贡献 ( I − P ( t ) ) g ( t ) ,总和为 g ( t ) ;其余项合成 A ( t ) x ∗ ( t ) 。因此式(4)确实解方程。
全线有界为什么消去所有齐次自由度
设 z 是一个全线有界齐次解。固定 t ,对 s < t ,由式(1)–(2)得
‖ P ( t ) z ( t ) ‖ = ‖ U ( t , s ) P ( s ) z ( s ) ‖ ≤ K e − α ( t − s ) ‖ z ‖ ∞ . 令 s → − ∞ ,得到 P ( t ) z ( t ) = 0 。再对 s > t 使用式(3),令 s → + ∞ ,得到 ( I − P ( t ) ) z ( t ) = 0 。于是 z = 0 。任意两个有界受迫解的差是这种齐次解,故唯一。
例子与边界
旋转中的增长、衰减方向
令
J = ( 0 − 1 1 0 ) , R ( θ ) = ( cos θ − sin θ sin θ cos θ ) , D = diag ( 1 , − 3 ) , A ( t ) = 3 J + R ( 3 t ) D R ( 3 t ) T . 演化矩阵和稳定投影为
U ( t , s ) = R ( 3 t ) diag ( e t − s , e − 3 ( t − s ) ) R ( 3 s ) T , P ( t ) = R ( 3 t ) diag ( 0 , 1 ) R ( 3 t ) T . 旋转保持欧氏范数,所以稳定方向向前以速率3衰减,另一方向向后以速率1衰减。式(1)–(3)可取 K = 1 , α = 1 。固定某时刻的坐标轴不能代替这里随时间旋转的投影。
对 g ( t ) = R ( 3 t ) ( 1 , 1 ) T ,在旋转坐标中式(4)给
y 1 = − ∫ t ∞ e t − s d s = − 1 , y 2 = ∫ − ∞ t e − 3 ( t − s ) d s = 1 3 , 所以 x ∗ ( t ) = R ( 3 t ) ( − 1 , 1 / 3 ) T 。任何多出非零第一齐次分量的解都向前指数偏离它。因此“全线唯一有界”不意味着“所有邻近初值向前趋近”。
中心方向和缓慢衰减都不够
标量 x ′ = 0 有无穷多个非零全线有界齐次解,故不可能满足上述二分条件。加常外力1后,x ′ = 1 没有全线有界解;这正是缺少衰减方向带来的障碍。
标量 x ′ = − x / ( 1 + | t | ) 的每个解在向前无穷时趋零,但从零时刻出发有 U ( t , 0 ) = 1 / ( 1 + t ) ,不能被固定的 K e − α t 控制。标量连续投影只能恒为0或1;选择1违反向前指数界,选择0则因向后增长违反式(3)。普通渐近衰减因此不足以推出一致指数二分。
推论与应用
从残差和有限积分取得明确误差
若 z 是全线有界C1近似函数,残差 r = z ′ − A z − g 连续有界,则 z − x ∗ 是外力为 r 的唯一有界解,故
(6) ‖ z − x ∗ ‖ ∞ ≤ 2 K α ‖ r ‖ ∞ . 有界性是这里的真实条件:一个很大的齐次增长解可以有零残差,却不接近有界响应。
若只将式(4)分别积到 t − L 与 t + L ,得到 x L ,则
(7) ‖ x ∗ ( t ) − x L ( t ) ‖ ≤ 2 K α e − α L ‖ g ‖ ∞ . 这是丢掉两条无限尾的误差,不包含有限区间数值求积或浮点舍入。式(6)还要求全线残差上界,有限采样不能自动替代这个上界。
周期系数的判定与周期响应
对连续T周期的 A ,Floquet理论的单周期矩阵与幂估计 理路 Floquet理论与单周期传播 Floquet theory · Floquet theorem · Floquet multipliers · 弗洛凯理论 将周期线性系统拆成周期因子和指数传播,以单周期矩阵判定稳定性,区分冻结谱、平均谱与实际乘子,并检查复对数和单位圆Jordan边界。 给出一个等价判据:全实线指数二分当且仅当没有乘子的模等于1。
确切地说,将 M = U ( T , 0 ) 的广义特征空间按乘子模小于1、大于1分成两部分,取投影 P 0 。两部分的 M k 或 M − k 分别指数衰减。令 X ( t ) = U ( t , 0 ) 、P ( t ) = X ( t ) P 0 X ( t ) − 1 ;由于 P 0 M = M P 0 ,P 为T周期。把任意时差分成整数个周期和余下不足一个周期的两端,再使用紧区间上的演化矩阵界,就得到统一的 K , α 。
反之,若存在模为1的乘子,其特征向量生成一个非零全线有界齐次解:整数圈只乘单位模数,两圈之间由一个紧周期内的矩阵控制。对复特征向量取一个非零实部或虚部即可。这与上面的唯一性证明矛盾。
若 g 也T周期,平移 x ∗ ( t + T ) 仍是同一方程的全线有界解,唯一性迫使它等于 x ∗ ( t ) 。这给周期受迫问题 理路 周期受迫线性ODE的相容与共振 Periodic boundary-value problem for linear ODEs · Periodic forcing compatibility · 周期线性边值问题 把周期响应的存在唯一性化为有限维单周期相容条件,保留共振无解与多解,给出残差转真实误差的证书,并区分周期性和吸引性。 一个更强的全线解释;但有限周期唯一性只要求乘子不等于1,仍允许其他单位圆乘子,条件比指数二分弱。
参考资料