Московский физико-технический институт
Опубликован: 25.10.2007 | Доступ: свободный | Студентов: 3916 / 1197 | Оценка: 4.50 / 4.33 | Длительность: 24:00:00
ISBN: 978-5-9556-0065-9
Специальности: Программист, Математик

Лекция 11: Численное решение краевых задач для систем обыкновенных дифференциальных уравнений

< Лекция 10 || Лекция 11: 123456789

10.7. Краевые задачи на собственные значения для обыкновенных дифференциальных уравнений

Краевые задачи на собственные значения достаточно часто встречаются в физических приложениях. Например, это задача определения собственных колебаний струны, сводящаяся к ОДУ вида

$ \frac{d}{dt} [k(t)\frac{du}{dt}] + \lambda r(t)u = 0.  $

В приведенном уравнении краевые условия зависят от способа закрепления струны.

Это и задачи собственных колебаний упругого стержня (ОДУ четвертого порядка), нахождения энергетических уровней атома водорода, вычисление критических нагрузок в теории стержней и оболочек и др.

В задачах на собственные значения добавляется еще один пристрелочный параметр\lambda, поэтому эти задачи часто решаются методом стрельбы. Приведем простейший пример — одно дифференциальное уравнение первого порядка с двумя краевыми условиями (второе краевое условие появляется из - за присутствия неизвестного параметра \lambda ):

$ \frac{du}{dt} + f(u, t, \lambda ) = 0, u(0) = u_1 , $

Если не учитывать правое краевое условие, то получим задачу Коши. Ее численное интегрирование приводит к некому значению на правом конце, зависящему от \lambda и, вообще говоря, не равному u2. Варируя параметр \lambda, можно добиться выполнения правого краевого условия с некоторой заданной точностью. При этом, разумеется, используются методы численного нахождения корней алгебраического уравнения, обычно метод касательных.

Второй пример — краевая задача на собственные значения для ОДУ второго порядка с нулевыми краевыми условиями:

\begin{gather*}
\frac{d^2 u}{dt^2 } + {b^{\prime}}(t)\frac{du}{dt} + [{c^{\prime}}(t) + \lambda ]u = 0, \\ 
u(0) = u(L) = 0.
\end{gather*}

Поскольку это уравнение второго порядка с неизвестным параметром \lambda (собственное значение дифференциального оператора), то для его решения требуется третье условие. Однако в силу линейности и однородности задачи решение определяется с точностью до произвольного постоянного множителя, что и является неявным заданием третьего условия. Его можно задать, например, следующим образом:

$ \frac{du(a)}{dt} = 1.  $

Трудности при использовании метода стрельбы возникают, если соответствующая задача Коши плохо обусловлена, а также в случае жестких краевых задач. В этих случаях появляется сильная зависимость численного решения от пристрелочного параметра \lambda.

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

$ \frac{d^2 u}{dt^2 } + \lambda u = 0, u(0) = u(L) = 0, $

имеющую точное решение

$  \lambda_k = \left({\frac{{\pi  k}}{L}}\right)^2, u_k (t) = \sin \left({\frac{{\pi  k t}}{L}}\right), k = 1, 2, \ldots  $

Непосредственной подстановкой показывается, что решениями соответствующей разностной задачи на собственные значения

$ \frac{{u_{n - 1} - 2u_n + u_{n + 1}}}{{\tau ^2 }} + \lambda u_n = 0, n = 1, 2, \ldots , N - 1, \tau N = L,  u_0 = u_n = 0, $

являются собственные значения и собственные функции

$  \lambda_k = \frac{4}{\tau^2} \sin^2 \frac {\pi k\tau }{2N}, \quad 
u^{\tau}_j = \sin \frac{\pi {kt}_n}{L}, \quad k, n = 1, 2, \ldots , N - 1, $

откуда видно, что

$  \lim \limits_{\tau \to 0} \frac{4}{\tau^2} \sin \frac{\pi k\tau}{2L} = 
(\frac{\pi k}{L})^2 = \lambda_k, $

т.е. имеет место сходимость решения разностной задачи к решению дифференциальной задачи \mathop {\lambda_k^\tau   \to \lambda_k}\limits_{\tau  \to 0} .

< Лекция 10 || Лекция 11: 123456789