- •ПРИКЛАДНАЯ МАТЕМАТИКА
- •Е.А.КОГАН
- •Москва 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
37
Уравнение (2.26) содержит в правой части известную функцию от x и относится к первому из рассмотренных нами типов (см. (2.2)). Интегрируя его последовательно три раза, окончательно получим общее решение исходного уравнения, содержащее 4 произвольных постоянных :
|
|
y′′ = |
1 +C1 lnx +C2 , |
|
|
|
|
|
|
||||||
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
y′ = lnx +C1x (lnx −1)+C2 x +C3 , |
|
|
|
|
||||||||||
|
|
|
x2 |
3 |
|
|
|
x2 |
|
|
|
|
|
||
|
y = x(lnx −1) |
+ C |
|
|
lnx − |
|
|
+ C |
|
|
+ C |
|
x + C |
|
. |
|
|
|
2 |
|
2 |
|
|
||||||||
|
|
1 2 |
|
|
|
|
2 |
|
3 |
|
4 |
|
|||
Пример. Решить уравнение y′2 + 2yy′′ = 0 |
при следующих начальных |
||||||||||||||
условиях: |
y(1) =1, y (1) =1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Уравнение не содержит явно независимой переменной x и потому от-
носится к третьему из перечисленных типов. Принимаем y (x) = z(y), тогда с |
||||||
|
|
|
|
|
|
′ |
учетом (2.15) уравнение примет вид |
|
|
|
|
|
|
z2 + 2 yz dz |
= 0. |
|
|
|
|
|
|
dy |
|
|
|
|
|
Сокращаем на z ≠ 0 ( решение z = 0 , |
то есть y′ = 0, y = C нужно исследо- |
|||||
вать отдельно ). |
|
|
|
|
|
|
z + 2 y dz = 0 . |
|
|
|
|||
|
dy |
|
|
|
|
|
Это уравнение с разделяющимися переменными: |
|
|
|
|||
2y dz = −z, |
dz |
= − dy . |
|
|
|
|
dy |
z |
2 y |
|
|
|
|
Интегрируя, получим |
|
|
|
|
|
|
lnz = −1ln y + lnC = ln |
C1 |
|
. |
|||
|
|
|
||||
2 |
|
1 |
|
y |
||
|
|
|
Отсюда следует
z = C1y .
и с учетом (2.14)
38 |
|
|
|
|
|
|
|
y |
′ |
|
C1 |
|
(2.27) |
||
|
|
|
|
||||
= |
|
y . |
|||||
|
|
Определим сразу произвольную постоянную C1. Подставляя (2.27) во второе начальное условие, получим С1 =1. Тогда
dydx = 1y .
Разделяя переменные, имеем y dy = dx, следовательно, после интегрирования получим
2 y3/ 2 |
= x +C2 . |
|
|
|
3 |
|
|
|
|
Находим C2 из первого начального условия: |
2 |
1=1+C2 , |
следовательно |
|
|
|
3 |
|
|
C2 = −1/ 3. Поэтому
23 y3/ 2 = x − 13.
Окончательное решение задачи Коши запишется в виде
|
3x −1 2 / 3 |
|
|
y = |
2 |
. |
|
|
|
|
|
2.2. Решение линейных неоднородных |
|
||
дифференциальных уравнений n-го порядка с |
|
||
постоянными коэффициентами |
|
||
Линейным дифференциальным уравнением n - го порядка с |
постоян- |
||
ными коэффициентами называется уравнение вида |
|
||
y(n) + a1y(n−1) + a2 y(n−2) |
+ + an−1y′+ an y = f (x), |
(2.28) |
содержащее неизвестную функцию и все её производные до n-го порядка включительно в первой степени. В (2.28) коэффициенты
(i =1,2,...,n). В общем случае, если хотя бы один из коэффициентов ai является функцией от x, то уравнение (2.28) будет уравнением с переменными коэффициентами.
39
Левая часть уравнения (2.28) называется линейным дифференциальным оператором и обозначается через
L(y) = y(n) + a1y(n−1) + a2 y(n−2) + + an−1y′+ an y.
Поэтому в компактной форме уравнение (2.28) запишется так:
L(y) = f (x).
Если правая часть уравнения f (x) = 0 , то уравнение (2.28) принимает вид
L(y) = 0 |
(2.29) |
и называется линейным однородным дифференциальным уравнением, в противном случае – неоднородным.
Справедлива теорема (о структуре общего решения линейного неоднородного дифференциального уравнения n-го порядка): общее реше-
ние линейного неоднородного дифференциального уравнения n – го порядка представляется в виде суммы
y = y0 + y , |
(2.30) |
где y0 - общее решение соответствующего линейного однородного уравненияL(y) = 0, y -частное решение неоднородного уравнения L(y) = f (x).
Полезна также следующая теорема, называемая принципом суперпо-
зиции:
Общее решение линейного неоднородного дифференциального уравнения
m
L(y) = ∑αk fk (x)
k =1
может быть представлено в виде
m
y = y0 + ∑αk yk ,
k =1
где y0 - общее решение соответствующего однородного уравнения L(y) = 0, αk = const, a yk - частные решения неоднородных уравнений вида
L(y) = fk (x), (k =1, 2, , m).
40
Из принципа суперпозиции следует, что решение исходного уравнения можно свести к решению нескольких более простых уравнений L(y) = fk (x). С физической точки зрения это означает, что результат сложного внешнего воздействия на некоторую систему (объект), характеризуемого
|
m |
функцией |
f (x) = ∑αk fk (x), можно представить как суперпозицию резуль- |
|
k =1 |
татов отдельных элементарных воздействий.
2.2.1. Построение общего решения линейного однородного дифференциального уравнения n –го порядка
Теорема (о структуре общего решения линейного однородного дифференциального уравнения n-го порядка): Общее решение линейного
однородного |
дифференциального уравнения |
n –го порядка |
L(y) = 0 пред- |
|
ставляется в виде линейной комбинации |
n |
линейно независимых частных |
||
решений этого уравнения: |
|
|
|
|
|
n |
|
|
|
|
y = ∑Ck yk = C1y1 |
+C2 y2 + +Cn yn . |
(2.31) |
|
|
k =1 |
|
|
|
Здесь |
C1, C2 , , Cn - произвольные постоянные (заметим, что ли- |
нейной комбинацией функций называется сумма произведений функций на различные постоянные числа, то есть выражение вида
С1y1 +C2 y2 + +Cn yn ); y1, y2 , , yn - частные линейно независимые ре-
шения однородного уравнения L(y) = 0 - такие решения, для которых, составленный из них определитель не равен нулю:
|
y1 |
y2 ... |
yn |
|
|
′ |
′ |
′ |
|
|
y1 |
y2 ... |
yn |
|
W = |
y1′′ |
y2′′ ... |
yn′′ |
≠ 0. |
|
....................................... |
|
||
|
y(n−1) |
y(n−1)... |
y(n−1) |
|
|
1 |
2 |
n |
|
Определитель W называется определителем Вронского или врон-
скианом.