
- •Федеральное государственное бюджетное
- •Глава 1 посвящена изучению общих вопросов теории уравнений (1) с точки зрения известных свойств операторов.
- •§2. Исследование поведения решений однородного уравнения на бесконечности
- •§3. Ограниченность решений однородного уравнения
- •§4.Условия существования ограниченного решения у неоднородного уравнения
- •Глава 2. Численные методы приближенного решения задачи Коши для однородных дифференциальных уравнений в банаховом пространстве
- •§5. Общие сведения о дифференциальных уравнениях
- •§6. Численные методы решения задачи Коши для обыкновенного дифференциального уравнения первого порядка
- •§ 7.Метод Эйлера
- •Решение.
- •§8.Метод Эйлера-Коши
- •Решение.
- •§ 9. Метод Рунге-Кутта
- •Решение.
- •Реализация метода Рунге-Кутта с помощью системы Maxima
- •§10. Решение краевых задач для линейных дифференциальных уравнений второго порядка методом конечных разностей
- •Решение.
- •Реализация конечно-разностного метода
- •§ 11. Метод сеток для решения дифференциальных уравнений в частных производных
- •Решение.
- •Решение.
- •Выполним реализацию конечно-разностного метода в системе компьютерной математики Maxima
- •Заключение
- •Список литературы
§2. Исследование поведения решений однородного уравнения на бесконечности
Поведение
решений задачи
на бесконечности существенно зависит
от расположения спектра оператора
.
Предположим,
что спектрлежит
внутри левой полуплоскости. Тогда из
следует на основании теоремыI.1.4
[1] оценка
при
любых
и
некоторых положительных постоянных
.
Верно
и обратное: если оценка
выполняется
для всякого решения
уравнения
,
то спектр
лежит
внутри левой полуплоскости. Действительно,
из
следует,
что
и остается применить теорему I.1.4 [1]:
Поведение
решения уравнения
в предположении, что спектр
лежит в левой полуплоскости, можно
охарактеризовать более точно, если
ввести новую эквивалентную норму в
пространстве
,
пользуясь формулой
Оказывается,
в этой норме решения уравнения
стремятся при
к нулю монотонно. Действительно,
И, таким образом,
Рассмотрим
теперь случай, когдапричем спектральное множество
непусто.
Пусть
-
спектральные проекторы, соответствующие
этому разложению спектра, и
- соответствующее прямое разложение
на инвариантные подпространства
оператора
.
Так как
инвариантны относительно операторов
,
то решение
уравнения
,
начинающееся в каком - нибудь из них,
уже не выходит из соответствующего
подпространства.
Введем
в
индефинитную норму[8]
Произведя
вычисления, аналогичные проведенным
выше, получим
Поэтому
и,
следовательно, индефинитная норма
любого решения уравнения
убывает.
Рассмотрим два частных случая.
Пусть
, т.е.
. В этом подпространстве
и нормы
и
эквивалентны. Из
следует, что
монотонно стремится к нулю.
Таким
образом, решения
уравнения
с начальным вектором
из подпространства
стремятся к нулю.
Пусть
, т.е.
.В этом подпространстве
является обычной нормой, эквивалентной
. Из
следует, что
Проинтегрировав это неравенство, мы видим, что
т.
е решения с начальными значениями из
неограниченно возрастают при
.
Заметим,
что из разложения
следует, что любое решение, для которого
,
неограниченно возрастает. В частности,
если весь спектр лежит внутри правой
полуплоскости, то
и ненулевое решение
уходит на бесконечность при
,
причем норма
монотонно возрастает.
Определение
2.1 Оператор
со спектром, распадающимся на два
спектральным множества соответственно
внутри правой и левой полуплоскости
каждое из которых непусто, называетсяэ-дихотомическим.
Вместе с ним будем называть э -
дихотомическим и дифференциальное
уравнение
.
Очевидно,
что фазовое пространство
э-дихотомического
уравнения
распадается в прямую сумму
,
причем решения, начинающиеся в
,
экспоненциально убывают, а решения,
начинающиеся в
,
экспоненциально растут.
Если
предположить, фазовое пространство
гильбертово,
то все рассуждения естественно проводить,
рассматривая вместо норм
гильбертовы нормы
[12].
В этом случае проведенная выше
перенормировка имеет простой геометрический
смысл.
Пусть
.
Тогда для решения уравнения
справедливо
соотношение
из
которого следует, что для угла
между направлением радиус-вектора
и вектора
,
касательного к интегральной кривой
уравнения
,
имеет
место оценка
Таким
образом,
и векторное поле касательных к интегральным
кривым уравнения
в каждой точке направлено существенно
внутрь сферы с центром в нуле, проходящей
через эту точку.
Рассмотрим
теперь более общую ситуацию, когда
спектр
лежит внутри левой полуплоскости. В
этом случае по теоремеI.
5.1.1[1] существует
ограниченный равномерно положительный
оператор
,
такой, что
Приведенные
ранее соображения останутся справедливыми,
если оценки производить в новой метрикеэквивалентной прежней.
В
самом деле, ограниченность и равномерная
положительность оператора
обеспечивают топологическую эквивалентность
норм.
С
другой стороны, если
,
то для решения
уравнения
будем иметь
Напомним,
что требуемый оператор
можно
получить как решение уравнения
где
- произвольный, равномерно - положительный
оператор. Полагая, например,
,[5]
мы получим
т.е.
Система
сфер в этом случае заменяется системой
эллипсоидовс
центром в нуле, внутрь которых входят
интегральные кривые. Рассмотрение в
банаховом пространстве нормы
,
аналогичной
обобщает приведенные выше геометрические
соображения, связанные с так называемым
методом Ляпунова[12].
Роль
эллипсоидов при этом играют выпуклые
центрально-симметрические тела,
ограниченные поверхностями
Несколько
более сложная, но достаточно ясная
геометрическая картина получается и в
том случае, когда спектр оператора
имеет также компоненту, внутри правой
полуплоскости, т.е. когда уравнение
является э-дихотомическим.
Рассмотрим
случай, когда фазовое пространство
гильбертово.
Справедлива теорема [12]
Теорема
2.1 Для
того, чтобы уравнение
было э - дихотомическим,
необходимо и достаточно, чтобы оператор
был равномерно
диссипативным
по отношению к некоторому эрмитову
индефинитному оператору
:
При
любом выборе оператораудовлетворяющего
этому соотношению, инвариантное
подпространство
оператора
,
соответствующее части спектра
,
лежащей внутри правой (левой) полуплоскости,
равномерно
отрицательно
(
положительно).
Индефинитная
форма
порождает в гильбертовом пространстве
две системы гиперболоидов: плюс
гиперболоиды
определяются уравнением
и минус-гиперболоиды
– при
.
Нетрудно
проверить, что из равномерной
диссипативности
оператора
следует для решений уравнения
соотношение
показывающее,
что индефинитная форма
убывает для любого решения
этого уравнения.
Если
точка
находится
на минус
гиперболоиде, то дальнейшее уменьшение
величины
с ростом
означает,
что траектория
пересекает гиперболоиды со все большими
по абсолютной величине отрицательными
значениями
,
а следовательно, уходит на бесконечность.
В
случае, когда точка
находится
на плюс – гиперболоиде, в подпространстве
уменьшение формы
которая на
эквивалентна обычной, приводит к
неограниченному приближению траектории
к центру.
Если
же точка
находится
на плюс – гиперболоиде, вне подпространства
,
то пересекая гиперболоиды
с
уменьшающимися значениями
,
траектория выходит на конус
и переходит на систему минус
гиперболоидов, удаляясь затем на
бесконечность.
Это
следует из того, что любое решение, для
которого
неограниченно возрастает по норме.
Заметим
теперь, что подбирая для э
дихотомического оператора
эрмитов оператор
из соотношения
где
спектральные
проекторы оператора
соответствующие инвариантным
подпространствам
мы придем к форме
Эта
форма аналогична (при
)
индефинитной норме
в банаховом пространстве
.