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

Системы линейных уравнений с постоянными коэффициентами методы их решения План. Системы линейных уравнений с постоянными коэффициентами методы их решения. Автономные системы дифференциальных уравнений. Классификация точек покоя системы двух однородных линейных дифференциальных уравнений с постоянными коэффициентами.

2015-01-30

60.86 KB

8 чел.


Поделитесь работой в социальных сетях

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


лухов Ю.П. Конспект лекций по высшей математике. Лекция № 39                                     7

Лекция 39

ТЕМА: Дифференциальные уравнения.

Системы линейных уравнений с постоянными коэффициентами, методы их решения

План.

  1.  Системы линейных уравнений с постоянными коэффициентами, методы их решения.
  2.  Определение устойчивости по Ляпунову и асимптотической устойчивости.
  3.   Автономные системы дифференциальных уравнений.
  4.  Фазовое пространство (плоскость), фазовая траектория. Точки покоя. Классификация точек покоя системы двух однородных линейных дифференциальных уравнений с постоянными коэффициентами.
  5.  Условия устойчивости точки покоя.

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

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

В частности, система линейных уравнений первого порядка с постоянными коэффициентами имеет вид:

         .                                                                  (36.1)

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

.

Тогда системе (36.1) эквивалентно матричное уравнение

.                                                  (36.2)

Если же рассмотреть линейный оператор  , уравнение (36.2) примет вид:

.                                                           (36.3)

Так как оператор L обладает свойствами линейности:

  1.  L[cX] = cL[X];

  1.  L[X1 + X2] = L[X1] + L[X2],

то для решений линейной однородной системы (36.3) (при F = 0) справедливы те же свойства:

если Х1 и Х2 – решения однородного уравнения , то и их линейная комбинация будет решением того же уравнения.

Можно ввести понятие линейной зависимости решений Х1, Х2,…, Хп:

Определение. Векторы (столбцы) Х1, Х2,…, Хп , где

, называются линейно зависимыми при , если существуют числа α12,…, αп, не все равные нулю, что

α1Х1 + α2Х2 +…+ αпХп0                                                      (36.4)

при . Если же тождество (36.4) справедливо только при всех αi = 0, векторы называются линейно независимыми.

Замечание. Назовем определителем Вронского для уравнения (36.4) определитель вида

,                                                          (36.5)

являющийся определителем системы уравнений, получаемых при координатной записи равенства (36.4). Можно показать, что так же, как и в случае решения линейного однородного уравнения, при W = 0 решения Х1, Х2,…, Хп линейно зависимы на [a,b]. Тогда справедлива следующая теорема:

Теорема 42.1. Линейная комбинация  п линейно независимых решений линейной однородной системы является общим решением этой системы.

Будем искать фундаментальную систему решений линейной однородной системы с постоянными коэффициентами

                                             (36.6)

в виде:

,                                            (36.7)

где αi – постоянные. Подставив (36.7) в (36.6) и сократив на ekt, получим:

 .                                        (36.8)

Для того, чтобы эта система имела ненулевое решение, необходимо и достаточно, чтобы ее главный определитель был равен нулю:

,                                             (36.9)

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

Если все корни характеристического уравнения различны, то, подставляя их последовательно в систему (36.8), можно найти соответствующие им значения  и тем самым п различных решений системы (36.6). Эти решения линейно независимы. Действительно, если бы существовали числа β1, β2,…, βп такие, что

                        ,  то в силу линейной независимости функций отсюда следовало бы, что  для каждого i. Но поскольку хотя бы одно из  не равно нулю, получаем, что все . Следовательно, найденные решения (36.7) линейно независимы, и общее решение системы имеет вид: ,    (36.10)

где ci – произвольные постоянные.

Пример.

. Составим характеристическое уравнение:

k1 = 1, k2 =5. Для k = 1 получаем систему для определения  : , то есть

. Примем , тогда . При k = 5  ,

. Тогда . Следовательно, общее решение системы имеет вид: .

В случае кратных корней характеристического уравнения решение системы (36.6) имеет вид

 , где γ – кратность корня ks.

Пример.

. Характеристическое уравнение имеет вид:

k1 = k2 = 3. Пусть x = (c1 + c2 t)e3t, y = (c3 + c4 t)e3t. Выразим постоянные с3 и с4 через с1 и с2. Для этого подставим найденные решения в одно из уравнений системы и приравняем коэффициенты при e3t и te3t: (3c1 + c2 + 3c2t)e3t = (2c1 + c3)e3t + (2c2 + c4)te3t, c3 = c1 + c2,

c4 = c2. Итак, общее решение системы получено в форме: x = (c1 + c2 t)e3t, y = (c1+ с2 + c2t)e3t.  

Замечание. Для неоднородной системы (36.1) общим решением, так же как для неоднородного уравнения, будет сумма общего решения соответствующей однородной системы и частного решения неоднородной системы. При подборе частных решений справедлив принцип суперпозиции.

Пример.

. Найдем частное решение в виде: . При подстановке получим: , откуда А = 3, В = 1. Прибавив к полученному частному решению общее решение соответствующей однородной системы, запишем общее решение исходной системы: x = c1et + 2c2e4t + 3e5t, y = -c1et + c2e4t + e5t.

Устойчивость решений дифференциальных уравнений и их систем

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

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

                                        (37.1)

с начальными условиями yi(t0) = yi0 . 

Определение . Решение φi (t)  (ǐ = 1,2,…,n) называется устойчивым по Ляпунову, если     такое, что для всякого решения yi (t) той же системы, начальные условия которого удовлетворяют неравенствам  , для всех  справедливы неравенства

                                                      (37.2)

(то есть близкие по значениям решения остаются близкими для всех ).

Если хотя бы для одного решения yi (t) неравенства (37.2) не выполняются, решение φi (t)  называется неустойчивым.

Если решение φi (t)  не только устойчиво по Ляпунову, но и удовлетворяет условию

                                                        (37.3)

при , то это решение называется асимптотически устойчивым.

Замечание. Одно условие (37.3) не обеспечивает устойчивость решения.

Фазовая плоскость.

Дифференциальное уравнение второго порядка

                                                         (37.4)

равносильно системе уравнений первого порядка

.                                       (37.5)

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

                                         (37.6)

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

,                                                       (37.7)

которое каждой точке ставит в соответствие наклон проходящей через нее интегральной кривой.

Точки покоя.

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

Замечание. Особые точки классифицируются по характеру фазовых траекторий в их окрестности.

Исследование на устойчивость некоторого решения    системы дифференциальных уравнений можно свести к исследованию тривиального решения – точки покоя, расположенной в начале координат, преобразуя систему к новым переменным:  - отклонениям прежних неизвестных от решения, исследуемого на устойчивость. В новых переменных система (37.1) принимает вид:

,            (37.8)

Простейшие типы точек покоя.

Исследуем расположение траекторий в окрестности точки покоя х = 0, у = 0 системы двух линейных однородных уравнений с постоянными коэффициентами:

, где  .                                  (37.9)

Характеристическое уравнение при этом имеет вид:

.

Рассмотрим различные наборы корней этого уравнения:

  1.  k1 и k2 действительны и различны. Тогда общее решение системы (37.9) можно задать так: . При этом возможны следующие случаи:

а) если k1 < 0 и k2 < 0, то точка покоя асимптотически устойчива, так как , и все точки, находящиеся в начальный момент t = t0 в любой δ – окрестности начала координат, при достаточно большом t переходят в точки, лежащие в сколь угодно малой ε – окрестности начала координат, а при  стремятся к началу координат. Такая точка покоя называется устойчивым узлом.

б) если k1 > 0, k2 >0, можно свести исследование к предыдущему случаю заменой t на –t. При этом фазовые траектории имеют такой же вид, но направление движения меняется на противоположное, то есть при увеличении t точка удаляется от начала координат, поэтому подобная точка покоя – неустойчивый узел – неустойчива по Ляпунову.

в) при k1 > 0, k2 < 0 точка покоя тоже неустойчива, так как движущаяся по траектории

точка с возрастанием t выходит из ε – окрестности начала координат. Точка покоя рассматриваемого типа называется седлом.

  1.  k1,2 = p ± qi . Тогда общее решение системы (37.9) можно представить в виде

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

а) p < 0, q0. Тогда  при , а тригонометрические функции являются ограниченными. Поэтому фазовые траектории являются спиралями, асимптотически приближающимися при  к началу координат. Таким образом, точка покоя асимптотически устойчива. Она называется устойчивым фокусом.

б) p > 0, q0. Изменяется направление движения по фазовым траекториям, следовательно, точки удаляются от начала координат и точка покоя неустойчива – неустойчивый фокус.

в) р = 0. Траекториями являются замкнутые кривые, окружающие точку покоя, называемую в этом случае центром. Такая точка покоя устойчива, так как можно подобрать такое δ, что замкнутые траектории, начальные точки которых лежат в δ – окрестности начала координат, не выходят за пределы ε – окрестности начала координат (x² (t) + y² (t) < ε² ).

  1.  Корни кратны: k1 = k2.

а) k1 = k2 < 0. Тогда общее решение   стремится к нулю при , и точка покоя вновь называется устойчивым узлом. При  получаем частный случай устойчивого узла – так называемый дикритический узел.

б) k1 = k2 > 0. Направление движения по траекториям меняется -  неустойчивый узел.

афедра информатики и высшей математики КГПУ



 

Другие похожие работы, которые могут вас заинтересовать.
10314. Системы линейных уравнений 50.64 KB
  Система уравнений называется линейной если все уравнения входящие в систему являются линейными. Систему уравнений принято записывать с помощью фигурной скобки например: Определение: Пара значений переменных обращающая в верное равенство каждое уравнение с двумя переменными входящих в систему называется решением системы уравнений. При решении системы линейных уравнений возможны следующие три случая: система не имеет решений; система имеет ровно одно решение; система имеет бесконечно много решений.
19443. Методы решения обыкновенных дифференциальных уравнений 72.36 KB
  Для начала рассмотрим метод Эйлера так как является самым простым из существующих численных методов решения дифференциальных уравнений и в конце сравним результаты. Метод Эйлера является явным одношаговым методом первого порядка точности основанном на аппроксимации интегральной кривой кусочно-линейной функцией...
19450. Численные методы решения нелинейных уравнений 156.56 KB
  Видим что обе части не являются алгебраическими и содержат тригонометрические формулы значит это трансцендентное уравнение для решения которого не существует формул для отыскания корней. Построим график для того чтобы примерно определить промежутки содержащие корни см. Для этого реализуем метод половинного деления. Для того чтобы использовать эту функцию напишем скрипт который будет выводить первые пять корней отмечать их на графике а также использовать встроенную функцию для проверки решения и вычислять резонансные частоты стержня.
6215. Численные методы решения обыкновенных дифференциальных уравнений 1.42 MB
  Порядком обыкновенного дифференциального уравнения называется порядок старшей производной от искомой функции. Общим интегралом уравнения. неявным образом причем число постоянных интегрирования равно порядку уравнения. Общим решением обыкновенного дифференциального уравнения называется функция.
13533. Простейшие дифференциальные уравнения 113.18 KB
  Заметим что многие ученые увлекались расшифровками обычно неудачными чужих анаграмм. Константа имеет смысл начальной координаты тела поскольку из 1 следует что и ее принято обозначать символом . Произвольная постоянная имеет смысл начальной скорости поскольку из формулы 6 ясно что . 1 Рассмотрим колебания грузика массой m на пружинке с коэффициентом жесткости k который лежит на плоском горизонтальном столе предполагая что трение грузика об поверхности стола отсутствует.
6903. Обыкновенные дифференциальные уравнения. Курс высшей математики 1.04 MB
  Обыкновенные дифференциальные уравнения. Решение различных геометрических физических и инженерных задач часто приводят к уравнениям которые связывают независимые переменные характеризующие ту ил иную задачу с какой – либо функцией этих переменных и производными этой функции различных порядков. Наивысший порядок производных входящих в уравнение называется порядком дифференциального уравнения. Общим решением дифференциального уравнения называется такая...
6217. Методы нахождения корней системы нелинейных уравнений 284.94 KB
  Методы нахождения корней системы нелинейных уравнений. Для системы из 2 уравнений это можно сделать графически но для систем высоких порядков удовлетворительных методов отделения корней не существует. Проблема решения системы 1 возникает при решении многих прикладных задач например поиска безусловного экстремума функций многих переменных с помощью необходимых условий...
7509. Система линейных одновременных уравнений 40.11 KB
  Идентификация модели. Уравнения 1 называются структурной формой модели экономического процесса. Коэффициенты ij и bij называются структурными коэффициентами модели. На практике достаточно часто рассматриваются модели в которых такие члены имеются в наличии.
8653. Решение систем линейных уравнений 91.38 KB
  Если число уравнений системы совпадает с числом неизвестных m=n то система называется квадратной. Решением линейной системы 2.2 называется набор чисел которые при подстановке вместо неизвестных обращают каждое уравнение системы в верное равенство. Для системы линейных уравнений матрица А = называется матрицей системы а матрица А= называется расширенной матрицей системы Определение.
841. Теоретическая информатика. Решение систем линейных уравнений методом Крамара 90.58 KB
  Информация ее виды и свойства Единицы количества информации: вероятностный и объемный подходы. Теоретическая информатика математическая дисциплина использующая методы математики для построения и изучения моделей обработки передачи и использования информации. Но как правило эти модели наполнены конкретным содержанием связанным со спецификой информации того объекта который нас интересует. В них разрабатываются методы позволяющие использовать достижения логики для анализа процессов переработки информации с помощью...
© "REFLEADER" http://refleader.ru/
Все права на сайт и размещенные работы
защищены законом об авторском праве.