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

Возникает еше другой вопрос: как практически вычислять собственные значения и собственные функџии? Существует систематический способ, который при достаточном терпении приводит к желаемому результату. Существуют также простые или плохо обоснованные способы, спеџифичные для разных частных случаев.

Предположим, что собственные значения существуют и посмотрим, как можно их вычислить. Наша задача записывается так:
\[
\begin{aligned}
\frac{d}{d x}\left(p \frac{d \varphi}{d x}\right)+\lambda \bar{q} \varphi & =0, \\
\left(\bar{\alpha}_{1} \varphi-\bar{\alpha}_{2} \varphi^{\prime}\right)_{0} & =0, \\
\left(\bar{\beta}_{1} \varphi+\bar{\beta}_{2} \varphi^{\prime}\right)_{l} & =0,
\end{aligned}
\]

причем
\[
\bar{\alpha}_{1} \bar{\alpha}_{2}>0, \bar{\beta}_{1} \bar{\beta}_{2}>0 .
\]

Мы можем несколько упростить уравнение (1), введя вместо $x$ новую переменную:
\[
\xi=\int_{0}^{x} \frac{d x}{p(x)} \text {. }
\]

Такое преобразование всегда возможно, так как
\[
p(x)>0 .
\]

Имеем:
\[
\frac{d \varphi}{d x}=\frac{d \varphi}{d \xi} \cdot \frac{1}{p},
\]

и, следовательно, уравнение (1) принимает вид
\[
\frac{d^{2} \varphi}{d \xi^{2}}+\lambda q(\xi) \varphi=0,
\]

где
\[
q(\xi)=p(x) \bar{q}(x) .
\]

Таким образом, наша задача записывается теперь так (вместо $\xi$ мы снова пишем $x$ ):
\[
\begin{aligned}
\frac{d^{2} \varphi}{d x^{2}}+\lambda q(x) \varphi & =0, \\
\left(\alpha_{1} \varphi-\alpha_{2} \varphi^{\prime}\right)_{0} & =0, \\
\left(\beta_{1} \varphi+\beta_{2} \varphi^{\prime}\right)_{l} & =0,
\end{aligned}
\]

где
\[
\alpha_{2}=\overline{\alpha_{2}}, \beta_{2}=\overline{\beta_{2}},
\]

Следует помнить, что здесь $l$ имеет другой смысл, чем в (3). Мы обозначили через $l$ значение $\xi$, соответствуюе прежнему $x=l$. Итак, мы можем рассматривать вместо (1) уравнение, в котором $p(x)=1$.

В обычной классической задаче о решении дифференџиального уравнения задаются $\varphi$ и $\varphi^{\prime}$ для одной точки. В интересующей нас задаче заданы два условия для $\varphi$ и $\varphi^{\prime}$ в двух точках. Но мы будем опираться при решении исследуемой задачи на то, что известно относительно классической задачи.

Уравнение (5) линейно. Если $\varphi_{1}$-решение уравнения (5), то $C \varphi$, где $C$-постоянная, тоже решение.

Известно, что уравнение (5) имеет фундаментальную систему решений $\varphi_{1}$ и $\varphi_{2}$, такую, что любое решение, вся совокупность решений, может быть представлена в виде
\[
\varphi=C_{1} \varphi_{1}+C_{2} \varphi_{2},
\]

где $C_{1}$ и $C_{2}$ – постоянные.
Если заданы $\varphi$ и $\varphi^{\prime}$ для $x=0$, то, подобрав соответственно $C_{1}$ и $C_{2}$, всегда можно найти решение (9), удовлетворяющее этим условиям.

Возьмем для простоты в качестве $p_{1}$ решение, удовлетворяющее при $x=0$ условиям
\[
\varphi_{1}(0)=1, \quad \varphi_{1}^{\prime}(0)=0,
\]

а в. качестве $\varphi_{2}$ – решение, удовлетворяющее условиям
\[
\varphi_{2}(0)=0, \quad \varphi_{2}^{\prime}(0)=1 .
\]

Такие решения существуют в силу основной теоремы теории дифференџиальных уравнений. Выясним, образуют ли эти решения фундаментальную систему. Вспомним, что необходимое и достаточное условие того, что $\varphi_{1}$ и $\varphi_{2}$ образуют фундаментальную систему, таково:
\[
W(x)=\left|\begin{array}{cc}
\varphi_{1}(x) & \varphi_{2}(x) \\
\varphi_{1}^{\prime}(x) & \varphi_{2}^{\prime}(x)
\end{array}\right|
eq 0 .
\]

Для нашего уравнения (5) очень легко доказать теорему, утверждающую, что если детерминант $W(x)$ отличен от нуля для какого-то одного значения $x$, то он отличен от нуля для всех значений $x$. Но в нашем случае
\[
W(0)=\left|\begin{array}{cc}
\digamma_{1}(0) & \varphi_{2}(0) \\
\varphi_{1}^{\prime}(0) & \overbrace{2}^{\prime}(0)
\end{array}\right|=1 .
\]

Следовательно, для любого $x$
\[
W(x):
eq 0
\]

и решения $\rho_{1}$ и $\varphi_{2}$ образуют фундаментальную систему.
Решим теперь основной вопрос: существует ли решение, удовлетворяющее граничным условиям (6) и (7). Запишем их сокращенно в таком виде:
\[
\begin{array}{l}
l_{1}(0)=0 ; \\
l_{2}(v)=0 .
\end{array}
\]

Операторы $l_{1}$ и $l_{2}$ линейны:
\[
\begin{array}{c}
l_{1,2}\left(\varphi_{1}+\varrho_{12}\right)=l_{1,2}\left(\rho_{1}\right)+l_{1,2}\left(\varrho_{2}\right) ; \\
l_{1,2}(C)=C l_{1,2}(\%) .
\end{array}
\]

Если краевая задача (5) – (7) имеет решение, то оно содержится в (9), так как в этом семействе функций содержатся все решения уравнения (5). Таким образом, вопрос сводится к тому, можно ли подобрать $C_{1}$ и $C_{2}$ так, чтобы удовлетворить условиям (6) и (7).

Укажем явно, что наши решения $\varphi_{1}$ и $\varphi_{2}$ зависят от параметра $\lambda$, т. е. будем писать их в виде
\[
p_{1}=\varphi_{1}(x, \lambda), \quad \varphi_{2}=\varphi_{2}(x, \lambda) .
\]

Потребуем, чтобы решение (9) удовлетворяло граничным условням:
\[
\left.\begin{array}{l}
C_{1} l_{1}\left(\varphi_{1}\right)+C_{2} l_{1}\left(\varphi_{2}\right)=0 \\
C_{1} l_{2}\left(\varphi_{1}\right)+C_{2} l_{2}\left(\varphi_{2}\right)=0
\end{array}\right\}
\]

Эти уравнения линейны и однородны по отношению к $C_{1}$ и $C_{2}$. Они всегда имеют тривиальное решение:
\[
C_{1}=C_{2}=0,
\]

но оно нас ие интересует. Нам нужно, чтобы по крайней мере одна из констант $C_{1}$ и $C_{2}$ не равнялась нулю. Уравнения (14) имеют нетривиальное решение только тогда, когда
\[
\Delta=\left|\begin{array}{ll}
l_{1}\left[\varphi_{1}(\lambda)\right] & l_{1}\left[\varphi_{2}(\lambda)\right] \\
l_{2}\left[\varphi_{1}(\lambda)\right] & l_{2}\left[\varphi_{2}(\lambda)\right]
\end{array}\right|=0 .
\]

Левая часть есть некоторая трансцендентная функция от $\lambda$. Нетривиальные решения для $C_{1}$ и $C_{2}$ могут быть только при таких $\%$, которые удовлетворяют трансџендентному уравнению (15). Таким образом, сразу видно, что краевая задача имеет решение не при всяком $\lambda$.

Итак, мы получили трансџендентное уравнение для определения интересующих нас чисел $\lambda$. Нам нужно уметь его решать, хотя бы приближенным способом.

Вернемся к уже рассмотренному случаю однородного стержня и посмотрим, что там дает изложенный способ. Возьмем простейшие краевые условия
\[
\varphi(0)=0, \quad \varphi(l)=0
\]

и положим для простоты:
\[
q(x)=1 .
\]

Пойдем только что указанным „лобовым“, систематическим путем. Для случая (16) мы можем написать в явном виде фундаментальную систему $\varphi_{1}$ и $\varphi_{2}$ :
\[
\varphi_{1}=\cos \sqrt{\lambda} x, \quad \varphi_{2}=\frac{1}{\sqrt{\lambda}} \sin \sqrt{\lambda} x .
\]

Уравнения (14) здесь таковы:
\[
\begin{array}{c}
C_{1} \cdot 1+C_{2} \cdot 0=0, \\
C_{1} \cos \sqrt{\lambda} l+C_{2} \frac{1}{\sqrt{\lambda}} \sin \sqrt{\lambda} l=0 .
\end{array}
\]

Условие (15) имеет, следовательно, вид
\[
\left|\begin{array}{cc}
1 & 0 \\
\cos \sqrt{\lambda} l & \frac{1}{\sqrt{\lambda}} \sin \sqrt{\lambda} l
\end{array}\right|=0,
\]

откуда
\[
\sin \sqrt{\lambda} l=0 .
\]

Это уже известное нам трансџендентное уравнение для собственных значений однородного стержня.

Вся трудность общей задачи (о неоднородном стержне)в решении трансцендентного уравнения (15). В большинстве случаев оно находится графически, но мы займемся принџипиальной стороной вопроса. Мы докажем, что уравнение (15) всегда имеет бесчисленное множество решений, образуюших дискретную совокупность. Это утверждение связано с определенным типом краевых условий. Аналогично обстоит дело и в случае периодических краевых условий
\[
\varphi(0)=\varphi(l), \quad \varphi^{\prime}(0)=\varphi^{\prime}(l)
\]

или похожих условий
\[
\varphi(0)=-\varphi(l), \quad \varphi^{\prime}(0)=-\varphi^{\prime}(l) .
\]

Но если, например, мы ставим условия
\[
\vartheta(0)=0, \varphi^{\prime}(0)=0,
\]

то детерминант $\perp=1$ и нельзя подобрать $\lambda$, удовлетворяющее уравнению (15).

В случае краевых условий (17) имеем $\perp=0$, и, следовательно, краевым условиям можно удовлетворить при любом $\lambda$. В этом легко убедиться путем элементарного рассуждения. Перенесем начало координат в середину стержня. Тогда условия (17) примут вид
\[
?\left(-\frac{l}{2}\right)=\varphi\left(\frac{l}{2}\right), \quad q^{\prime}\left(-\frac{l}{2}\right)=-\eta^{\prime}\left(\frac{l}{2}\right) .
\]

Им удовлетворяет любая функџия вида $\cos \sqrt{\lambda} x$, так как косинус четная функџия, а его производная – нечетная.

Этот аппарат очень красив на бумаге, но его трудно применить практически, так как трудно найти фундаментальную систему уравнения (5) при $q$, зависящем от $x$. В первый момент может даже показаться, что этот способ практически ничего не дает.

Надо указать путь для вычисления с достаточным приближением $\varphi_{1}$ и $\varphi_{2}$ как функџий $\lambda$. Общая теория говорит, что $\varphi$ является непрерывной функцией от $\lambda$, а кроме того (для всякого заданного $x$ ) она является целой функџией от $\lambda$. (Џелая функџия означает-не обращающаяся в конечной области в бесконечность; $\cos \lambda, e^{\lambda}-$ џелые функџии $\lambda$ ). Итак, $\varphi$ – џелая трансџендентная функция $\lambda$, и, следовательно, ее можно искать в виде ряда по степеням $\lambda$ :
\[
\varphi(x, \lambda)=u_{0}(x)+\lambda u_{1}(x)+\lambda^{2} u_{2}(x)+\ldots
\]

Представив решение в таком виде, можно обрезать ряд на какомнибудь члене и подставить полученное таким путем приближенное выражение в трансдендентное уравнение (15). Можно указать редепт нахождения ряда (18) и доказать, что ряд сходится для всех конечных значений $\lambda$. То и другое легко сделать.
Продифферендируем дважды ряд (18):
\[
\varphi^{\prime \prime}=u_{0}{ }^{\prime \prime}+\lambda u_{1}{ }^{\prime \prime}+\lambda^{2} u_{2}{ }^{\prime \prime}+\ldots
\]

Подставим (18) и (19) в уравнение (5) и приравняем нулю коэффициенты при всех степенях $\lambda$. Мы получим:
\[
u_{0}{ }^{\prime \prime}=0, u_{1}{ }^{\prime \prime}+q u_{0}=0, \ldots u_{n}{ }^{\prime \prime}+q u_{n-1}=0 .
\]

Такова последовательность уравнений, которым удовлетворяют функции $u_{n}$.

Условия (10), которым должно удовлетворять одно из фундаментальных решений, будут выполнены, если потребовать
\[
u_{0}(0)=0, u_{0}{ }^{\prime}(0)=1 \text {, }
\]

а для всех $n>0$
\[
u_{n}(0)=0, u_{n}{ }^{\prime}(0)=0 .
\]

Мы получаем в результате:
\[
u_{0}=x .
\]

Каждая функџия $u_{n}$ определяется из предыдущей $u_{n-1}$ простой квадратурой:
\[
u_{n}=-\int_{0}^{x} d \xi \int_{0}^{\xi} q\left(n_{i}\right) u_{n-1}\left(r_{i}\right) d n
\]

Итак, нахождение $\varphi_{1}(x)$ [и аналогично $\varphi_{2}(x)$ ] сведено к последовательному ряду квадратур. Мы оперировали с рядом формально, но можно показать, что ряд (18) сходится.

Оборвав ряд (18) на некотором члене, мы заменим трансцендентное уравнение (15) алгебраическим уравнением той или иной степени. При достаточно большом числе членов его корни будут как угодно близки к корням трансџендентного уравнения. Практически для основного колебания часто достаточно взять $3-4$ члена разложения. Для обертонов приходится брать больше членов и вычисление становится хлопотливым.

Техника выработала хорошие способы приближенного вычисления, но пока еще мало оправданные.

Перейдем к доказательству того, что существует бесчисленное множество собственных значений. Доказательства этой теоремы были даны уже давно. При этом доказательства проводились для некоторых частных граничных условий, а затем к ним сводили остальные граничные условия. Недавно появилась работа Прюфера (она изложена в третьем издании книги Бибербаха ${ }^{1}$ ), в которой доказательство проводится сразу для самого общего случая задачи Штурма-Лиувилля. Изложим әто доказательство.
Введем новую переменную
\[
y=\frac{d \varphi}{d x}
\]

и нап̈ишем вместо уравнения (5) два уравнения:
\[
\frac{d \varphi}{d x}=\gamma, \quad \frac{d \varphi}{d x}+\gamma q ?=0 .
\]

Положим:
\[
\%=p \sin \theta, \quad \%=\rho \cos \theta \text {. }
\]

Мы вводим, таким образом, две новые функции: $\hat{\imath}$ и $\theta$. Имеем:
\[
\left.\begin{array}{l}
\frac{d \varphi}{d x}=\rho^{\prime} \sin \theta+\theta^{\prime} \cos \theta ; \\
\frac{d x}{d x}=\rho^{\prime} \cos \theta-p \theta^{\prime} \sin \theta
\end{array}\right\}
\]

Подставляя (22) в (20), получаем уравнения:
\[
\begin{array}{l}
\rho^{\prime} \sin \theta+\rho \theta^{\prime} \cos \theta-\rho \cos \theta=0 \\
\rho^{\prime} \cos \theta-\rho \theta^{\prime} \sin \theta+\lambda q \rho \sin \theta=0
\end{array}
\]

Мы можем их преобразовать так, чтобы получилось одно уравнение для $\theta^{\prime}$ и другое для $\rho^{\prime}$. Умножим (23) на $\cos \theta$, а (24)на $\sin \theta$, и вычтем. Это дает по сокращении на? (мы далее увидим, что $\rho$ нигде не обращается в нуль):
\[
\theta^{\prime}=\cos ^{2} \theta+\lambda p \sin ^{2} \theta \text {. }
\]
1 [L. Bieberbach. Theorie der Differentialgleichungen. 3-e переработ. издание, Берлин, 1930.]

Далее умножим (23) на $\sin \theta$, а (24) – на $\cos \theta$, и сложим. Мы получим:
\[
\hat{q}^{\prime}=p(1-\lambda p) \sin \theta \cos \theta \text {. }
\]

Для 9 получилось уравнение (25), не содержащее $р$. Уравнение (26) решается посредством квадратуры
\[
\ln |\rho|=\int_{0}^{x}(1-\lambda p) \sin \theta \cos \theta d x+c .
\]

Чтобы ее вычислить, нужно сначала решить уравнение для $\theta$, а потом подставить решение в (26). При этом получатся две постоянные интегрирования.

Первый же шаг продвинул нас довольно значительно вперед: мы получили для $\theta$ одно уравнение. Но это куплено той џеной, что уравнение нелинейное. Можно легко модифиџировать подход так, чтобы было удобно находить решение количественно для случая медленно или мало меняюџихся $q(x)$. Но нас интересует сейчас другое – доказательство существования неограниченной последовательности собственных значений. Сейчас мы увидим, в чем заключается сила излагаемого метода.
Перепишем граничные условия (6) и (7) в таком виде:
\[
\left.\begin{array}{l}
\left.\left(\alpha_{1} ?-\alpha_{i}\right\rangle\right)_{0}=0, \\
\left(\alpha_{1} p+\alpha_{2} \chi\right)_{i}=0
\end{array}\right\}
\]

и подставим в них, вместо $\%$ и функии $\rho$ и $\theta$ :
\[
\left(\alpha_{1} \sin \theta-\alpha_{2} \cos \theta\right)_{0}=0, \quad\left(\psi_{1} \sin \theta+\beta_{2} \cos \theta\right)_{l}=0 .
\]

Мы видим, что в граничные условия вошла только 0 . Нужно удовлетворить граничным условиям, куда входит только $\theta$, тогда будут удовлетворены и исходные граничные условия (28).

Так как $\alpha_{1}$ и $\alpha_{2}$ – одного знака, можно удовлетворить первому граничному условию, взяв для $\theta_{\text {, }}$ определенное значение $і$ в первом квадранте:
\[
0<\gamma_{1}<\frac{\pi}{2} .
\]

Но $\beta_{1}$ и $\beta_{2}$ тоже имеют одинаковые знаки и, следовательно, в первом квадранте заведомо нельзя удовлетворить второму граничному условию. Можно пойти во второй квадрант, там мы найдем то, что нам нужно. Мы удовлетворим второму граничному условию, если возьмем $\theta_{l}=\gamma_{2}$, причем
\[
\pi_{2}<i_{2}<\pi \text {. }
\]

Теперь вся задача прояснилась. Найдя $\theta(x)$, удовлетворяющее уравнению (25) и равное $\gamma_{1}$ для $x=0$ и $\gamma_{2}$ для $x=l$, мы решим нашу основную задачу.

Если $\theta(x)$ для $x=0$ равно $\gamma_{1}+\pi, \gamma_{1}+2 \pi, \gamma_{1}+3 \pi, \ldots$, а для $x=l$ соответственно: $\gamma_{2}+\pi, \gamma_{2}+2 \pi, \gamma_{2}+3 \pi, \ldots$, то это дает то же самое, поскольку тангенс периодичен с периодом $\pi$ и краевые условия удовлетворяются, а $\sin ^{2}$ и $\cos ^{2}$ тоже периодичны с периодом $\pi$; и, следовательно, уравнение (25) тоже удовлетворится. Мы получим другое решение для $\theta$, а именно: $\theta_{2}=\theta_{1}+\pi$, но это не даст ничего нового по отношению к первому решению.
Но можно поставить другую задачу: оставить первое граничное условие $\theta(0)=\gamma_{1}$ и задать новое граничное условие для $x=l$, а именно: $\theta(l)=$ $=\gamma_{2}+n \pi$. Если мы решим эmy задачу, то для $\theta(x)$ получится другая функџия, отличная от $\theta(x)$ при $n=0$.

Если при любом целом $n$ мы сможем найти такое $\lambda$, что соответствующее $\theta(x)$ удовлетворяет граничным условиям
\[
\theta(0)=\gamma_{1}, \quad \theta(l)=\gamma_{2}+n \pi,
\]

то этим будет доказано, что исходное уравнение (20) имеет бесконечное множество собственных значений.

Доказательство проводится вполне строго. Оно становится гораздо нагляднее, если поясняется геометрически.

Нелинейное уравнение (25) мы решить не можем. Но мы можем вывести из него интересующие нас свойства функџии $\theta(x)$. Из уравнения (25) видно, что при $\lambda$ положительном $\theta(x)$ – монотонно возрастающая функция. В самом деле, ее производная всюду положительна (она не может обрашаться в нуль, так как $\cos \theta$ и $\sin \theta$ не могут быть равны нулю одновременно).

Рассмотрим уравнение (25) при различных $\lambda$. Величину $\lambda q(x)$ обозначим через $\sigma(x)$ :
\[
\theta^{\prime}=\cos ^{2} \theta+\sigma(x) \sin ^{2} \theta .
\]

Обозначим индексами 1 и 2 две различные функџии $\sigma(x)$ и соответствующие им $\theta(x)$ :
\[
\begin{array}{l}
\theta_{1}^{\prime}=\cos ^{2} \theta_{1}+\sigma_{1}(x) \sin ^{2} \theta_{1}, \\
\theta_{2}^{\prime}=\cos ^{2} \theta_{2}+\sigma_{2}(x) \sin ^{2} \theta_{2} .
\end{array}
\]

Пусть всюду
\[
\sigma_{2}(x)>\sigma_{1}(x) \text {. }
\]

Кривые $\theta_{1}(x)$ и $\theta_{2}(x)$ обе монотонно поднимаются с ростом $x$. Они могут пересекаться. Если они пересекаются, то только так, что при увеличении $x$ за точку пересечения $\theta_{2}$ идет выше $\theta_{1}$ (рис. 169, a). Они не могут пересечься так, как на рис. 169 , б поскольку производная $\theta_{1}^{\prime}$ меньше производной $\theta_{2}^{\prime}$ ( $\sigma_{1}$ меньше, чем $\sigma_{2}$ ). Значит, в точке пересечения $\theta_{1}^{\prime}<\theta^{\prime}$. Но, может быть, они пересекаются в точке, где $\sin \theta=0$ ? Для этого значения $x$ обе производные одинаковы. Тогда кривые касаются. Я утверждаю, однако, что и здесь кривая $\theta_{2}$ пойдет выше. Действительно, кривые касаются в одной точке (они не могут касаться по џелому отрезку). Для соседних интегральных кривых $\sin \theta_{2}
eq 0$. Различные интегральные кривые $\theta_{2}$ пересекаться между собой не могут (особых точек нет). Значит, и в этом случае после пересечения (касания) кривая $\theta_{2}$ пойдет выше $\theta_{1}$. В итоге мы заключаем, что если кривые $\theta_{1}$ и $\theta_{2}$ один раз пересеклись, то они уже не могут пересечься второй раз.

Будем теперь сравнивать функци $\theta(x)$ для различных і при одном и том же $q(x)$. Будем рассматривать решение как функџию от $\lambda$. Функция $\sigma(x)=\lambda$ растет при любых $x$ вместе с $\lambda$.

Теперь ясно, как мы докажем, что существует бесконечная последовательность собственных значений.

Если $\theta_{1}$ и $\theta_{2}$, соответствуюшие $\lambda=\lambda_{1}$ и $\lambda=\lambda_{2}\left(\lambda_{2}>\lambda_{1}\right)$, выходят при $x=0$ из одной точки, то значит, при $x>0$ всюду $\theta_{2}>\theta_{1}$, так как кризые не могут пересекаться. Чем больше значение $\lambda$, тем выше проходит кривая $\theta$ и тем выше она пересекает прямую $x=l$. Для того, чтобы догазать то, что нам нужно, необходимо доказать, что при $\lambda=0$ кривая $0(x)$ пересекает прямую $x=l$ ниже точки $\gamma_{2}$, а при $\lambda \rightarrow \infty$ точка пересечения с прямой $x=l$ не стремится к конечному пределу. Первое нужно для того, чтобы доказать существование нацменьшего собственного значения, второе – для того, чтобы убедиться в существовании бесконечного числа собственных значений без точек сгущения.

Если эти два утверждения удастся доказать, то теорема будет доказана.