Лекции по обыкновенным дифференциальным уравнениям. В.И.Дмитриев (немного урезанные) (1114451), страница 4
Текст из файла (страница 4)
Значит ∃(9.8)∂ y=U .∂ µТеорема доказана.Без доказательства приведем теорему о разложении решения возмущенной задачи помалому параметру µ .Т е о р е м а 9.3. Пусть в областиD : (t0 < t < t0 + T , y − y0 ≤ a,| µ |≤ ε ) функция f (t , y, µ )обладает непрерывными иравномерно ограниченными частными производными по y и µ до порядка (n+1)включительно. Тогда существует сегмент [t0 ; t0 + T ] , на котором для решения y (t ,m)возмущенной задачи (9.6) cправедливо асимптотическое представление∂ y (t ,0)µ n ∂ n y (t ,0)+ ... ++ O( µ n+1 )(9.9)y (t , µ ) = y (t ,0) + µn∂ µn! ∂µНеравенство Чаплыгина.Если имеются две задачи Коши dy dz= f1 (t , y ), = f 2 (t , z ), dt dt y (t0 ) = y0 , z (t0 ) = z0 ,причем в D выполняются условияf1 (t , y ) ≤ f 2 (t , y ) и y0 ≤ z0 , то при t0 < t < t0 + T имеем y (t ) ≤ z (t ) .Сингулярное возмущение дифференциального уравнения dy = f ( y, t , µ ), t ∈ [ 0, T ] , dt y (t = 0) = y0,|µ |≤ εвозникает, если f ( y, t , µ ) при µ → 0 имеет нерегулярность, т.е.
ведет себя особым(сингулярным) образом. Это, например,1) f ( y, t , µ = 0) не удовлетворяет условиям теоремы ∃ и ! решения2) f → ∞ при µ → 0 и т.п.Наиболее частый и практически важный случай – это малый параметр при старшейпроизводнойµ y((tn)) = F (t , y, y′,... y ( n−1) )(9.10)или, соответственно, система с малым параметром при одной производной24dy1= F1 (t , y )dtdy2= f 2 (t , y );dt...................µ⇒dy1 1= F1 (t , y ) = f1 (t , y, µ ) → ∞dt µпри µ → 0;(9.11)dyn= f n (t , y ).dtп.10. Линейное дифференциальное уравнение n-го порядка и егосвойства. Сведение к нормальной системе первого порядка.Существование решения.nLn ( y ) = ∑ ak (t ) y ( n−k ) (t ) = f ( t ) , a0 (t ) ≠ 0, t ∈ [t0 , t0 + T ] .(10.1)k =0f (t ) ≡ 0 уравнение однородное,f (t ) ≠ 0 уравнение неоднородное.Т е о р е м а 10.1.
Линейность уравнения сохраняется при замене переменногоилинейномпреобразованиифункцииt = ϕ (τ ), ϕ ∈ Cn , ϕ ′ ≠ 0y (t ) = α (t ) z (t ) + β (t ), α ,β ∈ Cn , α ≠ 0.Д о к а з а т е л ь с т в о.dy dy dτ1dy d 2 y1d 1dy 1. t = ϕ (τ ) ⇒===;2ϕ ′ (τ ) dτ dtϕ ′ (τ ) dτ ϕ ′ (τ ) dτ dt dt dtdk y.dτ kСледовательно, сохраняется линейность уравнения.2. y = α z + β ⇒ y′ = α z′ + α ′z + β ′ ⇒ y′′ = α z′′ + 2α ′z′ + α ′′z + β ′′и т.д. все линейно. Приведя подобные члены, получим линейное уравнение.Т е о р е м а 10.2.
Для линейного дифференциального уравнения выполняетсяпринцип суперпозиции m mLn ∑ ck yk = ∑ Ck Ln ( yk )(10.2) k =1 k =1Применение принципа суперпозиции:и т.д. y ( k ) линейная комбинацияM1) Для суммы правых частей f = ∑ f mm =1MLn ( ym ) = f m ⇒ y = ∑ ym .m =1Это суммирование источников.2) Разделение задачи Коши на неоднородную с нулевыми начальными данными и наоднородную с начальными данными.25 Ln ( y ) = f( n −1)( n −1) y (t0 ) = y0 , y′(t0 ) = y0′ , ..., y (t0 ) = y0y = U (t ) + V (t ) L(U ) = f , L(V ) = 0,( n −1)( n −1)(t0 ) = 0.(t0 )=y0( n−1) .U (t0 ) = 0, ..., UV (t0 ) = y0 , ..., V3) Разделение начальных данных для однородного уравнения.n Ln ( y ) = 0⇒ y = ∑U m (t ) y0( m−1)( n −1)( n −1)= y0m =1 y (t0 ) = y0 , ..., y L(U m ) = 0 (k )U m (t0 ) = 0 k ∈ [ 0, n − 1] , k ≠ m (m)U m (t0 ) = 14) Комбинация решений однородного уравнения есть решение однородного уравнения.Т е о р е м а 10.3. ( ∃ и ! решения на всем интервале).Если коэффициенты α k (t ) и правая часть f (t ) есть непрерывные функции приt ∈ [t0 , t0 + T ] , то решение ∃ и ! на всем интервале [t0 , t0 + T ].
(т.к. условия теоремы ∃ и !выполняются на всем интервале).Линейное дифференциальное уравнение (10.1) сводится к нормальной линейнойсистеме дифференциальных уравнений. Введем вектор-функциюu (t ) = ( u1 = y (t ), u2 = y′(t ), ..., un = y ( n−1) (t ) ) ,для которой получим нормальную линейную систему уравненийприm ∈ [1, n − 1]um+1 (t ) n(10.3)um′ (t ) = ak (t )b(t)u(t)+f(t)приm=n,b=.n − k +1k∑ ka0 (t ) k =1В общем случае нормальная линейная система уравнений записывается в виде:ˆ +F,u′(t ) = Au(10.4)где матрица Aˆ = {α (t )} m, k ∈ [1, n] . В дальнейшем мы будем подробно рассматриватьmkлинейную систему дифференциальных уравнений, т.к. уравнение n-го порядка сводится кчастному случаю такой системы.п.11. Линейное дифференциальное уравнение 2-го порядка.
Понижениепорядка уравнения. Уравнение Риккати.Рассмотрим линейное уравнение 2-го порядкаy′′(t ) + a1 (t ) y′(t ) + a2 (t ) y (t ) = f (t ) .26(11.1)Оно сводится к системе второго порядка введением вектор-функцииu (t ) = ( u1 (t ) = y (t ), u2 (t ) = y′(t ) ) , для которой получаем системуu1′(t ) = u2 (t )u2′ (t ) = − a2 (t )u1 (t ) − a1 (t )u2 (t ) + f (t )или1 0ˆ + f , Aˆ = 0u′(t ) = Au(11.2) −a −a , f = .f 1 2У линейного однородного уравнения ( f = 0 ) можно понизить порядок, введя новуюфункциюy′(t ).(11.3)Z (t ) =y (t )Тогдаy′(t ) = Z (t ) y (t ), y′′ = Z ′y + Zy′ = Z ′y + Z 2 y .(11.4)Подставив (11.4) в (11.1) при f = 0 , получимZ ′(t ) + Z 2 + a1 (t ) Z + a2 (t ) = 0 .(11.5)Полученное уравнение является уравнением Риккати.Общий вид уравнения Риккати:y′(t ) = p (t ) y 2 + q (t ) y + r (t ) ,Z (t )приводится к виду (11.5).которое заменой искомой функции y = −p(t )p′ Z ′(t ) + Z 2 (t ) − q (t ) + Z + r (t ) p (t ) = 0 .pВ уравнении (11.5) можно убрать член, содержащий Z , не изменяя коэффициентапри Z 2 с помощью замены искомой функцииa (t )(11.6)Z (t ) = u (t ) − 1 .2Тогда из (11.5) получимa 2 + 2a1′(11.7)u′(t ) = −u 2 (t ) + R(t ); R(t ) = 1− a2 .4Если R(t ) = const = R0 , то переменные разделяются и мы имеемdu= dtR0 − u 2илиu (t ) = R0 ⋅271− e−2 R0 t1+ e−2 R0 t(11.8)Тогда, согласно (11.3) и (11.6), получим−2 R0 ty′(t )1− e= R0 ⋅−2y (t )1+ eR0 t−a1 (t )= q (t ) .2(11.9)Откудаty (t ) = y (t0 )e∫ q (τ ) dτt0.(11.10)п.12.
Общая теория однородных линейных системобыкновенных дифференциальных уравнений.Линейная однородная системаn / yi (t ) = ∑ aik (t ) yi (t ), i ∈ [1, n ] , t ∈ [t0 , t0 + T ] ,k =1 y (t ) = y 0 .i i 0(12.1) y1 Если обозначить матрицу Aˆ = {aik (t )} , а y (t ) = ... , то задача Коши yn ˆ (t ), t ∈ [t , t + T ] y (t ) = Ay0 0,(12.2)0 y (t0 ) = yL( y ) ≡ y′ − Ay – линейный оператор, следовательно, к нему применим принципсуперпозицииM ML ∑ Cm( m ) y = ∑ Cm L( ( m ) y ) .(12.3) m=1 m=1(m)Черезy – обозначаем m-ое решение, чтобы отличить от m-ой производной y ( m ) .MЕсли f (t ) = ∑ am f(m), тоm =1Т е о р е м а 12.1. ПустьMy (t ) = ∑ am( m ) y , гдеm =1(1)L ( (m) y ) = f (m) .y (t ),..., ( n ) y (t ) - "n" решений однородной системыˆ = 0.y′ − Ay(12.4)Тогда матрица (1) y1 ,..., ( n ) y1 Wˆ (t ) = ..............
(1)(n) yn ,..., yn удовлетворяет матричному уравнениюˆ ˆ (t ) = 0Wˆ ′(t ) − AW(12.5)и, обратно, если матрица Wˆ (t ) удовлетворяет уравнению (12.5), то ее столбцы естьвектора, являющиеся решением уравнения (12.4).28Доказательство проводится покомпонентным дифференцированием.ˆ (C - постоянныйТ е о р е м а 12.2. Если Wˆ (t ) - решение (12.5), то y = WCˆ ˆ (Cˆ - постоянная матрица)вектор) удовлетворяет системе (12.4), а Zˆ = WCудовлетворяет матричному уравнению (12.5).Доказательство следует из принципа суперпозиций.О п р е д е л е н и е . Векторные функции (1) y (t ),... ( n ) y (t ) – линейно зависимы наинтервале τ = {t0 , t0 + T } , если ∃ ненулевой постоянный вектор C такой, что выполняетсятождествоˆ ≡ 0 при ∀t ∈τ .WC(12.6)(1)(n)Если условие (12.6) выполняется только при C ≡ 0 , то y (t ),..., y (t ) являютсялинейно независимыми.{(i )О п р е д е л е н и е .
Определителем Вронского для системы вектор- функцийy (t )}, i ∈ [1, n ] называется∆(t ) = Det Wˆ (t ) .Т е о р е м а 12.3. Если решения{ y}(k )(12.7)ˆ =0k ∈ [1, n ] однородной системы y′ − Ayлинейно зависимы на t ∈τ , то определитель Вронского ∆(t ) = 0 для ∀t ∈τ .Д о к а з а т е л ь с т в о.ˆ = 0 . Это линейно однороднаяИз линейной зависимости следует ∃C ≠ 0 такое, что WCсистема для C , следовательно, Det Wˆ = ∆(t ) = 0.Т е о р е м а 12.4. Если ∆(t ) = 0 хотя бы для одного t ∈τ , то ∆(t ) = 0 и для ∀t ∈τ , и ,следовательно,{ y} линейно зависимы на τ .(k )Д о к а з а т е л ь с т в о.Пусть при t = t0 ∈τ имеем ∆(t0 ) = 0.
Тогда ∃C ≠ 0 , которые удовлетворяют системеуравнений Wˆ (t0 )C = 0. Возьмем y (t ) = Wˆ (t )C . Согласно теореме 12.2 y решение задачиКошиˆ = 0 t ∈τ y′ − Ay y (t0 ) = 0.Следовательно, y ≡ 0 ∀t ∈τ по теореме единственности решения задачи Коши.ˆ = 0 для ∀t ∈τ ⇒ DetWˆ = ∆ (t ) = 0 для ∀t ∈τ .Тогда WC{ y}(k )Т е о р е м а 12.5 (альтернатива). Определитель Вронского ∆(t ) для решенияk ∈ [1, n ] однородной системы дифференциальных уравнений или ∆(t ) ≡ 0 для∀t ∈τ , что означает линейную зависимостьозначает линейную независимость{ y }.(k )29{ y },(k )или ∆(t ) ≠ 0 для ∀t ∈τ , чтоп.13. Фундаментальная система решений и общее решение для линейнойсистемы дифференциальных уравнений.О п р е д е л е н и е. Фундаментальной системой решений (Ф.С.Р.) однороднойсистемы уравнений называется "n" линейно независимых решений { ( k ) y } , k ∈ [1, n ] этойсистемы, а соответственно матрица Wˆ = { (1) y , (2) y ,..., ( n ) y } называется фундаментальнойматрицей системы.Фундаментальная матрица является решением матричного уравненияˆ ˆ (t ) ,Wˆ ′(t ) = AWпричем DetWˆ ≠ 0.Т е о р е м а 13.1.