常延迟微分方程与分步构造

形式陈述

今天的变化率可以依赖一段时间以前的状态。固定延迟τ>0，考虑实向量函数x及方程

x'(t)=f(t,x(t),x(t-\tau)),\quad t\ge0, \qquad x(s)=\phi(s),\quad -\tau\le s\le0. \tag{1}

这里的初始输入φ是整个闭区间上的连续函数，不是只有一个x(0)。解在[−τ,T]上连续，在正时间具有连续导数；零点只要求右导数，负时间的历史不必可微。它逐点满足式(1)，并完整保留所给历史。将初始时刻移到非零值只需同时平移全部时间。

本页使用一个容易逐项检查的全局充分条件：f在[0,∞)×Rᵈ×Rᵈ上连续；对每个有限T存在常数Lₓ,L_z≥0，使

\|f(t,u,v)-f(t,\widetilde u,\widetilde v)\| \le L_x\|u-\widetilde u\|+L_z\|v-\widetilde v\|, \quad 0\le t\le T. \tag{2}

两个状态位置分别满足全局Lipschitz控制，常数允许依赖T。任取连续φ，式(1)便有唯一解到每个有限T，因此有全部正时间上的唯一解。“全局”不表示解有界或稳定，只表示有限时刻不会失去存在性。

为什么能按延迟长度逐段求解

在第一个窗口0≤t≤τ，x(t−τ)=φ(t−τ)已经知道。先解普通初值问题

u'(t)=f(t,u(t),\phi(t-\tau)),\qquad u(0)=\phi(0). \tag{3}

右端连续且对u为Lₓ-Lipschitz，Picard–Lindelöf定理给局部存在唯一性。到第二个窗口，延迟位置落在刚完成的第一段内；将那段已知函数代入，解同样形式的ODE。一般第k段是[kτ,(k+1)τ]与目标窗口的交，起值由前段末端给出。

还须证明每一段确实能走完。固定T，令

M=\max_{[-\tau,0]}\|\phi\|,\qquad c=\max_{0\le t\le T}\|f(t,0,0)\|, \qquad L=L_x+L_z.

连续性和紧集极值定理使这两个最大值有限。在任何已经构造的时间范围内，记V(t)=max_{−τ≤s≤t}‖x(s)‖。积分式与式(2)给

V(t)\le M+ct+L\int_0^t V(s)\,ds \le (M+cT)e^{Lt},\qquad t\le T. \tag{4}

最后一步使用Grönwall不等式。若某段在有限端点停止，当前状态和延迟状态都由式(4)有界，f也在对应紧集上有界。积分式于是给出有限末端极限，再由局部ODE定理续接，产生矛盾。重复有限次便到T；每段的唯一性依次保证全部解唯一。

这份证明依赖τ为固定正数，延迟位置总落在已完成的历史内。它没有证明含x′(t−τ)的中立型方程、含未来值的方程，或状态依赖延迟的存在理论。

直觉

一张即时照片一般不能预测下一刻的速度。第一步需要读取φ(−τ)，再往前则依次读取φ的其余部分。把当前窗口算完后，它又成为下一窗口可查询的历史；分步法维护的对象是连续曲线，不是一串没有插值规则的离散点。

左图取τ=2/3，蓝色常数历史和灰色仿射历史在零点相接却给出不同未来。右图取数值例的τ=3/5、h=1/4：时刻一需要查询2/5处的连续折线，时间位置与插值权重共同确定输入。

在时刻t可以把历史写成一个函数

x_t(\theta)=x(t+\theta),\qquad -\tau\le\theta\le0. \tag{5}

其大小取上确界范数‖x_t‖∞。自治方程的未来由这一整段决定，而非由x(t)单独决定。连续历史空间里的演化通常只向前定义；即使未来唯一，也不能据此声称历史能唯一倒推。

例如方程x′=0会在时间τ以后把所有具有相同φ(0)的初始历史送成同一常数窗口。过去曲线不同，后来的整段状态仍可能相同。这没有违反逐段向前的唯一性。

例子与边界

同一个当前值，两种不同的第一段

取标量方程x′(t)=x(t−τ)。历史φ₁(s)=1和φ₂(s)=1+s/τ都连续，并且φ₁(0)=φ₂(0)=1。但在0≤t≤τ上分别有

x_1(t)=1+t,\qquad x_2(t)=1+\frac{t^2}{2\tau}. \tag{6}

第一条的初始右导数为1，第二条为0，因为它们在−τ处的值不同。到τ时两值为1+τ和1+τ/2。只保存x(0)=1将无法区分这两个合法问题。

常数历史生成延迟多项式

继续取φ=1。每个窗口只是积分已经求出的多项式，得到

x(t)= \begin{cases} 1+t,&0\le t\le\tau,\\ 1+t+(t-\tau)^2/2,&\tau\le t\le2\tau,\\ 1+t+(t-\tau)^2/2+(t-2\tau)^3/6,&2\tau\le t\le3\tau. \end{cases} \tag{7}

一般对t≥0，可以写成有限和

x(t)=1+\sum_{k=0}^{\lfloor t/\tau\rfloor} \frac{(t-k\tau)^{k+1}}{(k+1)!}. \tag{8}

在节点恰好多出的项为零。对正时间求导，首项给1，其余项重新编号后就是x(t−τ)；第一窗口则直接使用历史1。因此式(8)不仅在网格点匹配，而是整个连续区间的解。

τ=2/3时，第三段末端给

x(2)=3+\frac{(4/3)^2}{2}+\frac{(2/3)^3}{6} =\frac{319}{81}.

接缝会把不光滑向后传

式(7)在全部正时间为C¹，但它与负时间常数历史在0处的导数不匹配：左导数0，右导数1。在τ处二阶导数由0跳到1，在2τ处三阶导数跳变。连续延迟值保证一阶右端仍连续，并没有让所有高阶导数自动平滑。

若φ本来为C¹，要让0处也成为两侧C¹接缝，至少需要

\phi'(0)=f(0,\phi(0),\phi(-\tau)).

更高阶光滑还要相应的兼容条件。数值方法若跨越这些高阶断点直接套全段高阶Taylor余项，必须重新检查其正则性假设。

式(2)是充分条件，不是必要条件。删去全局增长控制后，仍可能有局部解，但未必存在到任意终点：取f(t,u,v)=u²便回到会有限时间爆破的方程。这个例子合法地忽略延迟位置，已经足以否定“有延迟就能无限续接”。

推论与应用

历史误差怎样传播

设x和y使用同一个f，但初始历史差不超过M₀。令E(t)=max_{−τ≤s≤t}‖x(s)−y(s)‖。由式(2)与积分式，

E(t)\le M_0+L\int_0^t E(s)\,ds, \qquad E(t)\le M_0e^{Lt}. \tag{9}

特别地，x(t)和整个当前窗口x_t都受这个界控制。它给有限时间连续依赖，通常不适合预测长时间衰减；即使真实系统耗散，取绝对值后的L仍非负。保留当前阻尼与延迟增益的符号，才能使用Halanay型比较得到衰减率和扰动底线。

数值历史是方程输入的一部分

若只保存节点值，下一次t−τ未必落在节点上。连续Euler延迟方法规定同一条折线同时负责输出和延迟查询，再检查它在整个区间的微分残差。只在图上把节点连起来、计算时却使用另一种近邻取值，实际上改变了算法。

分步构造提供存在和唯一的逻辑，不承诺每段ODE有初等闭式解。实际成本取决于各段求解、历史表示与查询。上面的多项式例可精确计算，一般非线性问题还需独立数值误差预算。

参考资料

Andrei D. Polyanin、Vsevolod G. Sorokin、Alexei I. Zhurov，Delay Ordinary and Partial Differential Equations，CRC Press，2024，§1.1.2印页1–3、§1.1.5印页13–15：历史初值、分步法和延迟指数例。本文的全局充分条件、增长界及历史比较在正文另给完整证明。
Gerald Teschl，Ordinary Differential Equations and Dynamical Systems，AMS，2012，Chapter 2，作为每个已知历史窗口的普通初值问题工具；一般延迟方程不被视作有限维ODE的无条件特例。

形式陈述 ​

为什么能按延迟长度逐段求解 ​

同一个当前值，两种不同的第一段 ​

常数历史生成延迟多项式 ​

接缝会把不光滑向后传 ​

历史误差怎样传播 ​

数值历史是方程输入的一部分 ​

反向链接