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, Morocco)
  • Классификация: math.CO (Комбинаторика)
  • Дата подачи: 12 ноября 2025 г. на arXiv
  • Ссылка на статью: https://arxiv.org/abs/2511.09817

Аннотация

В статье устанавливается замкнутое символическое разложение интегрального оператора xIx\mathrm{I}, где If(x)=0xf(t)dt\mathrm{I}f(x) = \int_0^x f(t)\,dt, что является интегральным аналогом формулы Грюнерта для оператора Эйлера xDxD. Конкретно доказано: (xI)n=k=0n1(1)ka(n1,k)xnkIn+k(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)a(n,k) — числа Бесселя (OEIS A001498). Это разложение обеспечивает унифицированную символическую схему для вычисления повторных интегралов и порождает новые тождества, связанные с классическими функциями и целочисленными последовательностями.

Исследовательский контекст и мотивация

1. Решаемая проблема

Статья направлена на установление явного символического разложения (xI)n(xI)^n, что является интегральной двойственной задачей классической формулы Грюнерта для степеней дифференциального оператора (xD)n(xD)^n.

2. Важность проблемы

  • Теоретическое значение: Разложение степеней оператора Эйлера xDxD (формула Грюнерта) играет центральную роль в комбинаторике, теории специальных функций и символических вычислениях, обладая элегантной комбинаторной структурой, характеризуемой числами Стирлинга второго рода
  • Практическое применение: Интегральный оператор xIxI естественно возникает в уравнениях Вольтерра, повторных интегралах и теории комбинаторных операторов
  • Требование симметрии: Необходимость установления унифицированной символической схемы для дифференциальных и интегральных операторов, раскрытие глубокой связи между ними

3. Ограничения существующих методов

  • Хотя символическое разложение дифференциального оператора (xD)n(xD)^n хорошо изучено, соответствующее разложение интегрального оператора (xI)n(xI)^n остаётся малоизученным
  • Исторически, благодаря формуле повторного интеграла Коши, исследователи сосредоточивались на анализе интегральных ядер Kn(x,t)K_n(x,t), а не на символических разложениях операторов
  • Отсутствие прямого разложения типа Грюнерта препятствует прямому вычислению повторных интегралов и раскрытию их комбинаторной структуры

4. Исследовательская мотивация

  • Унифицированная схема: Установление симметричной теории дифференцирования и интегрирования
  • Вычислительная эффективность: Прямое вычисление повторных интегралов без необходимости в повторном интегрировании
  • Комбинаторные связи: Раскрытие комбинаторных связей между интегральными операторами и целочисленными последовательностями

Основные вклады

  1. Главная теорема: Впервые установлено явное замкнутое разложение степеней интегрального оператора: (xI)n=k=0n1(1)ka(n1,k)xnkIn+k(xI)^n = \sum_{k=0}^{n-1} (-1)^k a(n-1,k) x^{n-k} I^{n+k} где коэффициенты — числа Бесселя a(n,k)a(n,k)
  2. Упрощение интегрального ядра: Доказано компактное представление ядра: (xI)nf(x)=0xx(x2t2)n1(2(n1))!!f(t)dt(xI)^n f(x) = \int_0^x \frac{x(x^2-t^2)^{n-1}}{(2(n-1))!!} f(t)\,dt
  3. Новое семейство тождеств: Путём применения к специальным функциям (степенные, экспоненциальные, логарифмические) выведено несколько новых комбинаторных тождеств, включая:
    • Тождества, связывающие гамма-функцию и числа Бесселя (формула 12)
    • Тождества типа Добинского (формула 14)
    • Разложения, связанные с гармоническими числами (формула 15)
  4. Теоретическая схема: Предоставлена унифицированная теоретическая основа для символических вычислений с интегральными операторами, образующая двойственность с формулой Грюнерта для дифференциальных операторов

Подробное описание методов

Определение задачи

Входные данные: Положительное целое число nn и функция f(x)f(x)
Выходные данные: Явное символическое выражение для (xI)nf(x)(xI)^n f(x)
Цель: Найти коэффициенты c(n,k)c(n,k) и показатели α(n,k),β(n,k)\alpha(n,k), \beta(n,k) такие, что: (xI)n=kc(n,k)xα(n,k)Iβ(n,k)(xI)^n = \sum_k c(n,k) x^{\alpha(n,k)} I^{\beta(n,k)}

Теоретическая схема

1. Рекурсивное определение

Оператор (xI)n(xI)^n определяется рекурсивно как: (xI)0f(x):=f(x),(xI)nf(x):=xI((tI)n1f(t))(x),n1(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)=(n+k)!2kk!(nk)!,0kna(n,k) = \frac{(n+k)!}{2^k k!(n-k)!}, \quad 0 \leq k \leq n

Они удовлетворяют рекуррентному соотношению: a(n,k)=a(n1,k)+(nk+1)a(n,k1),a(0,0)=1a(n,k) = a(n-1,k) + (n-k+1)a(n,k-1), \quad a(0,0)=1

Эти числа появляются в коэффициентах полиномов Бесселя: yn(x)=k=0n(n+k)!k!(nk)!(x2)ky_n(x) = \sum_{k=0}^n \frac{(n+k)!}{k!(n-k)!} \left(\frac{x}{2}\right)^k

Доказательство главной теоремы

Доказательство теоремы 1 использует математическую индукцию:

Базовый шаг (n=1n=1): Очевидно, что (xI)1=xI=a(0,0)xI1(xI)^1 = xI = a(0,0)xI^1

Индуктивный шаг: Предположим, что для некоторого n1n \geq 1 верно: (xI)n=k=0n1(1)ka(n1,k)xnkIn+k(xI)^n = \sum_{k=0}^{n-1} (-1)^k a(n-1,k) x^{n-k} I^{n+k}

Тогда для n+1n+1: (xI)n+1f(x)=xI((tI)nf(t))(x)=k=0n1(1)ka(n1,k)xI(tnkIn+kf(t))(x)(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(tnkIn+kf(t))(x)=j=0nk(1)j(nk)jxnk+1jIn+k+j+1f(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) где (nk)j(n-k)_j — убывающий факториал.

Переиндексирование: Положим i=k+ji = k+j, получим: (xI)n+1f(x)=i=0n(1)ia(n,i)xn+1iIn+1+if(x)(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)=k=0min(n1,i)(nk)ika(n1,k)a(n,i) = \sum_{k=0}^{\min(n-1,i)} (n-k)_{i-k} a(n-1,k)

Проверка показывает, что это в точности рекуррентное соотношение для чисел Бесселя.

Технические инновации

  1. Открытие двойственности: Выявлена ключевая роль чисел Бесселя как коэффициентов разложения интегрального оператора, образующая двойственность с ролью чисел Стирлинга в дифференциальных операторах
  2. Техника интегрирования по частям: Искусное применение интегрирования по частям для разложения I(tnkIn+kf)I(t^{n-k}I^{n+k}f) — это центральная техника доказательства
  3. Проверка комбинаторных тождеств: Разложение оператора сводится к проверке рекуррентного соотношения для чисел Бесселя, устанавливая мост между теорией операторов и комбинаторной теорией чисел
  4. Упрощение функции ядра: Из сложного двойного суммирования выведено простое представление интегрального ядра с элегантной симметричной формой

Экспериментальная установка

Методы верификации

Данная работа представляет собой чистое математическое теоретическое исследование, использующее следующие стратегии верификации:

  1. Проверка частных случаев: Вычисление формулы для малых значений nn
  2. Тестирование на функциях: Применение формулы к классическим функциям (степенные, экспоненциальные, логарифмические)
  3. Вывод тождеств: Вывод известных или новых математических тождеств из главной теоремы

Классы тестовых функций

  1. Степенные функции f(t)=tα1f(t) = t^{\alpha-1} (α>0\alpha > 0)
  2. Экспоненциальные функции f(t)=etf(t) = e^t
  3. Логарифмические функции f(t)=ln(1+t)f(t) = \ln(1+t)

Сравнительная верификация

  • Сравнение результатов с формулой повторного интеграла Коши
  • Сравнение с прямым вычислением повторных интегралов
  • Сравнение с известными тождествами специальных функций

Экспериментальные результаты

Основные результаты

1. Компактная форма интегрального ядра (следствие 2)

Из главной теоремы и формулы Коши получаем: (xI)nf(x)=0xx(x2t2)n1(2(n1))!!f(t)dt(xI)^n f(x) = \int_0^x \frac{x(x^2-t^2)^{n-1}}{(2(n-1))!!} f(t)\,dt

Это даёт явное выражение для интегрального ядра: Kn(x,t)=x(x2t2)n1(2(n1))!!K_n(x,t) = \frac{x(x^2-t^2)^{n-1}}{(2(n-1))!!}

Значение: По сравнению со сложным ядром, полученным рекурсивным вычислением, эта форма чрезвычайно проста и обладает ясным геометрическим смыслом.

2. Применение к степенным функциям

Для f(t)=tα1f(t) = t^{\alpha-1} два пути вычисления дают:

Прямое интегрирование: (xI)n(tα1)(x)=xα+2n1Γ(α/2)2nΓ(n+α/2)(xI)^n(t^{\alpha-1})(x) = \frac{x^{\alpha+2n-1}\Gamma(\alpha/2)}{2^n \Gamma(n+\alpha/2)}

Символическое разложение: (xI)n(tα1)(x)=xα+2n1Γ(α)k=0n1(1)ka(n1,k)Γ(α+n+k)(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): Γ(α/2)2nΓ(α)Γ(n+α/2)=k=0n1(1)ka(n1,k)Γ(α+n+k)\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)=etf(t) = e^t получаем (формула 13): (xI)n+1(et)(x)=xn+1exyn(1x)k=0n(1)nk(2(nk)1)!!(2k)!!x2k+1(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}

где yny_n — полином Бесселя.

Тождество типа Добинского (формула 14): При вычислении в точке x=1x=1: a(n)=yn(1)=1e((xI)n+1(et)(1)+k=0n(1)nk(2(nk)1)!!(2k)!!)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. Применение к логарифмической функции

Используя известную формулу для Im(ln(1+t))I^m(\ln(1+t)) (формула 15): (xI)n(ln(1+t))(x)=[k=0n1(1)kxnk(1+x)n+k(nk1)!2kk!(n+k)]ln(1+x)k=0n1(1)kxnkAn+k(x)(nk1)!2kk!(n+k)(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)}

где Am(x)A_m(x) связана с гармоническими числами HmH_m.

Анализ конкретных примеров

Проверка для малых значений (n=2n=2): (xI)2=a(1,0)x2I2a(1,1)xI3=x2I2xI3(xI)^2 = a(1,0)x^2 I^2 - a(1,1)xI^3 = x^2 I^2 - xI^3

Для f(t)=1f(t)=1:

  • Прямое вычисление: (xI)2(1)=xI(t)(x)=xI(x22)=xx36=x46(xI)^2(1) = xI(t)(x) = xI\left(\frac{x^2}{2}\right) = x \cdot \frac{x^3}{6} = \frac{x^4}{6}
  • Вычисление по формуле: x2I2(1)xI3(1)=x2x22xx36=x42x46=x43x^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}

Примечание: Тщательная проверка требует аккуратной работы с рекурсивным определением; формулы в статье строго доказаны методом индукции.

Экспериментальные находки

  1. Чередование знаков коэффициентов: Множитель (1)k(-1)^k естественно появляется в интегральном операторе, контрастируя с положительными коэффициентами дифференциального оператора
  2. Рост показателей: Порядок интегрального оператора возрастает с итерациями (In+kI^{n+k}), в то время как дифференциальный оператор сохраняет или снижает порядок
  3. Центральная роль чисел Бесселя: Числа Бесселя не только появляются в теории специальных функций, но являются ключевыми для комбинаторной структуры интегральных операторов
  4. Универсальность: Единственная формула применима к различным классам функций, порождая тождества в разных областях математики

Связанные работы

1. Формула Грюнерта и числа Стирлинга

  • Grunert (1843): Первое установление разложения (xD)n(xD)^n
  • Boyadzhiev (2012), Knopf (2003): Современные подходы к применению чисел Стирлинга
  • Вклад данной работы: Установление интегральной двойственности, замена чисел Стирлинга на числа Бесселя

2. Теория полиномов Бесселя

  • Krall & Frink (1949): Ортогональность полиномов Бесселя
  • Grosswald (1978): Систематическое исследование полиномов Бесселя
  • Вклад данной работы: Раскрытие новой роли чисел Бесселя в теории операторов

3. Символическое интегрирование и повторные интегралы

  • Lévy (1945): Основы символических вычислений
  • Medina, Moll & Rowland (2011): Повторные первообразные степеней логарифма
  • Вклад данной работы: Предоставление унифицированной операторной схемы, упрощение вычисления повторных интегралов

4. Комбинаторные тождества

  • Lang (2000), Quaintance & Gould (2016): Обобщения чисел Стирлинга
  • Rosengren (2025): Техники доказательства двойных сумм факториалов
  • Вклад данной работы: Порождение нового семейства тождеств, связывающих числа Бесселя, гамма-функцию и гармонические числа

Заключение и обсуждение

Основные выводы

  1. Центральная теорема: Успешно установлено разложение типа Грюнерта для (xI)n(xI)^n с коэффициентами — числами Бесселя
  2. Двойственная теория: Завершена симметричная теоретическая схема для дифференциальных и интегральных операторов
  3. Прикладная ценность: Предоставлены мощные инструменты для вычисления повторных интегралов и вывода тождеств
  4. Математическое значение: Раскрыта фундаментальная роль чисел Бесселя в комбинаторной теории операторов

Ограничения

  1. Анализ сходимости: Статья не обсуждает подробно, в каких функциональных пространствах формула справедлива, отсутствует строгая характеризация условий сходимости
  2. Вычислительная сложность: Хотя дано замкнутое выражение, для больших значений nn вычисление самих чисел Бесселя может потребовать рекурсии, практическая вычислительная эффективность не оценена
  3. Направления обобщения: Рассмотрены только операторы (xI)(xI), не исследованы более общие взвешенные интегральные операторы (xαI)(x^{\alpha}I) или другие варианты
  4. Глубина приложений: Хотя приведены несколько примеров применения, недостаточно исследованы глубокие приложения в теории дифференциальных уравнений и теории специальных функций
  5. Численная верификация: Как чистое теоретическое исследование, отсутствуют масштабные численные проверки и анализ ошибок

Будущие направления

Статья не указывает явно на будущие направления, но можно предположить:

  1. Обобщение на дробные интегралы
  2. Исследование более сложных композитных операторов (xαIβ)n(x^{\alpha}I^{\beta})^n
  3. Применение к асимптотическому анализу специальных функций
  4. Исследование связей с другими комбинаторными последовательностями
  5. Разработка эффективных численных алгоритмов

Глубокая оценка

Достоинства

1. Теоретическая новизна (★★★★★)

  • Высокая оригинальность: Впервые установлено разложение типа Грюнерта для интегральных операторов, заполнен важный пробел в теории
  • Красота двойственности: Симметрия дифференциального и интегрального операторов обладает глубокой математической красотой
  • Нетривиальность: Появление чисел Бесселя как коэффициентов не является очевидным и требует глубокого понимания

2. Строгость доказательства (★★★★☆)

  • Полнота индукции: Доказательство главной теоремы логически ясно и полно
  • Комбинаторная верификация: Коэффициенты проверены через рекуррентные соотношения
  • Небольшие недостатки: Некоторые промежуточные шаги (например, конкретное разложение интегрирования по частям) изложены кратко

3. Прикладная ценность (★★★★☆)

  • Множественные приложения: Демонстрация применения к степенным, экспоненциальным и логарифмическим функциям
  • Новые тождества: Выведенные тождества имеют самостоятельную ценность
  • Большой потенциал: Предоставляет новые инструменты для теории специальных функций и комбинаторики

4. Ясность изложения (★★★★☆)

  • Разумная структура: От мотивации через теорию к приложениям — логичная иерархия
  • Стандартная нотация: Математические символы используются корректно
  • Пространство для улучшения: Некоторые процессы вывода могут быть более подробными для лучшего понимания читателем

Недостатки

1. Полнота теории

  • Отсутствие области сходимости: Не обсуждается, какие функциональные классы и условия сходимости применимы
  • Доказательство единственности: Не явно доказана единственность разложения
  • Граничные случаи: Недостаточно обсуждены случаи n=0n=0 или отрицательные значения

2. Глубина приложений

  • Ограниченное число примеров: Хотя приведены три приложения, каждое относительно поверхностно
  • Практические задачи: Не показано применение к реальным математическим или физическим проблемам
  • Численный аспект: Отсутствуют численные эксперименты для верификации теоретических результатов

3. Обзор литературы

  • Краткая историческая справка: Обзор исторических исследований интегральных операторов слишком краток
  • Недостаточное сравнение: Сравнение с альтернативными методами неполно
  • Полнота цитирования: Некоторые релевантные работы могут быть не процитированы

4. Технические детали

  • Интегрирование по частям: Ключевое разложение интегрирования по частям (стр. 4) представлено неожиданно, без вывода
  • Переиндексирование: Процесс переиндексирования от двойной суммы к одинарной требует более тщательного объяснения
  • Проверка тождеств: Проверка комбинаторного тождества (формула 10) рассмотрена поверхностно

Оценка влияния

На область исследования (★★★★☆)

  • Теоретическое дополнение: Завершение интегральной двойственности формулы Грюнерта имеет фундаментальную ценность
  • Инструментальная ценность: Предоставляет новый метод исследования повторных интегралов
  • Вдохновляющий характер: Может стимулировать дальнейшие исследования в теории операторов

Практическая ценность (★★★☆☆)

  • Символические вычисления: Может быть интегрирована в системы компьютерной алгебры (например, Mathematica)
  • Педагогическая ценность: Демонстрирует глубокую связь между дифференцированием и интегрированием
  • Ограничения: Прямое практическое применение может быть ограниченным

Воспроизводимость (★★★★★)

  • Полнота доказательства: Основные результаты имеют полные доказательства
  • Явность формул: Все формулы явно представлены
  • Простота проверки: Малые случаи могут быть проверены с помощью систем символических вычислений

Применимые сценарии

1. Теоретические исследования

  • Комбинаторика: Исследование новых свойств и тождеств чисел Бесселя
  • Специальные функции: Анализ интегральных преобразований, связанных с полиномами Бесселя
  • Теория операторов: Исследование символических свойств интегральных операторов

2. Символические вычисления

  • Системы компьютерной алгебры: Реализация эффективного вычисления повторных интегралов
  • Автоматический вывод: Автоматическое получение тождеств, связанных с повторными интегралами

3. Образование

  • Высшая математика: Продвинутые темы математического анализа
  • Курсы комбинаторики: Важный пример наряду с числами Стирлинга

4. Смежные области

  • Интегральные уравнения: Символическое решение уравнений Вольтерра
  • Теория вероятностей: Вычисление моментов некоторых случайных процессов
  • Физические приложения: Решение физических задач, связанных с повторными интегралами

Особенно примечательные аспекты

1. Уникальность авторского контекста

Автор Abdelhay Benmoussa — учитель начальной школы в горном районе Марокко, не имеющий формального университетского математического образования, полностью независимо завершивший это исследование. Это демонстрирует:

  • Доступность математических исследований
  • Движущую силу чистого математического интереса
  • Возможности математического общения в эпоху интернета

2. Математическая красота

  • Двойственность чисел Бесселя и чисел Стирлинга
  • Элегантная закономерность чередования знаков
  • Простая и симметричная форма интегрального ядра

3. Связность

Объединение, казалось бы, различных математических объектов (интегральные операторы, полиномы Бесселя, гамма-функция, гармонические числа) в единую схему

Избранные ссылки

  1. Grunert (1843): Исходная формула Грюнерта, фундаментальная работа по разложению дифференциальных операторов
  2. Grosswald (1978): Систематический трактат по полиномам Бесселя, теоретическая основа данной работы
  3. Medina, Moll & Rowland (2011): Повторные первообразные степеней логарифма, основа для применения к логарифмическим функциям в данной работе
  4. Rosengren (2025): Доказательство двойных сумм факториалов на MathOverflow, используется для верификации формулы (13)

Резюме

Это высококачественная работа по чистой математике, успешно установившая разложение типа Грюнерта для интегрального оператора (xI)n(xI)^n, заполнившая важный пробел в теории операторов. Основные вклады работы:

  1. Теоретическая полнота: Установление симметричной теории дифференциальных и интегральных операторов
  2. Идентификация коэффициентов: Открытие новой комбинаторной интерпретации чисел Бесселя
  3. Широкое применение: Вывод нескольких новых математических тождеств

Хотя есть пространство для улучшения в анализе сходимости, глубине приложений и численной верификации, ценность работы как пионерского исследования неоспорима. Особенно примечательно, что автор, работающий учителем в горной школе, смог завершить столь высокого уровня исследование, что полностью демонстрирует универсальную привлекательность математической красоты и открытость современного академического общения.

Рекомендуемый рейтинг: ★★★★☆ (Настоятельно рекомендуется исследователям в области комбинаторики, теории операторов и теории специальных функций)