
- •ПРИКЛАДНАЯ МАТЕМАТИКА
- •Е.А.КОГАН
- •Москва 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

58
При этом значении силы Pкр первоначальная прямолинейная форма равно-
весия сжатого стержня становится неустойчивой.
Пример. Найти собственные значения и собственные функции однородной краевой задачи:
|
y′′ = λy, |
y(0) = y(l) = 0, |
l > 0. |
|
Если λ = 0, |
то y =C1x + C2 , |
и граничным условиям удовлетворяет |
||
только тривиальное решение |
y ≡ 0. |
Следовательно, |
λ = 0 не является соб- |
|
ственным значением. |
|
|
|
|
При λ > 0 |
общее решение уравнения y′′ = λy |
будет: |
y =C1eλx + C2e−
λx .
Подставляя это решение в граничные условия, убеждаемся в том, что они удовлетворяются только при y ≡ 0.
Если же λ < 0, то общее решение уравнения имеет вид
y =C1 sin− λx + C2 cos
− λx.
Подставляя это выражение в граничные условия, получим
C2 = 0,
C1 sin− λl = 0.
Чтобы краевая задача имела ненулевые решения, необходимо принять
|
|
|
|
|
|
nπ |
2 |
|
|||
sin − λl = 0. Тогда |
− λl = nπ |
и |
, |
n =1,2,... |
|||||||
λ = − |
l |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
Эти значения λ и являются собственными. Соответствующие им собственные функции с точностью до множителя равны
y =sin |
nπx |
, |
n =1,2,... |
|
l |
||||
|
|
|
Пример. Найти собственные значения и собственные функции однородной краевой задачи для уравнения (2.51) при граничных условиях
y(0) − y(1) = 0, y′(1) = 0.

59
Общее решение уравнения согласно (2.55) имеет вид
Вычисляем |
y = C1 cos λx +C2 sinλx. |
|
y (x) = −C1λsinλx +C2λcos λx, |
||
|
′ |
|
|
y(0) = C1, |
y(1) = C1 cos λ +C2 sinλ, |
|
′ |
|
|
y (1) = −C1λsinλ +C2λcos λ. |
Подставляя эти величины в граничные условия, после элементарных преобразований получим систему линейных однородных алгебраических уравнений относительно C1 и C2 :
(1− cosλ) C1 −sinλ C2 = 0,
− λsinλ C1 + λcosλ C2 = 0.
Ненулевое решение ее существует тогда и только тогда, когда определитель системы равен нулю:
1− cosλ |
−sinλ |
|
= 0. |
||
− λsinλ |
λcosλ |
|
Раскрывая этот определитель, получим уравнение относительно параметра
λ :
λ(cosλ −1) = 0,
корни которого являются собственными значениями задачи:
λ = 0, |
|
|
λ = 2πn, |
n =1,2,... |
|
Легко убедиться, что при λ = 0 y = const. Подставляя λ = 2πn, |
n =1,2,... |
в систему линейных однородных алгебраических уравнений относительно C1 и C2 , находим C2 = 0. Соответствующие собственные функции с точностью до множителя будут y = cos 2πnx.
2.4. Дифференциальные уравнения с переменными коэффициентами
Если коэффициенты ai (x), (i =1,2,...,n) линейного неоднородного уравнения
60 |
|
y(n) + a1(x) y(n−1) + a2 (x) y(n−2) + an−1(x) y′+ an (x) y = f (x) |
(2.61) |
и его правая часть f(x) представляют собой функции, которые определены и непрерывны на заданном интервале, то рассмотренные выше теоремы о структуре общего решения линейного неоднородного уравнения и соответствующего ему однородного уравнения остаются справедливыми. Остаются в силе также принцип суперпозиции решений для неоднородного уравнения и метод вариации произвольных постоянных. Но при этом нельзя искать фундаментальную систему решений однородного уравнения рассмотренным выше методом.
Одним из простейших линейных дифференциальных уравнений с переменными коэффициентами является уравнение Эйлера.
2.4.1. Уравнение Эйлера
Уравнением Эйлера называется линейное уравнение вида
(ax +b)n y(n) + a (ax |
+b)n−1 y(n−1) + a |
2 |
(ax +b)n−2 |
y(n−2) + + |
|
1 |
|
|
|
(2.62) |
|
+ an−1(ax +b)y′+ an y |
= f (x), |
|
|||
|
|
где a, b, a1, a2 , , an−1, an - константы. В наиболее распространенном
случае, при a =1, b = 0 |
уравнение Эйлера имеет вид |
xn y(n) + a1xn−1y(n−1) |
+ a2 xn−2 y(n−2) + + an−1xy′+ an y = f (x). (2.63) |
Как видно, уравнение Эйлера является уравнением с переменными коэффициентами специального вида, но оно приводится к уравнению с постоянными коэффициентами заменой независимой переменной x.
Для уравнения (2.63), полагая
x = et , (x > 0) → t = ln x, |
(2.64) |
находим производные разного порядка от функции y по новой переменной t:
|
′ |
|
dy |
|
|
dy dt |
|
−t |
dy |
|
|
|
|
|
−2t d 2 y |
dy |
|
||||||||||
y |
= |
|
= |
|
|
|
|
|
|
= e |
|
|
|
|
|
, |
|
y |
′′ |
= e |
|
|
2 − |
|
, |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
dx |
dt dx |
|
|
dt |
|
|
|
dt |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|||||||||||
|
|
|
|
|
|
|
3 |
y |
|
|
d |
2 |
y |
|
dy |
|
|
|
|
|
|
|
|
||||
|
y |
′′′ |
= e |
−3t |
d |
|
−3 |
|
+ 2 |
|
|
|
|
|
|
|
(2.65) |
||||||||||
|
|
|
|
3 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
dt |
|
dt |
, |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
…………………………………

61
и подставляя (2.65) в (2.63), получим уравнение с постоянными коэффициентами относительно функции y(t).
В общем случае, для уравнения(2.62) формулы перехода (2.65) примут
вид:
|
′ |
|
−t dy |
|
|
′′ |
2 |
−2t d 2 y |
dy |
|
||
y |
== ae |
|
, |
y |
= a e |
|
|
2 − |
|
, |
||
|
dt |
|
|
dt |
|
|||||||
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
d |
3 |
y |
|
d |
2 |
y |
|
dy |
|
y |
′′′ |
= a |
3 |
e |
−3t |
|
−3 |
|
+ 2 |
|
|||||
|
|
|
3 |
|
|
2 |
|
||||||||
|
|
|
dt |
dt |
dt |
, |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
…………………………………..
Пример. Решить уравнение
x2 y′′+3xy′+ y = 2 + x. |
(2.66) |
Уравнение (2.66) есть уравнение Эйлера второго порядка. Применим замену независимой переменной и производных по формулам (2.64), (2.65). Тогда уравнение примет вид
|
2t |
|
−2t d 2 y |
dy |
t −t dy |
|
t |
||||
e |
|
e |
|
|
|
|
|
+3e e |
|
+ y = 2 + e |
|
|
|
2 − |
|
|
|||||||
|
|
dt |
|
dt |
|
||||||
|
|
|
|
|
|
dt |
|
|
|
и после упрощений получим линейное неоднородное уравнение с постоянными коэффициентами
d 2 y |
+ 2 dy |
+ y = 2 +et . |
(2.67) |
|
dt 2 |
||||
dt |
|
|
Его общее решение
y = y0 + y1 + y2 ,
где y0 - общее решение соответствующего однородного уравнения: y0 = C1e−t +C2 t e−t ,
y1 - частное решение неоднородного уравнения
y2 - частное решение неоднородного уравнения
d 22y + 2 dy + y = 2 , dt dt
d 22y + 2 dy + y = et . dt dt