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

|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
97 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
F(p) = |
|
|
|
p − 2 |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
p2 + 6 p +16 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
p − 2 |
= |
|
|
|
p − 2 |
= |
|
|
|
|
p + 3 |
|
|
|
|
|
− |
|
|
|
5 |
|
|
|
|
|
|
= |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(p + 3)2 + ( |
|
)2 |
|
|
|
|
|
( |
|
)2 |
|
|
||||||||||||||||||
p2 + 6 p +16 |
|
(p + 3)2 + ( |
|
)2 |
(p + 3)2 + |
||||||||||||||||||||||||||||||||||||||
|
7 |
7 |
7 |
||||||||||||||||||||||||||||||||||||||||
|
p + 3 |
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
• |
|
|
−3t |
|
|
|
|
|
7 |
|
|
−3t |
|
|
|
|
|||||||||
= |
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
=• |
e |
cos |
|
7 t − |
|
e |
sin 7 t. |
|||||||||||||||||||
(p + 3)2 + ( |
|
|
|
)2 |
|
|
|
|
(p + 3)2 + ( |
|
|
)2 |
|
|
|
|
|
|
7 |
|
|
|
|||||||||||||||||||||
|
|
|
|
7 |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
7 |
|
|
7 |
|
|
|
|
|
|
При практическом применении преобразования Лапласа для нахождения оригинала по его изображению часто пользуются справочными таблицами соответствия между оригиналами и их изображениями [8]. Если необходимая функция в таблице отсутствует, то используют свойства преобразования Лапласа для приведения заданной функции к табличному виду.
3.3. Решение обыкновенных дифференциальных уравнений операционным методом
Пусть дано неоднородное линейное дифференциальное уравнение второго порядка с постоянными коэффициентами
y′′+ a1y′+ a2 y = f (t) |
(3.27) |
||
и заданы начальные условия |
, y (0) |
= y1, |
(3.28) |
y(0) = y0 |
|||
|
′ |
|
|
то есть сформулирована задача Коши.
Операционный метод решения дифференциальных уравнений базируется на том, что искомая функция y(t) и правая часть f (t) рассматриваются как оригиналы, и к уравнению (3.27) применяются теоремы дифференцирования оригинала и линейности.
Пусть y(t)• =• Y (p), |
|
f (t)• =• F(p) , |
тогда |
|
|||
y′(t)• =• pY (p) − y0 , |
|
y′′(t)• =• |
p2Y (p) − py0 − y1. |
(3.29) |
|||
Применяя к уравнению (3.27) теорему линейности, с учетом |
соотно- |
||||||
шений (3.29) получим |
|
|
|
|
|
|
|
p2Y (p)+ a pY (p) |
+ a |
Y (p)−(p + a ) y |
0 |
− y = F(p). |
(3.30) |
||
1 |
|
2 |
1 |
|
1 |
|
98
Алгебраическое уравнение ( 3.30) называется изображающим или операторным уравнением. Разрешая его относительно Y (p), найдем изображение искомого решения
Y (p) = |
F(p) |
|
+ |
(p + a1)y0 + y1 |
. |
|
|
p2 + a p + a |
2 |
|
p2 + a p + a |
2 |
|
|
1 |
|
1 |
|
Далее следует выполнить обратное преобразование Лапласа методами, указанными выше, и найти соответствующий изображению Y (p) оригинал y(t), который и будет решением задачи Коши (3.27)- (3.28).
Замечание. Достоинство операционного метода решения по сравнению с классическим методом решения неоднородных дифференциальных уравнений состоит в том, что начальные условия автоматически (естественным образом в процессе преобразований) входят в изображающее уравнение. Поэтому после выполнения обратного преобразования Лапласа сразу получается частное решение уравнения, удовлетворяющее заданным начальным условиям. Следовательно, при операционном методе не надо искать произвольные постоянные.
Недостаток метода - трудность обращения преобразования Лапласа, особенно в случае сложной правой части или уравнений высокого порядка.
Пример. Решить операционным методом уравнение
|
|
|
y′′+ y = 2cost |
|
y (+0) = −1. |
||||||||||
при заданных начальных условиях |
|
|
y(+0) = 0 , |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
Пусть |
y(t)• =• Y (p), |
cost • =• |
|
|
|
p |
. |
|
|
|
|
|
|||
p2 +1 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
С учетом начальных условий y′(t)• =• |
pY (p), |
y′′(t)• =• p2Y (p) +1. |
|||||||||||||
Изображающее уравнение примет вид |
|
|
|
|
|
|
|
||||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
2 p |
|||
|
p Y (p) +1 |
+Y (p) = |
|
|
. |
||||||||||
|
|
p2 +1 |
|||||||||||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y (p) = |
|
2 p |
|
− |
|
|
1 |
|
= Y (p)−Y (p). |
|||||
|
|
(p2 +1)2 |
|
p2 +1 |
|||||||||||
|
|
|
|
|
1 |
2 |
где

99
Y |
(p) = |
2 p |
|
, |
Y |
(p) = |
1 |
. |
|
|
|
||||||
1 |
|
(p2 +1) |
2 |
|
2 |
|
p2 +1 |
|
|
|
|
|
|
Для оригинала, соответствующего изображению Y1(p), воспользуемся теоремой дифференцирования изображения:
Y1(p) = (p22+p1)2 = − p21+1 ′• =• t sint.
Изображение Y2 (p), является табличным. Ему соответствует оригинал sint . Следовательно, решение задачи Коши примет вид
y(t) =t sint −sint = (t −1)sint.
100
ПРИЛОЖЕНИЕ 1
ВАРИАНТЫ РАСЧЕТНО – ГРАФИЧЕСКОЙ РАБОТЫ
ПО ОБЫКНОВЕННЫМ ДИФФЕРЕНЦИАЛЬНЫМ
УРАВНЕНИЯМ