Ядро интегрального уравнения: типы, свойства, методы

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

Что такое ядро интегрального уравнения

Линейное интегральное уравнение второго рода записывают так:

$$u(x) = f(x) + \lambda \int_a^b K(x, t)\, u(t)\, dt.$$

Здесь $u(x)$ - искомая функция, $f(x)$ - заданная правая часть, $\lambda$ - числовой параметр, а $K(x, t)$ - ядро интегрального уравнения: функция двух переменных, которая под интегралом «взвешивает» значения неизвестной функции $u(t)$ и собирает из них вклад в точке $x$.

Интуитивно ядро - это правило связи: $K(x, t)$ говорит, насколько сильно значение $u$ в точке $t$ влияет на результат в точке $x$. Если рассматривать интеграл как «непрерывную матрицу», то ядро - это её элементы, а само уравнение - непрерывный аналог системы линейных уравнений $u = f + \lambda K u$, где роль матрицы играет интегральный оператор с ядром $K$.

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

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

Ядро и пределы: Фредгольм против Вольтерра

Самое первое, что говорит ядро вместе с пределами интегрирования, - к какому классу относится уравнение.

• Если оба предела постоянны ($a$ и $b$ фиксированы), это уравнение Фредгольма. Ядро действует на всю область сразу, и задача может быть как разрешимой, так и нет - здесь работает альтернатива Фредгольма.
• Если верхний предел переменный ($x$ вместо $b$), это уравнение Вольтерра. Ядро «видит» только прошлое - значения $u$ от $a$ до текущей точки $x$. Формально это ядро Фредгольма, обращающееся в ноль при $t > x$.

Эта развилка важнее, чем кажется: для Вольтерра ряд приближений сходится при любом $\lambda$, а у Фредгольма сходимость требует малости параметра. Одно и то же по виду ядро в этих двух случаях ведёт себя совершенно по-разному.

Вырожденное ядро

Самый удобный для расчётов случай - вырожденное (сепарабельное) ядро, которое раскладывается в конечную сумму произведений функций только от $x$ и только от $t$:

$$K(x, t) = \sum_{i=1}^{n} a_i(x)\, b_i(t).$$

Чем хорошо такое ядро: интегральное уравнение с ним точно сводится к системе линейных алгебраических уравнений. Подставив разложение, выносим $a_i(x)$ за интеграл, а оставшиеся $\int_a^b b_i(t) u(t)\, dt$ обозначаем неизвестными числами $c_i$. Решение ищется в виде $u(x) = f(x) + \lambda \sum_i c_i a_i(x)$, и для коэффициентов $c_i$ получается обычная линейная система $n \times n$.

Например, ядро $K(x,t) = xt + x^2 t^2$ вырождено с $n = 2$: $a_1 = x,\, b_1 = t,\, a_2 = x^2,\, b_2 = t^2$. Любое уравнение с таким ядром решается школьной линейной алгеброй - без рядов и итераций.

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

Многие методы строятся на аппроксимации произвольного ядра вырожденным: усекают разложение в ряд (например, по полиномам), и задача становится конечномерной.

Разностное ядро и свёртка

Если ядро зависит только от разности аргументов,

$$K(x, t) = K(x - t),$$

его называют разностным, а интеграл в уравнении превращается в свёртку функций. Это резко упрощает решение: к свёртке применимо преобразование Лапласа (для Вольтерра с $a = 0$) или преобразование Фурье (для уравнений на всей оси). Преобразование переводит свёртку в произведение, и уравнение становится алгебраическим:

$$U(p) = F(p) + \lambda\, \mathcal{K}(p)\, U(p) \quad\Rightarrow\quad U(p) = \frac{F(p)}{1 - \lambda\, \mathcal{K}(p)}.$$

Остаётся выполнить обратное преобразование. Для уравнений с разностным ядром это самый короткий путь к ответу - никаких итераций и систем. Разностные ядра естественно возникают в задачах с инвариантностью к сдвигу: фильтрация сигналов, процессы с памятью, теория восстановления.

Симметричное ядро

Ядро называют симметричным, если

$$K(x, t) = K(t, x).$$

Это вещественный аналог самосопряжённого оператора, и у него богатая спектральная теория. Для симметричного непрерывного ядра справедлива теорема Гильберта-Шмидта: существует счётный набор вещественных собственных значений $\lambda_n$ и ортонормированная система собственных функций $\varphi_n(x)$, таких что

$$\varphi_n(x) = \lambda_n \int_a^b K(x, t)\, \varphi_n(t)\, dt.$$

Само ядро при этом раскладывается в ряд по собственным функциям (билинейное разложение), а решение уравнения выписывается как ряд по той же системе. Симметричность - это «хороший» случай: спектр вещественный, собственные функции ортогональны, а решение и условие разрешимости формулируются через собственные значения. Несимметричное ядро такой гарантии не даёт, и спектр может быть комплексным.

Слабо особое ядро

Не всякое ядро непрерывно. Слабо особое (полярное) ядро имеет интегрируемую особенность вида

$$K(x, t) = \frac{H(x, t)}{|x - t|^{\alpha}}, \qquad 0 < \alpha < 1,$$

где $H(x, t)$ непрерывна. В точке $t = x$ ядро обращается в бесконечность, но из-за условия $\alpha < 1$ интеграл всё равно сходится. Классический пример - ядро $\dfrac{1}{\sqrt{x - t}}$ в уравнении Абеля о таутохроне.

Для слабо особых ядер стандартные оценки непрерывного случая не работают напрямую, но теория всё равно строится: повторные ядра постепенно «сглаживают» особенность, и начиная с некоторого номера итерированное ядро становится непрерывным. Это позволяет применять обычные методы, лишь немного их доработав.

Повторные ядра и резольвента

Когда ядро не вырождено и не разностное, общий метод - построить решение через резольвенту $R(x, t; \lambda)$, которая зависит только от ядра:

$$u(x) = f(x) + \lambda \int_a^b R(x, t; \lambda)\, f(t)\, dt.$$

Резольвента собирается как ряд по повторным (итерированным) ядрам:

$$R(x, t; \lambda) = \sum_{n=0}^{\infty} \lambda^{n} K_{n+1}(x, t),$$

где первое повторное ядро $K_1 = K$, а каждое следующее получается сцеплением предыдущего с исходным:

$$K_{n+1}(x, t) = \int_a^b K(x, s)\, K_n(s, t)\, ds.$$

Этот ряд (ряд Неймана) для непрерывного ядра Фредгольма сходится при $|\lambda| < 1 / \|K\|$, а для ядра Вольтерра - при любом $\lambda$. Главная ценность резольвенты: найдя её один раз для данного ядра, можно решать уравнение с любой правой частью $f(x)$ - повторно работать с ядром не нужно.

Норма ядра и оценки сходимости

Свойства сходимости упираются в «размер» ядра. Для непрерывного ядра на квадрате $[a, b] \times [a, b]$ вводят оценку $M = \max |K(x, t)|$ и норму

$$\|K\| = \left( \int_a^b \int_a^b |K(x, t)|^2\, dx\, dt \right)^{1/2},$$

(норму Гильберта-Шмидта). От неё зависит, при каких $\lambda$ сходится ряд Неймана: для Фредгольма условие малости $|\lambda|\,\|K\| < 1$ как раз гарантирует сходимость и единственность решения. У Вольтерра переменный предел даёт для $n$-го повторного ядра оценку с факториалом в знаменателе вида $\dfrac{(M(x-a))^{n}}{n!}$, которая бьёт любую степень $\lambda$ - поэтому ряд там сходится всегда. Так одно число - норма ядра - отвечает на вопрос об устойчивости всей задачи.

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

• Считают ядром всю подынтегральную функцию. Ядро - это только множитель $K(x, t)$ при неизвестной $u(t)$, без самой $u$ и без правой части $f(x)$.
• Путают вырожденное и разностное ядро. Вырожденное - конечная сумма $\sum a_i(x) b_i(t)$ (метод: линейная система); разностное - $K(x - t)$ (метод: преобразование Лапласа или Фурье). Это разные структуры и разные методы.
• Ждут вещественного спектра у несимметричного ядра. Теорема Гильберта-Шмидта (вещественные собственные значения, ортогональные собственные функции) работает только для симметричного ядра $K(x,t) = K(t,x)$.
• Применяют оценки непрерывного ядра к слабо особому. При полярном ядре $|x-t|^{-\alpha}$ обычные оценки $\max|K|$ бессмысленны - нужна теория слабо особых уравнений с итерированными ядрами.
• Требуют малости $\lambda$ для Вольтерра. Условие $|\lambda|\,\|K\| < 1$ - это про Фредгольм; у Вольтерра ряд по повторным ядрам сходится при любом параметре.

FAQ

Чем ядро интегрального уравнения отличается от подынтегральной функции? Ядро - это только функция $K(x, t)$, стоящая множителем при неизвестной $u(t)$ под интегралом. Сама неизвестная $u(t)$ и правая часть $f(x)$ в ядро не входят. Подынтегральная функция целиком это произведение $K(x, t)\, u(t)$, а ядро - её фиксированная, заданная часть.

Как тип ядра влияет на выбор метода решения? Напрямую. Вырожденное ядро сводит уравнение к линейной системе; разностное ($K(x-t)$) решается преобразованием Лапласа или Фурье через свёртку; симметричное даёт разложение по собственным функциям; для произвольного непрерывного ядра строят резольвенту через повторные ядра. Поэтому классификация ядра - первый шаг.

Может ли интегральное уравнение не иметь решения из-за ядра? Да, для уравнений Фредгольма. Если $\lambda$ совпадает с собственным значением ядра, однородное уравнение имеет нетривиальное решение, и тогда по альтернативе Фредгольма неоднородное уравнение разрешимо лишь при условии ортогональности правой части собственным функциям. У Вольтерра такого не происходит - оно почти всегда корректно разрешимо.

Коротко

Ядро интегрального уравнения - функция двух переменных $K(x, t)$, которая под интегралом связывает значения неизвестной функции и определяет всю задачу. Постоянные пределы дают Фредгольма, переменный верхний предел - Вольтерра. По структуре ядра выбирают метод: вырожденное (сумма произведений функций от $x$ и $t$) сводит задачу к линейной системе; разностное $K(x-t)$ решается преобразованием Лапласа или Фурье; симметричное $K(x,t)=K(t,x)$ даёт вещественный спектр и разложение по собственным функциям; слабо особое полярное ядро требует отдельной теории. Для общего случая строят резольвенту через ряд повторных ядер, а норма ядра отвечает за сходимость этого ряда и устойчивость решения.