- •ПРИКЛАДНАЯ МАТЕМАТИКА
- •Е.А.КОГАН
- •Москва 2007
- •Кафедра “Прикладная и вычислительная математика"
- •Е.А. Коган
- •Москва 2007
- •1.1. Основные понятия
- •Дифференциальным уравнением называется уравнение, в которое неизвестная функция входит под знаком производной или дифференциала.
- •Обыкновенное дифференциальное уравнение n-го порядка в самом общем виде записывается так:
- •Для дифференциальных уравнений первого порядка различают общее, частное и особое решения, а также общий, частный и особый интегралы.
- •Таким образом, общее решение дифференциального уравнения можно определить как множество всех частных решений уравнения.
- •Особым решением дифференциального уравнения называется решение, которое не может быть получено из общего решения ни при одном частном значении произвольной постоянной.
- •Часто при интегрировании уравнения первого порядка не удается найти общее решение в явном виде, а получается конечное (не дифференциальное) соотношение вида
- •Знание изоклин позволяет во многих случаях даже для не интегрируемых явно дифференциальных уравнений получить графическое решение задачи Коши и выявить характер интегральных кривых.
- •Пример. Построить методом изоклин интегральную кривую уравнения
- •Очевидно, это уравнение с разделенными переменными. Интегрируя его, получим
- •Следовательно, общий интеграл уравнения будет
- •Интегрируем полученное уравнение с разделенными переменными
- •Тогда
- •Полагаем
- •Тогда уравнение примет вид
- •Дифференциальное уравнение вида
- •В результате уравнение приводится к однородному
- •Возвращаясь к старой переменной, получим
- •Подставляя (1.14) в (1.12), получаем дифференциальное уравнение с разделяющимися переменными относительно варьируемой функции v(x), интегрируя которое, находим эту функцию.
- •В результате, общее решение уравнения (1.12) может быть представлено в виде
- •Общее решение уравнения примет вид
- •Находим произвольную постоянную C из начального условия:
- •После разделения переменных получим
- •Тогда уравнение (1.18) примет вид
- •Следовательно,
- •Уравнением Бернулли называется уравнение вида
- •Уравнение вида
- •Следовательно, его общий интеграл, а значит, и общий интеграл уравнения (1.24) имеет вид
- •Это уравнение интегрируется непосредственно n раз. При каждом интегрировании порядок уравнения понижается на единицу, и появляется произвольная постоянная. В результате общее решение уравнения будет иметь вид
- •После интегрирования получим
- •Следовательно, общее решение уравнения (2.22) будет
- •Далее ищем решение уравнения (2.21) в форме, аналогичной по структуре выражению (2.23), но произвольную постоянную в (2.23) заменяем неизвестной функцией
- •Подставляя (2.24) в (2.21), получим
- •Отсюда следует
- •и с учетом (2.14)
- •N линейно независимых частных решений линейного однородного дифференциального уравнения n-го порядка образуют фундаментальную систему решений этого уравнения.
- •а общее решение запишется в виде
- •Рассмотрим его реализацию для линейных дифференциальных уравнений второго порядка:
- •Подстановка (2.42) в (2.40) приводит к следующей системе дифференциальных уравнений
- •Пример. Решить уравнение
- •Определитель этой системы
- •Поэтому для определения варьируемых функций согласно (2.44) получаем дифференциальные уравнения вида
- •Интегрируя уравнения (2.49), находим
- •Подставляя (2.50) в (2.48), получим общее решение уравнения в виде
- •Рассмотрим линейное однородное дифференциальное уравнение второго порядка
- •Пример. Найти собственные значения и собственные функции однородной краевой задачи:
- •Подставляя это решение в граничные условия, убеждаемся в том, что они удовлетворяются только при
- •Подставляя это выражение в граничные условия, получим
- •Чтобы краевая задача имела ненулевые решения, необходимо принять
- •Пример. Найти собственные значения и собственные функции однородной краевой задачи для уравнения (2.51) при граничных условиях
- •Общее решение уравнения согласно (2.55) имеет вид
- •Вычисляем
- •Ненулевое решение ее существует тогда и только тогда, когда определитель системы равен нулю:
- •Раскрывая этот определитель, получим уравнение относительно параметра
- •корни которого являются собственными значениями задачи:
- •Уравнением Эйлера называется линейное уравнение вида
- •Пример. Решить уравнение
- •Уравнение (2.66) есть уравнение Эйлера второго порядка. Применим замену независимой переменной и производных по формулам (2.64), (2.65). Тогда уравнение примет вид
- •Его общее решение
- •Эти решения находятся методом подбора (см. выше) и имеют вид
- •Поэтому общее решение уравнения (2.67) будет
- •Пример. Решить задачу Коши для уравнения
- •Общее решение уравнения имеет вид
- •Будем искать частное решение однородной системы в виде
- •Подставляя (2.99) в (2.98), получим систему
- •Характеристическое уравнение системы будет
- •Общее решение однородной системы (2.98) запишется в виде
- •Этот метод применим к решению систем неоднородных линейных уравнений n-го порядка. Ограничимся для простоты нормальной системой двух линейных уравнений с постоянными коэффициентами
- •Пусть общее решение однородной системы уравнений известно:
- •Пример. Решить систему
- •Общее решение однородной системы, согласно (2.105), имеет вид
- •Принимаем частное решение системы (2.111) в виде
- •Интегрируя эти уравнения, получим
- •Общее решение системы запишется в виде:
- •Пусть требуется решить краевую задачу для уравнения
- •Будем искать решение задачи в виде
- •Задача сводится к решению дифференциального уравнения
- •Ограничиваясь двучленным приближением, будем искать решение уравнения в виде
- •Выбранные функции удовлетворяют всем перечисленным выше требованиям. Они линейно независимы, непрерывно дифференцируемы и удовлетворяют граничным условиям (2.125).
- •Решение методом Бубнова
- •Решение методом наименьших квадратов
- •Решение методом коллокаций
- •В результате находим
- •Операционное исчисление представляет собой своеобразный и эффективный метод решения различных математических задач, прежде всего, дифференциальных уравнений. В основе операционного исчисления лежит понятие преобразования Лапласа.
- •Согласно формулам Эйлера [7]
- •Учитывая формулу (3.3) и применяя теорему линейности, получим
- •Аналогично
- •Таблица 4
- •Применяя формулу (3.25), находим
- •Поэтому
- •то применяя теорему линейности, окончательно найдем
- •Разложение данной дроби на простейшие имеет вид
- •После приведения к общему знаменателю получим
- •Пусть дано неоднородное линейное дифференциальное уравнение второго порядка с постоянными коэффициентами
- •Пример. Решить операционным методом уравнение
- •ПРИЛОЖЕНИЕ 1
- •Вариант № 1
- •Вариант № 2
- •Вариант № 3
- •Вариант № 4
- •Вариант № 5
- •Вариант № 6
- •Вариант № 7
- •Вариант № 8
- •Вариант № 9
- •Вариант № 11
- •Вариант № 12
- •Вариант № 13
- •Вариант № 14
- •Вариант № 15
- •Вариант № 16
- •Вариант № 17
- •Вариант № 18
- •Вариант № 19
- •Вариант № 20
- •Вариант № 21
- •Вариант № 22
- •Вариант № 23
- •Вариант № 24
- •Вариант № 25
- •Вариант № 26
- •Вариант № 27
- •Вариант № 29
- •ОПЕРАЦИОННОМУ ИСЧИСЛЕНИЮ
- •Таблица 2
- •Таблица 3
- •ЛИТЕРАТУРА
- •ПРИЛОЖЕНИЕ 1. ВАРИАНТЫ РАСЧЕТНО-ГРАФИЧЕСКОЙ
- •РАБОТЫ ПО ОБЫКНОВЕННЫМ ДИФФЕРЕНЦИАЛЬНЫМ
- •УРАВНЕНИЯМ……………………………………………………… 100
- •ПРИЛОЖЕНИЕ 2. ВАРИАНТЫ РАСЧЕТНО-ГРАФИЧЕСКОЙ
- •РАБОТЫ ПО ОПЕРАЦИОННОМУ ИСЧИСЛЕНИЮ…………… 131
41
N линейно независимых частных решений линейного однородного дифференциального уравнения n-го порядка образуют фундаментальную систему решений этого уравнения.
Для построения фундаментальной системы решений уравнения L(y) = 0 с постоянными коэффициентами его частные решения ищутся в виде показательных функций
y = ek x , |
(2.32) |
где k – неизвестные постоянные числа. Подстановка (2.32) в дифференциальное уравнение (2.29) приводит к алгебраическому уравнению вида
k n + a k n−1 |
+ a |
2 |
k n−2 |
+ + a |
n−1 |
k + a |
n |
= 0 . |
(2.33) |
1 |
|
|
|
|
|
|
Алгебраическое уравнение (2.33) той же степени, что и порядок дифференциального уравнения (2.29), с теми же коэффициентами на соответст-
вующих местах называется характеристическим уравнением.
Заметим, что характеристическое уравнение (2.33) получается из дифференциального уравнения (2.29) формальной заменой i – ой производной y(i) числом k i .
Характеристическое уравнение (2.33) имеет n корней (с учетом их кратности). В зависимости от вида корней соответствующие частные линейно независимые решения будут иметь различный вид (см. таблицу 1).
Составляя линейную комбинацию этих решений, получим общее ре-шение линейного однородного уравнения L(y) = 0 в виде (2.31).
Таким образом, построение общего решения линейных однородных дифференциальных уравнений n-го порядка с постоянными коэффициентами сводится к чисто алгебраической проблеме решения соответствующих характеристических уравнений.
Построение решения однородных линейных дифференциальных уравнений является обязательным первым этапом решения более общих неоднородных уравнений (2.28). Поэтому примеры построения фундаментальной системы решений для однородных уравнений приведены ниже.
42
Таблица 1
Вид частных решений линейного однородного уравнения L(y)=0 в зависимости от вида корней характеристического уравнения
|
Вид корней |
Вид частных решений |
||
|
|
|
|
|
1 |
Корни ki |
y1 = exp(k1x), |
|
|
|
(i =1,2, ,n) |
y2 = exp(k2 x), |
||
|
действительны |
……………… |
||
|
и различны |
yn = exp(kn x). |
||
|
|
|
|
|
2a |
Комплексные корни |
|
|
|
|
k1 =α + β i, |
y1 = exp(αx) cos β x , |
||
|
k2 =α − β i |
y2 = exp(αx) sinβ x. |
||
|
|
|
|
|
2б |
Мнимые корни |
y1 = cos βx, |
|
|
|
k1 = β i , |
y2 =sinβx. |
|
|
|
k2 = −β i |
|
|
|
3а |
Кратные |
y1 = exp(kx), |
|
|
|
действительные |
y2 = x exp(kx), |
||
|
корни |
y3 = x2 exp(kx), |
||
|
k1 = k2 = = km = k |
|||
|
……………….… |
|||
|
(m- кратность кор- |
|||
|
ym = xm−1 exp(kx), |
|||
|
ня) |
|||
|
|
|
|
|
|
|
|
|
|
3б |
Кратные |
y1 = exp(αx)cos βx , |
|
y2 = exp(αx) sin βx , |
|
комплексные |
y3 = xexp(αx)cos βx. |
|
y4 = xexp(αx) sinβx , |
|
корни: |
y5 = x2 exp(αx)cos βx , |
|
y6 = x2 exp(αx) sinβx. |
|
k1 =α + β i |
|
||
|
…………………………… |
|
………………………… |
|
|
кратности m, |
|
||
|
……………… |
|
…………… |
|
|
k2 =α − β i |
|
||
|
y2m−1 = xm−1 exp(αx)cos β x |
|
y2m = xm−1 exp(αx)sin βx |
|
|
кратности m |
|
||
|
|
|
|
|
|
|
|
|
|
43
Частное решение y неоднородного линейного дифференциального уравнения n-го порядка (2.28) может быть получено или методом подбора (методом неопределенных коэффициентов) или методом вариации произвольных постоянных.
2.2.2. Метод подбора частного решения
Этот метод называют ещё методом неопределенных коэффициентов.
Он не является универсальным и применим, если правая часть уравнения (2.28) в общем случае имеет вид
f (x) = P(x) eα x cos β x +Q(x)eα x sinβ x, |
(2.34) |
где P(x) и Q(x) - одночлены или многочлены (в общем случае различных степеней от x). Пусть при этом n - наивысшая степень одного из многочленов P(x) или Q(x).
Алгоритм построения частного решения неоднородного линейного дифференциального уравнения (2.28) следующий:
1.Находим корни характеристического уравнения (2.33).
2.Сравниваем заданную правую часть уравнения (2.28) с общим выражением (2.34), при котором применим метод подбора, и находим из этого сопоставления три числа: n, α, β.
3.Сравниваем "контрольное" комплексное число (α + β i) с корнями ха-
рактеристического уравнения и находим число m корней, совпавших с комплексным числом α + β i (если таких корней нет, то m=0).
4. Принимаем частное решение неоднородного уравнения (2.28) в виде
y |
|
= xm[R |
(x)eα x cos β x +T (x)eα x sinβ x], |
(2.35) |
|
n |
n |
|
|
где Rn (x), Tn (x) |
|
- многочлены одной и той же n-ой степени, но с неоп- |
ределёнными и различными коэффициентами.
5. Записываем решение (2.35) в развернутой форме в зависимости от n.
Так, если |
n = 0, то y = xm (Aeax cos β x + Beax sinβ x), |
||
если |
n =1, то |
y = xm [R1(x)eax cos β x +T1(x)eax sinβ x]= |
|
= xm [(Ax + B)eax cos β x + (Cx + D)eax sinβ x)]. |
|||
|
|
44
6. Подставляем y в исходное уравнение (2.28) и получаем систему алгебраических уравнений относительно неопределенных коэффициентов A, B,
C, D,...
Замечания.
1. Если правая часть уравнения (2.28) имеет более простой вид, например, содержит произведение степенной функции на показательную
f (x) = Pn (x)eax
(в частности, возможны случаи n=0 или (и) α = 0), или содержит только линейную комбинацию тригонометрических функций вида
f (x) = M cos βx + N sinβx,
где M и N - постоянные числа, то частные решения неоднородного уравнения можно искать в форме, указанной в таблице 2 (в неё для полноты включен также общий случай).
2. Правая часть уравнения L(y) = f (x) |
может содержать только функ- |
цию вида f (x) = M cos βx или функцию вида |
f (x) = N sinβx . Но частное ре- |
шение методом подбора следует искать в полной форме, содержащей и cos βx и sin βx (см. таблицу 2).
Процедура подбора неопределённых коэффициентов показана ниже на примерах.
Пример. Решить уравнение y′′ − 2 y′ = x3 −1.
1). Решаем сначала соответствующее однородное уравнение
y′′− 2 y′ = 0.
Составляем характеристическое уравнение, отыскивая частные решения уравнения в виде y = ek x . Получаем
|
|
|
|
|
|
|
k 2 − 2k = 0 . |
|
|
|
|||||
Корни этого уравнения k1 = 0, |
k2 |
= 2 - действительны и различны. Соответ- |
|||||||||||||
ствующие им частные линейно независимые решения (см. таблицу 1) |
|
||||||||||||||
y = e0 x =1, |
y |
2 |
= e2x . Поэтому общее решение однородного уравнения запи- |
||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
шется в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
0 |
=C y |
+ C |
2 |
y |
2 |
=C |
+ C |
e2x . |
(2.36) |
||
|
|
|
|
|
1 |
1 |
|
|
1 |
2 |
|
|
45
Таблица 2
Структура частного решения уравненияL(y) = f (x)
в зависимости от вида правой части
1 |
f (x) = Pn (x)eax , |
A. Если число α не совпадает ни с одним |
||||||||||
|
где Pn (x) -многочлен (или |
из корней характеристического уравнения: |
||||||||||
|
α ≠ k j |
( j =1,2, ,n), то частное решение |
||||||||||
|
одночлен) n-ой степени |
|||||||||||
|
следует принимать в форме |
|||||||||||
|
от x |
|||||||||||
|
|
|
|
|
|
|
|
y = Qn (x)eax , |
||||
|
|
где Qn (x) - многочлен n-ой степени с не- |
||||||||||
|
|
определенными коэффициентами. |
||||||||||
|
|
Б. Если α = k j (k j -корень кратности m), то |
||||||||||
|
|
|
|
|
|
|
|
y |
|
= xmQ (x)eax . |
||
|
|
|
|
|
|
|
|
|
|
n |
|
|
2 |
f (x) = M cos β x + N sinβ x, |
A. Если мнимое число β i не совпадает ни |
||||||||||
|
где M и N –заданные |
с одним из корней характеристического |
||||||||||
|
постоянные числа |
уравнения: β i ≠ k j |
(( j =1,2, ,n), то част- |
|||||||||
|
|
ное решение следует принимать в форме |
||||||||||
|
|
|
|
|
|
|
y = Acos βx + Bsinβx, |
|||||
|
|
где А и В–неопределенные коэффициенты. |
||||||||||
|
|
Б. Если β i = k j (k j - корень кратности m), |
||||||||||
|
|
то |
|
|
y = xm (Acos βx + Bsinβx). |
|||||||
3 |
f (x) = P(x)eax cos βx + |
А. Если комплексное число α + β i не сов- |
||||||||||
|
+ Q(x)eax sin βx, |
падает ни с одним из корней: α + β i ≠ k j |
||||||||||
|
где P(x) и Q(x)- много- |
( j =1,2, ,n), то частное решение следует |
||||||||||
|
члены (или одночлены) в |
принимать в форме |
|
|
||||||||
|
общем случае различных |
|
|
y |
|
= R |
(x)eα x cos βx |
+T (x)eα x sinβx, |
||||
|
степеней |
|
|
|
n |
|
|
|
|
n |
||
|
где Rn (x), Tn (x)- многочлены n-ой степени |
|||||||||||
|
|
(равной наивысшей степени одного из |
||||||||||
|
|
многочленов P(x) или Q(x)), но с неопреде- |
||||||||||
|
|
ленными и различными коэффициентами. |
||||||||||
|
|
Б. Если α + β i = k j |
(k j |
- корень кратности |
||||||||
|
|
m, то |
|
|
|
|
|
|
|
|||
|
|
y |
|
= xm [R |
n |
(x)eα x cos β |
x +T (x)eα x sinβ x]. |
|||||
|
|
|
|
|
|
|
|
|
|
n |
46
2). Находим частное решение заданного неоднородного уравнения ме-
тодом подбора, так как правая часть f (x) = x3 −1= P (x)- многочлен |
третьей |
||||
|
|
|
|
3 |
|
степени относится к первому из указанных в таблице 2 случаев. |
|
||||
Сравнивая функцию |
f (x) = x3 −1 с выражением f (x) = P (x)ea x , за- |
||||
|
|
|
|
n |
|
ключаем, что n = 3, α = 0. |
|
|
|||
Сравниваем |
α = 0 с корнями характеристического уравнения. |
Так как |
|||
α = k1, то m =1, |
и частное решение принимаем в виде |
|
|||
|
y |
|
= xQ (x) = x (Ax3 + Bx2 +Cx + D), |
(2.37) |
|
|
|
|
3 |
|
где A, B, C, D - неопределенные коэффициенты, подлежащие подбору. Подставляем (2.37) в исходное уравнение:
|
0 |
|
y = Ax4 + Bx3 +Cx2 + Dx, |
||
|
|
||||
|
−2 |
|
y′ = 4Ax3 +3Bx2 + 2Cx + D, |
||
|
1 |
|
y′′ =12Ax2 +6Bx + 2C. |
||
|
|
|
|
||
|
y′′ − 2y′ =12Ax2 + 6Bx + 2C −8Ax3 − 6Bx2 − 4Cx − 2D = x3 −1 |
||||
или |
−8Ax3 +(12A −6B)x2 +(6B −4C)x + 2C −2D = x3 −1. |
Приравниваем коэффициенты при одинаковых степенях x в левой и правой частях получающегося равенства многочленов и приходим к системе алгебраических уравнений относительно неизвестных коэффициентов:
x3 |
−8A =1, |
x2 |
12A −6B = 0, |
x1 |
6B −4C = 0, |
x0 |
2C −2D = −1. |
Решая эту систему, находим коэффициенты: A = -1/8, B = -1/4, C = -3/8,
D=1/8.
В результате, частное решение уравнения примет вид
y = − |
1x4 |
− |
1 x3 |
− |
3 x2 |
+ |
1x. |
(2.38) |
|
8 |
|
4 |
|
8 |
|
8 |
|
Складывая (2.36) и (2.38), получим общее решение уравнения в виде