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

グルネルト公式の積分類似

基本情報

論文ID: 2511.09817
タイトル: The Integral Analogue of Grunert's Formula
著者: Abdelhay Benmoussa（モロッコ、Anssis School Complex）
分類: math.CO（組合論）
提出日時: 2025年11月12日（arXivへ）
論文リンク: https://arxiv.org/abs/2511.09817

要約

本論文は、積分作用素 $x\mathrm{I}$ （ただし $\mathrm{I}f(x) = \int_0^x f(t)\,dt$ ）の閉じた記号展開式を確立する。これはオイラー作用素 $xD$ のグルネルト公式の積分類似である。具体的には以下を証明した： $(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)$ はベッセル数（OEIS A001498）である。この展開式は反復積分の計算に統一的な記号フレームワークを提供し、古典関数と整数列に関する新しい恒等式をもたらす。

研究背景と動機

1. 解決すべき問題

本論文は、積分作用素 $(xI)^n$ の明示的な記号展開を確立することを目指している。これは微分作用素 $(xD)^n$ の古典的なグルネルト公式の積分双対問題である。

2. 問題の重要性

理論的意義：オイラー作用素 $xD$ の冪展開（グルネルト公式）は、組合論、特殊関数、記号計算において中心的な役割を果たし、第2種スターリング数によって特徴付けられる優美な組合構造を持つ
実用的応用：積分作用素 $xI$ はボルテラ方程式、反復積分、組合作用素論に自然に現れる
対称性の必要性：微分と積分作用素の統一的な記号フレームワークを確立し、両者間の深い関連性を明らかにする

3. 既存方法の限界

微分作用素 $(xD)^n$ の記号展開は十分に研究されているが、積分作用素 $(xI)^n$ の対応する展開はほとんど注目されていない
歴史的には、コーシーの反復積分公式のため、研究者はより記号作用素展開よりも積分核 $K_n(x,t)$ の解析に関心を持ってきた
グルネルト型展開式の欠如は、反復積分の直接計算と組合構造の発見を妨げている

4. 研究動機

統一フレームワーク：微分と積分の対称理論を確立する
計算効率：反復積分を直接計算でき、重複積分が不要
組合的関連性：積分作用素と整数列間の組合的つながりを明らかにする

核心的貢献

主定理：積分作用素の冪の明示的な閉形式展開を初めて確立： $(xI)^n = \sum_{k=0}^{n-1} (-1)^k a(n-1,k) x^{n-k} I^{n+k}$ ここで係数はベッセル数 $a(n,k)$ である
積分核の簡略化：コンパクトな核表現を証明： $(xI)^n f(x) = \int_0^x \frac{x(x^2-t^2)^{n-1}}{(2(n-1))!!} f(t)\,dt$
新しい恒等式族：特殊関数（べき関数、指数関数、対数関数）への適用を通じて、複数の新しい組合恒等式を導出。以下を含む：
- ガンマ関数とベッセル数を含む恒等式（式12）
- ドビンスキー型恒等式（式14）
- 調和数に関連する展開（式15）
理論的フレームワーク：積分作用素の記号計算に統一的な理論基礎を提供し、微分作用素のグルネルト公式と双対を形成

方法の詳細説明

タスク定義

入力：正整数 $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. ベッセル数の性質

ベッセル数は以下のように定義される： $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$

これらの数はベッセル多項式の係数に現れる： $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)$

検証により、これはベッセル数の再帰関係そのものである。

技術的革新点

双対性の発見：ベッセル数が積分作用素展開の係数の重要な役割を担うことを認識。これは微分作用素におけるスターリング数の役割と双対をなす
部分積分技巧： $I(t^{n-k}I^{n+k}f)$ の部分積分展開を巧妙に利用。これが証明の中核技術
組合恒等式の検証：作用素展開をベッセル数の再帰関係の検証に帰着させ、作用素論と組合数論の橋渡しを確立
核関数の簡略化：複雑な二重和から簡潔な積分核表現を導出。優美な対称形式を持つ

実験設定

検証方法

本論文は純粋数学の理論研究であり、以下の検証戦略を採用：

特殊ケースの検証：小さい $n$ 値の場合を計算して公式を検証
関数適用テスト：古典関数（べき関数、指数関数、対数関数）への公式適用
恒等式導出：主定理から既知または新しい数学恒等式を導出

テスト関数クラス

べき関数 $f(t) = t^{\alpha-1}$ ( $\alpha > 0$ )
指数関数 $f(t) = e^t$
対数関数 $f(t) = \ln(1+t)$

比較検証

コーシーの反復積分公式の結果との比較
直接計算した反復積分との比較
既知の特殊関数恒等式との比較

実験結果

主要結果

1. 積分核のコンパクト形式（系2）

主定理とコーシー公式から、以下を得る： $(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}$ に対して、2つの計算経路が以下を与える：

直接積分： $(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)}$

これはベッセル数とガンマ関数を含む非自明な恒等式である。

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$ はベッセル多項式である。

ドビンスキー型恒等式（式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}$