ОДУ - 1 (1086549), страница 16

Файл №1086549 ОДУ - 1 (А.М. Денисов, А.В. Разгулин - Обыкновенные дифференциальные уравнения) 16 страницаОДУ - 1 (1086549) страница 162019-05-08СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 16)

, y n (t) – фундаментальная система решений, а Y (t) – фундаментальная матрица.Замечание 4.3.1. Фундаментальная матрица неединственна. Полагая в задаче Коши (4.9) начальное условие Y (t0 ) = B, det B 6= 0, мыполучим другую фундаментальную матрицу.Замечание 4.3.2.

Так как элементы aij (t) матрицы системы вещественны, то и фундаментальная матрица может быть выбранавещественной.4.3.2. Общее решение линейной однородной системыОпределение 4.3.2. Общим решением линейной однородной системы дифференциальных уравнений n-го порядка называется зависящееот n произвольных постоянных решение этого уравнения такое, чтолюбое другое решение системы может быть получено из него в результате выбора некоторых значений этих постоянных.98Глава 4.

Общая теория линейных системТеорема 4.3.2. Пусть Y (t) = (y 1 (t), y 2 (t), . . . , y n (t)) – фундаментальная матрица для линейной однородной системыdy(t)= A(t)y(t)dtна отрезке [a, b]. Тогда ее общее решение представимо в видеy OO (t) = c1 y 1 (t) + c2 y 2 (t) + · · · + cn y n (t) = Y (t)c,(4.10)где c1 , c2 , . . . , cn – произвольные постоянные, c = (c1 , c2 , . . . , cn ).Доказательство. По теореме 4.1.2 вектор-функция Y (t)c является решением однородной системы для любых c ∈ Cn . Согласно определениюобщего решения осталось показать, что для любого наперед заданногорешения y(t) линейной однородной системы найдется вектор константec ∈ Cn такой, что на отрезке [a, b] выполнено равенствоy(t) = Y (t)ec.(4.11)Для построения ec зафиксируем произвольное t0 ∈ [a, b] и вычислимy 0 = y(t0 ). Рассмотрим систему линейных алгебраических уравненийотносительно ec = (ec1 , ec2 , .

. . , ecn )> :Y (t0 )ec = y0 .(4.12)В силу невырожденности матрицы Y (t0 ) c определителем det Y (t0 ) 6=0 эта система имеет единственное решение ec = (ec1 , ec2 , . . . , ecn )> . Тогдафункции ye(t) = Y (t)ec и y(t) являются решениями одной и той же задачиКошиdy(t)(4.13)= A(t)y(t), y(t0 ) = y 0 ,dtи по теореме единственности обязаны совпадать, что доказывает (4.11).Отметим, что для фиксированного решения y(t) вектор констант ec ∈ Cnв представлении (4.11) определен однозначно.Следствие 4.3.1. В ходе доказательства теоремы 4.3.2 была фактически выведена формула для решения задачи Коши (4.13) с произвольным начальным вектором y 0 .

Действительно, из (4.12) имеемec = Y −1 (t0 )y 0 и после использования (4.11) получаемy(t) = Z(t, t0 )y 0 ,Z(t, t0 ) = Y (t)Y −1 (t0 ).(4.14)4.3. Фундаментальная система решений и общее решение системы99Функциональная матрица Z(t, t0 ) называется матрицантом. Какматричная функция переменной t она является решением следующейзадачи КошиdZ(t, t0 )= A(t)Z(t, t0 ),dtZ(t0 , t0 ) = Y (t0 )Y −1 (t0 ) = E.Замечание 4.3.3.

Так как элементы aij (t) матрицы системы вещественны, то и общее решение естественно искать в классе вещественнозначных функций. Тогда при выборе вещественной фундаментальной матрицы (это всегда возможно в рассматриваемом случае)формула (4.10) при c ∈ Rn дает общее вещественнозначное решениелинейной однородной системы дифференциальных уравнений.4.3.3. Общее решение линейной неоднородной системы, методвариации постоянныхРассмотрим линейную неоднородную систему с непрерывным вектором f (t) = (f1 (t), f2 (t), . . .

, fn (t))> :dy(t)= A(t)y(t) + f (t),dtt ∈ [a, b].(4.15)Как и в предыдущем пункте, Y (t) обозначает фундаментальную матрицу соответствующей (4.15) однородной системы dy(t)/dt = A(t)y(t) стой же самой матрицей коэффициентов A(t).Определение 4.3.3. Общим решением линейной неоднородной системы дифференциальных уравнений n-го порядка (4.15) называетсязависящее от n произвольных постоянных решение этой системы такое, что любое другое решение системы (4.15) может быть полученоиз него в результате выбора некоторых значений этих постоянных.Теорема 4.3.3.

Общее решение y OH (t) линейной неоднородной системы дифференциальных уравнений (4.15) представимо в видеy OH (t) = Y (t)c + y H (t),∀c = (c1 , c2 , . . . , cn )> ∈ Cn ,(4.16)где y H (t) – некоторое (частное) решение неоднородной системы (4.15).Доказательство. В силу линейности системы (4.15) вектор-функцияy OH (t) является решением (4.15) для любого вектора констант c ∈ Cn .100Глава 4.

Общая теория линейных системСогласно определению общего решения, осталось показать, что для любого наперед заданного решения ye(t) системы (4.15) найдется векторконстант ec ∈ Cn такой, что на отрезке [a, b] будет выполнено равенствоye(t) = Y (t)ec + y H (t).(4.17)Пусть ye(t) – решение (4.15). Разность y(t) = ye(t) − y H (t) двух решений неоднородной системы является решением однородной системы,dy(t)/dt = A(t)y(t). Тогда по теореме 4.3.2 об общем решении линейной однородной системы найдется такой вектор констант ec ∈ Cn , чтона рассматриваемом отрезке выполнено равенство y(t) = Y (t)ec, котороеприводит к (4.17).Построение одного из частных решений неоднородной системы может быть проведено методом вариации постоянных и выражено с помощью введенного в (4.14) матрицанта Z(t, τ ).Теорема 4.3.4.

Для любого t0 ∈ [a, b] формулаZty H (t) =Z(t, τ )f (t)dτ,t ∈ [a, b],(4.18)t0задает частное решение неоднородной системы (4.15), удовлетворяющее условию y H (t0 ) = 0.Доказательство. Воспользуемся методом вариации постоянных, согласно которому частное решение неоднородной системы ищется в виде,повторяющем структуру (4.10) общего решения однородной системы, вкотором вектор констант c заменен на пока произвольную непрерывно дифференцируемую вектор-функцию c(t) = (c1 (t), c2 (t), .

. . , cn (t))> ,а именно:y(t) = Y (t)c(t).(4.19)Поскольку фундаментальная матрица удовлетворяет однородномууравнению dY (t)/dt = A(t)Y (t), тоdy(t)dY (t)dc(t)dc(t)=c(t) + Y (t)= A(t)Y (t)c(t) + Y (t).dtdtdtdt(4.20)Подставляя выражения (4.19) и (4.20) в уравнение (4.15), получаемуравнение для определения вектор-функции c(t):Y (t)dc(t)= f (t).dt(4.21)4.4. Фундаментальная система решений: постоянная матрица101В силу невырожденности фундаментальной матрицы это уравнениеможно переписать в виде dc(t)/dt = Y −1 (t)f (t) и проинтегрировать отt0 до t. Полагая по определению, что интеграл от вектор-функции естьвектор, составленный из интегралов координатных функций, имеемZtc(t) =Y −1 (τ )f (τ )dτ.t0После подстановки в (4.19) окончательно получаемZty(t) = Y (t)c(t) = Y (t)Y−1Zt(τ )f (τ )dτ =t0Z(t, τ )f (τ )dτ.t0Следствие 4.3.2.

Решение y(t) = y(t; y 0 ) задачи Коши для линейнойнеоднородной системыdy(t)= A(t)y(t) + f (t),dtt ∈ [a, b]с заданным в точке t0 ∈ [a, b] начальным условиемy(t0 ) = y 0имеет видZty(t; y 0 ) = Z(t, t0 )y 0 +Z(t, τ )f (τ )dτ.(4.22)t04.4. Построение фундаментальной системы решенийдля линейной однородной системы дифференциальных уравнений с постоянной матрицейРассмотрим однородную систему обыкновенных дифференциальных уравнений с постоянной матрицей коэффициентов A(t) ≡ A =(ai,j ), ai,j ∈ R, i, j = 1, .

. . , n:dy(t)= Ay(t).dt(4.23)102Глава 4. Общая теория линейных системПо аналогии со скалярным уравнением y 0 (t) = ay(t), которое имеет решение y(t) = h exp{at} для любого h ∈ C, будем искать нетривиальныерешения системы (4.23) в видеy(t) = h exp{λt},h = (h1 , h2 , . . . , hn )> ∈ Cn ,λ ∈ C.(4.24)Подстановка вектор-функции (4.24) в систему (4.23) приводит к задаченахождения таких λ ∈ C, при которых система линейных алгебраических уравнений(4.25)(A − λE)h = θимеет нетривиальное решение h. Как известно из курса линейной алгебры, такие λ называются собственными значениями матрицы A, а отвечающие им векторы h – собственными векторами матрицы A. Собственные значения и только они являются корнями характеристического многочлена M (λ):M (λ) = det(A − λE) = 0.(4.26)4.4.1.

Построение фундаментальной системы решений, когдасуществует базис из собственных векторовПоскольку характеристический многочлен имеет степень n, то поосновной теореме алгебры у него имеется ровно n корней (собственныхзначений), с учетом их кратности λ1 , . . . , λn , λj ∈ C. Из курса линейнойалгебры известно, что существует не более, чем n линейно независимыхсобственных векторов матрицы A. Остановимся сначала на случае, когда количество линейно независимых собственных векторов в точностиравно n. Заметим, что в этом случае собственные векторы составляютбазис пространства Cn .Теорема 4.4.1. Пусть у матрицы A имеется ровно n линейно независимых собственных векторовh1 ,h2 ,...,hn ,отвечающих соответствующим собственным значениямλ1 ,λ2 ,...,λn .4.4. Фундаментальная система решений: постоянная матрица103Тогда вектор-функцииy 1 (t) = h1 exp{λ1 t}, y 2 (t) = h2 exp{λ2 t}, .

. . y n (t) = hn exp{λn t} (4.27)образуют фундаментальную систему решений (4.23) на произвольномотрезке [a, b].Доказательство. Рассмотрим произвольный отрезок [a, b]. Для любогоj = 1, . . . , n собственное значение λj и соответствующий собственныйвектор hj удовлетворяют уравнению (4.25), и тогда каждая из векторфункций y j (t) = hj exp{λj t} является решением системы (4.23) на [a, b]по построению.Докажем линейную независимость на отрезке [a, b] построенной системы функций. Для этого, согласно теореме 4.2.3 об альтернативе дляопределителя Вронского, достаточно убедиться, что det Y (t) 6= 0 длянекоторого t ∈ [a, b], где Y (t) = (y 1 (t), y 2 (t), .

. . , y n (t)). Рассмотрим отрезок [c, d], включающий в себя исходный отрезок [a, b] и точку t = 0:[a, b] ⊆ [c, d],0 ∈ [c, d].Вектор-функции из (4.27) являются решениями системы (4.23) на отрезке [c, d]. В принадлежащей этому отрезку точке t = 0 определительВронскогоdet Y (0) = det(h1 , h2 , . . . , hn ) 6= 0,так как в противном случае составляющие Y (0) столбцы – собственныевекторы h1 , h2 , . . . , hn – были бы линейно зависимыми. Согласно теореме 4.2.3 об альтернативе для определителя Вронского det Y (t) 6= 0 навсем отрезке [c, d], а значит и на его части [a, b].4.4.2. Построение фундаментальной системы решений, когдане существует базиса из собственных векторовРассмотрим случай, когда количество существующих у матрицы Aлинейно независимых собственных векторов строго меньше, чем порядок системы n.

Выпишем все попарно различные собственные значенияλj с соответствующими кратностями kj :λ1 ,k1 ,λ2 ,k2 ,...,...,λ`k`,,λi =6 λj при i 6= j,kj > 1, k1 + k2 + · · · + k` = n.104Глава 4. Общая теория линейных системПусть далее λ ∈ {λ1 , . . . , λ` } обозначает одно из собственных значений с соответствующей кратностью k. Покажем, что каждому такомусобственному значению можно сопоставить ровно k вектор-функций,являющихся решениями однородной системы (4.23). Если размерностьs = dim Ker(A − λE) собственного подпространства, определяющая количество линейно независимых собственных векторов для данного собственного значения, равна кратности собственного значения, s = k, тоискомые функции строятся согласно (4.27).Если размерность собственного подпространства меньше кратности собственного значения, s < k, то, как известно из курса линей111ной алгебры, можно выбрать собственные векторы h1 , h2 , .

. . , hs так,1что состоящая ровно из k векторов система собственных векторов hjmи присоединенных векторов hj , m = 2, . . . , pj , j = 1, . . . , s, pj > 1,p1 + p2 + · · · + ps = k, которую запишем в виде1h1 ,2h1 ,.........1hj ,2hj ,.........1hs ,2hs ,.........pjps. . . hj , . . . hs ,p1h1 ,удовлетворяет уравнениям1Ahj2AhjmAhjpjAhj=1λhj ,21= λhj + hj ,...mm−1= λhj + hj ,...pjpj −1= λhj + hj.(4.28)С помощью собственных и присоединенных векторов построим семейство из следующих k функцийy 1j (t)y 2j (t)...1= hj exp{λt},t 12=hj + hj exp{λt},1!(4.29)4.4. Фундаментальная система решений: постоянная матрица105 mt m−1 t2 m−2tq m−qym(t)=h+h+h+···+h+ ...jj1! j2! jq! jtm−11··· +hj exp{λt},(m − 1)!... pt pj −1 t2 pj −2tq pj −qjpy j j (t) = hj + hj+ hj+ · · · + hj+ ...1!2!q!tpj −1 1 ··· +h exp{λt},(pj − 1)! jj = 1, .

. . , s.Докажем, что все функции из построенного семейства являются решениями линейной однородной системы (4.23). Рассмотрим функциюmymj (t), вычислим ее производную dy j (t)/dt и сгруппируем результаттак, чтобы удобно было воспользоваться соотношениями (4.28). Имеемdy mj (t)=dt m−1t m−2 t2 m−3tq m−q−1tm−21= hj+ hj+ hj+ · · · + hj+ ··· +h +1!2!q!(m − 2)! jt m−1 t2 m−2tq m−qtm−2m2+ λhj+ · · · + λhj+ ··· +λh ++ λhj + λhj1!2!q!(m − 2)! jtm−11+λhj exp{λt} =(m − 1)! mtt2 m−2tqm−1m−q= Ahj + Ahj+ Ahj+ · · · + Ahj+ ...1!2!q!tm−11··· +Ahj exp{λt} = Ay mj (t),(m − 1)!m = 1, . .

Характеристики

Тип файла
PDF-файл
Размер
1,02 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Почему делать на заказ в разы дороже, чем купить готовую учебную работу на СтудИзбе? Наши учебные работы продаются каждый год, тогда как большинство заказов выполняются с нуля. Найдите подходящий учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6549
Авторов
на СтудИзбе
299
Средний доход
с одного платного файла
Обучение Подробнее