Skip to content

定理Theorem

指数二分性与全线有界响应

Exponential dichotomy · Green operator for bounded ODE solutions · 指数二分的Green核

用随时间移动的稳定与不稳定投影给双向指数界,构造全实线唯一有界受迫解,并提供残差、积分截断及周期系统的判定接口。

向前增长的方向不适合从遥远过去累积外力,但反向传播时它会衰减。指数二分把两类方向分别处理:稳定部分从过去积分,不稳定部分从未来积分。这样即使零解不吸引,也能找到全实线上唯一有界的受迫响应。

形式陈述 ​

设 A:R→Rn×n 连续,U(t,s) 是线性系统 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)‖≤Ke−α(t−s)(t≥s),(3)‖U(t,s)(I−P(s))‖≤Ke−α(s−t)(t≤s).

P(t) 是投影,但不要求正交。式(1)表示稳定方向随动力学一致搬运;式(2)控制它向前衰减;式(3)控制其补方向向后衰减。允许 P=I 或 P=0,并非强制两部分都非零。

有界响应定理。 若上述条件成立,g:R→Rn 连续有界,则 x′=A(t)x+g(t) 恰有一个全线有界解:

(4)x∗(t)=∫−∞tU(t,s)P(s)g(s)ds−∫t∞U(t,s)(I−P(s))g(s)ds.

两项绝对收敛,且

(5)‖x∗‖∞≤2Kα‖g‖∞.

若某部分为零,相应积分和对应的一个 K/α 可直接省去。这里的有界是同时控制过去和未来,不只是某个初值的向前有界。

直觉

两个积分的符号从哪里来 ​

稳定部分在时间 s 注入后,可向后来的时刻 t 传播并指数减弱,因此累积整个过去。补方向若这样做会放大过去误差;改用终端条件,沿未来时刻 s 向回传播到 t。常数变易中将终点移回当前时刻,正好带来式(4)的负号。

用式(2)–(3)积分 Ke−α|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)‖≤Ke−α(t−s)‖z‖∞.

令 s→−∞,得到 P(t)z(t)=0。再对 s>t 使用式(3),令 s→+∞,得到 (I−P(t))z(t)=0。于是 z=0。任意两个有界受迫解的差是这种齐次解,故唯一。

例子与边界

旋转中的增长、衰减方向 ​

令

J=(0−110),R(θ)=(cos⁡θ−sin⁡θsin⁡θcos⁡θ),D=diag(1,−3),A(t)=3J+R(3t)DR(3t)T.

演化矩阵和稳定投影为

U(t,s)=R(3t)diag(et−s,e−3(t−s))R(3s)T,P(t)=R(3t)diag(0,1)R(3t)T.

旋转保持欧氏范数,所以稳定方向向前以速率3衰减,另一方向向后以速率1衰减。式(1)–(3)可取 K=1,α=1。固定某时刻的坐标轴不能代替这里随时间旋转的投影。

对 g(t)=R(3t)(1,1)T,在旋转坐标中式(4)给

y1=−∫t∞et−sds=−1,y2=∫−∞te−3(t−s)ds=13,

所以 x∗(t)=R(3t)(−1,1/3)T。任何多出非零第一齐次分量的解都向前指数偏离它。因此“全线唯一有界”不意味着“所有邻近初值向前趋近”。

中心方向和缓慢衰减都不够 ​

标量 x′=0 有无穷多个非零全线有界齐次解,故不可能满足上述二分条件。加常外力1后,x′=1 没有全线有界解;这正是缺少衰减方向带来的障碍。

标量 x′=−x/(1+|t|) 的每个解在向前无穷时趋零,但从零时刻出发有 U(t,0)=1/(1+t),不能被固定的 Ke−αt 控制。标量连续投影只能恒为0或1;选择1违反向前指数界,选择0则因向后增长违反式(3)。普通渐近衰减因此不足以推出一致指数二分。

推论与应用

从残差和有限积分取得明确误差 ​

若 z 是全线有界C1近似函数,残差 r=z′−Az−g 连续有界,则 z−x∗ 是外力为 r 的唯一有界解,故

(6)‖z−x∗‖∞≤2Kα‖r‖∞.

有界性是这里的真实条件:一个很大的齐次增长解可以有零残差,却不接近有界响应。

若只将式(4)分别积到 t−L 与 t+L,得到 xL,则

(7)‖x∗(t)−xL(t)‖≤2Kαe−αL‖g‖∞.

这是丢掉两条无限尾的误差,不包含有限区间数值求积或浮点舍入。式(6)还要求全线残差上界,有限采样不能自动替代这个上界。

周期系数的判定与周期响应 ​

对连续T周期的 A,Floquet理论的单周期矩阵与幂估计给出一个等价判据:全实线指数二分当且仅当没有乘子的模等于1。

确切地说,将 M=U(T,0) 的广义特征空间按乘子模小于1、大于1分成两部分,取投影 P0。两部分的 Mk 或 M−k 分别指数衰减。令 X(t)=U(t,0)、P(t)=X(t)P0X(t)−1;由于 P0M=MP0,P 为T周期。把任意时差分成整数个周期和余下不足一个周期的两端,再使用紧区间上的演化矩阵界,就得到统一的 K,α。

反之,若存在模为1的乘子,其特征向量生成一个非零全线有界齐次解:整数圈只乘单位模数,两圈之间由一个紧周期内的矩阵控制。对复特征向量取一个非零实部或虚部即可。这与上面的唯一性证明矛盾。

若 g 也T周期,平移 x∗(t+T) 仍是同一方程的全线有界解,唯一性迫使它等于 x∗(t)。这给周期受迫问题一个更强的全线解释;但有限周期唯一性只要求乘子不等于1,仍允许其他单位圆乘子,条件比指数二分弱。

参考资料
关系图谱7 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系