где p,q – вещественные числа. Это же уравнение с нулевой правой частью, т.е. когда f(x)=0, называется линейным однородным ДУ. Вначале решают линейное однородное ДУ, а затем находят общее решение линейного неоднородного ДУ.
а) Линейные однородные уравнения дифференциальные уравнения с
постоянными коэффициентами.
y// +py/+qy = 0 . (14)
Они попадают в класс ДУ из раздела 3.3. Однако для них имеется более простой метод решения.
Пусть y1(x) и y2(x) - частные решения этого уравнения.
Определение. Два решения называются линейно зависимыми, если
они пропорциональны друг другу: y2(x) = α · y1(x)или y2(x) = b · y1(x) . Понятия линейной зависимости и линейной независимости применимы к любой паре функций. Например, и при линейно независимы, так как .Две функции ekxи z(x)ekxтакже линейно независимы, если z(x) – многочлен степени ≥ 1.
Для двух линейно независимых частных решений y1и y2ДУ(14) и произвольных постоянных C1и C2 линейная комбинация
y=C1 y1+ C2 y2 (15)
является решением ДУ (14). Эта формула в действительности даёт все решения линейного однородного ДУ (14),так как число (два) произвольных постоянных совпадает с порядком ДУ, и разным значениям произвольных постоянных соответствуют разные решения. Если решение дифференциального уравнения есть комплексная функция y(x)=u(x)+iv(x), то отдельно действительная часть u(x) и мнимая часть v(x) являются решениями.
Будем искать решение уравнения (14) в виде y(x)= ekx, где k - пока неизвестное число. Тогда y/=kekx, y//= k2ekx ,и уравнение примет вид:
k 2 ekx+ pkekx+qekx= 0.
Поделив обе части уравнения на ekx, получим алгебраическое, а именно, квадратное уравнение относительно неизвестного k:
k2+pk+q = 0
(это уравнение называется характеристическим), решение которого есть
.
Числа k1 и k2 называются корнями характеристического уравнения.
Рассмотрим три возможности ( >0, =0, <0) для подкоренного выражения. (Напомним, что p2 - 4 q называется дискриминантом характеристического уравнения.)
1. >0, тогда , и - два линейно независимые решения уравнения (14). Общее решение имеет вид
y(x)=C1+C2. (16)
2. =0, тогда . Формула y(x)= ekx дает только одно из двух линейно независимых решений: y1(x)= ekx. Второе решение y2(x) ищем в виде y2(x)= z(x)· ekx. Подставляя это в уравнение (14), получаем : z//=0, одним из решений которого является z = x. В результате получаем пару линейно независимых решений y1(x)= ekx, y2(x) =xekx иобщее решение ДУ (14):
Согласно формулам Эйлера eβix=cos βx+isin βx, e- βix= cos βx-isin βx.
Как указывалось выше, в случае комплексного решения отдельно eαx cos βx и eαx sin βx являются решениями. Следовательно, общее решение (вещественное) ДУ (14) имеет вид
y(x)= eαx(С1cos βx+С2 sin βx). (18)
Пример 11. y// -2y/-8y = 0.
Решение: k2-2k-8=0, k1=4, k2=-2. Согласно формуле (16) пишем
Ответ: y=C1e4x+C2e-2x.
Пример 12. y// - 6y/ + 9y = 0.
Решение: k0-6k+9=0, k1= k2=k =3. Согласно формуле (16)
Ответ: y = C1e3x+C2xe3x.
Пример13. y// - 6y /+ 13y = 0
Решение: k2-6k+13=0, k1=3-2i, k2=3+2i. Согласно формуле (18) имеем
Ответ:y = e3x(C1cos2x+C2sin2x).
б) Линейные неоднородные дифференциальные уравнения с постоянными коэффициентами.
Общее решение неоднородного линейного ДУ (13) равно сумме
y=C1y1+C2y2+y*
общего решения C1y1+C2y2 соответствующего однородного уравнения (14) и частного решения y* данного неоднородного уравнения (13).Частное решение y* определяется видом правой части f(x) в (13),и его можно явно определить для функций f(x) специального вида – произведений функции еlх и многочленов Pn(x) некоторой степени ≤ n . (При комплексном l берется вещественная часть такой такого произведения; см. ниже).
При вещественном l ,когда
f(x) = Pn(x)·elx ,
частное решение y* можно взять в виде
y* = xm·Qn(x)·elx,
где показатель m означает, со сколькими корнями k1 и k2 характеристического уравнения совпадает l ( показатель m может быть 0, или 1, или 2).
Если l=λ+μi- комплексное число, то правая часть f(x) cпециального вида есть вещественная часть произведения еlхна многочлен, т.е.
f(x) = eλx ·(Pn(x)·cos μx+ n(x)·sin μx).
Тогда частное решение ищется в виде
y* = xm·eλx ·( Qn(x)·cos μx + n(x)·sin μx),
где m=0. если λ+μi не совпадает ни с одним из корней k1 и k2 ,
и m=1, если λ+μi совпадает с одним из корней k1 и k2 .
Перечислим типичные случаи.
(а) f(x)=Pn(x) (Pn(x) - многочлен n -й степени), тогда y* ищем в виде:
y*=Qn(x), если k1≠ k 2, k1≠0, k 2≠ 0 ;
y*=xQn(x), если k1=0, k 2 ≠0;
y*=x2Qn(x), если k1= k 2=0.
(b) f(x)=aelxс вещественным l, тогда y* ищем в виде:
y*=Aelx, если l ≠k1, l≠k2;
y*=Axelx, если k1≠ k 2, l=k1либо l=k2;
y*=Ax2elx, если l=k1=k2.
(c) f(x)= eλx(a ·sinµ+b·cosµx) с комплексным l= λ+iμ, тогда y* ищем в виде:
y*= eλx· (A ·sin µx+B ·cos µx),если λ+iμ≠k 1 и λ+iμ≠k2;
y*=x·eλx (A ·sin µx+B · cosµx), если λ+iμ=k 1 или λ+iμ=k2.
Для нахождения частного решения y* применяется метод неопределённых коэффициентов : подставляем y* в уравнение (13) и находим неизвестные коэффициенты многочлена Qn(x) или неизвестные числа A,B.
Пример 14. y// -2y/-8y =2ex.
Решение: общее решение однородного уравнения есть С1e4x+С2e-2x, Частное ищем в виде y*=Аex. Подставляя y* в исходное уравнение, получим
A ex-2A ex-8A ex=2 ex. откуда находим А=-2/9 и y*=- 2/9 ex.
Ответ: y(x) = С1e4x+С2e-2x - 2/9 ex.
Пример 15. y// -2y/-8y = e4x .
Решение: общее решение однородного уравнения: С1e4x+С2e -2x; частное решение ищем в виде y*=Aхe4x. Подставляя в исходное уравнение y* и приравнивая коэффициенты при e4x , получим A = 1/6
и y*=1/6 e4x .
Ответ: y(x) = С1e4x+С2e-2x +1/6 xe4x.
Пример 16. y// -6y/+9y = e3x.
Решение: общее решение однородного уравнения есть С1e3x+С2xe3x, тогда частное решение ищем в виде y*=Aх2e3x. Подставляя y* в исходное уравнение и приравнивая коэффициенты при е3х ,получим: A =1/2,
y*=1/2· x2e3x
Ответ: y(x) = С1e3x+С2xe3x+ 1/2· x2e3x.
Пример 17. y// - 6y /+ 25y =50 x.
Решение: общее решение однородного уравнения есть e3x(C1 ·cos 4x+ C2·sin 4x); частное решение ищем в виде y*=Ax+B. Подставляя это в исходное уравнение и приравнивая коэффициенты при одинаковых степенях х,получим: A=2, B=12/25 , y*=2x+12/25.
Ответ: y = e3x(C1 cos 4x+ C2sin 4x)+2x+12/25.
Пример 18. y// -6y /+ 13y =sin x.
Решение: общее решение однородного уравнения естьe3x(C1 cos 2x+
C2sin 2x).Частное решение ищем в виде y*=A sin x +B cos x. Подставляя y* в исходное и приравнивая коэффициенты при sin x и cos x, получим А=1/15, В=1/30 , y*=1/15·sin x +1/30 ·cos x.
Ответ: y = e3x(C1cos2x+ C2sin2x)+ 1/15·sin x +1/30 ·cos x.