Задача Штурма-Лиувилля собственные значения: как искать

Задача Штурма-Лиувилля - это краевая задача для линейного дифференциального уравнения второго порядка, в которой ищут не просто решение, а те особые значения параметра, при которых нетривиальное решение вообще существует. Эти значения называют собственными, а отвечающие им решения - собственными функциями. Конструкция лежит в основе метода разделения переменных для уравнений теплопроводности, колебаний струны и Шрёдингера: ряд Фурье - это частный случай разложения по собственным функциям задачи Штурма-Лиувилля. Ниже разберём канонический вид оператора, как краевые условия порождают дискретный спектр, почему собственные функции ортогональны с весом и где чаще всего ошибаются при поиске собственных значений.

Канонический вид и оператор

Регулярная задача Штурма-Лиувилля на отрезке $[a, b]$ записывается в самосопряжённой форме:

$-\frac{d}{dx}\!\left( p(x)\,\frac{dy}{dx} \right) + q(x)\, y = \lambda\, w(x)\, y,$

где $p(x) > 0$ и весовая функция $w(x) > 0$ на всём отрезке, а $p$, $p'$, $q$, $w$ непрерывны. Левую часть удобно обозначить как действие дифференциального оператора $\mathcal{L}[y] = -(p y')' + q y$, и тогда уравнение приобретает вид задачи на собственные значения $\mathcal{L}[y] = \lambda w y$, в точности как $A\mathbf{x} = \lambda \mathbf{x}$ в линейной алгебре. Параметр $\lambda$ и есть искомое собственное значение, а функция $y$ - собственная функция. Любое линейное уравнение $a_2(x) y'' + a_1(x) y' + a_0(x) y + \lambda y = 0$ приводится к этой форме умножением на интегрирующий множитель $\mu(x) = \exp\!\big(\int (a_1/a_2)\,dx\big)/a_2$ - после чего $p = \mu a_2$, и оператор становится самосопряжённым.

Tool: постановка и поиск собственных значений

Свести своё уравнение к самосопряжённой форме, выписать характеристическое условие из краевых условий и найти первые собственные значения - рутина, в которой легко потерять знак или вес. Опишите ниже коэффициенты $p$, $q$, $w$, отрезок и тип краевых условий - и получите разбор: приведение оператора, уравнение на спектр, первые $\lambda_n$ и собственные функции.

Краевые условия и постановка задачи

Само уравнение задаёт лишь оператор; задачу замыкают граничные условия на концах. В регулярном случае берут разделённые (Робина) условия:

$\alpha_1 y(a) + \alpha_2 y'(a) = 0, \qquad \beta_1 y(b) + \beta_2 y'(b) = 0,$

где пары коэффициентов не обнуляются одновременно. Частные случаи - условия Дирихле ($y = 0$ на конце) и Неймана ($y' = 0$). Однородность условий принципиальна: правые части равны нулю, поэтому функция $y \equiv 0$ всегда формально подходит. Задача Штурма-Лиувилля ищет те $\lambda$, при которых, кроме тривиального, есть ещё и нетривиальное решение. Если концы отрезка особые (например, $p$ обращается в нуль) или условия периодические - задача называется сингулярной и ведёт себя тоньше; здесь рассматриваем регулярный случай.

Спектр: почему собственные значения дискретны

Ключевая теорема Штурма-Лиувилля утверждает, что у регулярной задачи существует бесконечная последовательность действительных собственных значений

$\lambda_1 < \lambda_2 < \lambda_3 < \cdots, \qquad \lambda_n \to +\infty,$

каждое простое (кратности один). Дискретность возникает из-за краевых условий: общее решение уравнения зависит от двух констант, граничные условия дают на них однородную систему $2\times 2$, и нетривиальное решение есть только тогда, когда определитель системы обращается в нуль. Это условие - трансцендентное уравнение на $\lambda$, его корни и образуют спектр. Например, для $y'' + \lambda y = 0$ с условиями Дирихле $y(0) = y(\ell) = 0$ получается $\sin(\sqrt{\lambda}\,\ell) = 0$, откуда $\lambda_n = (n\pi/\ell)^2$, а собственные функции $y_n = \sin(n\pi x/\ell)$ - те самые синусы ряда Фурье.

Собственные функции и число нулей

Собственная функция $y_n$, отвечающая $\lambda_n$, имеет ровно $n-1$ нулей внутри открытого интервала $(a, b)$ - это осцилляционная теорема Штурма. Поэтому основная мода $y_1$ знакопостоянна, вторая $y_2$ меняет знак один раз, и так далее: чем больше собственное значение, тем сильнее «дрожит» функция. Это даёт удобный контроль при численном счёте - если найденная функция для $\lambda_3$ пересекает ось не дважды, в расчёте ошибка. Нули соседних собственных функций к тому же чередуются (теорема о разделении нулей), что напоминает поведение синусов и косинусов с разными частотами.

Ортогональность с весом и нормировка

Самосопряжённость оператора $\mathcal{L}$ влечёт, что собственные функции, отвечающие разным собственным значениям, ортогональны по скалярному произведению с весом $w$:

$\langle y_m, y_n \rangle_w = \int_a^b y_m(x)\, y_n(x)\, w(x)\, dx = 0, \qquad m \neq n.$

Именно весовая функция $w$ из правой части уравнения задаёт «правильное» скалярное произведение - забыть её здесь нельзя. Нормируя каждую функцию на единичную норму $\|y_n\|_w = 1$, получают ортонормированную систему. Её полнота означает, что любую достаточно гладкую функцию $f$, согласованную с краевыми условиями, можно разложить в обобщённый ряд Фурье

$f(x) = \sum_{n=1}^{\infty} c_n\, y_n(x), \qquad c_n = \frac{\langle f, y_n \rangle_w}{\langle y_n, y_n \rangle_w},$

который сходится в среднеквадратичном. Это и есть рабочая лошадка метода разделения переменных: начальные и граничные данные раскладывают по $y_n$, а каждая мода эволюционирует независимо. Сам спектральный приём похож на анализ системы через её собственные значения - как и в случае линейной системы ОДУ с постоянными коэффициентами, где общее решение собирается из мод $e^{\lambda_k x}$ по спектру матрицы.

Знак собственных значений и оценка снизу

Часто важно знать не точные $\lambda_n$, а их знак - например, для устойчивости решения уравнения теплопроводности. Умножив уравнение на $y$ и проинтегрировав по частям, получают отношение Рэлея:

$\lambda = \frac{\int_a^b \big( p\,(y')^2 + q\,y^2 \big)\, dx + \big[\text{граничные члены}\big]}{\int_a^b w\, y^2\, dx}.$

Если $q \geq 0$ и граничные члены неотрицательны, числитель неотрицателен, а знаменатель положителен - значит, все $\lambda_n \geq 0$. Минимальное собственное значение $\lambda_1$ равно минимуму отношения Рэлея по всем допустимым функциям (вариационный принцип), что даёт практичный способ оценивать $\lambda_1$ сверху, подставляя пробные функции, и понимать структуру спектра без решения трансцендентного уравнения.

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

• Забывают весовую функцию $w$ в условии ортогональности. Интеграл $\int y_m y_n\, dx$ без множителя $w$ в общем случае не равен нулю; ортогональность есть только по скалярному произведению с весом.
• Не приводят уравнение к самосопряжённой форме. Без интегрирующего множителя оператор не самосопряжён, и теоремы об ортогональности и действительности спектра формально не применимы.
• Берут неоднородные краевые условия. Условия обязаны быть однородными ($=0$ справа), иначе $y \equiv 0$ не является решением и сама постановка задачи на собственные значения теряет смысл.
• Теряют знак в операторе $-(py')'$. Из-за минуса спектр обычно ограничен снизу и стремится к $+\infty$; если записать $+(py')'$, знаки $\lambda_n$ перевернутся, и анализ устойчивости даст обратный вывод.
• Считают $\lambda = 0$ всегда допустимым. Является ли нуль собственным значением, зависит от краевых условий: для Дирихле - нет, для Неймана - да (с постоянной собственной функцией).

FAQ

Чем собственное значение задачи Штурма-Лиувилля отличается от собственного значения матрицы? Идея та же - нетривиальное решение уравнения $\mathcal{L}[y] = \lambda w y$ при особых $\lambda$. Но оператор $\mathcal{L}$ действует в бесконечномерном пространстве функций, поэтому собственных значений счётно много, они уходят в бесконечность, а собственные «векторы» - это функции.

Зачем приводить уравнение к самосопряжённому виду? Самосопряжённость гарантирует действительность спектра, простоту собственных значений и ортогональность собственных функций с весом. Без неё разложение в ряд по $y_n$ не обосновано, и метод разделения переменных перестаёт работать.

Как связаны задача Штурма-Лиувилля и ряд Фурье? Классический ряд Фурье - частный случай: для $y'' + \lambda y = 0$ с условиями Дирихле или периодическими собственными функциями оказываются синусы и косинусы. Общая теория обобщает разложение на любой самосопряжённый оператор второго порядка.

Коротко

Задача Штурма-Лиувилля - это краевая задача $-(p y')' + q y = \lambda w y$ с однородными разделёнными условиями, в которой ищут собственные значения $\lambda$ и собственные функции $y$. У регулярной задачи спектр действителен, дискретен и неограничен: $\lambda_1 < \lambda_2 < \cdots \to +\infty$, причём $y_n$ имеет ровно $n-1$ нулей. Самосопряжённость оператора влечёт ортогональность собственных функций по скалярному произведению с весом $w$, а их полнота позволяет раскладывать функции в обобщённый ряд Фурье. Знак и оценку $\lambda_n$ даёт отношение Рэлея. На этом фундаменте стоят метод разделения переменных и спектральный анализ уравнений математической физики.