Volterra第二类方程与预解核

有些反馈读取过去所有时刻，却从不读取未来。把昨天的反馈再展开，会出现前天；继续展开时，时间必须严格排队。这种顺序把一个看似会不断放大的迭代压成带阶乘分母的级数。因此，即使一次迭代不是压缩映射，整个方程仍可以唯一求解。

形式陈述

给定的是三角核，不是一个方阵

固定有限的T>0，令 $\Delta_T=\{(t,s):0\le s\le t\le T\}$ 。设 $k:\Delta_T\to\mathbb C$ 、 $f:[0,T]\to\mathbb C$ 连续，参数 $\lambda\in\mathbb C$ 。求连续函数u，使

u(t)=f(t)+\lambda\int_0^t k(t,s)u(s)\,ds. \tag{1}

积分在每个有限区间上按Riemann积分理解，连续性保证其存在。称它为第二类，是因为未知量既在积分外，也在积分内；上限为当前t，使其为Volterra型。没有要求k只依赖t−s。

记 $(Ku)(t)=\int_0^t k(t,s)u(s)ds$ 、 $M=\max_{\Delta_T}|k|$ 、 $B=\|f\|_\infty$ 。那么对每个λ，式(1)恰有一个连续解，且

u=\sum_{n=0}^{\infty}\lambda^nK^nf, \qquad \|u\|_\infty\le Be^{|\lambda|MT}. \tag{2}

这里 $K^0$ 是恒等算子。级数在整个[0,T]上一致绝对收敛；不需要 $|\lambda|MT<1$ 。

预解核把无限反馈集中起来

定义

k_1(t,s)=k(t,s),\qquad k_{n+1}(t,s)=\int_s^t k(t,r)k_n(r,s)\,dr.

则

|k_n(t,s)|\le\frac{M^n(t-s)^{n-1}}{(n-1)!},\qquad R_\lambda(t,s)=\sum_{n=1}^{\infty}\lambda^{n-1}k_n(t,s) \tag{3}

在三角域上一致收敛，并满足

u(t)=f(t)+\lambda\int_0^tR_\lambda(t,s)f(s)ds, \qquad |R_\lambda(t,s)|\le Me^{|\lambda|M(t-s)}. \tag{4}

本页把最外面的λ留在式(4)中；另一些文献把它吸收到预解核，比较公式时必须先核约定。λ=0时 $R_0=k$ ，但式(4)的积分整体被零乘掉，解仍为f。

有限近似带着什么保证

从 $u_0=f$ 开始，令 $u_{N+1}=f+\lambda Ku_N$ ，得到 $u_N=\sum_{n=0}^N\lambda^nK^nf$ 。记x=|λ|MT，则

\|u-u_N\|_\infty \le B\sum_{n=N+1}^\infty\frac{x^n}{n!} \le B e^x\frac{x^{N+1}}{(N+1)!}. \tag{5}

若N+2>x，还可直接使用便于有理复算的界

\|u-u_N\|_\infty \le B\frac{x^{N+1}}{(N+1)!}\frac1{1-x/(N+2)}. \tag{6}

它是充分上界，不能把某一N未通过这份预算解读成真实误差必定超标。

直觉

阶乘来自有序时间的体积

先证 $|K^nf(t)|\le B M^nt^n/n!$ 。n=0显然成立。若第n步成立，则

|K^{n+1}f(t)| \le M\int_0^t\frac{BM^ns^n}{n!}ds =\frac{BM^{n+1}t^{n+1}}{(n+1)!}.

等价地，n重积分只遍历 $0<s_n<\cdots<s_1<t$ ，其体积为 $t^n/n!$ 。一般有界算子的粗界 $\|K^n\|\le\|K\|^n$ 看不见这个时间顺序，因而会错过所有参数下的收敛。

由M判别法，式(2)一致绝对收敛。每一项连续，一致极限保连续使u连续。又有 $\|Kg\|_\infty\le MT\|g\|_\infty$ ，故K可传入一致极限，直接得到u=f+λKu。这同时证明存在性和公式的合法性，不只是对未知逆算子作形式展开。

若u、v都是解，d=u−v满足d=λKd。反复代入给

|d(t)|\le\|d\|_\infty\frac{(|\lambda|Mt)^n}{n!}

对每个n成立。令n趋于无穷，右端趋零，便有d=0。这个唯一性证明也不要求短时间或小参数。

两种核复合必须指向同一个顺序

对连续核在紧三角域上的绝对值积分有限，Fubini换序允许把 $K^nf$ 写成 $\int_0^t k_n(t,s)f(s)ds$ 。式(3)的界同样逐次积分归纳。级数于是给两条核恒等式

R_\lambda(t,s)=k(t,s)+\lambda\int_s^t k(t,r)R_\lambda(r,s)dr =k(t,s)+\lambda\int_s^t R_\lambda(t,r)k(r,s)dr. \tag{7}

第二条不是假定两个任意积分算子可交换，而是同一个K的幂自然结合。改成两个不同核，交换次序通常会改变结果。

式(5)最后一步可用 $(N+1+j)!\ge(N+1)!j!$ 逐项放大尾和。式(6)则从第一份遗漏项开始：后继项与当前项的比最多x/(N+2)<1，于是由几何级数包住整个尾部。必须保留这个严格小于一的检验。

例子与边界

不是压缩，也能逐项得到正确答案

取T=1、λ=1、k(t,s)=2s、f=1。此时M=2，粗范数界为2，甚至 $\|K\|=1$ ，所以K不是严格压缩。直接作用于幂函数给

K(t^{2n})=\frac{t^{2n+2}}{n+1},\qquad K^n1=\frac{t^{2n}}{n!},\qquad u(t)=e^{t^2}.

用微分再核一次：u′=2tu、u(0)=1，恰好恢复原积分式。核只看s，不能误写成 $2(t-s)$ ；后者给出的是另一个方程。

另一份卷积核k(t,s)=t−s给 $K^n1=t^{2n}/(2n)!$ ，故λ=1时u=cosh t。两份反馈都能用积分转成ODE，但其方程分别是一阶变系数和二阶常系数；“积分上限相同”并不意味着动力学相同。

一般多项式核仍可以输出完整日志

取 $k(t,s)=1+t-s$ 、λ=1、f=1、T=1。前三轮为

\begin{aligned} u_0&=1,\\ u_1&=1+t+t^2/2,\\ u_2&=1+t+t^2+ t^3/3+t^4/24. \end{aligned}

这里 $K1=t+t^2/2$ 、 $K^21=t^2/2+t^3/3+t^4/24$ 。每次提交整份多项式后，检查器可独立积分恢复下一轮，而不必相信一张自报误差表。x=2、B=1，可直接将式(6)写成精确有理数。终点任务还改变λ、f与核，检验程序没有写死这个例子。

全参数可解是有限窗、连续三角核的结论

若把上限改成固定1，方程 $u(t)=1+\int_0^1u(s)ds$ 无解：积分两侧得到c=1+c。这属于固定区间反馈，其共振需要Fredholm兼容条件，不能沿用三角体积证明。

若k含 $(t-s)^{\alpha-1}$ 、0<α<1，M已经无限，本页的阶乘估计不能直接套用。分数阶积分会用Beta恒等式代替三角域的常数核体积，随后得到Gamma分母。无穷时间上的统一有界性同样需要新条件；有限T上每次都有解，不代表 $\sup_{t\ge0}|u(t)|$ 有限。

推论与应用

残差要在整段时间上认证

给任意连续候选v，令 $r=v-f-\lambda Kv$ 。真误差e=v−u满足e=r+λKe，所以

|e(t)|\le |r(t)|+|\lambda|M\int_0^t|e(s)|ds.

若已认证 $\|r\|_\infty\le\eta$ ，Grönwall不等式给

\|v-u\|_\infty\le\eta e^{|\lambda|MT}. \tag{8}

有限网格残差小不自动给η；多项式可用全区间极值或保守系数和包围，数值求积还须把积分误差加入残差。对Picard候选有精确恒等式 $r_N=-\lambda^{N+1}K^{N+1}f$ ，因此可以同时交先验尾界与后验残差界，而不要求两种界数值相同。

从中途继续时，过去并没有消失

对a<t，将原方程拆成

u(t)=\underbrace{f(t)+\lambda\int_0^a k(t,s)u(s)ds}_{f_a(t)} +\lambda\int_a^t k(t,s)u(s)ds.

历史以整份函数 $f_a(t)$ 进入以后时刻，通常不能只保留u(a)。特殊的可分离核可以用有限个累计量表示历史；任意连续核没有这项保证。重启时把 $f_a$ 换回f，会解出另一个问题。

若核多项式含J个单项式，当前候选次数d，按单项式精确乘法积分需要O(J(d+1))次有理算术，另计合并同次项与输出；N轮时应累加各轮的实际次数。一般核没有有限多项式闭包，必须另选并分析求积方法。本页的级数收敛不等于已经给出任意黑箱核的有限算法。

因果积分与分数阶初值终点要求交出完整Picard日志、可核验尾预算和重启时的历史项，再把有界核迁移到弱奇异核。

参考资料

Hermann Brunner，On the History of Numerical Methods for Volterra Integral Equations，CWI Newsletter 11，1986，印页3–6，式(1.1)、(2.3)–(2.4)及两条预解核恒等式。本文补出有限窗阶乘界、明确的两种尾预算及多项式证书。
Francesco Mainardi，A Tutorial on the Basic Special Functions of Fractional Calculus，2020，§4，式(4.1)–(4.7)：弱奇异的第二类Abel方程。它不满足本页有界核假设，需用分数阶积分的独立证明。

形式陈述 ​

给定的是三角核，不是一个方阵 ​

预解核把无限反馈集中起来 ​

有限近似带着什么保证 ​

阶乘来自有序时间的体积 ​

两种核复合必须指向同一个顺序 ​

不是压缩，也能逐项得到正确答案 ​

一般多项式核仍可以输出完整日志 ​

全参数可解是有限窗、连续三角核的结论 ​

残差要在整段时间上认证 ​

从中途继续时，过去并没有消失 ​

反向链接