2025-11-28T17:31:18.745539

The Integral Analogue of Grunert's Formula

Benmoussa

We establish a closed symbolic expansion for the operator $x\mathrm{I}$, where $\mathrm{I}f(x) = \int_0^x f(t)\,dt$, analogous to the Grunert formula for the Euler operator $xD$. Specifically, we prove that \[ (x\mathrm{I})^n = \sum_{k=0}^{n-1} (-1)^k\, a(n-1,k)\, x^{\,n-k} \mathrm{I}^{\,n+k}, \] where $a(n,k)$ are the Bessel numbers (OEIS \seqnum{A001498}). This expansion provides a unified symbolic framework for computing iterated integrals and yields new identities involving classical functions and integer sequences.

academic

The Integral Analogue of Grunert's Formula

基本信息

论文ID: 2511.09817
标题: The Integral Analogue of Grunert's Formula
作者: Abdelhay Benmoussa (Anssis School Complex, Morocco)
分类: math.CO (Combinatorics)
发表时间: 2025年11月12日提交至arXiv
论文链接: https://arxiv.org/abs/2511.09817

摘要

本文建立了积分算子 $x\mathrm{I}$ 的闭符号展开式，其中 $\mathrm{I}f(x) = \int_0^x f(t)\,dt$ ，这是欧拉算子 $xD$ 的Grunert公式的积分类比。具体证明了： $(x\mathrm{I})^n = \sum_{k=0}^{n-1} (-1)^k\, a(n-1,k)\, x^{n-k} \mathrm{I}^{n+k}$ 其中 $a(n,k)$ 是Bessel数（OEIS A001498）。该展开式为计算迭代积分提供了统一的符号框架，并产生了涉及经典函数和整数序列的新恒等式。

研究背景与动机

1. 要解决的问题

本文旨在建立积分算子 $(xI)^n$ 的显式符号展开，这是微分算子 $(xD)^n$ 经典Grunert公式的积分对偶问题。

2. 问题的重要性

理论意义：欧拉算子 $xD$ 的幂次展开（Grunert公式）在组合数学、特殊函数和符号计算中扮演核心角色，具有优美的组合结构，由第二类Stirling数表征
实际应用：积分算子 $xI$ 在Volterra方程、迭代积分和组合算子理论中自然出现
对称性需求：建立微分与积分算子的统一符号框架，揭示两者间的深层联系

3. 现有方法的局限性

虽然微分算子 $(xD)^n$ 的符号展开已被充分研究，但积分算子 $(xI)^n$ 的对应展开却鲜有关注
历史上，由于Cauchy重复积分公式，研究者更关注积分核 $K_n(x,t)$ 的分析，而非符号算子展开
缺乏直接的Grunert型展开式阻碍了迭代积分的直接计算和组合结构的揭示

4. 研究动机

统一框架：建立微分和积分的对称理论
计算效率：无需重复积分即可直接计算迭代积分
组合联系：揭示积分算子与整数序列间的组合连接

核心贡献

主要定理：首次建立了积分算子幂次的显式闭形式展开： $(xI)^n = \sum_{k=0}^{n-1} (-1)^k a(n-1,k) x^{n-k} I^{n+k}$ 其中系数为Bessel数 $a(n,k)$
积分核简化：证明了紧凑的核表示： $(xI)^n f(x) = \int_0^x \frac{x(x^2-t^2)^{n-1}}{(2(n-1))!!} f(t)\,dt$
新恒等式族：通过应用于特殊函数（幂函数、指数函数、对数函数），导出了多个新的组合恒等式，包括：
- 涉及Gamma函数和Bessel数的恒等式（式12）
- Dobinski型恒等式（式14）
- 调和数相关的展开式（式15）
理论框架：提供了积分算子符号计算的统一理论基础，与微分算子的Grunert公式形成对偶

方法详解

任务定义

输入：正整数 $n$ 和函数 $f(x)$
输出： $(xI)^n f(x)$ 的显式符号表达式
目标：找到系数 $c(n,k)$ 和指数 $\alpha(n,k), \beta(n,k)$ 使得： $(xI)^n = \sum_k c(n,k) x^{\alpha(n,k)} I^{\beta(n,k)}$

理论框架

1. 递归定义

算子 $(xI)^n$ 递归定义为： $(xI)^0 f(x) := f(x), \quad (xI)^n f(x) := x \cdot I((tI)^{n-1}f(t))(x), \quad n \geq 1$

2. Bessel数的性质

Bessel数定义为： $a(n,k) = \frac{(n+k)!}{2^k k!(n-k)!}, \quad 0 \leq k \leq n$

满足递推关系： $a(n,k) = a(n-1,k) + (n-k+1)a(n,k-1), \quad a(0,0)=1$

这些数出现在Bessel多项式的系数中： $y_n(x) = \sum_{k=0}^n \frac{(n+k)!}{k!(n-k)!} \left(\frac{x}{2}\right)^k$

主要定理的证明

定理1的证明采用数学归纳法：

基础步骤 ( $n=1$ )：显然 $(xI)^1 = xI = a(0,0)xI^1$

归纳步骤：假设对某个 $n \geq 1$ 成立： $(xI)^n = \sum_{k=0}^{n-1} (-1)^k a(n-1,k) x^{n-k} I^{n+k}$

则对 $n+1$ ： $(xI)^{n+1}f(x) = xI((tI)^n f(t))(x) = \sum_{k=0}^{n-1} (-1)^k a(n-1,k) xI(t^{n-k}I^{n+k}f(t))(x)$

关键步骤：对每一项应用分部积分： $I(t^{n-k}I^{n+k}f(t))(x) = \sum_{j=0}^{n-k} (-1)^j (n-k)_j x^{n-k+1-j} I^{n+k+j+1}f(x)$ 其中 $(n-k)_j$ 是下降阶乘。

重新索引：令 $i = k+j$ ，得到： $(xI)^{n+1}f(x) = \sum_{i=0}^n (-1)^i a(n,i) x^{n+1-i} I^{n+1+i}f(x)$

其中系数满足组合恒等式： $a(n,i) = \sum_{k=0}^{\min(n-1,i)} (n-k)_{i-k} a(n-1,k)$

通过验证，这正是Bessel数的递推关系。

技术创新点

对偶性发现：识别出Bessel数作为积分算子展开系数的关键角色，与Stirling数在微分算子中的角色形成对偶
分部积分技巧：巧妙运用分部积分展开 $I(t^{n-k}I^{n+k}f)$ ，这是证明的核心技术
组合恒等式验证：将算子展开归结为Bessel数递推关系的验证，建立了算子理论与组合数论的桥梁
核函数简化：从复杂的双重求和推导出简洁的积分核表示，具有优美的对称形式

实验设置

验证方法

本文作为纯数学理论研究，采用以下验证策略：

特殊情况验证：通过计算小 $n$ 值的情况验证公式
函数应用测试：将公式应用于经典函数（幂函数、指数函数、对数函数）
恒等式推导：从主定理导出已知或新的数学恒等式

测试函数类

幂函数 $f(t) = t^{\alpha-1}$ ( $\alpha > 0$ )
指数函数 $f(t) = e^t$
对数函数 $f(t) = \ln(1+t)$

对比验证

与Cauchy重复积分公式的结果对比
与直接计算的迭代积分对比
与已知的特殊函数恒等式对比

实验结果

主要结果

1. 积分核的紧凑形式（推论2）

从主定理和Cauchy公式，得到： $(xI)^n f(x) = \int_0^x \frac{x(x^2-t^2)^{n-1}}{(2(n-1))!!} f(t)\,dt$

这给出了积分核的显式表达： $K_n(x,t) = \frac{x(x^2-t^2)^{n-1}}{(2(n-1))!!}$

意义：相比通过递归计算得到的复杂核，此形式极其简洁且具有清晰的几何意义。

2. 幂函数应用

对 $f(t) = t^{\alpha-1}$ ，两种计算路径给出：

直接积分： $(xI)^n(t^{\alpha-1})(x) = \frac{x^{\alpha+2n-1}\Gamma(\alpha/2)}{2^n \Gamma(n+\alpha/2)}$

符号展开： $(xI)^n(t^{\alpha-1})(x) = x^{\alpha+2n-1}\Gamma(\alpha) \sum_{k=0}^{n-1} \frac{(-1)^k a(n-1,k)}{\Gamma(\alpha+n+k)}$

新恒等式（式12）： $\frac{\Gamma(\alpha/2)}{2^n \Gamma(\alpha)\Gamma(n+\alpha/2)} = \sum_{k=0}^{n-1} \frac{(-1)^k a(n-1,k)}{\Gamma(\alpha+n+k)}$

这是涉及Bessel数和Gamma函数的非平凡恒等式。

3. 指数函数应用

对 $f(t) = e^t$ ，得到（式13）： $(xI)^{n+1}(e^t)(x) = x^{n+1}e^x y_n\left(-\frac{1}{x}\right) - \sum_{k=0}^n \frac{(-1)^{n-k}(2(n-k)-1)!!}{(2k)!!} x^{2k+1}$

其中 $y_n$ 是Bessel多项式。

Dobinski型恒等式（式14）：在 $x=1$ 处求值： $a(n) = y_n(-1) = \frac{1}{e}\left((xI)^{n+1}(e^t)(1) + \sum_{k=0}^n \frac{(-1)^{n-k}(2(n-k)-1)!!}{(2k)!!}\right)$

这给出了OEIS序列A000806的新表示。

4. 对数函数应用

利用已知的 $I^m(\ln(1+t))$ 公式（式15）： $(xI)^n(\ln(1+t))(x) = \left[\sum_{k=0}^{n-1} \frac{(-1)^k x^{n-k}(1+x)^{n+k}}{(n-k-1)!2^k k!(n+k)}\right]\ln(1+x) - \sum_{k=0}^{n-1} \frac{(-1)^k x^{n-k}A_{n+k}(x)}{(n-k-1)!2^k k!(n+k)}$

其中 $A_m(x)$ 涉及调和数 $H_m$ 。

案例分析

小规模验证（ $n=2$ ）： $(xI)^2 = a(1,0)x^2 I^2 - a(1,1)xI^3 = x^2 I^2 - xI^3$

对 $f(t)=1$ ：

直接计算： $(xI)^2(1) = xI(t)(x) = xI\left(\frac{x^2}{2}\right) = x \cdot \frac{x^3}{6} = \frac{x^4}{6}$
公式计算： $x^2 I^2(1) - xI^3(1) = x^2 \cdot \frac{x^2}{2} - x \cdot \frac{x^3}{6} = \frac{x^4}{2} - \frac{x^4}{6} = \frac{x^4}{3}$