ODEs的主化：解的比较

ODE系统的主化：通过右端函数比较解

在动力系统分析中，常微分方程（ODE）的精确解往往难以获得。主化定理允许我们通过一个具有已知性质的“简单”系统的解来估计“复杂”系统的行为。如果一个系统的右端函数和初始条件在分量上均大于另一个系统，则解轨迹将保持这种不等式关系。

对于向量 $x=(x^1,\ldots,x^m)$, $y=(y^1,\ldots,y^m) \in \mathbb{R}^m$，记 $x \le y$ 若对所有 $i=1,\dots,m$ 有 $x^i \le y^i$。立方体 $[x,y]:=\{z\in\mathbb{R}^m\mid x\le z\le y\}$。

考虑一个区域 $D \subset \mathbb{R}^m$，使得对于 $x,y\in D$，$x\le y$ 蕴含 $[x,y]\subset D$。函数 $f_1,f_2\in C(I\times D,\mathbb{R}^m)$ 具有连续雅可比矩阵 $\frac{\partial f_i}{\partial x}\in C(I\times D,\mathbb{R}^{m\times m})$，其中 $I=[t_0,t_1]$。

两个柯西问题：

$\dot x = f_i(t,x)$, $x(t_0)=\hat x_i$, $i=1,2$。

主要定理的表述

设 $f_1(t,x)\le f_2(t,x)$ 对所有 $(t,x)\in I\times D$ 成立，且 $\hat x_1\le\hat x_2$。存在具有非负元素的对角矩阵 $U_i(t,x)$ 和具有非负元素的矩阵 $B_i(t,x)$：

$\frac{\partial f_i}{\partial x}=U_i + B_i$, $i=1,2$。

陈述： 若 $x_i(t)$ 是柯西问题在 $I$ 上的解，则对所有 $t\in I$ 有 $x_1(t)\le x_2(t)$。

这一性质允许我们用一个性质已知的简单系统来主化一个不可解系统。

证明结构

引入同伦 $y=y(t,s)$，其中 $s\in[0,1]$：

$y_t=(1-s)f_1(t,y)+s f_2(t,y)$, $y(t_0,s)=(1-s)\hat x_1 + s\hat x_2$。

在 $s=0$ 时，得到 $x_1(t)$；在 $s=1$ 时，得到 $x_2(t)$。只需证明对于固定的 $t$，$s\mapsto y(t,s)$ 是非递减的，即 $y_s\ge 0$。

对参数的微分

对 $s$ 微分：

$$\begin{aligned} y_{st} &= f_2 - f_1 + (u + b)y_s, \\ u(t,s) &= sU_2(t, y(t,s)) + (1-s)U_1(t, y(t,s)), \\ b(t,s) &= sB_2(t, y(t,s)) + (1-s)B_1(t, y(t,s)), \\ y_s|_{t=t_0} &= \hat x_2 - \hat x_1 \ge 0.\end{aligned}$$

单调性区域

$M=A=\{t\in I\mid y_s(\xi,s)\ge0,\ \xi\in[t_0,t],\ s\in[0,1]\}$。该集合非空（$t_0\in A$）且封闭。令 $\tau=\sup A$。

假设 $\tau<t_1$。积分形式：

$$y_s(t,s)=e^{\int_\tau^t u\,d\lambda}y_s(\tau,s)+\int_\tau^t e^{\int_\xi^t u\,d\lambda}(f_2-f_1+by_s)\,d\xi=:\Phi(y_s).$$

压缩方法

根据压缩映射原理，方程 $\Phi(y_s)=y_s$ 对于小的 $t-\tau>0$ 可解。零近似 $Y_0=e^{\int_\tau^t u\,d\lambda}y_s(\tau,s)\ge0$。归纳地 $Y_{k+1}=\Phi(Y_k)\ge0$，因为算子 $\Phi$ 保持非负性。极限 $y_s\ge0$ 在 $[\tau,\tau+\varepsilon]\times[0,1]$ 上成立，与 $\tau$ 的定义矛盾。

应用：Gronwall-Bellman不等式

对于非负函数 $w(t)\le C+\int_{t_0}^t v(\xi)w(\xi)\,d\xi$，$v\ge0$。

辅助函数 $\phi(t)=C+\int_{t_0}^t v w\,d\xi$，$w\le\phi$，$\dot\phi=vw\le v\phi$。

系统 1：$\dot x=f_1(t,x):=vx-(v\phi-\dot\phi)$，$x(t_0)=C$。括号内非负。

主化：$\dot x=f_2(t,x):=vx$，$x(t_0)=C$，解为 $x(t)=Ce^{\int_{t_0}^tv\,d\xi}$。

根据定理 $w\le\phi\le Ce^{\int_{t_0}^tv\,d\xi}$。

关键点：

• 当 $f_1\le f_2$ 且 $\hat x_1\le\hat x_2$ 时，主化有效。
• 雅可比矩阵分解为对角矩阵 $U$（非负）和矩阵 $B$（非负）。
• 通过同伦和压缩方法证明。
• 与Gronwall不等式直接关联，用于积分估计。
• 应用：用简单指数函数估计非线性ODE的解。

结论

该定理为无需显式积分的动力学定性分析提供了工具。雅可比矩阵的条件确保了同伦的单调性，这对于保持不等式至关重要。

— Editorial Team