Вышка - дифуры (1021455), страница 2
Текст из файла (страница 2)
Утв2: (решение ЛОДУ первого порядка): y′+p(x)y = 0 y′= -p(x)y dy/dx= -p(x)y Поделим на y≠0 => ∫dy/y = - ∫p(x)dx ln|y| = - ∫p(x)dx + C |y| = e-∫p(x)dx+ C |y|=eCe-∫p(x)dx y = ±k*e-∫p(x)dx k≠0
Утв3: (связь между ЛОДУ и ЛНДУ) Пусть y′ + p(x)y = q(x) y1(x) – его частное решение. Общее решение ЛНДУ: y=y1 +k*e-∫p(x)dx
Док-во: y′+p(x)y=q(x) замена: y=z+y1 z-? y′=z′+y1′ - подставим в уравнение. z′+y1′ + p(x)*(z+y1) = q(x) z′+y1′+p(x)z + p(x)y1 =q(x) т.к. y1 – частное решение => y1′ + p(x)y1 = q(x) => ура-ние имеет вид. z′ +p(x) =0 - ЛОДУ, из утв2 => z = k*e-∫p(x)dx y=y1+k*e-∫p(x)dx
Решение ЛНДУ: (1) Способ Бернулли: y′+p(x)y=q(x) ищем y в виде: y=UV U, V-? => ура-ние U′ V +V′ U + p(UV)=q(x) U′ V +U(V′ +pV)=q(x) 1) Выберем V(x) так чтобы V′+pV=0 т.е. найдем частное решение ЛОДУ V=e-∫p(x)dx при k=1 => ура-ние U′ V + U(V′ +pV)=q(x) => U′ C=q(x) U′=q(x)/V = q(x)*e∫p(x)dx => U′=q(x)*e∫p(x)dx => U= ∫q(x)*e∫p(x)dx dx +C => y=UV= e-∫p(x)dx (∫q(x)*e∫p(x)dx +C) – общее решение ЛНДУ.
(2) Способ Лагранжа: (метод вариаций произвольной постоянной) y′+p(x)y=q(x) (**) 1) решаем ЛОДУ y′+p(x)y=0 => y = k*e-∫p(x)dx 2) решаем ЛНДУ: y=k(x)*e-∫p(x)dx (*) Надо найти k(x)-? => Подставим (*) и (**): y′=k(x)′ *e-∫p(x)dx + k*e-∫p(x)dx *(-∫p(x)dx)′ = k′*e-∫p(x)dx - p(x)*k*e-∫p(x)dx => k′*e∫pdx =qx (обязательно должно произойти сокращения!!!) => k′ *e-∫pdx =q(x) => k′=q(x)*e-∫pdx => k(x) =∫q(x)*e∫pdx dx +C => общее решение ЛНДУ: y=(∫q(x)*e∫pdx dx +C)*e-∫pdx
(3) Уравнения сводящиеся к линейным: 1) Уравнение Бернулли: y′+p(x)y=q(x)ym m≠0 m≠1 Делим на ym (y=0 отдельно) => y-my′+p(x)y1-m=q(x) Замена: y1-m=z z-? Дифф. это равенство. z′ =(1-m)*y-m*y′ => y-m*y′= z′/(1-m) => уравнение имеет вид: z′/(1-m) + p(x)z =q(x) - линейное уравнение.
2) Уравнение вида: f ′(y)*y′+p(x)*f(y) =q(x) Замена: z=f(y) z-? => z′=f ′(y)*y′ => ура-ние сводится к виду: z′+p(x)z=q(x) – линейное уравнение
3) Уравнение вида: y′+p(x) = q(x)*eny Поделим на eny => e-ny*y′+e-ny*p(x) = q(x) => z = e-ny z′= -e-ny*n*y′ => -z′/n +p(x)*z =q(x) –линейное уравнение.
7. Уравнения не разрешенные относительно производных:
1) F(x,y,y′) =0 – уравнение не разрешенное
2) y′=f(x,y) - уравнение разрешенное относительно производной.
(1) Если удается выразить y′=fk(x,y) => проинтегрировать фk(x,y,C)=0 => общий интеграл исходного уравнения: ∏{k=1; m} Фk(x,y,C) =0 = Ф(x,y,C)
Условие для поиска особых решений: Нужно составить такую систему уравнений { Ф(x,y,C)=0 Фс′(x,y,C)=0 } => исключить из нее C
Пример: (y′)2+xy=y2+xy′ (y′)2-xy′+xy-y2=0 - квадратное уравнение. a) y1′=x-y z=x-y z′=1-y′ dz/(1-z) =dx z=1-ke-x => y=x+1+ke-x б) y2′=y - раздел. переменные. y=kex
(2) Метод введения параметра: Общее решение может выглядеть 1) y=φ(t,C) 2) Ф(y,x,C) – неявно задана.
(2.1) Рассмотрим x= φ(y′) введем параметр y′=t => y=φ(t) осталось выразить через t. x=φ(t) y′=t y′=dy/dx =t => dy=t*dx но x=φ(t) => dx=φ(t) dt => dy=tφ(t)dt => y=∫tφ(t)dt +C Общее решение x=φ(t) y=∫tφ(t)dt+C
Зам (условие для поиска особого решения): Пусть F(x,y,y′)=0 – исходное уравнение t=y′ {F(x,y,t)=0 Ft′(x,y,t)=0 -исключить t }
(3) y=ψ(y′) Пусть y′=t – параметр => y=ψ(t) => осталось выразить x(t) y′=dy/dx => dy=tdx но: y=ψ(t) => но: dy= ψ(t)dt => ψ(t)dt=tdx – ура-ние с разделю переменными. => x=∫(ψ(t)dt/t)dt+C y=ψ(t) -общее решение.
(4) Уравнение Лагранжа:
Опр: y=φ(y′)x+ψ(y′) – уравнение Лагранжа. φ(y′)≠y′
Решение: y′=t y=φ(t)x+ψ(t) => dy=(φ(t)x+ψ(t))dt + φ(t)dx y′=t=dy/dx dy=tdx => tdx=(φx+ψ)dt +φdx Поделим на dt => x=x(t) tx= φx+ψ+φx x(t-φ) = φx + ψ Поделим на (t-φ(t)) Решение t=φ(t) – находится отдельно. x-φx/(t-φ) =ψ/(t-φ) -φ/(t-φ)=f(t) ψ/(t-φ)=g(t) => x+f(t)x=g(t) - ЛНДУ после его решения получим x=α(t,C) => y= φ(t)*α(t,C)+ψ(t)=β(t,C) => x=α(t,C) y=β(t,C) -общее решение.
(5) уравнение Клеро: (частный случай (4)): y=y′x+ψ(y′) y′=t y=tx+ ψ(t) => dy=xdt +tdx + ψ′(t)dt dy=tdx => xdt + ψ′(t)dt=0 => dt(x+ ψ′(t))=0 1 случай: dt=0 t=C y=xC+ ψ(C) –общее решение. 2 случай: x+ ψ′(t)=0 x=- ψ(t) => y= -tψ(t)+ψ(t) - особое решение.
8. Уравнение высших порядков:
Общий вид: F(x,y,y′….y(n))=0 (*) – уравнение n-ого порядка. Выразим старшую производную через ост. производные: y(n) =f(x,y,y′…y(n-1)) - уравнение в нормальной форме (разреш. относит. старшей производной)
Опр: (Задача Коши): Найти решение ура-ния (*), удовлетворяет начальным данным: y=φ(x) - ? при x=x0 φ(x0) =y0 φ′(x0)=y0′ φ′′(x0)=y0′′ …. φ(n-1)(x0)=y0(n-1)
Пример: y′′′=1 => y′′=x+C1 => y′=x2/2+xC1+C2 => y=x3/6 +C1x2/2 +xC2+C3 Если имеются начальные условия => найдем C1 C2 C3
Теорема: Рассмотрим y(n) =f(x,y,y′…y(n-1)) Рассмотрим задачу Коши: (x0 y0 y0′ ….y0(n-1)) =M0 є Rn (точка n-мерного про-ва) Пусть f –непрер. в нек. окрестности точки M0 и f0′ fy′′ fy^(n-1)′ -непрер в этой окр-ти. => задача Коши имеет единственное решение. y=φ(x) определ. в нек. окр-ти точки x.
Опр: (общее решение): y=φ(x,С1….Сn) –общее решение в области DcRn если: 1) при любом C1…Cn :
y=φ(x) является решением (подходит в ура-ние)
2) при любом M0єRn (x0, y0′…y0(n-1)) Существует константы С10 … Cn0 : y=φ(x,C10 , …, Cn0) –явл. решением. y(x0)=y0 y′(x0)=y0′ … y(n-1)=y0(n-1)
Уравнение допускающие понижение порядка.
(1) Ура-ние вида: y(n) =f(x) -решается последовательным интегрированием. Введем z=y(n-1)(x)-? z′=f(x) y(n-1)(x)=∫f(x)dx+C1 y(n-2)(x)=∫(∫f(x)dx) +xC1+C2 и т.д.
(2) Уравнение не содержащее y,y′….y(n-1) Общий вид такого уравнения: F(x,y(k)..y(n))=0 Обозначим z=y(k) => z′=y(k+1) …. z(n-k)=y(n) => уравнение примет вид: F(x,z,z′…z(n-k))=0 Оно решаемо если: z=φ(x,C1 …Cn-k) обратная замена: y(k)=φ(x,C1 …Cn-k) => получен случай (1) => k раз интегрировать.
Пример: xy′′′=y′′ z=y′′-? => xz′=z dz/z=dx/x z=k1x y′′=k1x => y′=k1x2/2 +k2 => y=k1x3/6 +k2x+k3
(3) Уравнение не содержащее x: Общий вид: F(y,y′…y(n))=0 Замена: y′=z(y) z-? Пусть n=2 F(y,y′,y′′)=0 Выразим y′′ : y′′=(y′)′=dz/dx=(dz/dy)*(dy/dx)=z(*)*z => F(y, z, z(*)z) =0 => Если его удается решить, то получим z(y)=φ(y,C1) y′=φ(y,C1) - ура-ние разд. переменными => dy/φ(y,C1)=dx x+C2=∫(dy/φ(y,C1))
Пример: yy′′+(y′)2=1 dy ′/dx=dz/dx=(dz/dy)*(dy/dx)=z(*)*z => yz(*)z+z2=1 zyz(*)=1-z2 => dz/(1-z2) =dy/y (z= ±1) => ∫(2zdz/(1-z2) =∫2dy/y => ln|z2-1|=-2ln|y| +C z2-1=k/y2 => z2=k1/y2 +1 => (y′)2 =(k1+y2)/y2 y′= ± √(k1+y2)/y dy*y/√(k1+y2) = ±dx => x+k2= ±√(y2+k1) => (x+k2)2 =y2+k1 - Ответ – это общее решение в неявном виде.
9. Линейные уравнения n-ого порядка
Опр: y(n)+p1(x)y(n-1) +p2(x)y(n-2) +…+p(n-1)(x)y′ +pn(x)y=f(x) (*) D(p1)=…=D(pn)=D(f) =(a, b) , все функции непрер. => уравнения такого вида наз. ЛУ n-ого порядка.
Опр: L(y)=левая часть уравнения (*) – называется линейным дифф. оператором, примененным к функции y => (*): L(y)=f(x) –лин. неоднородное – ЛНДУ. На ряду с этим рассмотрим ура-ние (**) L(y)=0 - ЛОДУ
Утв1: (сво-ва лин. оператора): 1) L(ky)=k*L(y) 2) L(y1+y2)= L(y1)+L(y2) -это очевидно, т.к. оператор явл комбинацией произв.
Утв2: (единств. решения задачи Коши для ЛНДУ): Пусть L(y)=f(x) (x0, y0, y0′, ….y0(n-1)) – имеет единственное решение y=φ(x), опред. на (a, b) при любом y0, y0′, ….y0(n-1) – этот факт следует из: y(n)= f(x)-p1y(n-1)-….-pny – удовлетворяет условию теор о существовнии (параграф 8)
Утв3: (свойства решений ЛОДУ):
1) Пусть y1(x) – решение L(y)=0 , тогда: С1y1 –тоже решение.
2) y1 и y2 – решение (**) , тогда y=y1+y2 – тоже решение. Док-во: (2) – L(y1)=0 L(y2)=0 по условию => По утв1: L(y1+y2)= L(y1)+ L(y2) =0 => y=y1+y2 – тоже решение ЛОДУ.
Следствие: Пусть y1 …. ym - решения (**) тогда: y=C1y1+ …+ Cmym - тоже решение.
Утв4: (связь между ЛНДУ и ЛОДУ): Пусть L(y)=f(x) L(y)=0 y=y1(x) –частное решение (*) y=φ(x, C1..Cn) –общее решение (**), тогда: y=y1(x) + φ(x, C1..Cn) –общее решение (*)
Док-во: рассмотрим L(y) =f(x) замена: y=z+y1 z-? => ура-ние примет вид: L(z+y1)=f(x) => L(z)+L(y1)=f(x) но y1 –частное решение ЛНДУ => L(y1)=f(x) => исх. ура-ние примет вид: L(z)=0 – ЛОДУ По условию z= φ(x, C1..Cn) => y= φ(x, C1..Cn) +y1
Пример: y′′-y= -x ЛОДУ: y′′-y=0 y′′=y общее решение будет выглядеть так: y=C1ex+ C2e-x а частное решение: y1=x => Ответ: y=C1ex+C2e-x+ x
10. Решение ЛОДУ (L(y)=0 (**)):
Опр 1: y1(x); … ym(x) – наз. линейно незав. на (a, b) если верно: α1y1(x)+ ….+ αmym(x) ≡0 => α1=…= αm=0 иначе функция наз линейно завис (ЛЗ).
Пример: 1) ex e-x -ЛНЗ на R т.к. α1ex+ α2e-x ≠0 при (y1)2 +(y2)2≠0 2) 2ex 3ex => α1=3 α2= -2 => α1y1+ α2y2 =0 => ЛЗ
Утв: (свойство ЛЗ систем): Пусть y1 … ym – ЛЗ , то одна из них явл. линейной комбинацией других.
Док-во: y1 …. ym - ЛЗ => существуют α1 ….αm – не все равны 0 : α1y1+ ….+ αmym ≡0 Пусть αm≠0 => αmym = -α1y1 - … - αm-1 ym-1 ym = (-α1/αm)*y1 - …. - (-αm-1/αm)*ym-1
Опр: Пусть y1 … yn – решение ЛОДУ (**) Пусть y1 … yn – ЛНЗ на (a, b) тогда это наз. фундаментальной системой решений (ФСР)
Опр: Составим определитель (n×n) W(x) = │ y1 ….. yn ; y1′ ….. yn′ ; ……; y1(n-1) ….. yn(n-1) │ - это определитель Вронского для y1 ….. yn