Вышка - дифуры (1021455), страница 3
Текст из файла (страница 3)
Теорема: (условие ЛНЗ функций): y1 ….. yn – решение ЛОДУ будет ЛНЗ <=> при любом xє(a, b) W(x)≠0
Док-во: “<=” Дано: при любом xє(a, b) W(x)≠0 Но пусть при этом y1 … yn ЛЗ на (a, b) => α1 ….αn – не все равны 0 , т.ч. α1y1+ ….+ αnyn =0 по утв. о ЛЗ функциях yn – лин. комб. y1 … yn-1 => yn=β1y1+ …. + βn-1yn-1 βi = -αi/αn => W(x) = │y1 ….. yn-1 _ β1y1 + ...+ βn-1yn-1 ; y1′ ….. y′n-1 _ β1y1′ + ...+ βn-1y′n-1 ;......... ; y1(n-1) ….. yn-1(n-1) _ β1y1(n-1) + ...+ βn-1yn-1(n-1)│ => последний столбец – это комбинация первых (n-1) => W(x)=0 => это противоречие.
“=>” Дано: y1 … yn ЛНЗ Пусть сущ. x0 є (a, b) : W(x0)=0 => │y1(x0) ….. yn(x0); y1′(x0) ….. yn′(x0); ……; y1(n-1)(x0) ….. yn(n-1)(x0)│=0 Рассмотрим СЛЧ относительно C1 …Cn: Определитель системы равен W(x0)=0: { C1y1(x0) ….. Cnyn(x0)=0 ; C1y1′(x0) ….. Cnyn′(x0)=0 ; ……; C1y1(n-1)(x0) ….. Cnyn(n-1)(x0) =0 } => она имеет бесконечно много решений. => Пусть C1 …. Cn – ненулевое решение. (не все равны 0) Пусть y=C1y1 + … + Cnyn – т.к. C1 …. Cn –решение той системы => { y(x0)=0 ; y′(x0)=0 ; y(n-1)(x0)=0 } => т.к. в ЛОДУ р, р1,…, pn – непрер. на (a, b) => то любая задача Коши имеет единственное решение => имеем: y=C1y1 + …. + Cnyn - это решение задачи Коши (x0; 0,…, 0), но очевидно, что φ(x)≡0 – тоже решение этой задачи Коши => y(x)≡φ(x) => C1y1 + …. + Cnyn ≡0 , на (a, b) => не все Ci =0 => y1 .. yn – ЛЗ. => получено противоречие.
Теорема: (общее решение ЛОДУ): L(y)=0 (**) Пусть y1 … yn -ФСР , => y= C1y1 + …. + Cnyn –общее решение уравнения (**)
Док-во: Проверим 1) y= C1y1 + …. + Cnyn - очевидно (см. сво-ва решений) 2) Любая задача Коши x0є (a, b) y0=y0′ =…= y0(n-1) сущ. (C1 …Cn ) т.ч. y= C1y1 + …. + Cnyn - решение задачи Коши: Система: { C1y1(x0) ….. Cnyn(x0)= y0 ; C1y1′(x0) ….. Cnyn′(x0)= y0′ ; ……; C1y1(n-1)(x0) ….. Cnyn(n-1)(x0) = y0(n-1)} - ЛСУ относительно C1 …Cn. её det A = │ y1(x0) ….. yn(x0); y1′(x0) ….. yn′(x0); ……; y1(n-1)(x0) ….. yn(n-1)(x0)│= W(x0) y1 … yn – ФСР => W(x0)≠0 => по т. Кромера существ решен. C1 … Cn
11. Линейные уравнения с постоянными коэфф.:
y(n) +a1y(n-1) +…+ an-1y′ +any = f(x) L(y)=f(x) (*) L(y)=0 (**) Рассмотрим уравнение второго порядка (n=2): y′′ +py′+qy=0 (**)
Метод Эйлера: Решение ищется в виде y=eλx λ-? => подставим это в (**) => y′=λeλx y′′= λ2eλx eλx(λ2+pλ+q)=0 λ2+pλ+q=0 ($)
Опр: Ура-ние: (алгебраическое): λn +a1λ(n-1)+….. + anλ0 =0 – называется характеристическим уравнением.
(1) ($) имеет 2 вещ. корня: λ1 и λ2 => y1= eλ1x y2= eλ2x -решения. (**) Проверим, что y1 и y2 - ЛНЗ (=> ФСР y=C1y1 + C2y2 - общее решение). Для проверки сост. определитель Вронского: W(x) =│y1 , y2 ; y1′ , y2′ │=│ eλ1x , eλ2x ; λ1eλ1x , λ2eλ2x │= λ2e(λ1+λ2)x – λ1e(λ1+λ2)x = e(λ1+λ2)x (λ2 - λ1)≠0 => y1 ; y2 – ЛНЗ => Ответ: y=C1eλ1x + C2eλ2x
Зам: (о ЛНЗ двух функций): y1 y2 - ЛЗ <=> y1=αy2 при любом xє(a, b)
Док-во: => y1 y2 –ЛЗ => Пусть k1 k2 – не оба равны 0 k1y1+k2y2=0 => y1 – k2y2/ k1
<= y1=αy2 y1 – αy2=0 k1=1 k2= -α => ЛЗ
Пример: y′′-y=0 λ2-λ0=0 => λ2=1 λ1=1 λ2= -1 y1=ex y2=e -x =>Ответ: y=C1e x + C2e -x
(2) ($) имеет 2 компл. корня: (его D<0)
//Отступление: Про C-числа:
Опр: Пусть α =b+ai є C eα=ea *(cos b +i*sin b)
Опр: Пусть xєR Я функцию R→C (комплексная функция веществ. перем.)
y(x)=U(x)+i*V(x) U(x) и V(x) є R => y(k)(x) =U(k)(x) +i*V(k)(x) y=eα(x) (αєC) => y′=αeαx
Зам: L(y(x))=0 <=> {L(U)=0 и L(V)=0} - очевидно. //
Эти корни будут выглядеть след. образом: λ1,2 =a ± bi => eλ1x и eλ2x - решения (это комплексные функции). eλ1x = e(a+bi)x =eax+bxi =eax(cos bx +i*sin bx)= eaxcos bx +i*eax *sin bx => получаем решение. y1=eax *cos bx y2= eax *sin bx -ЛНЗ (т.к. y1 ≠αy2 ) => ФСР => Ответ: y= eax(C1cos bx +C2sin bx)
(3) ($) имеет 2 одинаковых веществ. корня: (D=0) => q=p2/4 => λ2+pλ+q=0, (λ+p/2)2=0 λ1,2 = -p/2 => y1 = e-px/2 – решение (**) y2 =x*e-px/2 - докажем y2 – решение (**) L(y2) =0 – проверим:
y′′+py′+p2y/4 =0 y′= e-px/2 – (xp/2)*e-px/2 y ′′= (-p/2)*e-px/2 - (p/2)*e-px/2 +(xp2/4)*e-px/2 = -pe-px/2 +(xp2/4)e-px/2 => подставим e-px/2 : e-px/2(-p+xp2/4 +p-xp2/2 +p2x/4) =0 => y2 –решение (**) y1 y2 – очевидно ЛНЗ => общее решение в 3 случае: y= e-px/2(C1 +C2x)
Решение методом Эйлера в общем случае (для ЛОДУ n-ого порядка с пост. вещ. коэфф.): y(n) +a1y(n-1) +…+ an-1y′ +any = 0
1) Составим характеристическое ура-ние: λn +a1λ(n-1) +an-1λ+ an=0 2) Решить его (n- корней с учетом кратностей) 3) Составить ФСР по правилам: ПРАВИЛО1: вещественному корню λ, кратности 1 соот 1 решение, кот. выглядит так: eλx ПРАВИЛО2: веществ. корню λ, кратности k, будет соот. k решений: eλx , xeλx , x2eλx , …. , x(k-1)eλx ПРАВИЛО3: паре комплексных сопряженных корней a±bi , соот. два решения: eaxcos bx ; eaxsin bx ПРАВИЛО4: паре комплексных сопряженных корней кратности k, соот. 2k решений: eaxcos bx ; xeaxcos bx ; x2eaxcos bx ; ….; xk-1eaxcos bx и eaxsin bx ; xeaxsin bx ;……; x(k-1)eaxsin bx
4) написать ответ: общее решение: (сложить все решения найденные в 3 шаге и домножить их на коэфф)
y=C1y1 +…+ Cnyn где y1 ….yn –ФСР
12. Решение ЛНДУ с пост. коэфф.:
Утв: 1) нужно подобрать 1 частное решение φ; 2) решить ЛОДУ: y=C1y1 + …. + Cnyn => Ответ: y= φ(x) + C1y1 + …. + Cnyn
Метод неопределенных коэфф. для подбора частных:
(I) L(y)=f(x) (*) f(x)= eaxPm(x) Pm- многочлен степени m. a: 1 случай: a – не явл. корнем характ. уравнения. Частное решение: φ(x)=eaxQm(x) Q-находится подстановкой. 2 случай: a – явл. корнем кратности k: φ(x)=xkeaxQm(x)
Пример: 1) y′′-5y′+6y=6x2-10x+2 y′′-5y′+6y=0 λ2-5 λ+6=0 λ1=2 λ23 y=C1e2x + C2e3x - решение ЛОДУ. Правая часть: e0xP2(x) a=0 - не корень => φ(x) = Ax2 +Bx+C A,B,C - ? y′=2A+B
y′′=2A 2A-10Ax-5B+6Ax2+6Bx+6C=6x2-10x+2 A=1 B=0 C=0
2) y′′-4y′+4y=2e2x λ1,2= 2 y= e2x(C1+C2x) Подбор частного решения: a=2 – корень кратности 2 => y=x2(e2x*b)=bx2e2x b-? находим производные и подставляем в уравнение => b=1 y1=x2e2x Ответ: у=у2x(C1+C2x +x2)
(II) L(y)=f(x) f(x)=eax(Pmcos bx +Qmsin bx) α=a+bi – контрольное число: 1 случай: α-не корень. y1=eax(Rs cos bx + Rs′ sin bx) Rs , Rs′ - неизвестные многочлены степени S ; S=max(m,n) 2 случай: корень кратности k (характ. ура-ния) => y1=xkeax(Rscos bx +Rs′sin bx)
Пример: 1) y′′+y′-2y=ex(cos x – 7sin x) => λ2+λ-2=0 λ1= -2 λ2=1 => y= C1ex +C2e-2x α=1+i - не корень => y1=ex(a*cos x+s*sin x) a, b-? Берем производные => a=2 b=1 y1= ex(2cos x+sin x) Ответ: y= C1ex +C2e-2x + ex(2cos x+sin x)
2) y′′+y=2sin x => λ= ±i y= C1cos x + C2sin x α=i (корень кратности 1) 2sin x=e0(2sin x+0cos x) y1=x(a*cos x +b*sin x) a, b-? b=0 a= -1 y1= -x*cos x Ответ: н=C1cos x +C2 sin x - x*cos x
Утв: (принцип наложения): L(y)=f1+f2 (ЛНДУ в общем случае) => y=φ(x, C1, … ,Сn) – общее решение ЛОДУ. Пусть y1 – частное решение: L(y)=f1 y2- частное решение: L(y)=f2 , тогда φ=y1+y2 –частное решение исходного уравнения. Док-во: Проверим факт: L(φ) =f1+f2 : L(y1+y2)=L(y1) +L(y2) =f1+f2 Пример: 1) y′′- 4y+8y=sin 2x α=2i – не корень y2=a*sin 2x +b*cos 2x a, b-? a=1/10 b=1/20 y2=(sin 2x)/20 +(cos 2x)/10 Ответ: y=e2x(C1cos 2x +C2sin 2x) + e2x/4 + y2
13. Метод вариации произвольных постоянных:
Рассмотрим частную ситуацию: Пусть n=2 y′′+p(x)y′+q(x)y=f(x) (*)
1) Решить ЛОДУ: y′′+p(x)y′+q(x)y=0 Пусть y1 и y2 – ФСР => │y1 , y2 ; y1′ , y2′ │≠0 при любом xє(a, b) y=C1y1 +C2y2
2) Считаем С1 и C2 – функциями от x , т.е. решение ищется в виде: y=C1(x)y1 +C2(x)y2 (3) C1(x) и С2(x) - ? подставим: y′=C1′y1 +C1y1′+C2′y2 +C2y2′ Дополнит. условия: С1′y1 +C2′y2 =0 (предположение) => y′=C1y1′ + C2y2′ (2) => y′′C1′y1′ +C1y1′′ +C2′y2′ +C2y2′′ (1) Подставим (1) (2) (3) в ура-ние (*): С1′y1′′+C1y1′′+C2y1′+C2y2′′+p(x)C1y1′ + p(x)C2y2′ +q(x)C1y1 +q(x)C2y2 =f(x) => C1′y1′+C2′y2′+C1(y1′′+py1′+q) + C2(y2′′+py2′+q)=f(x) y1′′+py1′+q=0 и y2′′+py2′+q=0 т.к. y1 и y2 – решения ЛОДУ => C1′y1′+C2′y2′ = f(x) имеем систему: { C1′y1+C2′y2 = 0 и C1′y1′+C2′y2′ = f(x) } Определитель: │y1 , y2 ; y1′ , y2′│ ≠0 => система имеет единственное решение. С1′=φ1(x) C2′=φ2(x) => C1=∫φ1dx +k1 C2=∫φ2dx +k2 Ответ: y=y1C1 + y2C2
В общем случае, при любом n: L(y)=f(x) - порядок n