- •ПРИКЛАДНАЯ МАТЕМАТИКА
- •Е.А.КОГАН
- •Москва 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
53
Это явление совпадения частоты возмущающей силы с частотой собственных колебаний системы называется резонансом.
2.2.3. Метод вариации произвольных постоянных
Если известно общее решение однородного линейного дифференциального уравнения, соответствующего заданному неоднородному, то его частное решение может быть найдено методом вариации произвольных посто-
янных (методом Лагранжа).
Рассмотрим его реализацию для линейных дифференциальных уравнений второго порядка:
y′′+ a1y′+ a2 y = f (x). |
(2.40) |
Сначала находим общее решение соответствующего линейного однородного уравнения L(y) = 0 в виде
y0 = C1y1 +C2 y2 , |
(2.41) |
где C1 и C2 - произвольные постоянные, y1, y2 - частные линейно независимые решения однородного уравнения.
Далее ищем решение неоднородного уравнения (2.40), аналогичное по структуре (2.41), но произвольные постоянные в (2.41) заменя ем неизвестными функциями, а именно принимаем
y =C1(x)y1 + C2 (x)y2 . |
(2.42) |
Подстановка (2.42) в (2.40) приводит к следующей системе дифференциальных уравнений
C′(x)y |
+ C′ |
(x)y |
2 |
= 0, |
(2.43) |
||
1 |
1 |
2 |
|
= f |
(x), |
||
C1′(x)y1′ + C2′(x)y2′ |
|
в которой первое уравнение вводится произвольно (как дополнительное условие, налагаемое на функции C1(x) и C2 (x)). Определитель этой системы - определитель Вронского
W = |
|
y1 |
y2 |
|
≠ 0, |
|
|
||||
|
|
y1′ |
y2′ |
|
|
|
|
|
|
|
|
|
|
|
54 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
так как функции |
y1 |
и |
|
y2 |
линейно независимы. |
|
Поэтому система (2.43), |
||||||||||||||||||
рассматриваемая |
как |
система линейных алгебраических уравнений относи- |
|||||||||||||||||||||||
тельно С1′(x), С2′(x), имеет |
решение и притом единственное. Оно представ- |
||||||||||||||||||||||||
ляется в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C1′(x) = − |
y2 f (x) |
, |
C2′(x) = |
y1 f (x) |
. |
|
|
|
(2.44) |
|||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||
|
|
W |
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W |
|
|
|
|
|
|||||
Интегрируя |
дифференциальные уравнения |
|
|
|
первого порядка (2.44), |
||||||||||||||||||||
находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
С (x) = − |
∫ |
|
y2 f (x) |
dx + C , |
C |
|
(x) = |
|
|
|
y1 f (x) |
dx +C . |
(2.45) |
||||||||||||
|
W |
|
|
∫ |
|
|
|||||||||||||||||||
1 |
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
W |
|
|
|
2 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Подставляя (2.45) |
|
|
в (2.42), получим общее решение неоднородного |
||||||||||||||||||||||
уравнения в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y =C1y1 |
+ C2 y2 − y1 |
∫ |
y2 |
f (x) |
|
dx + y2 |
∫ |
|
y1 f |
(x) |
dx. |
|
(2.46) |
||||||||||||
|
W |
|
W |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Последние два слагаемых в правой части формулы (2.46) определяют частное решение неоднородного уравнения (2.40).
Таким образом, если известна фундаментальная система решений соответствующего линейного однородного уравнения y1, y2 , то общее решение неоднородного уравнения может быть найдено с помощью квадратур (то есть неопредёленных интегралов от этих решений).
Замечание. Метод вариации произвольных постоянных является общим методом, пригодным для построения решения линейных неоднородных дифференциальных уравнений (как с переменными, так и с постоянными коэффициентами) при произвольной непрерывной правой части. Но он является принципиально более сложным, чем метод подбора, так как его реализация связана с интегрированием дифференциальных уравнений первого порядка. Поэтому для линейных неоднородных уравнений с постоянными коэффициентами и с правыми частями специального вида, указанными в таблице 2, проще применять метод подбора частного решения.
Пример. Решить уравнение
|
′′ |
|
′ |
e−2x |
|
|
y |
+ 4y |
+ 4y = x . |
(2.47) |
|||
|
|
55
1). Находим общее решение соответствующего однородного уравнения y′′+ 4y′+ 4y = 0.
Составляем характеристическое уравнение |
k 2 + 4k + 4 = 0. Его корни |
||||||||
k1 = k2 = −2. Следовательно, |
|
частные линейно независимые решения равны |
|||||||
(см. таблицу 1, случай 3а) y |
|
|
= e−2x , |
y |
2 |
= xe−2x , |
а общее решение |
||
1 |
|
|
|
|
|
|
|
||
y |
0 |
= C e−2x +C |
2 |
xe−2x . |
|
||||
|
1 |
|
|
|
|
|
2). Так как правая часть неоднородного уравнения (2.47) не относится ни к одному из рассмотренных в таблице 2 случаев, то частное решение находим методом вариации произвольных постоянных. Принимаем
y =C1(x)e−2x + C2 (x)xe−2x ,
(2.48) где C1(x) и C2 (x) -неизвестные (варьируемые) функции. Тогда их
производные C1′ (x) и C2′(x) могут быть найдены из решения сист |
емы |
(2.43): |
|
C1′(x)e−2x + C2′(x)− 2C1′(x)e−2x + C2′
Определитель этой системы
e−2x W = − 2e−2x
xe−2x |
= 0, |
|
|
|
|
|
||
(x)(e |
−2x |
− 2xe |
−2x |
) = |
e−2x |
. |
||
|
|
|
x |
|||||
|
|
|
|
|
|
|
|
|
xe−2x |
|
|
= e−4x . |
|
||||
|
|
|
||||||
e−2x − 2xe−2x |
|
|
||||||
|
|
|
|
|
Поэтому для определения варьируемых функций согласно (2.44) получаем дифференциальные уравнения вида
C1′(x) = − |
xe−2x e−2x |
= −1, |
C2′(x) = |
e−2xe−2x |
= |
1 |
. |
(2.49) |
|
xe−4x |
xe−4x |
x |
|||||||
|
|
|
|
|
|
||||
Интегрируя уравнения (2.49), находим |
|
|
|
|
|
|
|||
C1(x) = −x +C1, |
C2 (x) = lnx +C2 . |
|
|
(2.50) |
Подставляя (2.50) в (2.48), получим общее решение уравнения в виде
y = e−2x [C1 +C2 x + x(lnx −1)].
56
2.3. Задачи на собственные значения
Рассмотрим линейное однородное дифференциальное уравнение второго порядка
y′′+ λ2 y = 0 |
(2.51) |
c однородными граничными условиями |
(2.52) |
y(0) = 0, y (L) = 0 . |
|
′ |
|
Здесь предполагается, что 0 ≤ x ≤ L, λ −некоторый параметр, имеющий определенный физический смысл.
Задача (2.51), (2.52) является однородной краевой задачей. Особенность ее в том, что она имеет "тривиальное" решение y=0. Но кроме того, имеются еще определенные значения параметра λ , при кото рых задача имеет не равные тождественно нулю (нетривиальные) решения. Такие значения параметра λ называются собственными значениями, а соответствующие им ненулевые решения y(x) называются собственными функциями. Сама же задача отыскания собственных значений и собственных функ-
ций называется задачей на собственные значения или задачей Штурма-
Лиувилля. Эта задача представляет большой интерес для физики и для технических приложений. К ней, например, сводятся многие задачи устойчивости и колебаний упругих систем.
Рассмотрим в качестве примера классическую задачу об устойчивости консольно закрепленного стержня, сжатого продольной силой P. С математической точки зрения эта задача сводится именно к задаче на собственные значения (2.51), (2.52) [15]. При этом параметр λ определяется следующим
образом: |
|
λ2 = P / EI, |
(2.53) |
где EI-изгибная жесткость стержня длиной L. Тривиальное решение y=0 соответствует прямолинейной форме равновесия сжатого стержня. Для определения собственных значений параметра λ следует, во-первых, найти общее решение уравнения (2.51). Характеристическое уравнение, соответствующее уравнению (2.51):
k 2 + λ2 = 0, |
(2.54) |
57
имеет мнимые корни k1 = λi , k2 = −λi . Поэтому общее решение уравнения (2.51) принимает вид
y =C1 cosλx + C2 sinλx. |
(2.55) |
Из первого граничного условия (2.52) следует, что C1 = 0 . Удовлетворяя второму граничному условию, получаем уравнение
C2 λcosλL = 0 . |
|
Из условия существования нетривиального решения C2 |
≠ 0, λ ≠ 0 . |
Следовательно, необходимо принять |
|
cosλL = 0. |
(2.56) |
Из этого тригонометрического уравнения следует λ L =π / 2 + nπ, |
n = 0,1,2,... |
Поэтому собственные значения параметра λ равны |
|
λ = π(1+ 2n). |
(2.57) |
2L |
|
Соответствующие им собственные функции с точностью до постоянного множителя C2 определяются по формуле
y = sin |
π(1+ 2n)x |
. |
(2.58) |
|
|||
|
2L |
|
|
Формулы (2.57), (2.58) дают решение задачи (2.51), |
(2.52) на собст- |
венные значения. Как видно из (2.57), существует бесконечное множество
собственных значений параметра λ . |
Подставляя (2.57) в выражение (2.53), |
||||
получим |
|
|
|
|
|
P = |
π 2 |
(1 |
+ 2n)2 |
EI. |
(2.59) |
|
4L2 |
||||
|
|
|
|
Из формулы (2.59) следует, что существует множество значений сжимающей силы P, при которых возможны различные искривленные формы равновесия продольно сжатого стержня. Наименьшее значение сжимающей силы, соответствующее n=0, называется критическим и равно
P |
= |
π 2 EI |
. |
(2.60) |
|
4L2 |
|||||
кр |
|
|
|