- •ПРИКЛАДНАЯ МАТЕМАТИКА
- •Е.А.КОГАН
- •Москва 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
33 |
|
|
|
Например, для уравнения второго порядка |
y |
′′ |
= f (x, y, y ) при |
|
|
′ |
0 ≤ x ≤ L граничные условия могут иметь вид: y(0) = y(L) = 0 или y(0) = 0, y′(L) = 0.
В отличие от задачи Коши, решение которой существует и единс т- венно (при некоторых весьма общих условиях, налагаемых на правую часть уравнения (2.1)), краевая задача может не иметь решения или решение может быть не единственным.
2.1.Интегрирование дифференциальных уравнений n – го порядка методом понижения порядка
Если правая часть уравнения (2.1) является известной непрерывной
функцией |
f (x) |
только одной переменной x |
или не содержит искомую |
|||||
функцию y: |
|
′ ′′ |
y |
(n−1) |
), |
или не содержит явно независимую пере- |
||
f (x, y , y , , |
|
|
||||||
менную x: |
′ |
′′ |
(n−1) |
), |
то для решения уравнения (2.1) может быть |
|||
f (y, y , y , , y |
|
|
||||||
применен метод понижения порядка. |
|
|||||||
1. |
|
|
|
|
|
y(n) (x) = f (x). |
(2.2) |
Это уравнение интегрируется непосредственно n раз. При каждом интегрировании порядок уравнения понижается на единицу, и появляется произвольная постоянная. В результате общее решение уравнения будет иметь вид
y = ∫∫ ∫ f (x)dx +C1xn−1 +C2 xn−2 + +Cn−1x +Cn , (2.3)
где в правой части - n-кратный интеграл от функции f(x) и многочлен (n-1)-ой степени от x, коэффициентами которого являются n произвольных постоянных.
2. |
y |
(n) |
′ |
′′ |
, , y |
(n−1) |
). |
(2.4) |
|
= f (x, y , y |
|
|
Уравнение (2.4) не содержит искомой функции y(x). Рассмотрим процедуру интегрирования уравнения данного типа на примере уравнения второго порядка
y |
′′ |
= f (x, y ). |
(2.5) |
|
′ |
|
Понижение порядка достигается подстановкой
|
34 |
|
y (x) = z(x). |
(2.6) |
|
′ |
|
|
Тогда |
|
(2.7) |
y (x) = z (x), |
||
′′ |
′ |
|
и уравнение (2.5) приводится к уравнению первого порядка относительно функции z(x):
z (x) = f (x,z). |
(2.8) |
′ |
|
Интегрируя уравнение (2.8), находим его общий интеграл в виде
z(x) =ϕ(x,C1), |
(2.9) |
где С1 - произвольная постоянная. Далее в (2.9) заменяем левую часть согласно (2.6) и вновь получаем уравнение первого порядка относительно искомой функции y:
y (x) =ϕ(x,C1), |
(2.10) |
′ |
|
Интегрируя уравнение (2.10), находим общее решение исходного уравнения (2.4) в виде
y(x) =ξ(x,C1,C2 ). |
(2.11) |
Замечание. Если уравнение (2.4) не содержит ни искомой функции y, ни её производных до ( k-1) - го порядка включительно, то есть имеет вид
y(n) = f (x, y(k ), y(k +1), , y(n−1) ),
то его порядок может быть понижен сразу на k единиц подстановкой y(k ) = z(x).
3. |
y |
(n) |
′ |
′′ |
(n−1) |
). |
(2.12) |
|
= f (y, y , y |
, , y |
|
Это уравнение не содержит явно независимой переменной x. В частном случае уравнение второго порядка данного типа будет
y |
′′ |
= f (y, y ). |
(2.13) |
|
′ |
|
Понижение порядка достигается подстановкой
y (x) = z(y). |
(2.14) |
′ |
|
Тогда по правилу вычисления производной от сложной функции
|
|
35 |
|
|
|
′′ |
dz dy |
dz |
|
y |
(x) = dy dx |
= dy z . |
(2.15) |
Поэтому уравнение (2.13) при водится к уравнению первого порядка относительно функции z(y):
z dz = f (y,z). |
(2.16) |
dy |
|
Интегрируя уравнение (2.16), находим его общее решение в виде |
|
z(y) =ϕ(y,C1), |
(2.17) |
где C1 - произвольная постоянная. Далее в (2.17) заменяем левую часть согласно (2.14) и вновь приходим к уравнению первого порядка с разделяющимися переменными относительно искомой функции y : y′(x) =ϕ (y,C1), следовательно,
dy |
= dx. |
(2.18) |
|
ϕ(y,C ) |
|||
|
|
||
1 |
|
|
Интегрируя уравнение (2.18), окончательно получим общий интеграл исходного уравнения (2.13) в виде
|
|
|
y(x) =ξ(x,C1,C2 ). |
(2.19) |
|||
Пример. |
Решить уравнение x3 y I V |
+ x2 y′′′ =1. |
|
||||
Это уравнение |
не содержит явно искомой функции y(x) и её первых |
||||||
производных y , y |
′′ |
и относится ко второму из рассмотренных нами типов. |
|||||
′ |
|
|
|
|
|
|
|
Применяя подстановку |
|
|
|
(2.20) |
|||
|
|
|
y (x) = z(x), |
||||
|
|
|
′′′ |
|
|
|
|
получаем линейное неоднородное уравнение первого порядка |
|
||||||
|
|
|
x3 z′+ x2 z =1 |
|
|||
или |
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
|
|
|
|
|
z′ + x z = |
|
. |
(2.21) |
|
|
|
|
x3 |
Интегрируем его методом вариации произвольной постоянной. Сначала решаем соответствующее однородное уравнение
|
36 |
|
z′ + |
1 |
|
x z = 0. |
(2.22) |
Разделяем в нем переменные:
dzz = − dxx .
После интегрирования получим
lnz = −lnx + lnC1 = lnCx1 .
Следовательно, общее решение уравнения (2.22) будет
z = |
C1 |
. |
(2.23) |
|
|||
|
x |
|
Далее ищем решение уравнения (2.21) в форме, аналогичной по структуре выражению (2.23), но произвольную постоянную в (2.23) заменяем неизвестной функцией v(x):
|
|
|
|
|
z = |
v(x) |
. |
|
|
(2.24) |
||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
x |
|
|
|||
Подставляя (2.24) в (2.21), получим |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
′ |
+ v(x) = |
1 |
, |
||||||
|
|
|
|
|
v (x)x − v(x) |
|||||||||
|
|
|
|
|
x2 |
|
|
x2 |
x3 |
|
||||
откуда следует |
|
|
|
|
|
|
|
|
||||||
|
dv(x) |
= |
1 |
|
и, после разделения переменных, |
dv(x) = dx . |
||||||||
|
dx |
x2 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
x2 |
|||
|
|
Интегрируя это уравнение, находим v(x) = −1/ x + C1. Поэтому согласно |
||||||||||||
(2.24) имеем |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
z = − |
|
1 |
+ |
C1 |
. |
|
(2.25) |
||
|
|
|
|
|
|
x2 |
|
|
||||||
|
|
|
|
|
|
|
|
x |
|
|
Заменяя в выражении (2.25) z по формуле (2.20), приходим к уравнению третьего порядка относительно искомой функции y(x):
y |
′′′ |
|
1 |
|
C1 |
|
|
|
= − x2 |
+ |
x . |
(2.26) |
|||||
|