- •1.1 Дифференциальные уравнения первого порядка
- •1.2 Интегрируемые типы уравнений первого порядка
- •1.3 Формулировка теоремы существования и единственности для нормальной системы
- •1.4 Сведение общей системы дифференциальных уравнений к нормальной
- •1.5 Комплексные дифференциальные уравнения
- •2.1 Линейное однородное уравнение n-го порядка с постоянными коэффициентами
- •2.2 Линейное неоднородное уравнение с постоянными коэффициентами
- •Часть 2
- •1.1 Сведения из линейной алгебры и анализа
- •1.2 Основная теорема теории линейных систем уравнений с постоянными коэффициентами
- •1.3 Практическое вычисление матрицы экспоненты
- •2.1 Нормальная система линейных уравнений переменными коэффициентами
- •2.2 Линейное уравнение n-го порядка с переменными коэффициентами
- •1.1 Автономные системы дифференциальных уравнений и их фазовые пространства
- •1.2 Фазовая плоскость линейной однородной системы с постоянными коэффициентами
- •2.1 Основные понятия
- •2.2 Устойчивость положений равновесия
- •2.3 Предельные циклы
- •1.1 Сжатые отображения
- •1.2 Доказательство теорем существования и единственности для одного уравнения
- •1.2.1 Редукция к интегральному уравнению
Часть 2
1.1 Сведения из линейной алгебры и анализа
Матрица линейного оператора Пусть Xn - конечномерное линейное пространство размерности n над полем K. Линейным преобразованием A пространства Xnназывается отoбражение
обладающее свойствами:
Такие линейные преобразования мы будем назвать линейными операторами в Xn.
Всякий
линейный оператор A в Xn может
быть задан своей матрицей в некотором
базисе. Делается это следующим образом:
выберем в пространстве Xn какой-либо
базис
в
силу которых возникает
матрица Обратно, с помощью любой матрицы размерности n×n можно задать линейный оператор в Xn, матричное представление которого в некотором базисе совпадает с исходной матрицей . Действительно, выбрав в Xn некоторый базис , мы положим
Построим
теперь линейное преобразование A пространства Xn,
переводящее векторы
соответственно
в
Затем вектору сопоставим вектор
Преобразование,
переводящее
в Ясно, что матрица оператора A в базисе e = ( ) есть матрица .
Формула
разложения (1.4) позволяет
установить взаимно - однозначное
соответствие между
векторами
Xn вещественного
(комплексного) пространства и точками
координатного пространства
Легко видеть, что это отображение является изоморфным. Иными словами, любое конечномерное пространство можно отождествить с изоморфным ему координатным пространством соответствующей размерности. Пусть выбранная система координат есть и вектор имеет в ней координаты x1,,xn. Обозначим элементы матрицы Aпреобразования A, вычисленные в той же системе координат, через aij (i, j = 1,,n). Мы имеем
Согласно определению матрицы оператора A
Подставляя эти значения в предыдущее равенство, получим
Следовательно,
координаты нового вектора
Таким образом, координатный столбец нового вектора равен произведению матрицы линейного преобразования на координатный столбец старого вектора
Преобразование координат В предыдущем разделе установлено взаимно-однозначное соответствие между линейными операторами, действующими в n-мерном линейном пространстве и квадратными матрицами порядка n. Однако для этого потребовалось сначала выбрать в Xn определенную координатную систему. Изменив ее, мы изменим соответствие. В результате одному и тому же линейному оператору A в старой и новой координатных системах будут отвечать различные матрицы. Приведем связывающую их формулу.
Пусть
в Xn имеются
две координатные системы:
и
Матрица
называется матрицей
перехода от
координатной системы
к
системе
.
Возьмем произвольный вектор
Подставляя
во второе из этих равенств выражения
для векторов
Таким образом, приходим к формулам
Эти равенства и являются искомыми формулами преобразования координат. Матричная запись формул (1.8):
Пусть A -
матричное представление линейного
оператора A в
базисе
,
а
а по правилу преобразования координат (1.9) имеем
Из соотношений (1.10), (1.11) следует, что
Так
как вектор
Заметим, что с точки зрения теории линейных преобразований матрица перехода есть матрица линейного преобразования, переводящая старую координатную систему в новую, вычисленная в старой системе координат.
Пусть A: Xn Xn -
линейный оператор и A -
матрица, соответствующая оператору A в
некоторой системе координат. Отличный
от нуля вектор
Оказывается, что число тогда и только тогда является собственным значением оператора A, когда оно есть корень многочлена
Многочлен D()
называют характеристическим
многочленом матрицы A.
Известно, что коэффициенты многочлена D()
не зависят от выбора системы координат
в Xn,
а полностью определяются самим
оператором A.
Поэтому многочлен D()
называется тожехарактеристическим
многочленом оператора A,
а уравнение (1.14) характеристическим
уравнением этого оператора. Координатный
столбец
Овеществление
Через
Пусть
-
базис в
.
Каждый вектор
где
Таким
образом, легко видеть, что овеществление
пространства
,
будет 2n-мерным
вещественным линейным пространством С точки зрения оперирования с координатами овеществление означает следующее: Фиксируя в некоторую координатную систему, мы каждому вектору x можем сопоставить n-мерный координатный столбец [x] комплексных чисел (состоящий из комплексных координат данного вектора). После овеществления , этому же элементу сопоставляется 2n-мерный координатный столбец действительных чисел; первыми n элементами его являются действительные части комплексных координат, а затем в том же порядке перечисляются мнимые части.
Как
и в вещественном линейном пространстве
любой линейный оператор A:
допускает
матричное представление. Пусть
в
фиксирована
какая-либо система координат, скажем,
то по определению
Пусть
Таким образом,
Переписав формулы (1.15) в виде
и заметив, что
можно сделать следующий вывод:
Линейный
оператор A:
порождает
линейный оператор
Оператор
Комплексификация
Рассмотрим
теперь так называемую
операцию комплексификации. Пусть
Элементы
пространства
-
это пары
Операция
умножения элемента
на
комплексное число
Легко
проверить, что
-
комплексное линейное пространство.
Пусть (
)
- базис в
и,
следовательно, любой вектор
Пусть
Тогда элемент записывается в виде:
Другими
словами, векторы
Пусть A:
-
линейный оператор. Комплексификацией оператора A называют
линейный оператор
Матричная экспонента
Скалярную
функцию
где E обозначает единицу. Отождествляя вещественное n-мерное линейное пространство Xn (произвольной природы) с изоморфным ему координатным пространством , мы часто само пространство Xn будем обозначать через . Итак, пусть A: - линейный оператор. Чтобы определить экспоненту eA, потребуется понятие предела последовательности линейных операторов.
Говорят,
что в вещественном линейном
пространстве E задано
скалярное произведение, если каждой
паре векторов
1.
2.
3.
4. Зафиксируем в скалярное произведение , например,
и положим
Формула (1.24) задает
так называемую норму вектора
каковы
бы ни были векторы Норма, связанная со скалярным произведением по формуле (1.24), называется нормой, порожденной данным скалярным произведением. В конечномерном пространстве все нормы эквивалентны и поэтому можно пользоваться скалярным произведением (1.23) и порождаемой им евклидовой нормой
Пусть A: - линейный оператор. Нормой оператора A называется число
Заметим, что норма A оператора A в конечномерном пространстве конечна. Действительно, достаточно проверить этот факт в том случае, когда норма в задается формулой (1.26). Тогда
Отсюда
видим, что функция
поскольку S - компактное множество. Вернемся к формуле (1.28). Из неравенства
следует, что
Из неравенства (1.29) получаем, что
С другой стороны, обозначим через j0 такой номер, что
и положим
Согласно формуле (1.28)
Кроме
того, (так как
Таким образом, из (1.30), (1.32), (1.33) вытекает двухсторонняя оценка нормы оператора A в терминах его матрицы:
Множество
Метрическим пространством называют множество X произвольной природы, каждой паре элементов x, y которого поставлено в соответствие число (x,y) (называемое расстоянием между точками x и y), обладающее свойствами:
Данные соотношения называют аксиомами метрики, а свойство 3) - неравенством треугольника. Последовательность точек xn, n = 1,2, метрического пространства X называется фундаментальной последовательностью (илипоследовательностью Коши,) если для любого > 0 существует номер N() такой, что (xn,xm) < , для любых n, m N. Последовательность xn X, n = 1,2, называется сходящейся к точке x X, если для любого > 0, существует номер N() такой, что(xn,x) < , для любого n N. Метрическое пространство называется полным, если всякая фундаментальная последовательность его точек сходится к некоторой точке этого пространства.
На
множестве линейных операторов
Свойства 1) и 2) метрики очевидным образом следуют из определения (1.27) нормы оператора. Доказательство неравенства треугольника равносильно проверке неравенства
для произвольных операторов A, B . Докажем это неравенство. Пусть
Мы имеем:
Таким образом, для каждого справедливо неравенство
Отсюда, в силу (1.27), следует неравенство (1.36). Теорема 1.1.1 Пространство линейных операторов с метрикой (1.35) является полным метрическим пространством. Доказательство теоремы. Пусть Ak, k = 1,2, - фундаментальная последовательность точек пространства , так что для каждого > 0 найдется номер N() > 0 такой, что
Пусть
и
Таким
образом,
Пространство полное. Поэтому существует предел
Переходя в (1.37) к пределу при m , получаем, что
причем N()
- то же, что и (1.37) число,
не зависящее от
.
Точка
В силу (1.39) мы имеем, что
Таким образом,
Иначе говоря, Ak A, т.е. пространство полное. Пусть дано метрическое пространство X с метрикой и пусть при этом множество X является линейным пространством, а метрика удовлетворяет условиям:
Точки пространства X называют векторами. Функция ( ,0) называется нормой вектора и обозначается через . Нетрудно убедиться, что так определенная норма удовлетворяет свойствам (1.25) и при этом справедливо равенство
Такое пространство X называется нормированным, а функция называется нормой в X. Описанное выше пространство операторов является, таким образом, полным нормированным пространством; норма в нем задается формулой (1.27). Пусть X - полное нормированное пространство. Поскольку элементы X можно складывать и умножать на числа, то имеет смысл рассматривать ряды вида
Теория таких рядов буквально повторяет теорию числовых рядов. Теория функциональных рядов также может быть обобщена на функции со значениями в линейном пространстве X. Говоря точнее, пусть fi, i = 1,2 - функции, заданные на некотором множестве (вообще говоря, произвольной природы) T и принимающие значения в полном нормированном пространстве X:
Функциональный
ряд
Другими словами, функциональный ряд сходится в точке t T к функции f(t), если для произвольного > 0 найдется число N =N(,t) (зависящее от и, вообще говоря, от t) такое, что
В силу полноты пространства X имеет место критерий Коши сходимости ряда - для каждого > 0 cуществует номер N(,t) такой, что:
Ряд называют равномерно сходящимся на множестве T, если номер N в (1.43) (или в (1.44), что одно и то же) может быть выбран не зависящим от точек t T.
Ряд
называют абсолютно
сходящимся,
если сходится числовой ряд Признак Вейерштрасса. Если ряд из функций fi: T X мажорируется сходящимся числовым рядом:
то он сходится абсолютно и равномерно на T. Пусть f: T X - функция, определенная на некотором открытом множестве T вещественной прямой R со значениями в полном нормированном пространстве X. Функция f называется дифференцируемой в точке t = t0 T, если в X существует предел
то есть
Элемент
Дифференцирование
ряда.
Если ряд
функций fi: R X сходится
и ряд из производных Для случая X = R доказательство обоих теорем имеются в курсе анализа. На описанный выше случай эти доказательства переносится дословно. Определение экспоненты eA Пусть A: - линейный оператор. Экспонентой eA оператора A называется линейный оператор из в , определяемый формулой
где E - тождественный оператор в . Теорема 1.1.2 Ряд eA сходится в пространстве при любом A равномерно на каждом множестве
Доказательство. Пусть
оператор A
такой,
что A a.
Замечая, что
мажорирует ряд (1.45). Поскольку ряд (1.46) сходится, то по признаку Вейерштрасса ряд (1.45) равномерно сходится при A Sa. Упражнение. Вычислить матрицу eAt, если матрица A имеет вид:
Свойства экспоненты
Пусть T:
-
линейный оператор. Семейство
Если t может
изменяться на всей числовой
прямой < t <
+,
то Теорема 1.1.3 Семейство линейных операторов etA: , t R, является однопараметрической группой линейных операторов в . Доказательство теоремы. Мы уже знаем, что etA - линейный оператор, действующий в . Проверим, что имеют место свойства:
Обозначим через Sm(A) частичную сумму для ряда eA:
Заметим, что Sm -многочлен с неотрицательными коэффициентами относительно A. Мы должны доказать, что разность
стремится к нулю при m по норме пространства . Заметим, что m - это многочлен с неотрицательными коэффициентами относительно sA и tA. Действительно, преобразуем выражение
Так как
то предыдущее выражение принимает вид
Отсюда видно, что коэффициенты этого многочлена неотрицательны и, следовательно,
Положим
Предыдущие выкладки показывают, что выражение в правой части неравенства (1.50) равно разности
которая
при любых значениях и стремится
к нулю в силу формулы:
Формула (1.47) доказана. Для доказательства формулы (1.48) продифференцируем ряд
формально по t. Получим ряд из производных
Этот ряд сходится абсолютно и равномерно в любой области A a, t T и поэтому производная суммы ряда существует и равна сумме ряда из производных, т.е. AetA. Теорема доказана. |
