Интегральное уравнение Абеля: ядро, обращение, решение

Интегральное уравнение Абеля - это уравнение первого рода, в котором неизвестная функция стоит под интегралом, а ядро в точке $t = x$ обращается в бесконечность. Такая особенность кажется препятствием, но именно она делает уравнение разрешимым в замкнутом виде: для него есть явная формула обращения. Исторически Абель пришёл к нему, решая механическую задачу о кривой, по которой тело скатывается за одно и то же время из любой точки. Ниже разберём, как устроено слабо особое ядро, как выводится формула обращения, что даёт обобщение с произвольным показателем и как всё это связано с дробным интегрированием. Если нужно обратить конкретное уравнение, соберите условие в форме ниже.

Что такое уравнение Абеля

Классическое интегральное уравнение Абеля записывают так:

$$f(x) = \int_0^x \frac{u(t)}{\sqrt{x - t}}\, dt.$$

Здесь $f(x)$ - заданная функция (правая часть), $u(t)$ - искомая функция под интегралом, а множитель $1/\sqrt{x - t}$ - это ядро. Верхний предел переменный, как у уравнения Вольтерра, поэтому формально уравнение Абеля - частный случай Вольтерра первого рода. Отличие в том, что ядро не ограничено: при $t \to x$ знаменатель стремится к нулю и подынтегральное выражение уходит в бесконечность.

Несмотря на это, интеграл сходится: особенность типа $(x - t)^{-1/2}$ интегрируема, потому что показатель меньше единицы. Такое ядро называют слабо особым (weakly singular). Именно слабая, а не сильная особенность позволяет построить корректную теорию решения.

Рассчитать для своей задачи

Задача о таутохроне - откуда оно взялось

Абель вывел уравнение в 1823 году из чисто механической задачи. Представим гладкую кривую в вертикальной плоскости, по которой без трения скатывается тяжёлая бусина. Время спуска от точки старта на высоте $y$ до самой нижней точки зависит от формы кривой. Вопрос: какой должна быть кривая, чтобы это время было одним и тем же для любой стартовой высоты?

Кривая с таким свойством называется таутохроной (от греческого «равное время»). Пусть $T(y)$ - заданное время спуска как функция стартовой высоты, а $s(y)$ - длина дуги. Закон сохранения энергии даёт скорость $v = \sqrt{2g(y - \eta)}$ на промежуточной высоте $\eta$, и суммирование времени по дуге приводит ровно к уравнению Абеля:

$$T(y) = \int_0^y \frac{\varphi(\eta)}{\sqrt{2g\,(y - \eta)}}\, d\eta,$$

где $\varphi(\eta) = ds/d\eta$ - искомая зависимость длины дуги от высоты. Если потребовать $T(y) = \text{const}$, обращение уравнения даёт циклоиду. Так задача о кривой равного времени свелась к обращению интегрального оператора со слабо особым ядром.

Формула обращения Абеля

Главный результат - явная формула обращения: уравнение Абеля разрешимо в квадратурах. Для классического показателя $1/2$ решение записывается так:

$$u(x) = \frac{1}{\pi}\, \frac{d}{dx} \int_0^x \frac{f(t)}{\sqrt{x - t}}\, dt.$$

Обратите внимание на симметрию: чтобы обратить интеграл с ядром $1/\sqrt{x - t}$, нужно применить интеграл с тем же ядром к правой части, а затем продифференцировать и поделить на $\pi$. Это не случайность - она объясняется свойством свёртки степенных функций (см. раздел про дробное интегрирование ниже).

Вывод формулы опирается на тождество для бета-функции:

$$\int_t^x \frac{dx'}{\sqrt{(x' - t)(x - x')}} = \pi.$$

Подставив исходное уравнение под двойной интеграл и поменяв порядок интегрирования, получают $\int_0^x u(t)\, dt$, а дифференцирование по $x$ возвращает саму $u(x)$. Формула требует, чтобы $f(0) = 0$ и $f$ была достаточно гладкой; иначе обращение даёт обобщённую функцию.

Рассчитать для своей задачи

Обобщённое уравнение с показателем альфа

Классическое ядро - частный случай. Обобщённое уравнение Абеля имеет ядро с произвольным показателем $\alpha \in (0, 1)$:

$$f(x) = \int_0^x \frac{u(t)}{(x - t)^{\alpha}}\, dt.$$

При $\alpha = 1/2$ это исходное уравнение. Условие $0 < \alpha < 1$ обязательно: при $\alpha \ge 1$ особенность становится неинтегрируемой, а при $\alpha \le 0$ ядро вообще не особое и теория меняется. Формула обращения обобщается естественно:

$$u(x) = \frac{\sin(\pi \alpha)}{\pi}\, \frac{d}{dx} \int_0^x \frac{f(t)}{(x - t)^{1 - \alpha}}\, dt.$$

При $\alpha = 1/2$ множитель $\sin(\pi/2)/\pi = 1/\pi$, и показатель $1 - \alpha = 1/2$ - формула сводится к классической. Видно общее правило: прямое ядро имеет показатель $\alpha$, обратное - дополнительный показатель $1 - \alpha$, а нормировка задаётся синусом.

Связь с дробным интегрированием

Самый ёмкий способ понять уравнение Абеля - через дробное исчисление. Дробный интеграл Римана–Лиувилля порядка $\mu > 0$ определяется именно через ядро Абеля:

$$\bigl(I^{\mu} u\bigr)(x) = \frac{1}{\Gamma(\mu)} \int_0^x \frac{u(t)}{(x - t)^{1 - \mu}}\, dt.$$

Сравнив с обобщённым уравнением, видим: правая часть $f(x)$ - это с точностью до множителя $\Gamma(1 - \alpha)$ дробный интеграл порядка $1 - \alpha$ от функции $u$. Тогда обратить уравнение Абеля - значит применить дробную производную того же порядка. Полугрупповое свойство $I^{\mu} I^{\nu} = I^{\mu + \nu}$ объясняет, почему обратный оператор имеет дополнительный показатель: $I^{1-\alpha}$ обращается оператором $D^{1-\alpha} = \frac{d}{dx} I^{\alpha}$.

Этот взгляд переводит уравнение Абеля из разряда «исторического курьёза» в фундамент: оно задаёт само определение дробного интегрирования и порядок появления гамма-функции в нормировках. С таким же ядром строятся резольвента интегрального уравнения и теория уравнений со слабо особым ядром.

Численное решение и устойчивость

Уравнение Абеля - уравнение первого рода, поэтому задача его обращения некорректна по Адамару: производная в формуле обращения усиливает шум в данных $f(x)$. Малое возмущение правой части (например, ошибка измерения) может дать большую ошибку в $u(x)$. На практике это особенно важно в томографии и спектроскопии, где из проекций восстанавливают плотность.

Поэтому численно уравнение Абеля решают не прямым дифференцированием, а через регуляризацию: интегрируют формулу обращения по частям, чтобы перенести производную на гладкое ядро, или применяют метод произведения трапеций со специальной обработкой особенности. Сравните с уравнением Фредгольма второго рода - там задача корректна, и квадратурные схемы работают без регуляризации.

Частые ошибки

• Решать как обычное Вольтерра. Метод последовательных приближений для уравнения Абеля сходится плохо или не сходится: ядро не ограничено, оценки Вольтерра неприменимы. Нужна формула обращения.
• Забыть условие $f(0) = 0$. Без него формула обращения даёт сингулярный член; правая часть классического уравнения обязана обращаться в ноль в начальной точке.
• Перепутать показатели прямого и обратного ядра. Прямое ядро - показатель $\alpha$, обратное - $1 - \alpha$. Их легко поменять местами и получить неверную нормировку.
• Считать особенность сильной. Ядро $(x - t)^{-\alpha}$ при $0 < \alpha < 1$ интегрируемо - особенность слабая. Сильная (например $\alpha \ge 1$) дала бы расходящийся интеграл и совсем другую теорию.
• Дифференцировать численно без регуляризации. Прямое численное дифференцирование зашумлённой $f$ разрушает решение - задача некорректна, нужна регуляризация.

FAQ

Чем уравнение Абеля отличается от уравнения Вольтерра? Уравнение Абеля - частный случай Вольтерра первого рода с переменным верхним пределом, но с особым отличием: его ядро $(x - t)^{-\alpha}$ не ограничено и обращается в бесконечность при $t = x$. Обычное уравнение Вольтерра предполагает непрерывное (ограниченное) ядро. Из-за слабой особенности у Абеля есть явная формула обращения, а методы для ограниченного ядра к нему напрямую не применяются.

Почему задача обращения некорректна? Уравнение Абеля - первого рода, а такие задачи обычно неустойчивы: оператор интегрирования сглаживает, поэтому обратный оператор содержит дифференцирование, которое усиливает шум. Малое изменение правой части $f(x)$ может сильно изменить решение $u(x)$. На практике это лечат регуляризацией.

Где уравнение Абеля применяется кроме механики? В абелевой томографии и спектроскопии плазмы: по измеренным интегральным проекциям (например, яркости вдоль лучей) восстанавливают радиальное распределение плотности или излучательной способности - это обращение уравнения Абеля. Также оно лежит в основе дробного интегрирования и встречается в задачах с памятью.

Коротко

Интегральное уравнение Абеля - это уравнение Вольтерра первого рода со слабо особым ядром $(x - t)^{-\alpha}$, $0 < \alpha < 1$. Оно появилось из задачи о таутохроне и разрешимо в замкнутом виде: формула обращения применяет интеграл с дополнительным показателем $1 - \alpha$, дифференцирует результат и нормирует на $\sin(\pi\alpha)/\pi$. По сути это операция дробного дифференцирования, а само уравнение задаёт определение дробного интеграла. Задача обращения некорректна, поэтому численно её решают с регуляризацией.