Линейные дифференциальные уравнения

Содержание



Введение

1. Линейные системы

1.1 Предварительные определения и обозначения

1.2 Линейные однородные системы

1.3 Неоднородные линейные системы

1.4 Линейные системы с постоянными коэффициентами

1.5 Линейные системы с периодическими коэффициентами

2. Линейные дифференциальные уравнения

2.1 Линейные дифференциальные уравнения порядка n

2.2 Линейные уравнения с аналитическими коэффициентами

2.3 Асимптотическое поведение решений некоторых линейных систем

3. Решение задач

Заключение

Список литературы




Введение



При изучении физических явлений часто не удается непосредственно найти законы, связывающие физические величины, но сравнительно легко устанавливается зависимость между теми же величинами, их производными или дифференциалами. 

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

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

Решением дифференциального уравнения называется n раз дифференцируемая функция , удовлетворяющая уравнению во всех точках своей области определения. 

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

Задача нахождения решения обыкновенного дифференциального уравнения или системы обыкновенных дифференциальных уравнений, удовлетворяющего некоторым начальным условиям, называется задачей Коши.

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

Цель дипломной работы – изучить понятие линейных дифференциальных уравнений.

В связи с поставленной целью необходимо выполнить следующие задачи:

1) Рассмотреть понятие линейных систем;

2) Изучить линейные дифференциальные уравнения различных порядков, в том числе с аналитическими коэффициентами;

3) Решить предложенные практические задания.




1. Линейные системы



1.1 Предварительные определения и обозначения



Пусть А = (aij) – квадратная матрица порядка n, где aij – комплексные числа. Определим норму А следующим образом:



. 									(1.1)



Если n-мерный вектор х представлять как матрицу с n строками и одним столбцом, то норма вектора совпадает с нормой x, определенной по формуле (1). Легко видеть, что норма обладает следующими свойствами:



(I) |A+B|  |A|+|B|,

(II) |AB|  |A|*|B|,

(III) |Ax|  |A|*|x|,



где А и В – матрицы, х – n-мерный вектор.

По определению, расстояние между двумя матрицами А и В равно |A-B|, и это расстояние удовлетворяет обычным свойствам метрики.

Нулевая матрица будет обозначаться через О, единичная – через Е. В случае опасности смешения размерностей эти квадратные матрицы порядка n будут обозначаться соответственно через Оn и Еn. 

Заметим, что | Оn | = 0 и | Еn | = n, а не 1.

Комплексно сопряженной матрицей для А = (aij) называется матрица , где  - комплексно сопряженные числа для aij. Транспонированная матрица обозначается через и определяется так: . Сопряженная матрица для А определяется так: . Заметим, что |A*|=||=||=|A|. Далее, (АВ)*=В*А*. Определитель матрицы А обозначается как det А.

Если det А = 0, то матрица А называется особой. Не особая матрица имеет обратную матрицу А-1, которая удовлетворяет соотношениям

А А-1 = А-1А = Е.

Многочлен det (λЕ-А) степени n от λ называется характеристическим многочленом для матрицы А, а его корни – характеристическими корнями А. Если эти корни обозначены λi (i = 1, …, n), то 



det (λЕ-А) = 



Две квадратные матрицы А и В порядка n называются подобными, если существует Неособая квадратная матрица Р порядка n, такая что



В = РАР-1.



Если А и В подобны, то они имеют один и тот же характеристический многочлен, ибо



det (λЕ-В) = det (Р(λЕ-А)Р-1)= det Р* det (λЕ-А)* det Р-1= det (λЕ-А).



В частности, коэффициенты многочлена det (λЕ-А) при степенях λ инвариантны относительно преобразования подобия. Два наиболее важных инварианта - det А и sp A – определитель и след А соответственно.

Приведем следующий фундаментальный результат о канонической форме матрицы.

Теорема 1.1 Каждая квадратная матрица А порядка n и подобная матрица вида





где J0 – диагональная матрица с элементами λ1, λ2,…, λq и 



 (i = 1, …, s).



Здесь λj , j = 1, …, q+s, - характеристические корни А, не обязательно различные. Если λj – простой корень, то он встречается в J0 и поэтому, если все корни различны, А подобна диагональной матрице







Из теоремы 1.1 непосредственно следует, что 



det А = , sp A = 



где произведение и сумма распространены на все корни, причем каждый корень считается столько раз, каков а его кратность. Матрицы Ji имеют вид

Ji = λq+iЕri+Zi ,



где Ji – квадратная матрица порядка ri и 







Матрицы Ji можно представить также в виде λq+iЕri+γZi, где γ – любая постоянная, отличная от нуля. 

Последовательность матриц {Аm} имеет своим пределом А, если для любого ε > 0 существует такое целое число N, что при p, q > N



|Aq - Ap| <ε.



Очевидно, что последовательность {Аm} сходится в том и только в том случае, когда сходится каждая из последовательностей компонент, а отсюда следует, что {Аm} сходится в том и только в том случае, когда существует предельная матрица, к которой и сходится эта последовательность.

Бесконечный ряд 







называется сходящимся, если сходится последовательность частных сумм, а суммой ряда называется предельная матрица для частных сумм. Важное значение при изучении линейных уравнений имеет специальный ряд, который называется экспонентной матрицей А, а именно:

 									(1.2)



где Аm есть m-я степень А. Ряд, определяющий еА, сходится для всех А, июо для любых положительных целых p и q







а последнее выражение есть разность Коши для ряда еА, сходящегося для всех конечных |A|. Далее,



 |еА|(n-1) + е|А|. 									(1.3)



Для матриц, вообще говоря, равенство еА+В = еА еВ неверно. Это равенство верно, если А и В коммутируют. Далее будет показано, что 



 det еА = еspА, 									(1.4)



и поэтому еА есть неособая матрица для всех А. Так как –А коммутирует с А, то е-А = (еА)-1.

Каждая матрица А удовлетворяет своему характеристическому уравнению det (λЕ-А) = 0, и это замечание часто бывает полезно для эффективного вычисления еА. 

Пусть В – неособая матрица. Покажем, что существует матрица А (называемая логарифмом В), такая, что еА = В. В самом деле, если в имеет каноническую форму J теоремы 1, то А, очевидно, можно представить в виде





при условии, что еАi = Jj, j = 0, 1, …, s. Легко также проверить, что А0 можно представить в виде







Далее, 







где Zj – нильпотентная матрица, определенная в теореме 1.1. так как высшие степени Zj равны нулю, то ряд 







содержит лишь конечное число членов и поэтому сходится. Положим, по определению, сумму этого ряда, который на самом деле является многочленом от, равной





Таким образом, 







есть многочлен от . С другой стороны. Из тождества





(|x| < 1)



следует после приведения справа подобных членов, коэффициенты при хk, k2, равны нулю, а коэффициент при х равен единице. Отсюда следует тот же результат для F, и поэтому 







Отсюда легко получаем, что Аj можно представить в виде







Пользуясь тем, что для каждой матрицы М



(PMP-1)k = PM k P-1 (k = 1, 2, …),



нетрудно видеть, что 





Отсюда следует, что результат, полученный для канонической матрицы В, переносится на произвольную неособую матрицу В. В самом деле, если J = eA и B = PJP-1, то В = , где = PАP-1. естественно, что матрица А не единственна.

Если Ф – произвольная квадратная матрица порядка n из функций, определенная на действительном i-интревале I (элементы матрицы могут быть действительными или комплексными функциями), то Ф называется непрерывной, дифференцируемой ли аналитической на I, если все элементы Ф соответственно непрерывны, дифференцируемы или аналитичны на I. Если Ф на I дифференцируема, то через обозначается произвольная матрица. Заметим, что если матрицы Ф, Ψ дифференцируемы, то 



  								(1.5)



и, вообще говоря, .

Если в точке t производная матрица  (t) существует и матрица Ф – неособая, то матрица Ф-1 в точке t дифференцируема. Это следует из равенства



  



где , а  - алгебраические дополнения элементов . Из равенств (1.5) и Ф Ф-1=Е следует, что 



   						(1.6)



Если матрица А на t-интервале I непрерывна и Ф удовлетворяет уравнению  (t) = А(t)Ф(t), то



  								(1.7)



а в интегральной форме



   					(1.8)



1.2 Линейные однородные системы



Пусть А – непрерывная квадратная матрица порядка n, элементами которой служат непрерывные комплексные функции, определенные на t-интервале I. Линейная система 



   (ЛО)



Называется линейной однородной системой порядка n. Для любого ξ и для τI существует единственное решение φ системы (ЛО) на интервале I, удовлетворяющее условию φ(τ) = ξ. Замечание: если каждый элемент матрицы А измерим на I и 



  , (*)



где m интегрируема по Лебегу на I, то существует единственное решение φ системы (ЛО), удовлетворяющее условию φ(τ) = ξ. В дальнейшем будем полагать, что для А выполняется по крайней мере условие (*).

Нулевая вектор-функция на I является решением системы (ЛО). Это решение называется тривиальным. Если решение системы (ЛО) равно нулю для некоторого , то в силу теоремы единственности оно равно нулю тождественно на I.

Теорема 2.1. Множество всех решений системы (ЛО) на интервале I образует n-мерное векторное пространство над полем комплексных чисел.

Доказательство. Если φ1 и φ2 – решения (ЛО) и с1 , с2 – комплексные числа, то с1φ1 + с2φ2 также является решением (ЛО). Это показывает, что решения образуют векторное пространство.

Чтобы доказать, что это пространство n-мерно, следует показать, что существует n линейно зависимых решений φ1 , φ2 , …, φn , таких, что каждое другое решение системы (ЛО) есть линейная комбинация (с комплексными коэффициентами) этих φi . Пусть ξi , i=1, 2, …, n – линейно независимые векторы n-мерного х-пространства. Например, за ξi можно взять вектор со всеми компонентами, равными нулю, кроме i-й, которая равна 1. Тогда, по теореме существования, если , то существуют решения φi, i=1, 2, …, n, системы (ЛО), для которых φi(τ) = ξi . Покажем, что эти решения удовлетворяют поставленному выше условию.

Если бы решения φi были линейно зависимы, то существовали бы n комплексных чисел , не равных одновременно нулю и таких, что



 .



Отсюда следует равенство 







противоречащее предположению о том, что векторы ξi линейно независимы.

Если φ – некоторое решение (ЛО) на I, такое, что φ(τ)=ξ , то можно найти (единственным образом определенные) постоянные сi , удовлетворяющие равенству



,



ибо векторы ξi образуют базис n-мерного х-пространства. Поэтому функция







есть решение (ЛО), принимающее при t = τ значение ξ, и, следовательно, в силу теоремы единственности 







Итак, каждое решение φ есть (единственная) линейная комбинация φi и теорема 2.1 полностью доказана.

Всякое множество φ1 , φ2 , …, φn линейно зависимых решений системы (ЛО) называется базисом или фундаментальным множеством решений системы (ЛО).

Если Ф – матрица, n столбцов которой являются n линейно независимыми решениями (ЛО) на I, то Ф называется фундаментальной матрицей системы (ЛО). Очевидно, Ф удовлетворяет матричному уравнению 



 . 							(2.1)



Под матричным дифференциальным уравнением, соответствующим системе (ЛО) на I, подразумевается задача отыскания квадратной матрицы Ф порядка n, столбцы которой являются решениями системы (ЛО) на I. Эта задача обозначается так:



 . 								(2.2)



Матрица Ф называется решением задачи (2.2) на I, и Ф удовлетворяет уравнению (2.1). Из теоремы 2.1. следует, что зная фундаментальную матрицу системы (ЛО), которая является, разумеется, частным решением уравнения (2.2), мы будем знать полную систему решений системы (ЛО).

 Теорема 2.2. Для того, чтобы решение-матрица уравнения (2.2) была фундаментальной матрицей, необходимо и достаточно, чтобы det Ф(t) 0 для .

Замечание. Если det Ф(t) 0 для некоторого , то в силу (1.8) det Ф(t) 0 для всех t.

Доказательство теоремы 2.2. Пусть Ф – фундаментальная матрица, столбцами которой являются векторы φj , и пусть φ – некоторое нетривиальное решение системы (ЛО). В силу теоремы 2.1 существуют единственным образом определенные постоянные с1 , с2 , …, сn , не равные все нулю и такие, что 







или, выражая при помощи матрицы Ф, 







где с – вектор-столбец с элементами с1 , с2 , …, сn . Это соотношение при каждом  есть система n линейных уравнений с неизвестными с1 , с2 , …, сn , имеющая единственное решение для каждого φ(τ). Поэтому det Ф(τ) 0 и, по сделанному выше замечанию, det Ф(t) 0 для каждого . Заметим, что это доказывает линейную независимость векторов-столбцов фундаментальной матрицы для каждого .

Наоборот, пусть Ф – матрица-решение уравнения (2.2) и пусть det Ф(t) 0 для каждого . Таким образом, векторы-столбцы матрицы Ф линейно независимы для каждого .

Матрица из векторов-столбцов может иметь определитель, тождественно равный нулю на интервале I, даже при линейно независимых векторах.

Например, пусть







для каждого действительного интервала I. Содержание теоремы 2.2 состоит в том, что этого не может случиться для векторов, которые являются решениями системы (ЛО).

Теорема 2.3. Если Ф – фундаментальная матрица для системы (ЛО) и С – (комплексная) постоянная неособая матрица, то ФС также является фундаментальной матрицей системы (ЛО). Каждая фундаментальная матрица системы (ЛО) может быть представлена в такой форме при помощи некоторой неособой матрицы С.

Доказательство. Из (2.1), если Ф – фундаментальная матрица, вытекает, что 



 ,


Или







и, следовательно, ФС есть матрица-решение уравнения (2.2). Так как



det (ФС)=( det Ф)( det С) 0,



то ФС – фундаментальная матрица.

Наоборот, если Ф1 и Ф2 – две фундаментальные матрицы , то Ф2 = Ф1С, где С – некоторая постоянная неособая матрица. Для доказательства этого положим Ф1-1Ф2 = Ψ(t). Тогда Ф2 = Ф1Ψ. Дифференцируя это равенство, получим, что . Отсюда и из (2.1) следует, что , или . Поэтому  и, следовательно, матрица ψ = С постоянна. Она неособая, так как этим свойством обладают Ф1 и Ф2.

Замечания. Если предполагать, что Ф2 – решение, то матрица С может быть особой.

Заметим, что если Ф – фундаментальная матрица системы (ЛО) и С – постоянная неособая матрица, то СФ, вообще говоря, не является фундаментальной матрицей. 

Две различные однородные системы не могут иметь одну и ту же фундаментальную матрицу, ибо из уравнения (ЛО) следует, что  Поэтому матрица Ф определяет матрицу А однозначно, хотя обратное утверждение и неверно.

Сопряженные системы. Если Ф - фундаментальная матрица для системы (ЛО), то 





или, переходя к сопряженным матрицам,







Поэтому - фундаментальная матрица для сопряженной системы (ЛО) и матричное уравнение 



  .								(2.3)



Система (2.3) называется сопряженной для системы (ЛО) и матричное уравнение 



   							(2.4)



называется сопряженным для уравнения (2.2.). Это соответствие симметрично, ибо (ЛО) и (2.2) сопряжены (2.3) и (2.4) соответственно.

Теорема 2.4. Если Ф - фундаментальная матрица для системы (ЛО), то Ψ есть фундаментальная матрица для сопряженной системы (2.3) в том и только в том случае, когда



 										 (2.5)



где С – постоянная неособая матрица.

Доказательство. Если Ф - фундаментальная матрица для системы (ЛО) и Ψ есть фундаментальная матрица системы (2.3), то так как - фундаментальная матрица частного вида уравнения (2.3), 





где D – некоторая постоянная неособая матрица (теорема 2.3). Поэтому







и можно принять С = D*.

Наоборот, если Ф - фундаментальная матрица для сопряженной системы (ЛО) и удовлетворяет (2.5), то  или  и, следовательно, в силу теоремы 2.3 Ψ - фундаментальная матрица сопряженной системы (2.3).

Если А = - А*, то , будучи фундаментальной матрицей для системы (2.3), является также фундаментальной матрицей для системы (ЛО). Поэтому в силу теоремы 2.3  или 



  										(2.6)



где С – постоянная неособая матрица. Из уравнения (2.6), в частности, следует, что эвклидова длина каждого вектора-решения системы (ЛО) постоянна.

Понижение порядка однородной системы. Если известно m (0<m<n) линейно зависимых решений системы (ЛО), то можно понизить порядок системы на m единиц, и следовательно, дело сведется к решению линейной системы порядка n-m.

Предположим, что φ1 , φ2 , …, φm - m линейно независимых векторов, которые являются решениями системы (ЛО) на интервале I. Пусть вектор φj имеет компоненты φij (i = 1, …, n). Тогда ранг прямоугольной матрицы с элементами φij (i = 1, …, n; j = 1, …, m) для каждого  равен m, так как ее столбцы линейно независимы. Это означает, что для каждого  в этой матрице найдется отличный от нуля определитель порядка m. Выберем некоторое  и предположим для определенности, что в точке t0 отличен от нуля определитель матрицы Фm с элементами φij (i = 1, …, m; j = 1, …, m). Тогда в силу непрерывной зависимости определителя от его элементов φij и непрерывности функций φij в окрестности t0 получим, что det Фm(t) 0 для t из некоторого интервала , содержащего t0. Пусть  - один из таких интервалов; процесс понижения проведем для . (Идея этого процесса является модификацией метода вариации произвольных постоянных.)

Пусть матрица U имеет своими первыми m столбцами векторы φ1 , φ2 , …, φm и своими n-m столбцами – векторы еm+1, …, en, где ej – вектор-столбец со всеми нулевыми элементами, кроме j-го, который равен 1. Очевидно, что U – неособая матрица на . Сделаем в (ЛО) подстановку



 x = Uy. 										(2.7)



Заметим, что решениям х = φj (j = 1, …, m) при подстановке (2.7) соответствуют решения y = ej (j = 1, …, m). Поэтому подстановку (2.7) можно рассматривать как систему относительно y, которая должна иметь решения ej (j = 1, …, m).

Подставляя (2.7) в систему (ЛО), получаем







или в координатах,



 (i = 1, …, m),

 (i = m+1, …, n).



Выражая то обстоятельство, что векторы φj с компонентами φij являются решениями системы (ЛО), получаем



 (i = 1, …, n; j = 1, …, m),



откуда следует, что 



 (i = 1, …, m),

  (i = m+1, …, n). 				(2.8)



Так как det Фm 0, то из первых m уравнений (2.8) можно выразить производные  (i = 1, …, m) через φij , aik и yk (k = m+1, …, n), и затем эти значения  подставить в остальные формулы (2.8). Мы получим совокупность уравнений первого порядка, которым удовлетворяют функции yi (i = m+1, …, n) вида



  (i = m+1, …, n), 					(2.9)



т.е. линейную систему порядка n-m.

Рассуждая в обратном порядке, предположим, что , …,  ( имеет компоненты  (i, j = m+1, …, n)) есть фундаментальная система решений на  для системы (2.9). Пусть  - матрица с элементами  (i, j = m+1, …, n). Очевидно, что det  0 на . Для каждого j = m+1, …, n пусть  (i = 1, …, m) определяется с помощью квадратур уравнений 

  						(2.10)



(i = 1, …, m; p = m+1, …, n).

Пусть  (p = m+1, …, n) обозначает с компонентами (i = 1, …, n) и пусть (p = 1, …, m). Так как  (p = m+1, …, n) удовлетворяют системе (2.9) и первым m уравнениям (2.8), то они должны также удовлетворять остальным n-m уравнениям (2.8), и поэтому  (p = m+1, …, n) являются решениями системы (2.8). таким образом, если Ψ – матрица со столбцами  (p = m+1, …, n) и 



Ф=U Ψ,



то Ф есть матрица-решение (ЛО) на I. U – неособая матрица.

Так как det =det  на , то Ф есть неособая матрица на  и, следовательно, является фундаментальным решением для системы (ЛО) на I.

Теорема 2.5. Пусть φ1 , φ2 , …, φm (m < n) – m известных линейно независимых решений системы (ЛО), причем φj (j = 1, …, m) имеют компоненты φij (i = 1, …, n). Предположим, что определитель матрицы с элементами φij (i, j = 1, …, m) на некотором подинтервале интервала i не обращается в нуль. Тогда с помощью подстановки (2.7) задачу определения n линейно независимых решений системы (ЛО) на  можно свети к решению системы (2.9)порядка n-m и к квадратурам (2.10).

Избавимся теперь от ограничения, что матрица Фm неособая на некотором интервале. Ясно, что прямоугольная матрица с элементами φij (i = 1, …, n; j = 1, …, m) имеет ранг m в силу независимости решений φj (j = 1, …, m). Таким образом, для каждого t = t0 существует неособая квадратная матрица порядка m, которую мы получим, выбирая m строк i1, …, im из прямоугольной матрицы с n строками и m столбцами. В силу непрерывности эта матрица неособая на некотором интервале . 

Хорошо известно и легко доказывается, что существует такая постоянная неособая матрица Т, применив которую к любому вектору х с n компонентами, получим матрицу Тх, имеющую своим m первыми компонентами компоненты вектора х с номерами i1, …, im. Полагая =Тх, мы заменим (ЛО) аналогичной системой, для которой выполняется первоначальное ограничение. Так как х=Т-1, то утверждение для х следует из доказанного уже утверждения для .



1.3 Неоднородные линейные системы



Пусть А – неособая квадратная матрица порядка n из непрерывных функций, определенных на действительном t-интервале I и b – непрерывный вектор на I, не равный тождественно нулю. Система уравнений 



 +b(t)  (ЛН)



называется линейной неоднородной системой порядка n. Если элементы А и b непрерывны и даже измеримы и мажорируются суммируемой функцией на I, то существует единственное решение φ системы (ЛН), для которого



φ(τ) = ξ ,



где и | ξ | < . Единственность решения следует из того, что если бы существовало два решения φ1 и φ2, то из разность φ = φ1 - φ2 была бы решением однородной системы (ЛО) на I при φ(τ) = 0. Но, по теореме единственности для (ЛО), разность φ должна равняться на I нулю тождественно и, следовательно, φ1 = φ2.

Если известна фундаментальная матрица Ф системы (ЛО), то легко найти решение системы (ЛН).

Теорема 3.1. Если Ф - фундаментальная матрица для системы (ЛО), то функция 



   						(3.1)



есть решение системы (ЛН), удовлетворяющее начальному условию 

φ(τ) = 0 ().

Доказательство получается непосредственно при помощи прямой проверки.

Интуитивные соображения, с помощью которых можно получить выражение (3.1), заключаются в следующем6 для каждого постоянного вектора с функция Фс является решением системы (ЛО). Метод состоит в том, что мы рассматриваем с как функцию, определенную на I, и находим, какой должна быть с (если она существует), для того, чтобы функция φ = Фс была решением неоднородной системы (ЛН).

Пусть φ = Фс – решение системы (ЛН). Тогда



= +Ф=АФс + Ф= А φ + Ф= А φ + b,



где последнее равенство следует из (ЛН). Поэтому Ф= b, или 



=b.



Последнее уравнение всегда разрешимо, причем если с(τ) = 0, то 


 .



Итак, φ определяется по формуле (3.1).

Легко видеть, что в условиях теоремы 3.1 решение системы (ЛН), удовлетворяющее условию φ(τ) = ξ (и | ξ | < ), дается в виде



  ,				 (3.2)



где  - решение системы (ЛО), удовлетворяющее условию



φh(τ) = ξ .



Формула (3.1) (или (3.2)) называется формулой вариации постоянных для системы (ЛН).

Заметим, что формулу (3.1) можно записать в виде



 ,



где Ψ – фундаментальная матрица системы



,



сопряженной системе (ЛО). Другая форма записи формулы (3.1) такова:


,



однако здесь необходимо ограничение .



1.4 Линейные системы с постоянными коэффициентами



Пусть А – постоянная квадратная матрица порядка n и рассмотрим соответствующую однородную систему



. 										(4.1)



Если n = 1, то (4.1) имеет очевидное решение еtА, и решение, которое при t = τ равно ξ , имеет вид е(t-τ)Аξ . Оказывается, что решение имеет эту форму и в том случае, когда х, ξ являются векторами произвольной конечной размерности и А – квадратная матрица порядка n.

Теорема 4.1. Фундаментальная матрица Ф системы (4.1) дается формулой



Ф(t) = еtА (|t| < ), 								(4.2)



и решение φ системы (4.1), удовлетворяющее условию 



 φ(τ) = ξ (|τ | < , | ξ | < ),



имеет вид 



 φ(t) = е(t-τ)Аξ (|t| < ). 								(4.3)




Доказательство. Так как е(t-Δt)А = еtА еΔtА, то из определения производной легко получаем, что







Поэтому Ф(t) = еtА есть решение системы (4.1). Так как Ф(0) = Е, то из (1.8) следует, что det Ф(t) = еtspА . Итак, Ф – фундаментальная матрица. Теперь формула (4.3) очевидна.

Замечание. Заметим, что выражение  не обязано быть решением системы , если матрицы А(t) и не коммутируют. Они коммутируют, когда матрица А(t) либо постоянная, либо диагональная.

Интересно исследовать структуру фундаментальной матрицы (4.2). пусть J – каноническая форма матрицы А, указанная в теореме 1.1, и предположим, что Р – неособая постоянная матрица, такая, что АР = PJ.

Тогда 



  							(4.4)



и J имеет вид



  							(4.5)



где J0 – диагональная матрица с элементами λ1, λ2,…, λq и 

  (i = 1, …, s). 		(4.6)



Далее, 



  						(4.7)



и легкое вычисление показывает, что 



  						(4.8)



Так как , то . Таким образом, 



  				(4.9)

где  - квадратная матрица порядка ri (n = q + r1 + … + rs). Поэтому, если известна каноническая форма (4.5), (4.6) матрицы А, то фундаментальная матрица еtА системы (4.1) дается в явном виде формулой (4.4), в которой еtJ может быть вычислена из (4.7), (4.8), (4.9).

Другая фундаментальная матрица системы (4.1) такова:



 Ψ(t) = еtАP = P еtJ. 								(4.10)



Пусть матрица Р имеет своими столбцами векторы р1, …, рn. Столбцы матрицы Ψ. Которые мы обозначаем через ψ1 , ψ 2 , …, ψ n , образуют совокупность n линейно независимых решений системы (4.1) и из (4.10) и вида матрицы J получаем



 , , …, ,

, 

,

,

,

.



Так как АР = PJ, то векторы р1, …, рn удовлетворяют соотношениям



Ар1 = λ1р1,…, Арq = λqрq,

Арq+1 = λ q+1рq+1,

Арq+2 = рq+1 + λ q+1рq+2,

Арn-rs+1 = λ q+sр n-rs+1,

Арn-rs+2 = р n-rs+1+λ q+sр n-rs+2,

Арn = р n-1+λ q+sр n.



Решения ψj выражаются посредством независимых векторов р1, …, рn из предыдущей последовательности уравнений.

Формула вариации постоянных (3.1) в применении к неоднородной системе



 +b(t) , 								(4.11)



где А – постоянная матрица, дает для решения φ системы (4.11), удовлетворяющего условию φ(τ) = 0, , выражение



 .



Решение φ системы (4.11), удовлетворяющее условию φ(τ) = ξ , где , | ξ |< , имеет вид



 .



1.5 Линейные системы с периодическими коэффициентами



Рассмотрим линейную однородную систему



  , 							(5.1)




где А – матрица элементами которой служат непрерывные комплексные функции, и 



  									(5.2)



для некоторой постоянной ω 0. В этом случае (5.1) называется периодической системой с ω-периодом А. Основной результат для таких систем состоит в том, что фундаментальную матрицу можно представить как произведение периодической матрицы с тем же периодом ω и матрицы-решения для системы с периодическими коэффициентами.

Теорема 5.1. Если Ф – фундаментальная матрица для системы (5.1), то тем же свойством обладает матрица



Ψ(t) = Ф(t+ω) .



Каждой такой матрице Ф соответствует периодическая неособая матрица Р с периодом ω и постоянная матрица R, такие, что 



 Ф(t) = P(t)etR. 									(5.3)



Доказательство. Так как



  , 



то в силу (5.2)



 .



Поэтому Ψ есть матрица-решение системы (5.1), и эта матрица фундаментальная, так как det Ψ(t) = det Ф(t+ω) 0 для .

Следовательно, существует постоянная матрица С, такая что



 Ф(t+ω) = Ф(t)С, 									(5.4)



и, сверх этого, существует постоянная матрица R, такая что



 С = еωR. 										(5.5)



Из (5.4) и (5.5) получаем



 Ф(t+ω) = Ф(t) еωR. 								(5.6)



Определим матрицу Р по формуле



Р(t) = Ф(t) е-tR. 									(5.7)



Тогда, используя (5.6), получаем



Р(t+ω) = Ф(t+ω) е-(t+ω)R = Ф(t) еωR е-(t+ω)R = Ф(t) е-tR = Р(t).



Так как матрицы Ф(t) и е-tR для  неособые, то Р(t) такая же, и это завершает доказательство.

Значение теоремы 5.1 состоит в том, что значение фундаментальной матрицы Ф на интервале длины ω, например , дает возможность определить Ф на всей числовой прямой. В самом деле, матрица С в (5.5) определяется как формула Ф-1(0)Ф(ω), а отсюда R определяется как (lnC)/ω. Теперь матрица Р(t) определена на интервале (0,ω) по формуле (5.7), а так как Р(t) имеет период ω, то она определяется на интервале . Теперь матрица Ф определена на интервале  по формуле (5.3).

Если Ф1 – некоторая другая фундаментальная матрица системы (5.1), для которой выполняется (5.2), то 



Ф = Ф1Т,



где Т – некоторая постоянная неособая матрица. Из (5.6) следует, что 



 Ф1(t+ω)Т = Ф1(t)ТеωR,



или 



 Ф1(t+ω) = Ф1(t)(ТеωRТ-1). 							(5.8)



Таким образом, в силу (5.8) каждая фундаментальна матрица Ф1 определяет матрицу ТеωRТ-1, подобную еωR . Наоборот, если Т – любая постоянная неособая матрица, то существует фундаментальная матрица Ф1 системы (5.1), такая, что выполняется (5.8). Следовательно, хотя Ф не определяет R однозначно, множество всех фундаментальных матриц системы (5.1), а следовательно матрица А, определяет однозначно все связанные с R величины, инвариантные относительно подобных преобразований. В частности, множество всех фундаментальных матриц системы (5.1) определяет однозначно множество характеристических корней, а именно характеристические корни матрицы С = еωR. Обозначим эти корни через λ1, λ2,…, λn и назовем их мультипликаторами, соответствующими матрице А. Ни один из мультипликаторов не равен нулю, ибо = det еωR 0. Характеристические корни матрицы R называются характеристическим показателями.

Интересно выяснить явный вид множества n линейно независимых векторов-решений системы (5.1). Пусть Т – постоянная неособая матрица, такая что матрица Т-1RТ = J имеет каноническую форму, указанную в теореме 1.1, и положим Ф1 = ФТ, Р1 = РТ. Тогда из (5.3) следует 



 Ф1(t) = Р1(t) еtJ , Р1(t+ω) = Р1(t).						 (5.9)



Поэтому, если ρi – характеристические корни R, то матрица еtJ имеет вид







где 







и 



 (i = 1, …, s; q+= n).



Очевидно, что λi = еωρi , и поэтому, хотя сами корни ρi определяются неоднозначно, но их действительные части определяются однозначно. Из (5.9) следует, что столбцы φ1 , φ2 , …, φn матрицы Ф1, которые образуют множество n линейно независимых решений системы (5.1), имеют вид:



,

 , 

,

, 

, 					(5.10)

,

,

.



В этих формулах р1 , р2 , …, рn - периодические векторы-столбцы матрицы Р1. 

Из (5.10) очевидно, что если Reρi < 0 или, что эквивалентно, | λi | < 1, то при  решения φi(t) экспоненциально убывают.

Из (5.6) следует, что Ф(ω) = Ф(0)еωR, и поэтому λi можно рассматривать как характеристические корни матрицы Ф-1(0)Ф(ω). В частности, если Ф(0) = Е, то еωR = Ф(ω) и λi являются характеристическими корнями матрицы Ф(ω). Так как



  						(5.11)



то n-й корень можно определить из (5.11), если известны n-1 корней λi.

Действительная неособая матрица С не обязана иметь действительный логарифм, т.е. не всегда существует действительная матрица В, такая что еВ = С. В самом деле, матрица с одной строкой и одним столбцом С = -1 доставляет соответствующий пример. Однако справедливо утверждение, что для действительной матрицы С существует действительная матрица В, такая, что С2 = еВ .

Используя это при доказательстве теоремы 5.1, нетрудно получить следующий результат6 если в системе (5.1) матрица А (t) действительная периодическая с периодом ω, то каждой действительной фундаментальной матрице Ф соответствует действительная матрица Р периода 2ω и действительная постоянная матрица R, такие, что 



Ф(t) = Р(t)еtR.




2. Линейные дифференциальные уравнения



2.1 Линейные дифференциальные уравнения порядка n



Предположим, что n+1 коэффициентов а0, а1, …, аn представляют собой непрерывные (комплексные) функции, определенные на действительном t-интервале I, и пусть Ln обозначает формальный дифференциальный оператор



;



это означает, что если функция g имеет n производных на I, то 



 



Далее предположим, что а0(t) 0 для . Тогда, по определению, уравнение



  



(в подробной записи  ()) есть дифференциальное уравнение



  ;



оно называется линейным однородным дифференциальным уравнением порядка n. Соответствующая этому уравнению система есть векторное уравнение 




  (6.1)



где 



  				(6.2) 



Так как (6.1) – линейная система с непрерывной на I матрицей коэффициентов А, то существует единственный вектор-решение φ на I системы (6.1), удовлетворяющий условию



 



где , |ξ| < . Таким образом, φ1 – первая компонента  - удовлетворяет условиям



  …,  					(6.3)



Так как φ1 – решение уравнения , то это решение удовлетворяет условиям (6.3).

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

Если φ1, …, φn – n решений уравнения , то матрица

  						(6.4)



есть матрица-решение для (6.1). Определитель этой матрицы называется вронскианом уравнения , соответствующим решениям φ1, …, φn, и обозначается через W(φ1, …, φn). При фиксированных φ1, …, φn он является функцией t на I и его значение в точке t обозначается W(φ1, …, φn)(t). Из того, что для линейной системы вида (6.1)



  ,



следует (замечая из (6.2), что spA = - а1/а0)



 W(φ1, …, φn)(t) = W(φ1, …, φn)(τ) . 		(6.5)



Теорема 6.1. Необходимое и достаточное условие того, чтобы n решений φ1, …, φn уравнения  на интервале I были линейно независимы, заключается в том, что 



W(φ1, …, φn)  0 .



Каждое решение уравнения  есть линейная комбинация с комплексными коэффициентами любых n линейно независимых решений.

Доказательство. Если решения φ1, …, φn на I линейно зависимы. То существуют постоянные с1, …, сn, не все равные нулю, такие, что





Отсюда следует тождество



 (k = 0, 1, …, n-1),



и поэтому векторы  с компонентами  (i = 1, 2, …, n) линейно зависимы на I. Наоборот, если векторы  линейно зависимы, то тем же свойством обладают решения φ1, …, φn уравнения. Из теоремы 2.2 следует, что необходимое и достаточное условие линейной независимости векторов  заключается в том, что det Ф(t) 0 на I, где Ф – матрица (6.4). Но это требование в точности совпадает с условием W(φ1, …, φn)(t) 0 на I. В силу (6.5), если W(φ1, …, φn)(τ) 0 для некоторого , то W(φ1, …, φn)(t) 0 для любого .

Так как каждый вектор-решение задачи (6.1), (6.2) есть линейная комбинация n линейно независимых векторов-решений, то каждое решение уравнения  есть линейная комбинация n линейно независимых решений этого уравнения. Это доказывает теорему.

Ввиду свойства, указанного в теореме 6.1, множество n линейно независимых решений уравнения  называется базисом или фундаментальным множеством для уравнения .

Теорема 6.2. Пусть φ1, …, φn - n функций, имеющих на интервале I непрерывные производные порядка n, и пусть W(φ1, …, φn)(t) 0 на I. Тогда существует единственное однородное дифференциальное уравнение порядка n (с коэффициентом при х(n), равным единице), для которого эти функции образуют фундаментальное множество, а именно:

  							(6.6)



Замечание. Вронскиан W(х, φ1, …, φn) представляет собой определитель матрицы, первая строка которой состоит из элементов х, φ1, …, φn, а последующие строки являются последовательными производными первой строки до порядка n включительно.

Доказательство. Очевидно, что W(φi, φ1, …, φn) = 0 (i = 1, …, n), ибо в этом определителе имеются два одинаковых столбца. Разложение числителя W(х, φ1, …, φn) левой части уравнения (6.6) по элементам первого столбца показывает, что (6.6) есть дифференциальное уравнение порядка n. Коэффициент при х(n) в W(х, φ1, …, φn) равен (-1)n W(φ1, …, φn), т.е. в (6.6) при х(n) равен единице. Так как W(φ1, …, φn) 0, то из теоремы 6.1 следует, что φ1, …, φn образуют для (6.6) фундаментальное множество.

Единственность уравнения (6.6) следует из того, что соответствующие векторы  с компонентами  определяют матрицу коэффициентов (6.2), соответствующие системе (6.1), однозначно. Так имеется взаимно однозначное соответствие между линейными уравнениями порядка n и линейными системами вида (6.1), (6.2), то доказательство завершено.

Если одно или более из решений уравнения  известно, то использование соответствующей системы (6.1) позволяет понизить порядок уравнения. Более прямо достигает цели следующий процесс, который является методом вариации постоянных применительно к уравнению . Пусть  и положим х = уφ1. Тогда уравнение  является для у линейным дифференциальным уравнением порядка n, которое имеет решение у = 1, ибо φ1 есть решение уравнения . Поэтому в новом уравнении коэффициент при у должен обращаться в нуль. Рассматривая это уравнение относительно переменной u = , получим уравнение (n-1)-го порядка. Если φ2 не зависит от φ1 и , то  есть решение (n-1)-го порядка для u, которое может быть аналогично сведено к уравнению (n-2)-го порядка, и т.д.

Сопряженные уравнения. С формальным оператором Ln тесно связан другой линейный оператор  порядка n, называемый сопряженным для Ln и определяемый следующим образом:







Иначе говоря, если g – функция на I, для которой произведение  (k = 0, 1, …, n) имеет n – k производных на I, то







Уравнение 



 



(в подробной записи 



),



называемое сопряженным для Lnх = 0 на I, определяется как задача отыскания функции φ (решения), на I, такой, что произведение  (k = 0, 1, …, n) имеет n – k производных на I,удовлетворяющей на I уравнению





Если  на I и φ – решение уравнения, имеющее n производных на I, то, используя правило дифференцирования произведения, получаем







разделив на , видим, что φ есть решение дифференциального уравнения порядка n рассмотренного выше типа.

Рассмотрим тот специальный случай оператора Lnх когда а0 = 1. Для системы (6.1), (6.2), ассоциированной с уравнением 



  						(6.7)



сопряженная система имеет вид



  , 								(6.8)



где в силу (6.2)



 						 (6.9)



Записывая (6.8) в компонентах, получаем в силу (6.9)



   (k = 2, …, n). 				(6.10)

Таким образом, если φ1, …, φn – решение системы (6.10), для которого  и 







существуют, то дифференцируя k-е равенство (6.10) k-1 раз и решая относительно, получаем







Поэтому φn удовлетворяет уравнению







которое является сопряженным к (6.7).

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

Теорема 6.3. (Тождество Лагранжа). Предположим, что в Ln  на I (k = 0, 1, …, n). Если u, v – произвольные (комплексные) функции на I, имеющие n производных, то 



  (), 							(6.11)



где [uv] – форма относительно величин (u, , …,  ) и (v, , …, ), задаваемая равенством

  					(6.12)



Доказательство. Пользуясь правилом дифференциального произведения, имеем для m = 0, 1, …, n







Таким образом, получаем









что доказывает формулу (6.11).

Следствие (Формула Грина). Если ак в Ln и u, v такие же как и в теореме 6.3, то для любых 



 					 (6.13)



где [u, v](ti) – значение [u, v] при t = ti.

Доказательство. Следует проинтегрировать тождество (6.11) в пределах от t1 до t2 .

Если ψ – известное решение уравнения  на I, то в силу (6.11) отыскание решения Lnх = 0 сводится к отысканию функции φ, удовлетворяющей уравнению (n-1)-го порядка





Неоднородное линейное уравнение порядка n. Предположим, что на действительном t-интервале I а0  0, а1, …, аn и b – непрерывные функции, и рассмотрим уравнение



,



которое совпадает с уравнением







Это уравнение (в случае b(t)  0) называется неоднородным линейным уравнением порядка n. Соответствующая этому уравнению система имеет вид



  , 							(6.14)



где А – матрица (6.2) и - вектор-столбец со всеми нулевыми элементами, кроме последнего, который равен b/а0. Таким образом, соответствующая уравнению Lnх = b(t) система (6.14) есть линейная неоднородная система; существование и единственность решения для системы (6.14) обеспечивает, как обычно, существование и единственность решения для уравнения Lnх = b(t).

Теорема 6.4. Если φ1, …, φn – фундаментальное множество для однородного уравнения 

 ( на I),



то решение ψ неоднородного уравнения 



Lnх = b(t) ( на I),



удовлетворяющее условию



  (, |ξ| < ),



имеет вид



  			(6.15)



где - решение уравнения Lnх = 0, для которого , и Wk(φ1, …, φn) –определитель, получаемый из W(φ1, …, φn) в результате замены k-го столбца на (0, …, 0, 1).

Доказательство. В силу (3.1) первая компонента ψ = ψ1 вектора-решения  системы (6.14), для которого = 0, имеет вид 







где  - элемент, находящийся на пересечении первой строки и n-го столбца матрицы Ф(t)Ф-1(s). Напомним, что на пересечении i-й строки и j-го столбца матрицы Ф(t) стоит элемент  и что det Ф(t) = W(φ1, …, φn)(t). Далее, на пересечении i-й строки и n-го столбца матрицы Ф-1 стоит элемент







где  - алгебраическое дополнение элемента  в Ф. Поэтому 







где Wk(φ1, …, φn)(s) определен в формуле теоремы. Таким образом, решение ψ уравнения Lnх = b(t) , удовлетворяющее условию = 0, имеет вид







и очевидно, что (6.15) дает решение, удовлетворяющее условию , если .

Линейное уравнение порядка n с постоянными коэффициентами. Рассмотрим тот случай, когда в Ln все коэффициенты а0 = 1, а1, …, аn – постоянные. В этом случае можно предполагать, что I есть вся числовая ось. Далее, уравнению



 						 (6.16)



соответствует система 



  										(6.17)

где А – постоянная матрица



  				(6.18)



Можно предполагать, что для (6.16) можно указать фундаментальное множество решений, и точный вид этих функций зависит от характеристического многочлена f(λ) = det (λE - A) постоянной матрицы А в (6.18).

Лемма. Характеристический многочлен для матрицы А в (6.18) имеет вид



 f(λ) = λn + а1 λn-1 +…+ аn. 							(6.19)



Заметим, что f(λ) может быть получено из Ln(х) формальной заменой х(k) на λk.

Доказательство проводится по индукции. Для n = 1 А = - а1; значит det(λE1 - A) = λ + а1 и, следовательно, (6.19) верно для n = 1. Предположим, что результат справедлив для n – 1. Разложим определитель



det(λEn - A) = 



по элементам первого столбца и заметим, что коэффициент при λ есть определитель (n-1)-го порядка, именно det(λEn-1 – A1), где







Поэтому λdet(λEn-1 – A1) = λn + а1 λn-1 +…+ аn--1λ. Единственный другой ненулевой элемент в первом столбце есть аn и его алгебраическое дополнение равно 1. Поэтому det(λEn – A) = λn + а1 λn-1 +…+ аn--1λ + аn, что и требовалось доказать.

Теорема 6.5 (без доказательства). Пусть λ1, …, λn - различные корни характеристического уравнения 



f(λ) = λn + а1 λn-1 +…+ аn = 0



и пусть кратность корня λi равна mi (i = 1, …, s ). Тогда фундаментальное множество для (6.16) дается n функциями



 tkeλi (k = 0, 1,…, mi – 1; i = 1, 2,…, s).					 (6.20)



2.2 Линейные уравнения с аналитическими коэффициентами



Предположим, что А – квадратная матрица порядка n и b – n-мерный вектор, определенные и аналитические в односвязной области D z-плоскости, и пусть . Используя метод последовательных приближений, нетрудно показать, что линейная система 



 								 (7.1)

при условии 



 



имеет в D единственное аналитическое решение .

В самом деле, пусть  и пусть С – дуга длины L, лежащая в D, соединяющая точки z0 и z1 и имеющая непрерывно вращающуюся касательную. Обозначим через s длину дуги вдоль С, начиная от точки z0. Выберем постоянную К настолько большой, чтобы было | A(z) | < K и || < K для . Пусть  и 







причем интеграл берется вдоль С, так что приближения  определены на С. Нетрудно получить оценки









Очевидно, эти оценки справедливы для всех точек z в D, достижимых из z0 другой длины L, на которой |A(z)| и || ограничены постоянной K. Отсюда следует, что эти оценки справедливы в каждом фиксированном множестве R, содержащемся в D. Так как каждая функция  аналитична в R, то из равномерной сходимости  следует, что предельная функция  также аналитична в R. Далее,





Это доказывает утверждение для R и, следовательно, для D.

Кроме того, все теоремы, доказанные в пп. 1.2 и 1.3, будучи существенно алгебраической природы, справедливы для системы (7.1).

Соответственно этому, если n+1 функций а1, …, аn, b аналитичны в D, то линейное уравнение порядка n 



  					(7.2)



имеет в D единственное решение, удовлетворяющее условиям 



, , …, ,



где w1, …, wn – n данных комплексных чисел. Наконец, все результаты п. 2.1 распространяются очевидным образом на случай (7.2). 



2.3 Асимптотическое поведение решений некоторых линейных систем



Если коэффициенты линейной системы дифференциальных уравнений при стремятся к постоянным, то иногда возможно охарактеризовать поведение решений.

Здесь рассматривается проблема для действительного переменного. Рассмотрим пример







где v – действительная дифференцируемая функция, для которой , r – интегрируемая функция и 



, 



для некоторого t0. (На самом деле достаточно, чтобы функция v имела в интервале  ограниченную вариацию.) Без ограничения общности можно в дальнейшем предполагать, что t0 = 0. Из доказанной ниже теоремы следует, что рассматриваемое уравнение имеет два решения φ и ψ, такие, что









при , а ψ имеет аналогичное поведение с заменой i на –i.

Этот результат показывает, что функция r нисколько не влияет на грубую асимптотику. Однако случай



 (0 < α < 1)



убеждает нас в том, что влияние v существенно. Эти асимптотические формулы показывают также, что если положить в уравнении функцию r(t) равной нулю, а 1 + v(t) – постоянной, то результат будет отличаться от точного только членом о(1) при .

В дальнейшем будет рассматриваться линейная система 



  							(8.1)

которая включает как частный случай предыдущий пример.

Теорема 8.1. Пусть А – постоянная матрица с различными характеристическими корнями μj, j = 1, …, n. Пусть матрица V дифференцируема и удовлетворяет условию



  									(8.2)



и пусть  при . Пусть матрица R интегрируема и 



 . 									(8.3)



Обозначим корни уравнения det (A + V(t) - λE) = 0 через λj(t), j = 1, …, n. Очевидно, что можно, если это необходимо, переставить μj так, чтобы . Для каждого k положим







Допустим, что все j, 1 j  n, попадают в один из двух классов I1 и I2, где



, если  при 



и 



  , 						(8.4)

 , если  ; 				(8.5) 



здесь k фиксировано и К – постоянная. Пусть pk – характеристический вектор А, соответствующий μk, так что



 Аpk = μk pk.										 (8.6)



Тогда существует решение φk системы (8.1) и число t0, 0  t0 , такие, что 

  						(8.7)



Доказательство. Если условия теоремы выполняются для всех k, 1  k  n, и Ф – матрица со столбцами φ1, …, φn, то Ф – фундаментальная матрица, так как det Ф(t) 0 для больших t, ибо pk линейно независимы.

Предположим в начале, что А + V(t) для t  t0 имеет диагональный вид А(t) причем t0 выбрано так, что 



  								(8.8)



Пусть Ψ(t) – диагональная матрица:



Ψ(t) = 



так что 


  									(8.9)



Пусть еК - вектор-столбец со всеми нулевыми элементами, за исключением k-го, который равен 1, и ψk – вектор, определенный равенством 







При фиксированном k и I1, I2, определенных согласно неравенствам (8.4), (8.5), положим



Ψ = Ψ1 + Ψ2,



где диагональные матрицы Ψ1 и Ψ2 содержат элементы Ψ, соответствующие столбцам с индексами j, принадлежащим соответственно I1 и I2. Тогда



  (j = 1, 2). 								(8.10)



Рассмотрим теперь уравнение





  							(8.11)



Можно непосредственно проверить, что если уравнение (8.11) имеет решение φ, то 




= (A + R) φ. 									(8.12)



Последнее уравнение имеет рассматриваемый нами вид (8.1)

Пусть φ0(t) = 0 и 	



  						(8.13)



Тогда φ1(t) = ψk(t) и для t  t0 



 						 (8.14)



Каждый элемент диагональной матрицы  имеет вид 



  



или равен нулю. Но для t0  τ t







Поэтому для t0  τ t





Точно также для τ  t получим







Используя эти неравенства, получаем из (8.13)





 



Из (8.8) и (8.14) теперь по индукции следует







Отсюда следует равномерная сходимость последовательности {φj} на каждом конечном подинтервале интервала [t0,). Так как φj непрерывно, то предельная функция φ также непрерывна и, очевидно, 



 							 (8.15)



Покажем теперь, что 



 					 (8.16)

Это будет установлено, если мы покажем, что при t→ ∞



  			(8.17)



И



 			 (8.18)



Доказательство соотношения (8.18) сразу получается из (8.15) и (8.5). Доказательство соотношения (8.17) основывается на равенстве



 						 (8.19)



которое является следствием (8.4). Каково бы ни было ε > 0, можно подобрать такое t1, что







Поэтому, обозначая левую часть (8.17) через J(t), получаем 







Из (8.19) следует, что 





Так как ε произвольно, то (8.17) доказано. Таким образом, теорема доказана для случая A + V(t) = A(t), если за φ взята φk.

Доказательство теоремы 8.1 вытекает из следующей леммы.

Лемма. Пусть A и V удовлетворяют условиям теоремы 8.1. Тогда существует матрица S(t), которая при t → ∞ стремится к постоянной неособой матрице Т, такая что



 S(A + V) = ΛS, 									(8.20)



где Λ(t) – диагональная матрица с диагональными элементами λj(t), j = 1, 2, …, n. При t → ∞ λj(t) → μj , где μj – характеристические корни матрицы А. Кроме того, для некоторого t0



  									(8.21)



Доказательство теоремы 8.1. Так как S(t) → T при t → ∞ и Т – неособая матрица, то S(t) – неособая матрица для всех достаточно больших t. Выберем t0 настолько большим, чтобы не только (8.21) выполнялось, но и S-1(t) существовала для tt0. Тогда, полагая в (8.1) у = S(t)х, получаем



  (tt0).						 (8.22)



Пусть = . Тогда из (8.3) и (8.21) следует, что норма  интегрируема. Таким образом, данное выше доказательство теоремы 8.1 для специального случая годится для уравнения (8.22), так что (8.22) имеет решение θk, для которого







Поэтому (8.1) имеет решение S-1 θk = φk. Так как S-1(t) → T-1, то Аpk = μkpk. Это завершает доказательство теоремы 8.1.




3. Решение задач



Задача 1. Пусть матрица А и вектор b – интегрируемые функции от t на интервале [a, c]. Пусть 



|A(t)|  k(t),

|b(t)|  k(t),



Пусть τ [a, c] и рассмотрим начальную задачу



, х(τ) = ξ.



Доказать, что существует решение φ на [a, c] в том смысле, что



 на [a, c].



Доказательство.

Используем последовательные приближения. Пусть φ0(t) = ξ и 



, j  0.



Тогда 


















Пусть  тогда 











Значит, 



 где 



Следовательно, последовательность {φj} сходится равномерно на [a, c].




Задача 2. Решить систему дифференциальных уравнений 2-го порядка методом Эйлера. 



Решение:



Система дифференциальных уравнений второго порядка имеет вид:







Приведем систему к системе дифференциальных уравнений первого порядка. Сделаем замену:







Пусть заданы начальные условия







и выбран шаг h по оси x.

Метод Эйлера для решения системы дифференциальных уравнений 2-го порядка в общем виде: 



,







где j – номер шага.


Заключение



В дипломной работе рассмотрены вопросы решения линейных дифференциальных уравнений и систем уравнений. 

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

Большое значение, которое имеют дифференциальные уравнения для математики и особенно для ее приложений, объясняются тем, что к решению таких уравнений сводится исследование многих физических и технических задач. 

Задача нахождения решения обыкновенного дифференциального уравнения или системы обыкновенных дифференциальных уравнений, удовлетворяющего некоторым начальным условиям, называется задачей Коши.

Решением будет функция, график которой касается каждой своей точкой соответствующего отрезка. Каждое отдельное решение называется частным решением дифференциального уравнения; если удается найти формулу, содержащую все частные решения (за исключением, быть может, нескольких особых), то говорят, что получено общее решение. Частное решение представляет собой одну функцию, в то время как общее – целое их семейство. Решить дифференциальное уравнение – это значит найти либо его частное, либо общее решение.

Линейные уравнения – это уравнения «первой степени» – неизвестная функция и ее производные входят в такие уравнения только в первой степени. Таким образом, линейное дифференциальное уравнение первого порядка имеет вид dy/dx + p(x) = q(x), где p(x) и q(x) – функции, зависящие только от x. Его решение всегда можно записать с помощью интегралов от известных функций. Многие другие типы дифференциальных уравнений первого порядка решаются с помощью специальных приемов. 

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

В работе также проведено решение конкретных заданий, связанных с нахождением решения дифференциальных уравнений и систем дифференциальных уравнений.

Таким образом, дифференциальные уравнения играют существенную роль и в других науках, таких, как биология, экономика и электротехника; в действительности, они возникают везде, где есть необходимость количественного (числового) описания явлений (коль скоро окружающий мир изменяется во времени, а условия изменяются от одного места к другому).




Список литературы



		Арнольд В.И. Обыкновенные дифференциальные уравнения. М.: Наука, 1966. – 384 с.

		Бабушка И., Витасек Э., Прагер М. Численные процессы решения дифференциальных уравнений. М.: Мир, 1969. – 428 с.

		Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. Приближение функций, дифференциальные и интегральные уравнения. М.: Наука, 1967. – 439 с.

		Дородницын А.А. Численные методы решения дифференциальных и интегральных уравнений и квадратурные формулы. Сборник статей. М.: Наука, 1964. – 386 с.

		Дьяченко В.Ф. Основные понятия вычислительной математики. М.: Наука, 1972. – 563 с.

		Ильин В.А., Позняк Э.Г. Основы математического анализа. – ч. 1. М.: Наука, 1973. – 591 с.

		Карташев А.П., Рождественский Б.Л. Обыкновенные дифференциальные уравнения и основы вариационного исчисления. М.: Наука, 1976. – 472 с.

		Коллатц Л. Функциональный анализ и вычислительная математика. М.: Мир, 1969. – 475 с.

		Марчук Г.И. Методы вычислительной математики. М.: Наука, 1977. – 623 с.

		Михлин С.Г., Смолицкий X.Л. Приближенные методы решения дифференциальных и интегральных уравнений. М.: Наука, 1965. – 352 с.

		Островский А.М. Решение уравнений и систем уравнений. М.: ИЛ, 1963. – 349 с.

		Степанов В.В. Курс дифференциальных уравнений. М.: Гостехиздат, 1963. – 461 с.

		Тихонов А.Н., Самарский А.А. Уравнения математической физики. М., 1977. – 522 с.

		Понтрягин Л.С. Обыкновенные дифференциальные уравнения. М., 1982. – 549 с.

		Полянин А. Д. Справочник по линейным уравнениям математической физики. М.: Физматлит, 2001. – 419 с.

		 Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М., 1984. – 463 с.

		Штеттер X. Анализ методов дискретизации для обыкновенных дифференциальных уравнений. М.: Мир, 1978. – 275 с.

		Эрроусмит Д., Плейс К. Обыкновенные дифференциальные уравнения. Качественная теория с приложениями. М., 1986. – 478 с.

		Интернет-источник: www.mathematics.ru.

		Интернет-источник: www.nsu.ru/matlab. 

		Интернет-источник: ru.wikipedia.org.

		Интернет-источник: www.matclub.ru.