- •ПРИКЛАДНАЯ МАТЕМАТИКА
- •Е.А.КОГАН
- •Москва 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
10
1.3. Дифференциальные уравнения с разделенными и разделяющимися переменными
Дифференциальное уравнение с разделенными переменными имеет
вид
M (x) dx + N (y)dy = 0 . |
(1.7) |
В уравнении с разделенными переменными перед дифференциалом dx стоит функция только одной переменной x, а перед дифференциалом dy стоит функция переменной y. Такие уравнения можно почленно интегрировать. В результате получим
∫M (x)dx + ∫N(y)dy = C
или
F(x)+Ф(y) = C, |
(1.8) |
где
F(x) = ∫M (x)dx, Ф(y) = ∫N(y)dy .
Конечное (не дифференциальное) соотношение (1.8) и является общим интегралом уравнения (1.7).
Пример. Решить уравнение ex dx + ln ydy = 0.
Очевидно, это уравнение с разделенными переменными. Интегрируя его, получим
∫ex dx + ∫ln ydy =C.
Следовательно, общий интеграл уравнения будет
ex + y(ln y −1) = C.
Дифференциальное уравнение вида
M1(x)N1(y)dx + M 2 (x)N2 (y)dy = 0 , |
(1.9) |
в котором коэффициенты при дифференциалах распадаются на множители, каждый из которых зависит только от одной переменной, называется диффе-
ренциальным уравнением с разделяющимися переменными. Уравнение
(1.9) делением обеих частей на произведение функций N1(y)M 2 (x) ≠ 0 приводится к уравнению с разделенными переменными
|
|
|
|
|
11 |
|
|
||
|
|
M1(x) |
dx + |
|
N2 (y) |
dy = 0, |
|||
|
|
|
|
|
|||||
|
|
M 2 (x) |
|
N1(y) |
|||||
общий интеграл которого |
|
|
|
|
|
|
|
|
|
∫ |
|
M1(x) |
dx + |
∫ |
N2 (y) |
dy =C. |
|||
|
|
|
|||||||
|
M 2 (x) |
|
|
N1(y) |
Пример. Решить равнение (1+ y2 )xdx + (1+ x2 ) ydy = 0 . Разделяем переменные делением на выражение (1+ y2 )(1+ x2 ) ≠ 0 :
1+xx2 dx +1+yy2 dy = 0.
Интегрируем полученное уравнение с разделенными переменными
∫1+xx2 dx + ∫1+yy2 dy = C1.
Тогда
12 ln(1+ x2 )+ 12 ln(1+ y2 ) = C1
Так как C1 - произвольная постоянная, принимая ее для упрощения полученного выражения в виде
C1 = 12 lnC,
представим общий интеграл уравнения в виде (1+ x2 )(1+ y2 ) = C.
1.4. Однородные дифференциальные уравнения
Если уравнение y′ = f (x, y) или M (x, y)dx + N (x, y)dy = 0 не меняется при замене x на kx, а y на ky, то оно называется однородным.
Однородное дифференциальное уравнение подстановкой t = xy
приводится к уравнению с разделенными переменными. Пример. Решить уравнение
(x + y)dx + (y − x)dy = 0.
|
|
|
|
|
|
|
|
|
|
|
|
12 |
|
|
|
|
|
|
|
||
Преобразуем уравнение к виду |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
y′ = |
|
x + y |
. |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
x − y |
|
|
|
|
|
||||||
Так как |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
kx + ky |
= |
x + y |
, то исходное уравнение однородное. |
|||||||||||||||||
|
kx −ky |
|
|||||||||||||||||||
|
|
x − y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Полагаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t = y / x, → y = tx и |
y |
′ |
|
′ |
|
|||||||||||
|
|
|
|
|
|
= t x + t. |
|
||||||||||||||
Тогда уравнение примет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
′ |
x + tx 1+ t |
|
|
dt |
|
1+t |
|
|
1+t 2 |
|
1+t 2 |
||||||||||
|
|
1−t |
или |
x dx |
= |
1−t −t = |
1−t |
, |
xdt = 1−t dx. |
||||||||||||
t x + t = x −tx = |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1+ t |
2 |
|
|
|
|
|
|
|||
Разделив обе части уравнения на |
|
|
|
|
приходим к уравнению с разде- |
||||||||||||||||
x |
|
1− t |
, |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ленными переменными |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
1−t |
|
|
|
dt = dx . |
|
|
|
|
|||||
|
|
|
|
|
|
|
1+t 2 |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
Интегрируя его, находим
arctg t − 12 ln (1+t 2 ) = lnx +C
или
arctg t −ln(x 1+t 2 ) = C.
Возвращаясь к старой переменной, получим общий интеграл исходного уравнения в виде
C = arctg xy − lnx2 + y2 .
1.5. Дифференциальные уравнения, приводящиеся к однородным
Дифференциальное уравнение вида
|
′ |
|
|
a1x +b1y + c1 |
|
||||
y |
= |
|
|
||||||
|
|||||||||
|
f a |
2 |
x +b y + c |
2 |
|
||||
|
|
|
|
|
2 |
|
13
называется приводящимся к однородному. В частности, |
к этому классу |
|||||||||||
относится уравнение вида |
|
|
|
|
|
|
|
|
|
|||
y′ = |
|
|
a1x +b1y +c1 |
(1.10) |
||||||||
|
|
|
|
|
|
|
|
. |
||||
|
a |
2 |
x +b y + c |
2 |
||||||||
|
|
|
|
|
|
2 |
|
|
||||
Некоторые из коэффициентов (но не одновременно c1 |
и c2 ) могут быть |
|||||||||||
равны нулю. |
|
|
|
|
|
|
|
|
|
|||
Следует различать два случая: |
|
|
|
|
|
|
|
|
||||
1). Если определитель δ = |
|
a1 |
b1 |
|
≠ 0 , |
|
то уравнение (1.10) приводит- |
|||||
|
|
|
||||||||||
|
|
|
a |
2 |
b |
|
|
|
|
|
||
|
|
|
|
|
2 |
|
|
|
|
|
||
ся к однородному подстановкой |
|
|
|
|
|
|
|
|
|
|||
x =u +α, |
|
y = v + β , |
(1.11) |
|||||||||
где постоянные a и β определяются из системы уравнений: |
|
a1α +b1β + c1 = 0,a2α +b2 β + c2 = 0.
Действительно, учитывая, что du = dx, dv = dy, следовательно, dydx = dudv , и
подставляя (1.11) в (1.10), получим однородное уравнение относительно новой функции v(u).
dv |
= |
a1u +b1v |
. |
|
du |
|
a |
u +b v |
|
|
|
2 |
2 |
|
Полагая далее
t = uv ,
приводим последнее уравнение к уравнению с разделяющимися переменными.
2) Если определитель
δ = |
|
a1 |
b1 |
|
= 0, |
|
|
|
|||||
|
|
a |
2 |
b |
|
|
|
|
|
2 |
|
|
то уравнение (1.10) сразу приводится к уравнению с разделенными переменными заменой