《Algorithms》：Optimal Methods for Approximating the Riemann–Liouville Fractional Integral in the Sobolev Space
编辑推荐：
研究人员针对索伯列夫空间 L2m(t,1) 中的右端黎曼-刘维尔分数阶积分，构造并分析了在 N+1 个等距节点处利用 φ,φ′,…,φ(m?1) 的最优求积公式。分析通过误差泛函的皮亚诺核（Peano kernel）展开，给出了该核的精确闭式表达；由此得到误差
研究人员针对索伯列夫空间 L2m(t,1) 中的右端黎曼-刘维尔分数阶积分，构造并分析了在 N+1 个等距节点处利用 φ,φ′,…,φ(m?1) 的最优求积公式。分析通过误差泛函的皮亚诺核（Peano kernel）展开，给出了该核的精确闭式表达；由此得到误差泛函范数、基于正交投影论证的最优系数存在唯一性，以及可计算的误差界。主要结果为结构性结论：依附于系数的函数构成网格上 m?1 次间断分段多项式空间的一组基，从而使最优性问题化为按面片解耦的 L2 投影。因此，全局最优（Sard）系数（而非仅序贯最优系数）对任意 m 均以闭式给出，计算代价为 O(m3N) 次运算，无全局线性系统，条件数与 N、h、α 及 t 无关。全局最优规则在 m 次多项式上精确，对光滑被积函数以 hm 收敛，而序贯规则在 m 次多项式上精确；研究人员量化了二者间隙。研究人员还证明了误差范数的锐渐近性：对 α>0 具 hm+min(α,1) 因子，α=1 处恰为对数因子，α<1 时损失半阶。在五个分数阶、三个求值点与八类给定索伯列夫正则性被积函数上的广泛数值研究验证了每一条理论陈述，公式与乘积、样条及高斯-雅可比（Gauss–Jacobi）求积比较，并应用于阿贝尔（Abel）积分方程。
研究背景方面，分数阶微积分因非局部性广泛用于流变学、扩散、控制与黏弹性建模，其中右端黎曼-刘维尔分数积分（right-sided Riemann–Liouville fractional integral）在逆向时间动力学中自然出现。传统求积规则在该弱奇异核下精度下降，而在索伯列夫空间（Sobolev space，记为 L2m(t,1)）中构造最优求积公式可最小化误差泛函范数，但已有工作多限于左端积分或固定局部多项式度，对任意阶 m 的右端积分系统序贯优化仍具挑战。为此，研究人员在《Algorithms》发表本文，旨在给出右端积分在 L2m(t,1) 中全局最优系数的闭式与结构化结论。
关键技术方法上，研究人员以误差泛函的皮亚诺核（Peano kernel）表示为核心框架，利用泰勒积分余项导出核的闭式；借助基引理证明系数依附函数构成网格上面片解耦的间断分段多项式基，从而将最优性化为 L2 正交投影；分别推导序贯（cascade）最优与全局 Sard 最优系数，前者逐阶固定前阶系数优化当前阶，后者通过面片独立 Gram 矩阵（在移位 Legendre 基下为单位阵）直接求解；辅以变量缩放引理、欧拉-麦克劳林（Euler–Maclaurin）余项与伯努利（Bernoulli）常数做渐近分析；数值部分采用 mpmath 40 位精度 Python 脚本，以 tanh–sinh 求积作参照，在五分数阶、三求值点、八类已知索伯列夫正则性被积函数（含解析、振荡、端点奇异等）上验证，并与乘积梯形/Simpson、Ciesielski–Grodzki 三次/五次样条乘积求积及 Gauss–Jacobi 求积比较，亦用于 Abel 方程 u(t)=∫t1(x?t)?1/2f(x)dx 反演。
研究结果各节浓缩如下。引言指出标准规则因核 (x?t)α?1 弱奇异而失效，综述乘积求积、样条求积、Jacobi 谱法与信息基复杂度理论，定位本文属 Hilbert/Sobolev 类最小最坏误差传统。问题陈述节给出求积公式（3）、误差泛函（4）（5）、极值函数（8）与格林函数表达，将 ∥R∥ 化为（9）（10）闭式。皮亚诺核框架节证定理2：误差可表为 ∫t1K(s)φ(m)(s)ds，核（12）（13）闭式，推论1给可算误差界；引理1证系数函数成 m?1 次间断分段多项式基；推论2由正交投影得存在唯一性； remark1 指出全局 Gram 面片解耦、条件数恒1。缩放引理2表明区间缩放下系数满足 Ck,m(a)=aα+mCk,m(1)。序贯最优节定理3、4给 m=1,2 系数即乘积梯形权；定理5证第 m 阶系数为面片均值之差，?m=O(hm+min(α,1))；定理6给锐渐近（17）（18），α=1 对数、α<1 半阶损失；命题2证序贯规则在 m 次多项式精确。全局最优节定理7证 Sard 系数由面片投影（19）（20）确定，推论4得 O(m3N) 无全局系统、条件数1；命题3证全局规则在 m 次多项式精确且光滑时 O(hm)；remark4 用表1量化序贯次优比（α=1 无损，m 奇/偶点态差距 hm/hm+1 对 hm+1/hm）。稳定性节命题4证函数值权 ≥0 且和为精确质量 (α)?1(1?t)α，Lebesgue 常数最小；数值扰动放大有 N 无关界。数值结果节表3–9 验 α=0.5,1,1.5 收敛阶、定理6指数、常数、点态阶、有限光滑饱和、同均匀网格下全局最优较乘积 Simpson 优9阶、较三次样条同阶但信息省3倍、Abel 方程6节点达 10?12。节点布置节命题5证内节点驻值条件（Peano 核两侧等幅），表10 显示最优/分级网格（22）对 α<1 恢复全阶 hm，均匀网格仅 hm+α。局限节指明需 φ∈Cm 于节点、α→0 首格降级可用（22）修、Gauss–Jacobi 自由节点更优、张量积可推广而非一般域。
讨论与结论翻译部分：研究人员讨论称，全局最优系数闭式、面片解耦 L2 投影结构解释了为何无全局线性系统且条件数与 N,h,α,t 无关；序贯方案虽丢失全局最优常数（表1 因子），但层次复用、便宜约3.4倍且在 α=1 恰优，故保留作自适应升阶。节点优化仅对 α<1 有界增益，均匀网格配合首格分级（22）已近最优。结论节原文总结：在 L2m(t,1) 中经皮亚诺核构造右端黎曼-刘维尔积分最优求积，系数函数成 m?1 次间断分段多项式基使投影按面片解耦，从而对任意 m 给全局 Sard 系数闭式、O(m3N) 代价、面片条件数1、光滑时 hm 收敛、序贯规则 m 次精确且次优比有界；证误差范数锐渐近含 α=1 对数过渡与 α<1 半阶损、权正值性致最小 Lebesgue 常数、五分数阶三求值点八被积函数数值全覆盖理论，并与乘积/样条/Gauss–Jacobi 及 Abel 方程比对；开放问题为 α→0 分级节点全分析与一般域上 cubature 面板解耦失效。
