- •ПРИКЛАДНАЯ МАТЕМАТИКА
- •Е.А.КОГАН
- •Москва 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
|
|
|
|
|
|
|
|
|
|
|
|
14 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u = a1x + b1y. |
|
|
|
|
|
|
|
||||||||
Пример. Решить уравнение y |
′ |
= |
|
x + y + |
2 |
. |
|
|
|
|
|
|||||||||||
|
|
|
x +1 |
|
|
|
|
|
|
|||||||||||||
В этом уравнении a1 = b1 =1, |
|
c1 = 2, |
a2 |
|
=1, b2 |
= 0, |
c2 |
=1. Поэтому |
||||||||||||||
δ = |
|
1 |
1 |
|
≠ 0 . Полагая |
x = u +α, |
|
|
y = v + β, |
находим α |
и |
β из системы |
||||||||||
|
|
|
|
|||||||||||||||||||
|
|
1 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
системы уравнений: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
α |
+ β + 2 = 0, |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
α +1= 0. |
|
|
|
|
|
|
|
|
|||||||
Следовательно, α = −1, |
|
β = −1, |
|
|
и формулы перехода от старых перемен- |
|||||||||||||||||
ных к новым и обратно примут вид: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
x = u −1, |
y = v −1, |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
u = x +1, |
v = y +1. |
|
|
|
|
||||||||||
|
|
|
В результате уравнение приводится к однородному |
|
|
|||||||||||||||||
|
|
|
|
|
dv |
= u −1+ v −1+ 2 |
= u + v =1+ |
v |
. |
|
|
|||||||||||
|
|
|
|
|
du |
|
u −1+1 |
|
|
v |
|
u |
u |
|
|
|
||||||
|
|
|
Полагая далее t = v / u, → |
|
v = tu, |
′ |
=t u + t, |
приходим к уравнению |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
с разделяющимися переменными относительно функции t: |
|
|
||||||||||||||||||||
|
|
|
|
|
′ |
|
|
|
du |
, |
|
t = lnu + lnC = lnCu. |
|
|||||||||
|
|
|
|
|
t u +t =1+t, d t = |
u |
|
|
||||||||||||||
Возвращаясь к старой переменной, получим |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
v |
= lnCu, |
|
|
y +1 |
= lnC(x +1), |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
x +1 |
|
|
|
|
|||||||||
|
|
|
|
|
|
u |
|
|
|
|
|
|
|
|
|
|
|
y = (x +1) lnC(x +1) −1.
1.6. Линейные дифференциальные уравнения
Линейным дифференциальным уравнением первого порядка называется уравнение вида
y′+ P(x) y = Q(x), |
(1.12) |
15
содержащее искомую функцию и ее производную в первой степени. Функции P(x) и Q(x) предполагаются непрерывными.
Если правая часть уравнения Q(x) = 0, то уравнение (1.12) называется линейным однородным уравнением, в противном случае - линейным неоднородным.
Рассмотрим интегрирование этого уравнения методом вариации про-
извольной постоянной (методом Лагранжа). В соответствии с этим мет о-
дом сначала ищется решение соответствующего линейного однородного
уравнения: y′+ P(x) y = 0 . |
Разделяя в нем переменные, получим его общее |
|||
решение в виде |
|
|
|
|
dy |
= −P(x) y, |
dy |
= −P(x)dx, ln y = −∫P(x)dx +lnC, |
|
dx |
|
y |
|
|
|
|
|
y =C e−∫P(x)dx . |
(1.13) |
Далее ищется общее решение исходного уравнения с ненулевой правой частью в виде (1.13), но произвольная постоянная в (1.13) заменяется неизвестной функцией:
y = v(x) e−∫P(x)dx . |
(1.14) |
Подставляя (1.14) в (1.12), получаем дифференциальное уравнение с разделяющимися переменными относительно варьируемой функции v(x), интегрируя которое, находим эту функцию.
В результате, общее решение уравнения (1.12) может быть представлено в виде
y = e−∫P(x)dx [ Q(x)e∫P(x)dx dx + C]. |
(1.15) |
∫ |
|
Формула (1.15) дает общее решение линейного дифференциального |
|
уравнения первого порядка в форме Коши. |
|
Пример. Решить задачу Коши для уравнения y′− 4y = x |
при началь- |
ном условии y(0)=1. |
|
Находим сначала общее решение линейного однородного уравнения y′− 4y = 0 .
Оно имеет вид
16
y = C e−∫P(x)dx = C e−∫−4dx = C e4x .
Заменяем произвольную постоянную C в этом решении неизвестной функцией v(x):
y = v(x)e4x .
Вычисляем производную y′
|
|
|
|
|
|
y |
′ |
|
|
|
|
′ |
|
|
4x |
+ 4v(x)e |
4x |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
= v (x)e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
и подставляя |
y и |
y′ |
в исходное уравнение, получим |
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
′ |
4x |
+ |
4v(x)e |
4x |
− 4v(x)e |
4x |
|
|
|
|
|
|
′ |
|
|
4x |
= x, |
|
||||||||||||||||
|
v (x)e |
|
|
|
|
|
= x, → v (x)e |
|
|
|
|||||||||||||||||||||||||
′ |
|
−4x |
, v(x) |
|
= ∫xe |
−4x |
dx + C |
= − |
1 |
e |
−4x |
|
+ |
1 |
|
+ C. |
|
||||||||||||||||||
v |
(x) = xe |
|
|
|
|
4 |
|
|
x |
4 |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Общее решение уравнения примет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
4x |
|
|
|
|
|
1 |
|
−4x |
|
|
|
1 |
|
|
4x |
|
|
4x |
|
|
|
x +1/ 4 |
|
|||||||
|
y =ν(x)e |
|
= |
− |
|
|
e |
|
x |
+ |
|
+ C e |
|
|
=Ce |
|
|
− |
|
|
. |
||||||||||||||
|
|
|
4 |
|
4 |
|
|
|
|
|
4 |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Находим произвольную постоянную C из начального условия: |
|||||||||||||||||||||||||||||||||||
при x = 0 |
|
|
|
|
C − |
|
|
1 |
|
=1, |
→ |
|
C = |
17. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
16 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Следовательно, решение задачи Коши будет
y =161 [17e4x − (4x +1)].
Решение линейного дифференциального уравнения (1.12) может быть также получено, если искомую функцию представить в виде произведения двух произвольных функций (метод Бернулли [4]):
|
|
y =u (x) v(x). |
(1.16) |
|||
Тогда |
|
y |
′ |
=u v + v u. |
(1.17) |
|
|
|
|
′ |
′ |
|
|
Подставляя (1.16) и (1.17) в (1.12), получим |
(1.18) |
|||||
[u |
′ |
+ P(x)u]v + v u = Q(x). |
||||
|
|
|
|
′ |
|
17
Функцию u(x) подбираем так, чтобы она была одним из решений урав-
нения
u′+ P(x)u = 0 . |
|
После разделения переменных получим |
|
u = e−∫ P(x)dx . |
(1.19) |
Тогда уравнение (1.18) примет вид
v′e−∫P(x)dx = Q(x).
Следовательно,
dv =Q(x)e∫P(x)dx dx.
Интегрируя это уравнение с разделенными переменными, находим функцию v:
v = ∫Q(x)e∫P(x)dx dx + C . |
(1.20) |
Подставляя (1.19) и (1.20) в (1.16), получим общее решение уравнения
(1.12) в виде (1.15).
1.7. Уравнение Бернулли
Уравнением Бернулли называется уравнение вида
y′+ P(x) y = Q(x) yn , |
(1.21) |
где n ≠ 0;1.
Уравнение Бернулли является нелинейным, но оно приводится к линейному следующим преобразованием:
1) Обе части уравнения умножаются на y−n , тогда
y−n y′+ P(x) y1−n = Q(x).
2) Далее применяется подстановка
z = y1−n .
Тогда по правилу дифференцирования сложной функции получим
dz |
= (1−n) y |
−n |
y′, следовательно, |
y |
−n |
y′ = |
1 dz |
|
dx |
|
|
|
dx . |
||||
|
|
1− n |
В результате уравнение становится линейным относительно функции z:
18 |
|
1 |
|
1− n z′ + P(x)z =Q(x). |
(1.22) |
Уравнение (1.22) может быть решено методом вариации произвольной постоянной или методом Бернулли.
Пример. Решить уравнение
|
|
|
|
|
|
|
|
|
|
|
|
y′+ |
y |
lnx |
y |
2 |
. |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
Умножим обе части уравнения на |
|
|
|
y−2 : |
|
|
|
y−2 y′+ |
y−1 |
= |
lnx |
. |
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
x |
x |
|
|
|||||||||||||||||||||||||||||||
Положим |
z = y |
−1 |
, |
|
→ |
|
|
z |
′ |
= −y |
−2 |
|
′ |
и |
уравнение преобразуется в линей- |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
y , |
|||||||||||||||||||||||||||||||
ное: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z′− |
z |
|
|
|
|
|
lnx |
. |
|
|
|
|
|
|
|
|
(1.23) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
= − x |
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
Находим сначала решение соответствующего линейного однородного |
|||||||||||||||||||||||||||||||||||||||||
уравнения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z′ |
− |
z |
= 0. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
dz |
= |
z |
, |
|
|
|
ln |
|
z |
|
= ln |
|
x |
|
+ ln C, |
|
|
z = Cx. |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
dx |
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Решение неоднородного уравнения (1.23) отыскиваем в виде |
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z = v(x) x, |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
тогда |
|
|
|
|
|
|
|
|
|
|
|
z |
′ |
|
= v (x) x + v(x). |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
dv(x) x + v(x) − v(x) x = −lnx , |
|
|
|
|
|
dv(x) x = −lnx |
, |
|
∫dv(x) = −∫ |
lnx |
dx +C. |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
2 |
|||||||||||||||||||||||||||||||||||
dx |
|
|
x |
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
x |
|||||||
После интегрирования получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
v(x) = |
lnx |
+ |
1 |
+ C, |
|
|
|
|
|
|
|
|
z = v(x) x |
lnx |
+ |
1 |
|
|
||||||||||||||||||||||
|
|
x |
x |
|
|
|
|
|
|
|
|
= |
x |
|
x |
+ C x, |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
поэтому общее решение исходного уравнения будет иметь вид