Поможем написать учебную работу
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Ларин Александр Александрович “Курс высшей математики. Часть 3.”
Содержание КВМ Часть 3.
mailto: aalar@yandex.ru
В Ы С Ш Е Й
М А Т Е М А Т И К И
Краткий курс лекций
ЧАСТЬ 3
2000
Обыкновенные дифференциальные уравнения.
Решение различных геометрических, физических и инженерных задач часто приводят к уравнениям, которые связывают независимые переменные, характеризующие ту ил иную задачу, с какой либо функцией этих переменных и производными этой функции различных порядков.
В качестве примера можно рассмотреть простейший случай равноускоренного движения материальной точки.
Известно, что перемещение материальной точки при равноускоренном движении является функцией времени и выражается по формуле:
В свою очередь ускорение a является производной по времени t от скорости V, которая также является производной по времени t от перемещения S. Т.е.
Тогда получаем: - уравнение связывает функцию f(t) с независимой переменной t и производной второго порядка функции f(t).
Определение. Дифференциальным уравнением называется уравнение, связывающее независимые переменные, их функции и производные (или дифференциалы) этой функции.
Определение. Если дифференциальное уравнение имеет одну независимую переменную, то оно называется обыкновенным дифференциальным уравнением, если же независимых переменных две или более, то такое дифференциальное уравнение называется дифференциальным уравнением в частных производных.
Определение. Наивысший порядок производных, входящих в уравнение, называется порядком дифференциального уравнения.
Пример.
- обыкновенное дифференциальное уравнение 1 го порядка. В общем виде записывается .
- обыкновенное дифференциальное уравнение 2 го порядка. В общем виде записывается
- дифференциальное уравнение в частных производных первого порядка.
Определение. Общим решением дифференциального уравнения называется такая дифференцируемая функция y = (x, C), которая при подстановке в исходное уравнение вместо неизвестной функции обращает уравнение в тождество.
Свойства общего решения.
1) Т.к. постоянная С произвольная величина, то вообще говоря дифференциальное уравнение имеет бесконечное множество решений.
2) При каких- либо начальных условиях х = х0, у(х0) = у0 существует такое значение С = С0, при котором решением дифференциального уравнения является функция у = (х, С0).
Определение. Решение вида у = (х, С0) называется частным решением дифференциального уравнения.
Определение. Задачей Коши (Огюстен Луи Коши (1789-1857)- французский математик) называется нахождение любого частного решения дифференциального уравнения вида у = (х, С0), удовлетворяющего начальным условиям у(х0) = у0.
Теорема Коши. (теорема о существовании и единственности решения дифференциального уравнения 1- го порядка)
Если функция f(x, y) непрерывна в некоторой области D в плоскости XOY и имеет в этой области непрерывную частную производную , то какова бы не была точка (х0, у0) в области D, существует единственное решение уравнения , определенное в некотором интервале, содержащем точку х0, принимающее при х = х0 значение (х0) = у0, т.е. существует единственное решение дифференциального уравнения.
Определение. Интегралом дифференциального уравнения называется любое уравнение, не содержащее производных, для которого данное дифференциальное уравнение является следствием.
Пример. Найти общее решение дифференциального уравнения .
Общее решение дифференциального уравнения ищется с помощью интегрирования левой и правой частей уравнения, которое предварительно преобразовано следующим образом:
Теперь интегрируем:
- это общее решение исходного дифференциального уравнения.
Допустим, заданы некоторые начальные условия: x0 = 1; y0 = 2, тогда имеем
При подстановке полученного значения постоянной в общее решение получаем частное решение при заданных начальных условиях (решение задачи Коши).
Определение. Интегральной кривой называется график y = (x) решения дифференциального уравнения на плоскости ХОY.
Определение. Особым решением дифференциального уравнения называется такое решение, во всех точках которого условие единственности Коши (см. Теорема Коши. ) не выполняется, т.е. в окрестности некоторой точки (х, у) существует не менее двух интегральных кривых.
Особые решения не зависят от постоянной С.
Особые решения нельзя получить из общего решения ни при каких значениях постоянной С. Если построить семейство интегральных кривых дифференциального уравнения, то особое решение будет изображаться линией, которая в каждой своей точке касается по крайней мере одной интегральной кривой.
Отметим, что не каждое дифференциальное уравнение имеет особые решения.
Пример. Найти общее решение дифференциального уравнения: Найти особое решение, если оно существует.
Данное дифференциальное уравнение имеет также особое решение у = 0. Это решение невозможно получить из общего, однако при подстановке в исходное уравнение получаем тождество. Мнение, что решение y = 0 можно получить из общего решения при С1 = 0 ошибочно, ведь C1 = eC 0.
Далее рассмотрим подробнее приемы и методы, которые используются при решении дифференциальных уравнений различных типов.
Дифференциальные уравнения первого порядка.
Определение. Дифференциальным уравнением первого порядка называется соотношение, связывающее функцию, ее первую производную и независимую переменную, т.е. соотношение вида:
Если такое соотношение преобразовать к виду то это дифференциальное уравнение первого порядка будет называться уравнением, разрешенным относительно производной.
Преобразуем такое выражение далее:
Функцию f(x,y) представим в виде: тогда при подстановке в полученное выше уравнение имеем:
Далее рассмотрим подробнее типы уравнений первого порядка и методы их решения.
Уравнения вида y = f(x).
Пусть функция f(x) определена и непрерывна на некотором интервале
a < x < b. В таком случае все решения данного дифференциального уравнения находятся как . Если заданы начальные условия х0 и у0, то можно определить постоянную С.
Определение. Дифференциальное уравнение называется уравнением с разделяющимися переменными, если его можно записать в виде
.
Такое уравнение можно представить также в виде:
Перейдем к новым обозначениям
Получаем:
После нахождения соответствующих интегралов получается общее решение дифференциального уравнения с разделяющимися переменными.
Если заданы начальные условия, то при их подстановке в общее решение находится постоянная величина С, а, соответственно, и частное решение.
Пример. Найти общее решение дифференциального уравнения:
Интеграл, стоящий в левой части, берется по частям (см. Интегрирование по частям.):
Чтобы проверить правильность полученного ответа продифференцируем его по переменной х.
- верно
Пример. Найти решение дифференциального уравнения при условии у(2) = 1.
при у(2) = 1 получаем
Итого: или - частное решение;
Проверка: , итого
- верно.
Пример. Решить уравнение
- общий интеграл
- общее решение
Пример. Решить уравнение
Пример. Решить уравнение при условии у(1) = 0.
Интеграл, стоящий в левой части будем брать по частям (см. Интегрирование по частям. ).
Если у(1) = 0, то
Итого, частный интеграл: .
Пример. Решить уравнение .
Для нахождения интеграла, стоящего в левой части уравнения см. Таблица основных интегралов. п.16. Получаем общий интеграл:
Пример. Решить уравнение
Преобразуем заданное уравнение:
Получили общий интеграл данного дифференциального уравнения. Если из этого соотношения выразить искомую функцию у, то получим общее решение.
Пример. Решить уравнение .
; ;
Допустим, заданы некоторые начальные условия х0 и у0. Тогда:
Получаем частное решение
Однородные уравнения.
Определение. Функция f(x, y) называется однородной n го измерения относительно своих аргументов х и у, если для любого значения параметра t (кроме нуля) выполняется тождество:
Пример. Является ли однородной функция
Таким образом, функция f(x, y) является однородной 3- го порядка.
Определение. Дифференциальное уравнение вида называется однородным, если его правая часть f(x, y) есть однородная функция нулевого измерения относительно своих аргументов.
Любое уравнение вида является однородным, если функции P(x, y) и Q(x, y) однородные функции одинакового измерения.
Решение любого однородного уравнения основано на приведении этого уравнения к уравнению с разделяющимися переменными.
Рассмотрим однородное уравнение
Т.к. функция f(x, y) однородная нулевого измерения, то можно записать:
Т.к. параметр t вообще говоря произвольный, предположим, что . Получаем:
Правая часть полученного равенства зависит фактически только от одного аргумента , т.е.
Исходное дифференциальное уравнение таким образом можно записать в виде:
Далее заменяем y = ux, .
таким образом, получили уравнение с разделяющимися переменными относительно неизвестной функции u.
Далее, заменив вспомогательную функцию u на ее выражение через х и у и найдя интегралы, получим общее решение однородного дифференциального уравнения.
Пример. Решить уравнение .
Введем вспомогательную функцию u.
.
Отметим, что введенная нами функция u всегда положительна, т.к. в противном случае теряет смысл исходное дифференциальное уравнение, содержащее .
Подставляем в исходное уравнение:
Разделяем переменные:
Интегрируя, получаем:
Переходя от вспомогательной функции обратно к функции у, получаем общее решение:
Уравнения, приводящиеся к однородным.
Кроме уравнений, описанных выше, существует класс уравнений, которые с помощью определенных подстановок могут приведены к однородным.
Это уравнения вида .
Если определитель то переменные могут быть разделены подстановкой
где и - решения системы уравнений
Пример. Решить уравнение
Получаем
Находим значение определителя .
Решаем систему уравнений
Применяем подстановку в исходное уравнение:
Заменяем переменную при подстановке в выражение, записанное выше, имеем:
Разделяем переменные:
Переходим теперь к первоначальной функции у и переменной х.
Итого, выражение является общим интегралом исходного дифференциального уравнения.
В случае если в исходном уравнении вида определитель то переменные могут быть разделены подстановкой
Пример. Решить уравнение
Получаем
Находим значение определителя
Применяем подстановку
Подставляем это выражение в исходное уравнение:
Разделяем переменные:
Далее возвращаемся к первоначальной функции у и переменной х.
таким образом, мы получили общий интеграл исходного дифференциального уравнения.
Линейные уравнения.
Определение. Дифференциальное уравнение называется линейным относительно неизвестной функции и ее производной, если оно может быть записано в виде:
при этом, если правая часть Q(x) равна нулю, то такое уравнение называется линейным однородным дифференциальным уравнением, если правая часть Q(x) не равна нулю, то такое уравнение называется линейным неоднородным дифференциальным уравнением.
P(x) и Q(x)- функции непрерывные на некотором промежутке a < x < b.
Линейные однородные дифференциальные уравнения.
Рассмотрим методы нахождения общего решения линейного однородного дифференциального уравнения первого порядка вида
.
Для этого типа дифференциальных уравнений разделение переменных не представляет сложностей.
Общее решение:
Линейные неоднородные дифференциальные уравнения.
Для интегрирования линейных неоднородных уравнений (Q(x)0) применяются в основном два метода: метод Бернулли и метод Лагранжа.
Метод Бернулли.
(Якоб Бернулли (1654-1705) швейцарский математик.)
Суть метода заключается в том, что искомая функция представляется в виде произведения двух функций .
При этом очевидно, что - дифференцирование по частям.
Подставляя в исходное уравнение, получаем:
Далее следует важное замечание т.к. первоначальная функция была представлена нами в виде произведения, то каждый из сомножителей, входящих в это произведение, может быть произвольным, выбранным по нашему усмотрению.
Например, функция может быть представлена как
и т.п.
Таким образом, можно одну из составляющих произведение функций выбрать так, что выражение .
Таким образом, возможно получить функцию u, проинтегрировав, полученное соотношение как однородное дифференциальное уравнение по описанной выше схеме:
Для нахождения второй неизвестной функции v подставим поученное выражение для функции u в исходное уравнение с учетом того, что выражение, стоящее в скобках, равно нулю.
Интегрируя, можем найти функцию v:
; ;
Т.е. была получена вторая составляющая произведения , которое и определяет искомую функцию.
Подставляя полученные значения, получаем:
Окончательно получаем формулу:
, С2 - произвольный коэффициент.
Это соотношение может считаться решением неоднородного линейного дифференциального уравнения в общем виде по способу Бернулли.
Метод Лагранжа.
( Ларганж Жозеф Луи (1736-1813) - французский математик, през. Берлинской АН,
поч. чл. Пет. АН (1776)).
Метод Лагранжа решения неоднородных линейных дифференциальных уравнений еще называют методом вариации произвольной постоянной.
Вернемся к поставленной задаче:
Первый шаг данного метода состоит в отбрасывании правой части уравнения и замене ее нулем.
Далее находится решение получившегося однородного дифференциального уравнения:
.
Для того, чтобы найти соответствующее решение неоднородного дифференциального уравнения, будем считать постоянную С1 некоторой функцией от х.
Тогда по правилам дифференцирования произведения функций получаем:
Подставляем полученное соотношение в исходное уравнение
Из этого уравнения определим переменную функцию С1(х):
Интегрируя, получаем:
Подставляя это значение в исходное уравнение, получаем:
.
Таким образом, мы получили результат, полностью совпадающий с результатом расчета по методу Бернулли.
При выборе метода решения линейных дифференциальных уравнений следует руководствоваться простотой интегрирования функций, входящих в исходный интеграл.
Далее рассмотрим примеры решения различных дифференциальных уравнений различными методами и сравним результаты.
Пример. Решить уравнение
Сначала приведем данное уравнение к стандартному виду:
Применим полученную выше формулу:
Уравнение Бернулли.
Определение. Уравнением Бернулли называется уравнение вида
где P и Q функции от х или постоянные числа, а n постоянное число, не равное 1.
Для решения уравнения Бернулли применяют подстановку , с помощью которой, уравнение Бернулли приводится к линейному.
Для этого разделим исходное уравнение на yn.
Применим подстановку, учтя, что .
Т.е. получилось линейное уравнение относительно неизвестной функции z.
Решение этого уравнения будем искать в виде:
Пример. Решить уравнение
Разделим уравнение на xy2:
Полагаем
.
Полагаем
Произведя обратную подстановку, получаем:
Пример. Решить уравнение
Разделим обе части уравнения на
Полагаем
Получили линейное неоднородное дифференциальное уравнение. Рассмотрим соответствующее ему линейное однородное уравнение:
Полагаем C = C(x) и подставляем полученный результат в линейное неоднородное уравнение, с учетом того, что:
Получаем:
Применяя обратную подстановку, получаем окончательный ответ:
Уравнения в полных дифференциалах (тотальные).
Определение. Дифференциальное уравнение первого порядка вида:
называется уравнением в полных дифференциалах, если левая часть этого уравнения представляет собой полный дифференциал некоторой функции
Интегрирование такого уравнения сводится к нахождению функции u, после чего решение легко находится в виде:
Таким образом, для решения надо определить:
1) в каком случае левая часть уравнения представляет собой полный дифференциал функции u;
2) как найти эту функцию.
Если дифференциальная форма является полным дифференциалом некоторой функции u, то можно записать:
Т.е. .
Найдем смешанные производные второго порядка, продифференцировав первое уравнение по у, а второе по х:
Приравнивая левые части уравнений, получаем необходимое и достаточное условие того, что левая часть дифференциального уравнения является полным дифференциалом. Это условие также называется условием тотальности.
Теперь рассмотрим вопрос о нахождении собственно функции u.
Проинтегрируем равенство :
Вследствие интегрирования получаем не постоянную величину С, а некоторую функцию С(у), т.к. при интегрировании переменная у полагается постоянным параметром.
Определим функцию С(у).
Продифференцируем полученное равенство по у.
Откуда получаем:
Для нахождения функции С(у) необходимо проинтегрировать приведенное выше равенство. Однако, перед интегрированием надо доказать, что функция С(у) не зависит от х. Это условие будет выполнено, если производная этой функции по х равна нулю.
Теперь определяем функцию С(у):
Подставляя этот результат в выражение для функции u, получаем:
Тогда общий интеграл исходного дифференциального уравнения будет иметь вид:
Следует отметить, что при решении уравнений в полных дифференциалах не обязательно использовать полученную формулу. Решение может получиться более компактным, если просто следовать методу, которым формула была получена.
Пример. Решить уравнение
Проверим условие тотальности:
Условие тотальности выполняется, следовательно, исходное дифференциальное уравнение является уравнением в полных дифференциалах.
Определим функцию u.
;
Итого,
Находим общий интеграл исходного дифференциального уравнения:
Уравнения вида y = f(y) и x = f(y).
Решение уравнений, не содержащих в одном случае аргумента х, а в другом функции у, ищем в параметрической форме, принимая за параметр производную неизвестной функции.
Для уравнения первого типа получаем:
Делая замену, получаем:
В результате этих преобразований имеем дифференциальное уравнение с разделяющимися переменными.
Общий интеграл в параметрической форме представляется системой уравнений:
Исключив из этой системы параметр р, получим общий интеграл и не в параметрической форме.
Для дифференциального уравнения вида x = f(y) с помощью той же самой подстановки и аналогичных рассуждений получаем результат:
Уравнения Лагранжа и Клеро.
( Алекси Клод Клеро (1713 1765) французский математик
ин. поч. член Петерб. АН )
Определение. Уравнением Лагранжа называется дифференциальное уравнение, линейное относительно х и у, коэффициенты которого являются функциями от y.
Для нахождения общего решение применяется подстановка p = y.
Дифференцируя это уравнение,c учетом того, что , получаем:
Если решение этого (линейного относительно х) уравнения есть то общее решение уравнения Лагранжа может быть записано в виде:
Определение. Уравнением Клеро называется уравнение первой степени (т.е. линейное) относительно функции и аргумента вида:
Вообще говоря, уравнение Клеро является частным случаем уравнения Лагранжа.
С учетом замены , уравнение принимает вид:
Это уравнение имеет два возможных решения:
или
В первом случае:
Видно, что общий интеграл уравнения Клеро представляет собой семейство прямых линий.
Во втором случае решение в параметрической форме выражается системой уравнений:
Исключая параметр р, получаем второе решение F(x, y) = 0. Это решение не содержит произвольной постоянной и не получено из общего решения, следовательно, не является частным решением.
Это решение будет являться особым интегралом. ( См. Особое решение. )
Далее рассмотрим примеры решения различных типов дифференциальных уравнений первого порядка.
Пример. Решить уравнение с заданными начальными условиями.
Это линейное неоднородное дифференциальное уравнение первого порядка.
Решим соответствующее ему однородное уравнение.
Для неоднородного уравнения общее решение имеет вид:
Дифференцируя, получаем:
Для нахождения функции С(х) подставляем полученное значение в исходное дифференциальное уравнение:
Итого, общее решение:
C учетом начального условия определяем постоянный коэффициент C.
Окончательно получаем:
Для проверки подставим полученный результат в исходное дифференциальное уравнение: верно
Ниже показан график интегральной кривой уравнения.
Пример. Найти общий интеграл уравнения .
Это уравнение с разделяющимися переменными.
Общий интеграл имеет вид:
Построим интегральные кривые дифференциального уравнения при различных значениях С.
С = - 0,5 С = -0,02 С = -1 С = -2
С = 0,02 С = 0,5 С = 1 С = 2
Пример. Найти решение дифференциального уравнения, удовлетворяющее заданным начальным условиям.
Это уравнение с разделяющимися переменными.
Общее решение имеет вид:
Найдем частное решение при заданном начальном условии у(0) = 0.
Окончательно получаем:
Пример. Решить предыдущий пример другим способом.
Действительно, уравнение может быть рассмотрено как линейное неоднородное дифференциальное уравнение.
Решим соответствующее ему линейное однородное уравнение.
Решение неоднородного уравнения будет иметь вид:
Тогда
Подставляя в исходное уравнение, получаем:
Итого
С учетом начального условия у(0) = 0 получаем
Как видно результаты, полученные при решении данного дифференциального уравнения различными способами, совпадают.
При решении дифференциальных уравнений бывает возможно выбирать метод решения, исходя из сложности преобразований.
Пример. Решить уравнение с начальным условием у(0) = 0.
Это линейное неоднородное уравнение. Решим соответствующее ему однородное уравнение.
Для линейного неоднородного уравнения общее решение будет иметь вид:
Для определения функции С(х) найдем производную функции у и подставим ее в исходное дифференциальное уравнение.
Итого
Проверим полученное общее решение подстановкой в исходное дифференциальное уравнение.
(верно)
Найдем частное решение при у(0) = 0.
Окончательно
Пример. Найти решение дифференциального уравнения
с начальным условием у(1) = 1.
Это уравнение может быть преобразовано и представлено как уравнение с разделенными переменными.
С учетом начального условия:
Окончательно
Пример. Решить дифференциальное уравнение с начальным условием у(1) = 0.
Это линейное неоднородное уравнение.
Решим соответствующее ему однородное уравнение.
Решение неоднородного уравнения будет иметь вид:
Подставим в исходное уравнение:
Общее решение будет иметь вид:
C учетом начального условия у(1) = 0:
Частное решение:
Пример. Найти решение дифференциального уравнения с начальным условием у(1) = е.
Это уравнение может быть приведено к виду уравнения с разделяющимися переменными с помощью замены переменных.
Обозначим:
Уравнение принимает вид:
Получили уравнение с разделяющимися переменными.
Сделаем обратную замену:
Общее решение:
C учетом начального условия у(1) = е:
Частное решение:
Второй способ решения.
Получили линейное неоднородное дифференциальное уравнение. Соответствующее однородное:
Решение исходного уравнения ищем в виде:
Тогда
Подставим полученные результаты в исходное уравнение:
Получаем общее решение:
Пример. Решить дифференциальное уравнение с начальным условием у(1)=0.
В этом уравнении также удобно применить замену переменных.
Уравнение принимает вид:
Делаем обратную подстановку:
Общее решение:
C учетом начального условия у(1) = 0:
Частное решение:
Второй способ решения.
Замена переменной:
Общее решение:
уравнений первого порядка.
у a
b
A S
x
Как уже говорилось выше (см. Интегральные кривые. ), линия S, которая задается функцией, являющейся каким- либо решением дифференциального уравнения, называется интегральной кривой уравнения
Производная y является угловым коэффициентом касательной к интегральной кривой.
В любой точке А(х, у) интегральной кривой этот угловой коэффициент касательной может быть найден еще до решения дифференциального уравнения.
Т.к. касательная указывает направление интегральной кривой еще до ее непосредственного построения, то при условии непрерывности функции f(x, y) и непрерывного перемещения точки А можно наглядно изобразить поле направлений кривых, которые получаются в результате интегрирования дифференциального уравнения, т.е. представляют собой его общее решение.
Определение. Множество касательных в каждой точке рассматриваемой области называется полем направлений.
С учетом сказанного выше можно привести следующее геометрическое истолкование дифференциального уравнения:
1) Задать дифференциальное уравнение первого порядка это значит задать поле направлений.
2) Решить или проинтегрировать дифференциальное уравнение это значит найти всевозможные кривые, у которых направление касательных в каждой точке совпадает с полем направлений.
Определение. Линии равного наклона в поле направлений называются изоклинами.
Численные методы решения дифференциальных уравнений.
Известные методы точного интегрирования дифференциальных уравнений позволяют найти решение в виде аналитической функции, однако эти методы применимы для очень ограниченного класса уравнений. Большинство уравнений, встречающихся при решении практических задач нельзя проинтегрировать с помощью этих методов.
В таких случаях используются численные методы решения, которые представляют решение дифференциального уравнения не в виде аналитической функции, а в виде таблиц значений искомой функции в зависимости от значения переменной.
Существует несколько методов численного интегрирования дифференциальных уравнений, которые отличаются друг от друга по сложности вычислений и точности результата.
Рассмотрим некоторые из них.
Метод Эйлера.
(Леонард Эйлер (1707 1783) швейцарский математик )
Известно, что уравнение задает в некоторой области поле направлений. Решение этого уравнения с некоторыми начальными условиями дает кривую, которая касается поля направлений в любой точке.
Если взять последовательность точек х0, х1, х2, …. и заменить на получившихся отрезках интегральную кривую на отрезки касательных к ней, то получим ломаную линию.
y
M2
M1 M3
M0
y0 M4
0 x0 x1 x2 x3 x4 x
При подстановке заданных начальных условий (х0, у0) в дифференциальное уравнение получаем угловой коэффициент касательной к интегральной кривой в начальной точке
Заменив на отрезке [x0, x1] интегральную кривую на касательную к ней, получаем значение
Производя аналогичную операцию для отрезка [x1, x2], получаем:
Продолжая подобные действия далее, получаем ломаную кривую, которая называется ломаной Эйлера.
Можно записать общую формулу вычислений:
Если последовательность точек хi выбрать так, чтобы они отстояли друг от друга на одинаковое расстояние h, называемое шагом вычисления, то получаем формулу:
Следует отметить, что точность метода Эйлера относительно невысока. Увеличить точность можно, конечно, уменьшив шаг вычислений, однако, это приведет к усложнению расчетов. Поэтому на практике применяется так называемый уточненный метод Эйлера или формула пересчета.
Суть метода состоит в том, что в формуле вместо значения
берется среднее арифметическое значений f(x0, y0) и f(x1, y1). Тогда уточненное значение:
Затем находится значение производной в точке . Заменяя f(x0, y0) средним арифметическим значений f(x0, y0) и , находят второе уточненное значение у1.
Затем третье:
и т.д. пока два последовательных уточненных значения не совпадут в пределах заданной степени точности. Тогда это значение принимается за ординату точки М1 ломаной Эйлера.
Аналогичная операция производится для остальных значений у.
Подобное уточнение позволяет существенно повысить точность результата.
При использовании компьютерной версии “Курса высшей математики” возможно запустить программу, которая решает любое дифференциальное уравнение первого порядка методом Эйлера и уточненным методом Эйлера. На каждом шаге вычислений подробно выводятся все указанные выше значения.
Для запуска программы дважды щелкните на значке
Примечание: Для запуска программы необходимо чтобы на компьютере была установлена программа Maple ( Waterloo Maple Inc.) любой версии, начиная с MapleV Release 4.
Метод Рунге Кутта.
Метод Рунге Кутта является более точным по сравнению с методом Эйлера.
Суть уточнения состоит в том, что искомое решение представляется в виде разложения в ряд Тейлора. (См. Формула Тейлора. )
Если в этой формуле ограничиться двумя первыми слагаемыми, то получим формулу метода Эйлера. Метод Рунге Кутта учитывает четыре первых члена разложения.
.
В методе Рунге Кутта приращения yi предлагается вычислять по формуле:
где коэффициенты ki вычисляются по формулам:
Пример. Решить методом Рунге Кутта дифференциальное уравнение при начальном условии у(0) = 1 на отрезке [0; 0,5] с шагом 0,1.
Для i = 0 вычислим коэффициенты ki.
Последующие вычисления приводить не будем, а результаты представим в виде таблицы.
i |
xi |
k |
yi |
yi |
|
0 |
0 |
1 |
0,1000 |
0,1104 |
1 |
2 |
0,1100 |
||||
3 |
0,1105 |
||||
4 |
0,1155 |
||||
1 |
0,1 |
1 |
0,1210 |
0,1325 |
1,1104 |
2 |
0,1321 |
||||
3 |
0,1326 |
||||
4 |
0,1443 |
||||
2 |
0,2 |
1 |
0,1443 |
0,1569 |
1,2429 |
2 |
0,1565 |
||||
3 |
0,1571 |
||||
4 |
0,1700 |
||||
3 |
0.3 |
1 |
0,1700 |
0,1840 |
1,3998 |
2 |
0,1835 |
||||
3 |
0,1842 |
||||
4 |
0,1984 |
||||
4 |
0,4 |
1 |
0,1984 |
0,2138 |
1,5838 |
2 |
0,2133 |
||||
3 |
0,2140 |
||||
4 |
0,2298 |
||||
5 |
0,5 |
1,7976 |
Решим этот же пример методом Эйлера.
Применяем формулу
Производя аналогичные вычисления далее, получаем таблицу значений:
i |
0 |
1 |
2 |
3 |
4 |
5 |
xi |
0,0 |
0,1 |
0,2 |
0,3 |
0,4 |
0,5 |
yi |
1 |
1,1 |
1,22 |
1,362 |
1,528 |
1,721 |
Применим теперь уточненный метод Эйлера.
i |
0 |
1 |
2 |
3 |
4 |
5 |
xi |
0,0 |
0,1 |
0,2 |
0,3 |
0,4 |
0,5 |
yi |
1 |
1,1 |
1,243 |
1,400 |
1,585 |
1,799 |
Для сравнения точности приведенных методов численного решение данного уравнения решим его аналитически и найдем точные значения функции у на заданном отрезке.
Уравнение является линейным неоднородным дифференциальным уравнением первого порядка. Решим соответствующее ему однородное уравнение.
Решение неоднородного уравнения имеет вид
Общее решение:
C учетом начального условия:
Частное решение:
Для сравнения полученных результатов составим таблицу.
i |
xi |
yi |
|||
Метод Эйлера |
Уточненный метод Эйлера |
Метод Рунге - Кутта |
Точное значение |
||
0 |
0 |
1 |
1 |
1 |
1 |
1 |
0,1 |
1,1 |
1,1 |
1,1104 |
1,1103 |
2 |
0,2 |
1,22 |
1,243 |
1,2429 |
1,2428 |
3 |
0,3 |
1,362 |
1,4 |
1,3998 |
1,3997 |
4 |
0,4 |
1,528 |
1,585 |
1,5838 |
1,5837 |
5 |
0,5 |
1,721 |
1,799 |
1,7976 |
1,7975 |
Как видно из полученных результатов метод Рунге Кутта дает наиболее точный ответ. Точность достигает 0,0001. Кроме того, следует обратить внимание на то, ошибка (расхождение между точным и приближенным значениями) увеличивается с каждым шагом вычислений. Это обусловлено тем, что во первых полученное приближенное значение округляется на каждом шаге, а во вторых тем, что в качестве основы вычисления принимается значение, полученное на предыдущем шаге, т.е. приближенное значение. Таким образом происходит накопление ошибки.
Это хорошо видно из таблицы. С каждым новым шагом приближенное значение все более отличается от точного.
При использовании кмпьютерной версии “Курса высшей математики” возможно запустить программу, которая решает любое дифференциальное уравнение первого порядка рассмотренным выше методом Рунге- Кутта. Программа подробно выводит результаты вычислений на каждом шаге.
Для запуска программы дважды щелкните на значке
Примечание: Для запуска программы необходимо чтобы на компьютере была установлена программа Maple ( Waterloo Maple Inc.) любой версии, начиная с MapleV Release 4.
Дифференциальные уравнения высших порядков.
Определение. Дифференциальным уравнением порядка n называется уравнение вида:
В некоторых случаях это уравнение можно разрешить относительно y(n):
Так же как и уравнение первого порядка, уравнения высших порядков имеют бесконечное количество решений.
Определение. Решение удовлетворяет начальным условиям , если
Определение. Нахождение решения уравнения , удовлетворяющего начальным условиям , называется решением задачи Коши.
Теорема Коши. (Теорема о необходимых и достаточных условиях существования решения задачи Коши).
Если функция (n-1) й переменных вида в некоторой области D (n-1)- мерного пространства непрерывна и имеет непрерывные частные производные по , то какова бы не была точка () в этой области, существует единственное решение уравнения , определенного в некотором интервале, содержащем точку х0, удовлетворяющее начальным условиям .
Дифференциальные уравнения высших порядков, решение которых может быть найдено аналитически, можно разделить на несколько основных типов.
Рассмотрим подробнее методы нахождения решений этих уравнений.
Уравнения, допускающие понижение порядка.
Понижение порядка дифференциального уравнения основной метод решения уравнений высших порядков. Этот метод дает возможность сравнительно легко находить решение, однако, он применим далеко не ко всем уравнениям. Рассмотрим случаи, когда возможно понижение порядка.
Уравнения вида y(n) = f(x).
Если f(x) функция непрерывная на некотором промежутке a < x < b, то решение может быть найдено последовательным интегрированием.
…………………………………………………………….
Пример. Решить уравнение с начальными условиями x0 = 0; y0 = 1;
Подставим начальные условия:
Получаем частное решение (решение задачи Коши): .
Ниже показана интегральная кривая данного дифференциального уравнения.
и ее производных до порядка k 1 включительно.
Это уравнения вида:
В уравнениях такого типа возможно понижение порядка на k единиц. Для этого производят замену переменной:
Тогда получаем:
Теперь допустим, что полученное дифференциальное уравнение проинтегрировано и совокупность его решений выражается соотношением:
Делая обратную подстановку, имеем:
Интегрируя полученное соотношение последовательно k раз, получаем окончательный ответ:
Пример. Найти общее решение уравнения .
Применяем подстановку
Произведя обратную замену, получаем:
Общее решение исходного дифференциального уравнения:
Отметим, что это соотношение является решением для всех значений переменной х кроме значения х =0.
Уравнения, не содержащие явно независимой переменной.
Это уравнения вида
Порядок таких уравнений может быть понижен на единицу с помощью замены переменных
и т.д.
Подставляя эти значения в исходное дифференциальное уравнение, получаем:
Если это уравнение проинтегрировать, и - совокупность его решений, то для решения данного дифференциального уравнения остается решить уравнение первого порядка:
Пример. Найти общее решение уравнения
Замена переменной:
1)
Для решения полученного дифференциального уравнения произведем замену переменной:
С учетом того, что , получаем:
Общий интеграл имеет вид:
2)
Таким образом, получили два общих решения.
Линейные дифференциальные уравнения высших порядков.
Определение. Линейным дифференциальным уравнением n го порядка называется любое уравнение первой степени относительно функции у и ее производных вида:
где p0, p1, …,pn функции от х или постоянные величины, причем p0 0.
Левую часть этого уравнения обозначим L(y).
Определение. Если f(x) = 0, то уравнение L(y) = 0 называется линейным однородным уравнением, если f(x) 0, то уравнение L(y) = f(x) называется линейным неоднородным уравнением, если все коэффициенты p0, p1, p2, … pn постоянные числа, то уравнение L(y) = f(x) называется линейным дифференциальным уравнением высшего порядка с постоянными коэффициентами.
Отметим одно важное свойство линейных уравнений высших порядков, которое отличает их от нелинейных. Для нелинейных уравнений частный интеграл находится из общего, а для линейных наоборот, общий интеграл составляется из частных. Линейные уравнения представляют собой наиболее изученный класс дифференциальных уравнений высших порядков. Это объясняется сравнительной простотой нахождения решения. Если при решении каких либо практических задач требуется решить нелинейное дифференциальное уравнение, то часто применяются приближенные методы, позволяющие заменить такое уравнение “близким” к нему линейным.
Рассмотрим способы интегрирования некоторых типов линейных дифференциальных уравнений высших порядков.
произвольными коэффициентами.
Рассмотрим уравнение вида
Определение. Выражение называется линейным дифференциальным оператором.
Линейный дифференциальный оператор обладает следующими свойствами:
1)
2)
Решения линейного однородного уравнения обладают следующими свойствами:
1) Если функция у1 является решением уравнения, то функция Су1, где С постоянное число, также является его решением.
2) Если функции у1 и у2 являются решениями уравнения, то у1 +у2 также является его решением.
Структура общего решения.
Определение. Фундаментальной системой решений линейного однородного дифференциального уравнения n го порядка на интервале (a, b) называется всякая система n линейно независимых на этом интервале решений уравнения.
Определение. Если из функций yi составить определитель n го порядка
,
то этот определитель называется определителем Вронского.
( Юзеф Вроньский (1776 1853) польский математик и философ - мистик)
Теорема. Если функции линейно зависимы, то составленный для них определитель Вронского равен нулю.
Теорема. Если функции линейно независимы, то составленный для них определитель Вронского не равен нулю ни в одной точке рассматриваемого интервала.
Теорема. Для того, чтобы система решений линейного однородного дифференциального уравнения была фундаментальной необходимо и достаточно, чтобы составленный для них определитель Вронского был не равен нулю.
Теорема. Если - фундаментальная система решений на интервале (a, b), то общее решение линейного однородного дифференциального уравнения является линейной комбинацией этих решений.
,
где Ci постоянные коэффициенты.
Применение приведенных выше свойств и теорем рассмотрим на примере линейных однородных дифференциальных уравнений второго порядка.
уравнения второго порядка.
Из вышеизложенного видно, что отыскание общего решения линейного однородного дифференциального уравнения сводится к нахождению его фундаментальной системы решений.
Однако, даже для уравнения второго порядка, если коэффициенты р зависят от х, эта задача не может быть решена в общем виде.
Тем не менее, если известно одно ненулевое частное решение, то задача может быть решена.
Теорема. Если задано уравнение вида и известно одно ненулевое решение у = у1, то общее решение может быть найдено по формуле:
Таким образом, для получения общего решения надо подобрать какое либо частное решение дифференциального уравнения, хотя это бывает часто довольно сложно.
постоянными коэффициентами.
Решение дифференциального уравнения вида или, короче, будем искать в виде , где k = const.
Т.к. то
При этом многочлен называется характеристическим многочленом дифференциального уравнения.
Для того, чтобы функция являлась решением исходного дифференциального уравнения, необходимо и достаточно, чтобы
т.е.
Т.к. ekx 0, то - это уравнение называется характеристическим уравнением.
Как и любое алгебраическое уравнение степени n, характеристическое уравнение имеет n корней. Каждому корню характеристического уравнения ki соответствует решение дифференциального уравнения.
В зависимости от коэффициентов k характеристическое уравнение может иметь либо n различных действительных корней, либо среди действительных корней могут быть кратные корни, могут быть комплексно сопряженные корни, как различные, так и кратные.
Не будем подробно рассматривать каждый случай, а сформулируем общее правило нахождения решения линейного однородного дифференциального уравнения с постоянными коэффициентами.
1) Составляем характеристическое уравнение и находим его корни.
2) Находим частные решения дифференциального уравнения, причем:
a) каждому действительному корню соответствует решение ekx;
б) каждому действительному корню кратности m ставится в соответствие m решений:
в) каждой паре комплексно сопряженных корней характеристического уравнение ставится в соответствие два решения:
и .
г) каждой паре m кратных комплексно сопряженных корней характеристического уравнения ставится в соответствие 2m решений:
3) Составляем линейную комбинацию найденных решений.
Эта линейная комбинация и будет являться общим решением исходного линейного однородного дифференциального уравнения с постоянными коэффициентами.
Пример. Решить уравнение .
Составим характеристическое уравнение:
Общее решение имеет вид:
Пример. Решить уравнение
Это линейное однородное дифференциальное уравнение с переменными коэффициентами второго порядка. Для нахождения общего решения необходимо отыскать какое - либо частное решение.
Таким частным решением будет являться функция
Исходное дифференциальное уравнение можно преобразовать:
Общее решение имеет вид:
Окончательно:
Пример. Решить уравнение
Составим характеристическое уравнение:
Общее решение:
Пример. Решить уравнение
Характеристическое уравнение:
Общее решение:
Пример. Решить уравнение
Характеристическое уравнение:
Общее решение:
Пример. Решить уравнение
Характеристическое уравнение:
Общее решение:
Пример. Решить уравнение
Характеристическое уравнение:
Общее решение:
Пример. Решить уравнение
Характеристическое уравнение:
Общее решение:
Пример. Решить уравнение
Это уравнение не является линейным, следовательно, приведенный выше метод решения к нему неприменим.
Понизим порядок уравнения с помощью подстановки
Тогда
Окончательно получаем:
Это выражение будет общим решением исходного дифференциального уравнения. Полученное выше решение у1 = С1 получается из общего решения при С = 0.
Пример. Решить уравнение
Производим замену переменной:
Общее решение:
с произвольными коэффициентами.
Рассмотрим уравнение вида
С учетом обозначения можно записать:
При этом будем полагать, что коэффициенты и правая часть этого уравнения непрерывны на некотором интервале ( конечном или бесконечном).
Теорема. Общее решение линейного неоднородного дифференциального уравнения в некоторой области есть сумма любого его решения и общего решения соответствующего линейного однородного дифференциального уравнения.
Доказательство. Пусть Y некоторое решение неоднородного уравнения.
Тогда при подстановке этого решения в исходное уравнение получаем тождество:
Пусть - фундаментальная система решений линейного однородного уравнения . Тогда общее решение однородного уравнения можно записать в виде:
Далее покажем, что сумма является общим решением неоднородного уравнения.
Вообще говоря, решение Y может быть получено из общего решения, т.к. является частным решением.
Таким образом, в соответствии с доказанной теоремой, для решения линейного неоднородного дифференциального уравнения необходимо найти общее решение соответствующего однородного уравнения и каким- то образом отыскать одно частное решение неоднородного уравнения. Обычно оно находится подбором.
На практике удобно применять метод вариации произвольных постоянных.
Для этого сначала находят общее решение соответствующего однородного уравнения в виде:
Затем, полагая коэффициенты Ci функциями от х, ищется решение неоднородного уравнения:
Можно доказать, что для нахождения функций Ci(x) надо решить систему уравнений:
Пример. Решить уравнение
Решаем линейное однородное уравнение
Решение неоднородного уравнения будет иметь вид:
Составляем систему уравнений:
Решим эту систему:
Из соотношения найдем функцию А(х).
Теперь находим В(х).
Подставляем полученные значения в формулу общего решения неоднородного уравнения:
Окончательный ответ:
Таким образом, удалось избежать нахождения частного решения неоднородного уравнения методом подбора.
Вообще говоря, метод вариации произвольных постоянных пригоден для нахождения решений любого линейного неоднородного уравнения. Но т.к. нахождение фундаментальной системы решений соответствующего однородного уравнения может быть достаточно сложной задачей, этот метод в основном применяется для неоднородных уравнений с постоянными коэффициентами.
коэффициентами.
Уравнения с правой частью специального вида.
Представляется возможным представить вид частного решения в зависимости от вида правой части неоднородного уравнения.
Различают следующие случаи:
I. Правая часть линейного неоднородного дифференциального уравнения имеет вид:
где - многочлен степени m.
Тогда частное решение ищется в виде:
Здесь Q(x)- многочлен той же степени, что и P(x), но с неопределенными коэффициентами, а r число, показывающее сколько раз число является корнем характеристического уравнения для соответствующего линейного однородного дифференциального уравнения.
Пример. Решить уравнение .
Решим соответствующее однородное уравнение:
Теперь найдем частное решение исходного неоднородного уравнения.
Сопоставим правую часть уравнения с видом правой части, рассмотренным выше.
Частное решение ищем в виде: , где
Т.е.
Теперь определим неизвестные коэффициенты А и В.
Подставим частное решение в общем виде в исходное неоднородное дифференциальное уравнение.
Итого, частное решение:
Тогда общее решение линейного неоднородного дифференциального уравнения:
II. Правая часть линейного неоднородного дифференциального уравнения имеет вид:
Здесь Р1(х) и Р2(х) многочлены степени m1 и m2 соответственно.
Тогда частное решение неоднородного уравнения будет иметь вид:
где число r показывает сколько раз число является корнем характеристического уравнения для соответствующего однородного уравнения, а Q1(x) и Q2(x) многочлены степени не выше m, где m- большая из степеней m1 и m2.
Заметим, что если правая часть уравнения является комбинацией выражений рассмотренного выше вида, то решение находится как комбинация решений вспомогательных уравнений, каждое из которых имеет правую часть, соответствующую выражению, входящему в комбинацию.
Т.е. если уравнение имеет вид: , то частное решение этого уравнения будет где у1 и у2 частные решения вспомогательных уравнений
и
Для иллюстрации решим рассмотренный выше пример другим способом.
Пример. Решить уравнение
Правую часть дифференциального уравнения представим в виде суммы двух функций f1(x) + f2(x) = x + (-sinx).
Составим и решим характеристическое уравнение:
Получаем: Т.е.
Итого:
Анализируя функцию f2(x), получаем:
Таким образом,
Итого:
Т.е. искомое частное решение имеет вид:
Общее решение неоднородного дифференциального уравнения:
Рассмотрим примеры применения описанных методов.
Пример. Решить уравнение
Составим характеристическое уравнение для соответствующего линейного однородного дифференциального уравнения:
Общее решение однородного уравнения:
Теперь найдем частное решение неоднородного уравнения в виде:
Воспользуемся методом неопределенных коэффициентов.
Подставляя в исходное уравнение, получаем:
Частное решение имеет вид:
Общее решение линейного неоднородного уравнения:
Пример. Решить уравнение
Характеристическое уравнение:
Общее решение однородного уравнения:
Частное решение неоднородного уравнения: .
Находим производные и подставляем их в исходное неоднородное уравнение:
Получаем общее решение неоднородного дифференциального уравнения:
Нормальные системы обыкновенных дифференциальных уравнений.
Определение. Совокупность соотношений вида:
где х- независимая переменная, у1, у2,…,уn искомые функции, называется системой дифференциальных уравнений первого порядка.
Определение. Система дифференциальных уравнений первого порядка, разрешенных относительно производных от неизвестных функций называется нормальной системой дифференциальных уравнений.
Такая система имеет вид:
(1)
Для примера можно сказать, что график решения системы двух дифференциальных уравнений представляет собой интегральную кривую в трехмерном пространстве.
Теорема. (Теорема Коши). Если в некоторой области (n-1) мерного пространства функции … непрерывны и имеют непрерывные частные производные по , то для любой точки этой области существует единственное решение
системы дифференциальных уравнений вида (1), определенное в некоторой окрестности точки х0 и удовлетворяющее начальным условиям
Определение. Общим решением системы дифференциальных уравнений вида (1) будет совокупность функций , , … , которые при подстановке в систему (1) обращают ее в тождество.
уравнений с постоянными коэффициентами.
При рассмотрении систем дифференциальных уравнений ограничимся случаем системы трех уравнений (n = 3). Все нижесказанное справедливо для систем произвольного порядка.
Определение. Нормальная система дифференциальных уравнений c постоянными коэффициентами называется линейной однородной, если ее можно записать в виде:
(2)
Решения системы (2) обладают следующими свойствами:
1) Если y, z, u решения системы, то Cy, Cz, Cu , где C = const тоже являются решениями этой системы.
2) Если y1, z1, u1 и y2, z2, u2 решения системы, то y1 + y2, z1 + z2, u1 + u2 тоже являются решениями системы.
Решения системы ищутся в виде:
Подставляя эти значения в систему (2) и перенеся все члены в одну сторону и сократив на ekx, получаем:
Для того, чтобы полученная система имела ненулевое решение необходимо и достаточно, чтобы определитель системы был равен нулю, т.е.:
В результате вычисления определителя получаем уравнение третьей степени относительно k. Это уравнение называется характеристическим уравнением и имеет три корня k1, k2, k3. Каждому из этих корней соответствует ненулевое решение системы (2):
Линейная комбинация этих решений с произвольными коэффициентами будет решением системы (2):
Пример. Найти общее решение системы уравнений:
Составим характеристическое уравнение:
Решим систему уравнений:
Для k1:
Полагая (принимается любое значение), получаем:
Для k2:
Полагая (принимается любое значение), получаем:
Общее решение системы:
Этот пример может быть решен другим способом:
Продифференцируем первое уравнение:
Подставим в это выражение производную у =2x + 2y из второго уравнения.
Подставим сюда у, выраженное из первого уравнения:
Обозначив , получаем решение системы:
Пример. Найти решение системы уравнений
Эта система дифференциальных уравнений не относится к рассмотренному выше типу, т.к. не является однородным (в уравнение входит независимая переменная х).
Для решения продифференцируем первое уравнение по х. Получаем:
Заменяя значение z из второго уравнения получаем: .
С учетом первого уравнения, получаем:
Решаем полученное дифференциальное уравнение второго порядка.
Общее решение однородного уравнения:
Теперь находим частное решение неоднородного дифференциального уравнения по формуле
Общее решение неоднородного уравнения:
Подставив полученное значение в первое уравнение системы, получаем:
Пример. Найти решение системы уравнений:
Составим характеристическое уравнение:
Если принять = 1, то решения в этом случае получаем:
Если принять = 1, то получаем:
Если принять = 3, то получаем:
Общее решение имеет вид:
Элементы теории устойчивости.
Теория устойчивости решений дифференциальных уравнений является одним из разделов качественной теории дифференциальных уравнений, которая посвящена не нахождению какого либо решения уравнения, а изучению характера поведения этого решения при изменении начальных условий или аргумента.
Этот метод особенно важен, т.к. позволяет делать вывод о характере решения без непосредственного нахождения этого решения. Т.е. даже в тех случаях, когда решение дифференциального уравнения вообще не может быть найдено аналитически.
Пусть имеется некоторое явление, описанное системой дифференциальных уравнений:
(1)
и начальные условия:
Для конкретного явления начальные условия определяются опытным путем и поэтому неточны.
Теорема. (о непрерывной зависимости решения от начальных условий)
Если правая часть дифференциального уравнения непрерывна и по переменной у имеет ограниченную частную производную на области прямоугольника, ограниченного , то решение
, удовлетворяющее начальным условиям , непрерывно зависит от начальных данных, т.е. для любого , при котором если
то при условии, что
где
Эта теорема справедлива как для одного дифференциального уравнения, так и для системы уравнений.
Определение. Если - решение системы дифференциальных уравнений, то это решение называется устойчивым по Ляпунову, если для любого , такое, что для любого решения той же системы, начальные условия которого удовлетворяют неравенствам
справедливы неравенства
(Ляпунов Александр Михайлович (1857 1918) академик Петерб. АН)
Т.е. можно сказать, что решение (t) устойчиво по Ляпунову, если близкие к нему по начальным условиям решения остаются близкими и при t t0.
Если , то решение (t) называется асимптотически устойчивым.
Исследование на устойчивость по Ляпунову произвольного решения системы можно свести к исследованию на устойчивость равного нулю решения некоторой другой системы, которая получена из данной заменой неизвестных функций:
Тогда:
(2)
Система (2) имеет тривиальное (равное нулю) решение
Теорема. Решение системы (1) устойчиво по Ляпунову тогда и только тогда, когда устойчиво по Ляпунову тривиальное решение системы (2).
Это тривиальное решение называется положением равновесия или точкой покоя.
Определение. Точка покоя системы (2) устойчива по Ляпунову, если для любого такое, что из неравенства
следует
.
Теорема. (Теорема Ляпунова). Пусть задана система
имеющая тривиальное решение .
Пусть существует дифференцируемая функция , удовлетворяющая условиям:
1) 0 и v = 0 только при у1 = у2 = … = уn =0, т.е. функция v имеет минимум в начале координат.
2) Полная производная функции v вдоль фазовой траектории (т.е. вдоль решения yi(t) системы (1)) удовлетворяет условию:
при
Тогда точка покоя устойчива по Ляпунову.
Если ввести дополнительное требование, чтобы вне сколь угодно малой окрестности начала координат выполнялось условие
где - постоянная величина, то точка покоя асимптотически устойчива.
Функция v называется функцией Ляпунова.
Классификация точек покоя.
Рассмотрим систему двух линейных дифференциальных уравнений с постоянными коэффициентами
Характеристическое уравнение этой системы имеет вид:
Рассмотрим следующие возможные случаи:
1) Корни характеристического уравнения действительные, отрицательные и различные.
Точка покоя будет устойчива. Такая точка покоя называется устойчивым узлом.
2) Корни характеристического уравнения действительны и
или .
В этом случае точка покоя также будет устойчива.
3) Хотя бы один из корней положителен.
В этом случае точка покоя неустойчива, и такую точку называют неустойчивым седлом.
4) Оба корня характеристического уравнения положительны .
В этом случае точка покоя неустойчива, и такую точку называют неустойчивым узлом.
Если полученного решения системы исключить параметр t, то полученная функция дает траекторию движения в системе координат XOY.
Возможны следующие случаи:
Устойчивый узел. Неустойчивый узел. Седло.
5) Корни характеристического уравнения комплексные .
Если р = 0, т.е. корни чисто мнимые, то точка покоя (0, 0) устойчива по Ляпунову.
Такая точка покоя называется центром.
Если p< 0, то точка покоя устойчива и называется устойчивым фокусом.
Если p > 0, то точка покоя неустойчива и называется неустойчивым фокусом.
Уравнения математической физики.
Уравнения в частных производных.
Определение. Дифференциальным уравнением в частных производных называется уравнение относительно неизвестной функции нескольких переменных, ее аргументов и ее частных производных различных порядков.
Порядком дифференциального уравнения в частных производных называется порядок старшей производной, входящей в это уравнение. Решением уравнения будет некоторая функция , которая обращает уравнение в тождество.
производных первого порядка.
Дифференциальное уравнение в частных производных первого порядка от функции можно в общем виде записать как
Линейное уравнение в частных производных имеет вид:
, (1)
где Xi некоторые заданные функции.
Очевидно, что одним из решений такого уравнения будет функция u = C.
Рассмотрим систему уравнений:
(2)
или - такая система называется нормальной.
Общее решение этой системы имеет вид:
Если разрешить эти уравнения относительно постоянных С, получим:
Каждая из функций является интегралом системы (2).
Теорема. Если - интеграл системы (2), то функция - решение уравнения (1).
физики.
1) Волновое уравнение. (Уравнение колебаний струны, электроколебания, крутильные колебания вала и др.) Это простейшее уравнение гиперболического типа.
2) Уравнение теплопроводности. (Уравнение Фурье) Это простейшее уравнение параболического типа. Описывает процессы теплопроводности, фильтрации жидкости и газа, некоторые вопросы теории вероятностей.
3) Уравнение Лапласа. Это простейшее уравнение эллиптического типа. Описывает магнитные и электрические поля, гидродинамику, диффузию и др.
В этих уравнениях функция u зависит от двух переменных, однако, задача может быть расширена для случая трех переменных:
1) Волновое уравнение:
2) Уравнение теплопроводности:
3) Уравнение Лапласа:
Рассмотрим подробнее каждое из этих уравнений.
Уравнение колебаний струны.
Определение. В математической физике струной называется тонкая нить, в которой возможно возникновение напряжений только в продольном, но не в поперечном направлении.
Пусть концы натянутой струны закреплены в точках х = а и x = b, возникающие в ней напряжения обозначим Т. Будем также считать, что плотность струны постоянна на всем ее протяжении.
Допустим, что в момент t0 = 0 струна выведена из состояния равновесия и совершает малые колебания.
Отклонение струны в каждой точке с координатой х в момент времени t обозначим как
u
C
B
A
D
0 a x x+x b x
На произвольный элемент длины нити (х, х + х) действуют две силы натяжения
и . При этом:
Если считать колебания малыми, то можно принять:
Тогда проекция силы на ось u:
Проекция силы на ось u:
Находим сумму этих проекций:
Выражение, стоящее в правой части равенства получено в результате применения теоремы Лагранжа ( см. Теорема Лагранжа ) к выражению, стоящему слева.
Произведение массы на ускорение рассматриваемого элемента струны равно:
где - плотность струны.
Приравнивая полученное выражение к значению проекции силы, получим:
Или
Для полного определения движения струны полученного уравнения недостаточно. Функция u(x, t) должна еще удовлетворять граничным условиям, описывающим состояние струны на концах (в точках x = a и x = b) и начальным условиям, описывающим состояние струны в момент времени t = 0.
Совокупность граничных и начальных условий называется краевыми условиями.
Таким образом, задача Коши состоит в нахождении решения линейного дифференциального уравнения с частными производными второго порядка при начальных условиях
и краевых условиях
.
Начальные условия показывают, в каком положении находится струна в начальный момент времени и скорость каждой ее точки в начальный момент времени.
Функции f(x) и F(x) заданы.
Краевые условия показывают, что концы струны закреплены в точках a = 0, b = l
Решение задачи Коши методом разделения переменных.
(Метод Фурье.)
Решение уравнения
будем искать в виде при граничных условиях:
Тогда X(0) = X(l) = 0.
Подставим решение в исходное уравнение:
Можно показать, что функции Х и Т имеют вид:
Все решения исходного дифференциального уравнения, удовлетворяющие граничным условиям, можно записать в виде:
Окончательно решение уравнения колебаний струны можно записать в виде:
где
Решение задачи Коши методом Даламбера.
( Жан Лерон ДЛамбер (1717 1783) французский математик)
В случае если длина струны очень велика, то на колебания, возникающие в середине струны, концы струны влияния практически не оказывают. Поэтому, рассматривая колебания бесконечной струны, уравнение
решается только при начальных условиях:
Для нахождения решения введем новые переменные:
Тогда исходное уравнение принимает вид:
Решением этого уравнения будет функция , где и - некоторые функции, которые будем считать дважды дифференцируемыми.
Получаем:
Если продифференцировать полученный ответ, получим:
Т.е. .
Далее с использованием начальных условий находим функции и .
Проинтегрировав последнее равенство на отрезке [0, x], получаем:
Тогда:
Решение задачи Коши получаем в виде:
Эта формула называется формулой Даламбера.
Уравнение теплопроводности.
Температуру физического тела в произвольной точке с координатами (x, y, z) в момент времени t можно представить в виде функции:
Составим дифференциальное уравнение:
Выражение называется оператором Лапласа.
Тогда составленное нами дифференциальное уравнение принимает вид:
и называется уравнением теплопроводности в пространстве.
В качестве частных случаев рассматривают:
- уравнение теплопроводности в стержне,
- уравнение теплопроводности на плоскости.
В случае рассмотрения уравнения теплопроводности в стержне искомая функция u(x, t) должна удовлетворять записанному выше дифференциальному уравнению, начальному условию и граничным условиям .
В результате решения дифференциального уравнения методом Фурье получим:
Отметим, что распространение тепла в теле называется стационарным, если функция u не зависит от времени t.
Уравнение Лапласа.
Определение. Функция называется гармонической на области , если она имеет непрерывные частные производные второго порядка на области и удовлетворяет условию
,
где - оператор Лапласа.
Уравнение называется уравнением Лапласа.
Если на некоторой границе Г тела поддерживать постоянную температуру , где f заданная функция, то внутри тела установится единственная постоянная температура. С физической точки зрения это утверждение очевидно, однако, данный факт может быть доказан математически.
Математическое доказательство этого факта называется задачей Дирихле.
(Петер Густав Дирихле (1805 1859) немецкий математик)
Решение задачи Дирихле для круга.
Пусть в плоскости XOY имеется круг радиуса R с центром в начале координат и на его окружности задана функция f(), где - полярный угол.
Требуется найти функцию , которая удовлетворяет уравнению Лапласа
и при
Запишем уравнение Лапласа в полярных координатах:
Полагаем Подставляя это соотношение в уравнение Лапласа, получаем:
Таким образом, имеем два уравнения:
Общее решение первого уравнения имеет вид:
Решение второго уравнения ищем в виде: . При подстановке получим:
Общее решение второго уравнения имеет вид: .
Подставляя полученные решения в уравнение , получим:
Эта функция будет решением уравнения Лапласа при любом k 0.
Если k = 0, то следовательно .
Решение должно быть периодическим, т.к. одно и то же значение будет повторяться через 2. (Тогда рассматривается одна и та же точка круга.) Поэтому В0 = 0.
Решение должно быть конечным и непрерывным, поэтому D0 = 0.
Окончательно получаем:
При этом:
Если подставить эти коэффициенты в полученную выше формулу и произвести упрощение, получаем окончательный результат решения задачи Дирихле, который называется интегралом Пуассона.
(Симеон Дени Пуассон (1781 1840) французский математик)
Ряды.
Основные определения.
Определение. Сумма членов бесконечной числовой последовательности называется числовым рядом.
При этом числа будем называть членами ряда, а un общим членом ряда.
Определение. Суммы , n = 1, 2, … называются частными (частичными) суммами ряда.
Таким образом, возможно рассматривать последовательности частичных сумм ряда S1, S2, …,Sn, …
Определение. Ряд называется сходящимся, если сходится последовательность его частных сумм. Сумма сходящегося ряда предел последовательности его частных сумм.
Определение. Если последовательность частных сумм ряда расходится, т.е. не имеет предела, или имеет бесконечный предел, то ряд называется расходящимся и ему не ставят в соответствие никакой суммы.
Свойства рядов.
1) Сходимость или расходимость ряда не нарушится если изменить, отбросить или добавить конечное число членов ряда.
2) Рассмотрим два ряда и , где С постоянное число.
Теорема. Если ряд сходится и его сумма равна S, то ряд тоже сходится, и его сумма равна СS. (C 0)
3) Рассмотрим два ряда и . Суммой или разностью этих рядов будет называться ряд , где элементы получены в результате сложения (вычитания) исходных элементов с одинаковыми номерами.
Теорема. Если ряды и сходятся и их суммы равны соответственно S и , то ряд тоже сходится и его сумма равна S + .
Разность двух сходящихся рядов также будет сходящимся рядом.
Сумма сходящегося и расходящегося рядов будет расходящимся рядом.
О сумме двух расходящихся рядов общего утверждения сделать нельзя.
При изучении рядов решают в основном две задачи: исследование на сходимость и нахождение суммы ряда.
Критерий Коши.
(необходимые и достаточные условия сходимости ряда)
Для того, чтобы последовательность была сходящейся, необходимо и достаточно, чтобы для любого существовал такой номер N, что при n > N и любом p > 0, где р целое число, выполнялось бы неравенство:
.
Доказательство. (необходимость)
Пусть , тогда для любого числа найдется номер N такой, что неравенство
выполняется при n>N. При n>N и любом целом p>0 выполняется также неравенство . Учитывая оба неравенства, получаем:
Необходимость доказана. Доказательство достаточности рассматривать не будем.
Сформулируем критерий Коши для ряда.
Для того, чтобы ряд был сходящимся необходимо и достаточно, чтобы для любого существовал номер N такой, что при n>N и любом p>0 выполнялось бы неравенство
.
Однако, на практике использовать непосредственно критерий Коши не очень удобно. Поэтому как правило используются более простые признаки сходимости:
1) Если ряд сходится, то необходимо, чтобы общий член un стремился к нулю. Однако, это условие не является достаточным. Можно говорить только о том, что если общий член не стремится к нулю, то ряд точно расходится. Например, так называемый гармонический ряд является расходящимся, хотя его общий член и стремится к нулю.
Пример. Исследовать сходимость ряда
Найдем - необходимый признак сходимости не выполняется, значит ряд расходится.
2) Если ряд сходится, то последовательность его частных сумм ограничена.
Однако, этот признак также не является достаточным.
Например, ряд 1-1+1-1+1-1+ … +(-1)n+1+… расходится, т.к. расходится последовательность его частных сумм в силу того, что
Однако, при этом последовательность частных сумм ограничена, т.к. при любом n.
Ряды с неотрицательными членами.
При изучении знакопостоянных рядов ограничимся рассмотрением рядов с неотрицательными членами, т.к. при простом умножении на 1 из этих рядов можно получить ряды с отрицательными членами.
Теорема. Для сходимости ряда с неотрицательными членами необходимо и достаточно, чтобы частные суммы ряда были ограничены.
Признак сравнения рядов с неотрицательными членами.
Пусть даны два ряда и при un, vn 0.
Теорема. Если un vn при любом n, то из сходимости ряда следует сходимость ряда , а из расходимости ряда следует расходимость ряда .
Доказательство. Обозначим через Sn и n частные суммы рядов и . Т.к. по условию теоремы ряд сходится, то его частные суммы ограничены, т.е. при всех n n M, где М некоторое число. Но т.к. un vn, то Sn n то частные суммы ряда тоже ограничены, а этого достаточно для сходимости.
Пример. Исследовать на сходимость ряд
Т.к. , а гармонический ряд расходится, то расходится и ряд .
Пример. Исследовать на сходимость ряд
Т.к. , а ряд сходится ( как убывающая геометрическая прогрессия), то ряд тоже сходится.
Также используется следующий признак сходимости:
Теорема. Если и существует предел , где h число, отличное от нуля, то ряды и ведут одинаково в смысле сходимости.
Признак Даламбера.
(Жан Лерон Даламбер (1717 1783) французский математик)
Если для ряда с положительными членами существует такое число q<1, что для всех достаточно больших n выполняется неравенство
то ряд сходится, если же для всех достаточно больших n выполняется условие
то ряд расходится.
Предельный признак Даламбера.
Предельный признак Даламбера является следствием из приведенного выше признака Даламбера.
Если существует предел , то при < 1 ряд сходится, а при > 1 расходится. Если = 1, то на вопрос о сходимости ответить нельзя.
Пример. Определить сходимость ряда .
Вывод: ряд сходится.
Пример. Определить сходимость ряда
Вывод: ряд сходится.
Признак Коши. (радикальный признак)
Если для ряда с неотрицательными членами существует такое число q<1, что для всех достаточно больших n выполняется неравенство
,
то ряд сходится, если же для всех достаточно больших n выполняется неравенство
то ряд расходится.
Следствие. Если существует предел , то при <1 ряд сходится, а при >1 ряд расходится.
Пример. Определить сходимость ряда .
Вывод: ряд сходится.
Пример. Определить сходимость ряда .
Т.е. признак Коши не дает ответа на вопрос о сходимости ряда. Проверим выполнение необходимых условий сходимости. Как было сказано выше, если ряд сходится, то общий член ряда стремится к нулю.
,
таким образом, необходимое условие сходимости не выполняется, значит, ряд расходится.
Интегральный признак Коши.
Если (х) непрерывная положительная функция, убывающая на промежутке [1;), то ряд (1) + (2) + …+ (n) + … = и несобственный интеграл одинаковы в смысле сходимости.
Пример. Ряд сходится при >1 и расходится 1 т.к. соответствующий несобственный интеграл сходится при >1 и расходится 1. Ряд называется общегармоническим рядом.
Следствие. Если f(x) и (х) непрерывные функции на интервале (a, b] и то интегралы и ведут себя одинаково в смысле сходимости.
При использовании компьютерной версии “Курса высшей математики” возможно запустить программу, исследующую на сходимость числовые ряды по всем рассмотренным выше признакам. Достаточно ввести общий член ряда и нажать Enter. Все признаки будут проверяться по очереди.
Для запуска программы дважды щелкните на значке:
Примечание: Для запуска программы необходимо чтобы на компьютере была установлена программа Maple ( Waterloo Maple Inc.) любой версии, начиная с Maple V Release 4.
Знакопеременные ряды.
Знакочередующиеся ряды.
Знакочередующийся ряд можно записать в виде:
где
Признак Лейбница.
Если у знакочередующегося ряда абсолютные величины ui убывают и общий член стремится к нулю , то ряд сходится.
Абсолютная и условная сходимость рядов.
Рассмотрим некоторый знакопеременный ряд (с членами произвольных знаков).
(1)
и ряд, составленный из абсолютных величин членов ряда (1):
(2)
Теорема. Из сходимости ряда (2) следует сходимость ряда (1).
Доказательство. Ряд (2) является рядом с неотрицательными членами. Если ряд (2) сходится, то по критерию Коши для любого >0 существует число N, такое, что при n>N и любом целом p>0 верно неравенство:
По свойству абсолютных величин:
То есть по критерию Коши из сходимости ряда (2) следует сходимость ряда (1).
Определение. Ряд называется абсолютно сходящимся, если сходится ряд .
Очевидно, что для знакопостоянных рядов понятия сходимости и абсолютной сходимости совпадают.
Определение. Ряд называется условно сходящимся, если он сходится, а ряд расходится.
Признаки Даламбера и Коши для знакопеременных рядов.
Пусть - знакопеременный ряд.
Признак Даламбера. Если существует предел , то при <1 ряд будет абсолютно сходящимся, а при >1 ряд будет расходящимся. При =1 признак не дает ответа о сходимости ряда.
Признак Коши. Если существует предел , то при <1 ряд будет абсолютно сходящимся, а при >1 ряд будет расходящимся. При =1 признак не дает ответа о сходимости ряда.
Свойства абсолютно сходящихся рядов.
1) Теорема. Для абсолютной сходимости ряда необходимо и достаточно, чтобы его можно было представить в виде разности двух сходящихся рядов с неотрицательными членами.
Следствие. Условно сходящийся ряд является разностью двух расходящихся рядов с неотрицательными стремящимися к нулю членами.
2) В сходящемся ряде любая группировка членов ряда, не изменяющая их порядка, сохраняет сходимость и величину ряда.
3) Если ряд сходится абсолютно, то ряд, полученный из него любой перестановкой членов, также абсолютно сходится и имеет ту же сумму.
Перестановкой членов условно сходящегося ряда можно получить условно сходящийся ряд, имеющий любую наперед заданную сумму, и даже расходящийся ряд.
4) Теорема. При любой группировке членов абсолютно сходящегося ряда (при этом число групп может быть как конечным, так и бесконечным и число членов в группе может быть как конечным, так и бесконечным) получается сходящийся ряд, сумма которого равна сумме исходного ряда.
5) Если ряды и сходятся абсолютно и их суммы равны соответственно S и , то ряд, составленный из всех произведений вида взятых в каком угодно порядке, также сходится абсолютно и его сумма равна S - произведению сумм перемножаемых рядов.
Если же производить перемножение условно сходящихся рядов, то в результате можно получить расходящийся ряд.
Функциональные последовательности.
Определение. Если членами ряда будут не числа, а функции от х, то ряд называется функциональным.
Исследование на сходимость функциональных рядов сложнее исследования числовых рядов. Один и тот же функциональный ряд может при одних значениях переменной х сходиться, а при других расходиться. Поэтому вопрос сходимости функциональных рядов сводится к определению тех значений переменной х, при которых ряд сходится.
Совокупность таких значений называется областью сходимости.
Так как пределом каждой функции, входящей в область сходимости ряда, является некоторое число, то пределом функциональной последовательности будет являться некоторая функция:
Определение. Последовательность {fn(x)} сходится к функции f(x) на отрезке [a,b], если для любого числа >0 и любой точки х из рассматриваемого отрезка существует номер N = N(, x), такой, что неравенство
выполняется при n>N.
При выбранном значении >0 каждой точке отрезка [a,b] соответствует свой номер и, следовательно, номеров, соответствующих всем точкам отрезка [a,b], будет бесчисленное множество. Если выбрать из всех этих номеров наибольший, то этот номер будет годиться для всех точек отрезка [a,b], т.е. будет общим для всех точек.
Определение. Последовательность {fn(x)} равномерно сходится к функции f(x) на отрезке [a,b], если для любого числа >0 существует номер N = N(), такой, что неравенство
выполняется при n>N для всех точек отрезка [a,b].
Пример. Рассмотрим последовательность
Данная последовательность сходится на всей числовой оси к функции f(x)=0, т.к.
Построим графики этой последовательности:
sinx
Как видно, при увеличении числа n график последовательности приближается к оси х.
При использовании компьютерной версии “Курса высшей математики” возможно запустить программу, которая исследует на сходимость знакочередующиеся ряды и определяет характер сходимости. Достаточно ввести общий член ряда и множитель, определяющий знак и нажать Enter. Все рассмотренные выше признаки будут проверены по очереди.
Для запуска программы дважды щелкните на значке:
Примечание: Для запуска программы необходимо чтобы на компьютере была установлена программа Maple ( Waterloo Maple Inc.) любой версии, начиная с Maple V Release 4.
Функциональные ряды.
Определение. Частными (частичными) суммами функционального ряда называются функции
Определение. Функциональный ряд называется сходящимся в точке (х=х0), если в этой точке сходится последовательность его частных сумм. Предел последовательности называется суммой ряда в точке х0.
Определение. Совокупность всех значений х, для которых сходится ряд называется областью сходимости ряда.
Определение. Ряд называется равномерно сходящимся на отрезке [a,b], если равномерно сходится на этом отрезке последовательность частных сумм этого ряда.
Теорема. (Критерий Коши равномерной сходимости ряда)
Для равномерной сходимости ряда необходимо и достаточно, чтобы для любого числа >0 существовал такой номер N(), что при n>N и любом целом p>0 неравенство
выполнялось бы для всех х на отрезке [a,b].
Теорема. (Признак равномерной сходимости Вейерштрасса)
(Карл Теодор Вильгельм Вейерштрасс (1815 1897) немецкий математик)
Ряд сходится равномерно и притом абсолютно на отрезке [a,b], если модули его членов на том же отрезке не превосходят соответствующих членов сходящегося числового ряда с положительными членами :
т.е. имеет место неравенство:
.
Еще говорят, что в этом случае функциональный ряд мажорируется числовым рядом .
Пример. Исследовать на сходимость ряд .
Так как всегда, то очевидно, что .
При этом известно, что общегармонический ряд при =3>1 сходится, то в соответствии с признаком Вейерштрасса исследуемый ряд равномерно сходится и притом в любом интервале.
Пример. Исследовать на сходимость ряд .
На отрезке [-1,1] выполняется неравенство т.е. по признаку Вейерштрасса на этом отрезке исследуемый ряд сходится, а на интервалах (-, -1) (1, ) расходится.
Свойства равномерно сходящихся рядов.
1) Теорема о непрерывности суммы ряда.
Если члены ряда - непрерывные на отрезке [a,b] функции и ряд сходится равномерно, то и его сумма S(x) есть непрерывная функция на отрезке [a,b].
2) Теорема о почленном интегрировании ряда.
Равномерно сходящийся на отрезке [a,b] ряд с непрерывными членами можно почленно интегрировать на этом отрезке, т.е. ряд, составленный из интегралов от его членов по отрезку [a,b] , сходится к интегралу от суммы ряда по этому отрезку.
3) Теорема о почленном дифференцировании ряда.
Если члены ряда сходящегося на отрезке [a,b] представляют собой непрерывные функции, имеющие непрерывные производные, и ряд, составленный из этих производных сходится на этом отрезке равномерно, то и данный ряд сходится равномерно и его можно дифференцировать почленно.
На основе того, что сумма ряда является некоторой функцией от переменной х, можно производить операцию представления какой либо функции в виде ряда (разложения функции в ряд), что имеет широкое применение при интегрировании, дифференцировании и других действиях с функциями.
На практике часто применяется разложение функций в степенной ряд.
При использовании компьютерной версии “Курса высшей математики” возможно запустить программу, которая определяет интервал сходимости для произвольного функционального ряда.
Для запуска программы дважды щелкните на значке
Примечание: Для запуска программы необходимо чтобы на компьютере была установлена программа Maple ( Waterloo Maple Inc.) любой версии, начиная с MapleV Release 4.
Степенные ряды.
Определение. Степенным рядом называется ряд вида
.
Для исследования на сходимость степенных рядов удобно использовать признак Даламбера.
Пример. Исследовать на сходимость ряд
Применяем признак Даламбера:
.
Получаем, что этот ряд сходится при и расходится при .
Теперь определим сходимость в граничных точках 1 и 1.
При х = -1: ряд сходится по признаку Лейбница (см. Признак Лейбница.).
При х = 1: ряд расходится (гармонический ряд).
Теоремы Абеля.
(Нильс Хенрик Абель (1802 1829) норвежский математик)
Теорема. Если степенной ряд сходится при x = x1 , то он сходится и притом абсолютно для всех .
Доказательство. По условию теоремы, так как члены ряда ограничены, то
где k- некоторое постоянное число. Справедливо следующее неравенство:
Из этого неравенства видно, что при x<x1 численные величины членов нашего ряда будут меньше ( во всяком случае не больше ) соответствующих членов ряда правой части записанного выше неравенства, которые образуют геометрическую прогрессию. Знаменатель этой прогрессии по условию теоремы меньше единицы, следовательно, эта прогрессия представляет собой сходящийся ряд.
Поэтому на основании признака сравнения делаем вывод, что ряд сходится, а значит ряд сходится абсолютно.
Таким образом, если степенной ряд сходится в точке х1, то он абсолютно сходится в любой точке интервала длины 2 с центром в точке х = 0.
Следствие. Если при х = х1 ряд расходится, то он расходится для всех .
Таким образом, для каждого степенного ряда существует такое положительное число R, что при всех х таких, что ряд абсолютно сходится, а при всех ряд расходится. При этом число R называется радиусом сходимости. Интервал (-R, R) называется интервалом сходимости.
Отметим, что этот интервал может быть как замкнутым с одной или двух сторон, так и не замкнутым.
Радиус сходимости может быть найден по формуле:
Пример. Найти область сходимости ряда
Находим радиус сходимости .
Следовательно, данный ряд сходится при любом значении х. Общий член этого ряда стремится к нулю.
Теорема. Если степенной ряд сходится для положительного значения х=х1 , то он сходится равномерно в любом промежутке внутри .
Действия со степенными рядами.
1) Интегрирование степенных рядов.
Если некоторая функция f(x) определяется степенным рядом: , то интеграл от этой функции можно записать в виде ряда:
2) Дифференцирование степенных рядов.
Производная функции, которая определяется степенным рядом, находится по формуле:
3) Сложение, вычитание, умножение и деление степенных рядов.
Сложение и вычитание степенных рядов сводится к соответствующим операциям с их членами:
Произведение двух степенных рядов выражается формулой:
Коэффициенты сi находятся по формуле:
Деление двух степенных рядов выражается формулой:
Для определения коэффициентов qn рассматриваем произведение , полученное из записанного выше равенства и решаем систему уравнений:
Разложение функций в степенные ряды.
Разложение функций в степенной ряд имеет большое значение для решения различных задач исследования функций, дифференцирования, интегрирования, решения дифференциальных уравнений, вычисления пределов, вычисления приближенных значений функции.
Возможны различные способы разложения функции в степенной ряд. Такие способы как разложение при помощи рядов Тейлора и Маклорена были рассмотрены ранее. (См. Формула Тейлора. )
Существует также способ разложения в степенной ряд при помощи алгебраического деления. Это самый простой способ разложения, однако, пригоден он только для разложения в ряд алгебраических дробей.
Пример. Разложить в ряд функцию .
Суть метода алгебраического деления состоит в применении общего правила деления многочленов:
1 1 - x
1 x 1 + x + x2 + x3 + …
x
x x2
x2
x2 x3
x3
……….
Если применить к той же функции формулу Маклорена
,
то получаем:
……………………………….
Итого, получаем:
Рассмотрим способ разложения функции в ряд при помощи интегрирования.
С помощью интегрирования можно разлагать в ряд такую функцию, для которой известно или может быть легко найдено разложение в ряд ее производной.
Находим дифференциал функции и интегрируем его в пределах от 0 до х.
Пример. Разложить в ряд функцию
Разложение в ряд этой функции по формуле Маклорена было рассмотрено выше.
(См. Функция y = ln(1 + x).) Теперь решим эту задачу при помощи интегрирования.
При получаем по приведенной выше формуле:
Разложение в ряд функции может быть легко найдено способом алгебраического деления аналогично рассмотренному выше примеру.
Тогда получаем:
Окончательно получим:
Пример. Разложить в степенной ряд функцию .
Применим разложение в ряд с помощью интегрирования.
Подинтегральная функция может быть разложена в ряд методом алгебраического деления:
1 1 + x2
1 + x2 1 x2 + x4- …
- x2
- x2 x4
x4
x4 + x6
………….
Тогда
Окончательно получаем:
степенных рядов.
С помощью степенных рядов возможно интегрировать дифференциальные уравнения.
Рассмотрим линейное дифференциальное уравнение вида:
Если все коэффициенты и правая часть этого уравнения разлагаются в сходящиеся в некотором интервале степенные ряды, то существует решение этого уравнения в некоторой малой окрестности нулевой точки, удовлетворяющее начальным условиям.
Это решение можно представить степенным рядом:
Для нахождения решения остается определить неизвестные постоянные ci.
Эта задача решается методом сравнения неопределенных коэффициентов. Записанное выражение для искомой функции подставляем в исходное дифференциальное уравнение, выполняя при этом все необходимые действия со степенными рядами (дифференцирование, сложение, вычитание, умножение и пр.)
Затем приравниваем коэффициенты при одинаковых степенях х в левой и правой частях уравнения. В результате с учетом начальных условий получим систему уравнений, из которой последовательно определяем коэффициенты ci.
Отметим, что этот метод применим и к нелинейным дифференциальным уравнениям.
Пример. Найти решение уравнения c начальными условиями y(0)=1, y(0)=0.
Решение уравнения будем искать в виде
Подставляем полученные выражения в исходное уравнение:
Отсюда получаем:
………………
Получаем, подставив начальные условия в выражения для искомой функции и ее первой производной:
Окончательно получим:
Итого:
Существует и другой метод решения дифференциальных уравнений с помощью рядов. Он носит название метод последовательного дифференцирования.
Рассмотрим тот же пример. Решение дифференциального уравнения будем искать в виде разложения неизвестной функции в ряд Маклорена.
Если заданные начальные условия y(0)=1, y(0)=0 подставить в исходное дифференциальное уравнение, получим, что
Далее запишем дифференциальное уравнение в виде и будем последовательно дифференцировать его по х.
После подстановки полученных значений получаем:
Ряды Фурье.
( Жан Батист Жозеф Фурье (1768 1830) французский математик)
Тригонометрический ряд.
Определение. Тригонометрическим рядом называется ряд вида:
или, короче,
Действительные числа ai, bi называются коэффициентами тригонометрического ряда.
Если ряд представленного выше типа сходится, то его сумма представляет собой периодическую функцию с периодом 2, т.к. функции sinnx и cosnx также периодические функции с периодом 2.
Пусть тригонометрический ряд равномерно сходится на отрезке [-; ], а следовательно, и на любом отрезке в силу периодичности, и его сумма равна f(x).
Определим коэффициенты этого ряда.
Для решения этой задачи воспользуемся следующими равенствами:
Справедливость этих равенств вытекает из применения к подынтегральному выражению тригонометрических формул. Подробнее см. Интегрирование тригонометрических функций.
Т.к. функция f(x) непрерывна на отрезке [-; ], то существует интеграл
Такой результат получается в результате того, что .
Получаем:
Далее умножаем выражение разложения функции в ряд на cosnx и интегрируем в пределах от - до .
Отсюда получаем:
Аналогично умножаем выражение разложения функции в ряд на sinnx и интегрируем в пределах от - до .
Получаем:
Выражение для коэффициента а0 является частным случаем для выражения коэффициентов an.
Таким образом, если функция f(x) любая периодическая функция периода 2, непрерывная на отрезке [-; ] или имеющая на этом отрезке конечное число точек разрыва первого рода, то коэффициенты
существуют и называются коэффициентами Фурье для функции f(x).
Определение. Рядом Фурье для функции f(x) называется тригонометрический ряд, коэффициенты которого являются коэффициентами Фурье. Если ряд Фурье функции f(x) сходится к ней во всех ее точках непрерывности, то говорят, что функция f(x) разлагается в ряд Фурье.
Достаточные признаки разложимости в ряд Фурье.
Теорема. (Теорема Дирихле) Если функция f(x) имеет период 2 и на отрезке
[-;] непрерывна или имеет конечное число точек разрыва первого рода, и отрезок
[-;] можно разбить на конечное число отрезков так, что внутри каждого из них функция f(x) монотонна, то ряд Фурье для функции f(x) сходится при всех значениях х, причем в точках непрерывности функции f(x) его сумма равна f(x), а в точках разрыва его сумма равна , т.е. среднему арифметическому предельных значений слева и справа. При этом ряд Фурье функции f(x) сходится равномерно на любом отрезке, который принадлежит интервалу непрерывности функции f(x).
Функция f(x), для которой выполняются условия теоремы Дирихле называется кусочно монотонной на отрезке [-;].
Теорема. Если функция f(x) имеет период 2, кроме того, f(x) и ее производная f(x) непрерывные функции на отрезке [-;] или имеют конечное число точек разрыва первого рода на этом отрезке, то ряд Фурье функции f(x) сходится при всех значениях х, причем в точках непрерывности его сумма равна f(x), а в точках разрыва она равна . При этом ряд Фурье функции f(x) сходится равномерно на любом отрезке, который принадлежит интервалу непрерывности функции f(x).
Функция, удовлетворяющая условиям этой теоремы, называется кусочно гладкой на отрезке [-;].
Разложение в ряд Фурье непериодической функции.
Задача разложения непериодической функции в ряд Фурье в принципе не отличается от разложения в ряд Фурье периодической функции.
Допустим, функция f(x) задана на отрезке [a, b] и является на этом отрезке кусочно монотонной. Рассмотрим произвольную периодическую кусочно монотонную функцию f1(x) c периодом 2Т b-a, совпадающую с функцией f(x) на отрезке [a, b].
y
f(x)
- 2T a b +2T + 4T x
Таким образом, функция f(x) была дополнена. Теперь функция f1(x) разлагается в ряд Фурье. Сумма этого ряда во всех точках отрезка [a, b] совпадает с функцией f(x), т.е. можно считать, что функция f(x) разложена в ряд Фурье на отрезке [a, b].
Таким образом, если функция f(x) задана на отрезке, равном 2 ничем не отличается от разложения в ряд периодической функции. Если же отрезок, на котором задана функция, меньше, чем 2, то функция продолжается на интервал (b, a + 2) так, что условия разложимости в ряд Фурье сохранялись.
Вообще говоря, в этом случае продолжение заданной функции на отрезок (интервал) длиной 2 может быть произведено бесконечным количеством способов, поэтому суммы получившихся рядов будут различны, но они будут совпадать с заданной функцией f(x) на отрезке [a,b].
Ряд Фурье для четных и нечетных функций.
Отметим следующие свойства четных и нечетных функций:
1)
2) Произведение двух четных и нечетных функций является четной функцией.
3) Произведение четной и нечетной функций нечетная функция.
Справедливость этих свойств может быть легко доказана исходя из определения четности и нечетности функций.
Если f(x) четная периодическая функция с периодом 2, удовлетворяющая условиям разложимости в ряд Фурье, то можно записать:
Таким образом, для четной функции ряд Фурье записывается:
Аналогично получаем разложение в ряд Фурье для нечетной функции:
Пример. Разложить в ряд Фурье периодическую функцию с периодом T = 2 на отрезке [-;].
Заданная функция является нечетной, следовательно, коэффициенты Фурье ищем в виде:
Получаем:
.
Построим графики заданной функции и ее разложения в ряд Фурье, ограничившись первыми четырьмя членами ряда.
Ряды Фурье для функций любого периода.
Ряд Фурье для функции f(x) периода Т = 2l, непрерывной или имеющей конечное число точек разрыва первого рода на отрезке [-l, l] имеет вид:
Для четной функции произвольного периода разложение в ряд Фурье имеет вид:
Для нечетной функции:
Ряд Фурье по ортогональной системе функций.
Определение. Функции (х) и (х), определенные на отрезке [a, b], называются ортогональными на этом отрезке, если
Определение. Последовательность функций 1(x), 2(x), …, n(x), непрерывных на отрезке [a, b], называется ортогональной системой функций на этом отрезке, если все функции попарно ортогональны.
Отметим, что ортогональность функций не подразумевает перпендикулярности графиков этих функций.
Определение. Система функций называется ортогональной и нормированной (ортонормированной), если
Определение. Рядом Фурье по ортогональной системе функций 1(x), 2(x), …,n(x) называется ряд вида:
коэффициенты которого определяются по формуле:
,
где f(x) = - сумма равномерно сходящегося на отрезке [a, b] ряда по ортогональной системе функций. f(x) любая функция, непрерывная или имеющая конечное число точек разрыва первого рода на отрезке [a, b].
В случае ортонормированной системы функций коэффициенты определяются:
При использовании компьютерной версии “Курса высшей математики” возможно запустить программу, которая разлагает в ряд Фурье произвольную функцию.
Для запуска программы дважды щелкните на значке
Примечание: Для запуска программы необходимо чтобы на компьютере была установлена программа Maple ( Waterloo Maple Inc.) любой версии, начиная с MapleV Release 4.
Интеграл Фурье.
Пусть функция f(x) на каждом отрезке [-l,l], где l любое число, кусочно гладкая или кусочно монотонная, кроме того, f(x) абсолютно интегрируемая функция, т.е. сходится несобственный интеграл
Тогда функция f(x) разлагается в ряд Фурье:
Если подставить коэффициенты в формулу для f(x), получим:
Переходя к пределу при l, можно доказать, что и
Обозначим
При l un 0.
Можно доказать, что предел суммы, стоящий в правой части равенства равен интегралу
Тогда - двойной интеграл Фурье.
Окончательно получаем:
- представление функции f(x) интегралом Фурье.
Двойной интеграл Фурье для функции f(x) можно представить в комплексной форме:
Преобразование Фурье.
Определение. Если f(x) любая абсолютно интегрируемая на всей числовой оси функция, непрерывная или имеющая конечное число точек разрыва первого рода на каждом отрезке, то функция
называется преобразованием Фурье функции f(x).
Функция F(u) называется также спектральной характеристикой функции f(x).
Если f(x) функция, представимая интегралом Фурье, то можно записать:
Это равенство называется обратным преобразованием Фурье
Интегралы и называются соответственно косинус - преобразование Фурье и синус преобразование Фурье.
Косинус преобразование Фурье будет преобразованием Фурье для четных функций, синус преобразование для нечетных.
Преобразование Фурье применяется в функциональном анализе, гармоническом анализе, операционном исчислении, теории линейных систем и др.
Элементы теории функций комплексного переменного.
Определение. Если каждому комплексному числу z из некоторого множества D по некоторому закону поставлено в соответствие определенное комплексное число w из множества G, то на этой области задана однозначная функция комплексного переменного, отображающая множество D на множество G.
w = f(z)
Множество D называется областью определения, множество G областью значений функции.
Комплексную функцию можно записать в виде:
u, v действительные функции от переменных х и у.
Если каждому z D соответствует несколько различных значений w, то функция w=f(z) называется многозначной.
Определение. Функция имеет предел в точке z0, равный числу А = a + ib, если
Свойства функций комплексного переменного.
Для функций комплексного переменного f(z) и g(z) справедливы следующие свойства:
1)
2)
3)
Определение. Функция называется непрерывной в точке z0, если выполняется равенство
Основные трансцендентные функции.
Определение. Трансцендентными называются аналитические функции, которые не являются алгебраическими.
Если аргументом показательной или тригонометрических функций является комплексное число, то определение этих функций, вводимое в элементарной алгебре теряет смысл.
Рассмотрим разложение в степенной ряд следующих функций:
См. Представление функций по формуле Тейлора.
Функции ez, cosz, sinz связаны между собой формулой Эйлера (см. Уравнение Эйлера.) Эта формула может быть очень легко получена сложением соотвествующих рядов.
Также справедливы равенства:
Для тригонометрических функций комплексного аргумента справедливы основные тригонометрические тождества (синус и косинус суммы, разности и т.д.), которые справедливы для функций действительного аргумента.
Определение. Гиперболическим синусом, косинусом, тангенсом и котангенсом называются соответственно функции:
Гиперболические функции могут быть выражены через тригонометрические:
Гиперболические функции sh z и ch z имеют период 2i, а функции th z и cth z период i.
Пример. Найти sin(1+2i).
Определение. Логарифмическая функция комплексного аргумента определяется как функция, обратная показательной.
Если w = u + iv, то и Arg ew = = v.
Тогда eu = .
Итого:
Для комплексного числа z = a + ib
Определение. Выражение называется главным значением логарифма.
Логарифмическая функция комплексного аргумента обладает следующими свойствами:
1)
2)
3)
4)
Обратные тригонометрические функции комплексного переменного имеют вид:
Производная функций комплексного переменного.
Определение. Производной от однозначной функции w = f(z) в точке z называется предел:
Определение. Функция f(z), имеющая непрерывную производную в любой точке области D называется аналитической функцией на этой области.
Правила дифференцирования функций комплексного аргумента не отличаются от правил дифференцирования функций действительной переменной.
Аналогично определяются производные основных функций таких как синус, косинус, тангенс и котангенс, степенная функция и т.д.
Производные гиперболических функций определяются по формулам:
Вывод правил интегрирования, значений производных основных функций ничем не отличается от аналогичных операций с функциями действительного аргумента, поэтому подробно рассматривать их не будем.
Условия Коши Римана.
(Бернхард Риман (1826 1866) немецкий математик)
Рассмотрим функцию комплексной переменной , определенную на некоторой области и имеющую в какой либо точке этой области производную
Стремление к нулю z0 может осуществляться в следующих случаях:
1)
2)
В первом случае:
Во втором случае:
Тогда должны выполняться равенства:
Эти равенства называются условиями Коши Римана, хотя еще раньше они были получены Эйлером и Даламбером.
Теорема. Если функция имеет производную в точке
z = x + iy, то ее действительные компоненты u и v имеют в точке (х, у) частные производные первого порядка, удовлетворяющие условию Коши Римана.
Также справедлива и обратная теорема.
На основании этих теорем можно сделать вывод, что из существования производной следует непрерывность функции.
Теорема. Для того, чтобы функция была аналитической на некоторой области необходимо и достаточно, чтобы частные производные первого прядка функций u и v были непрерывны на этой области и выполнялись условия Коши Римана.
Интегрирование функций комплексной переменной.
Пусть - непрерывная функция комплексного переменного z, определенная в некоторой области и L кривая, лежащая в этой области.
у
В
L
А
х
Кривая L задана уравнением
Определение. Интеграл от функции f(z) вдоль кривой L определяется следующим образом:
Если учесть, что , то
Теорема. (Теорема Коши) Если f(z) - аналитическая функция на некоторой области, то интеграл от f(z) по любому кусочно гладкому контуру, принадлежащему этой области равен нулю.
Интегральная формула Коши.
Если функция f(z) аналитическая в односвязной замкнутой области с кусочно гладкой границей L.
D
z0
Тогда справедлива формула Коши:
где z0 любая точка внутри контура L, интегрирование по контуру производится в положительном направлении (против часовой стрелки).
Эта формула также называется интегралом Коши.
Ряды Тейлора и Лорана.
(Пьер Альфонс Лоран (1813 1854) французский математик)
Функция f(z), аналитическая в круге , разлагается в сходящийся к ней степенной ряд по степеням (z z0).
Коэффициенты ряда вычисляются по формулам:
Степенной ряд с коэффициентами такого вида называется рядом Тейлора.
Рассмотрим теперь функцию f(z), аналитическую в кольце . Эта функция может быть представлена в виде сходящегося ряда:
Ряд такого вида называется рядом Лорана. При этом функция f(z) может быть представлена в виде суммы:
Ряд, определяющий функцию f1(x), называется правильной частью ряда Лорана, а ряд, определяющий функцию f2(x), называется главной частью ряда Лорана.
Если предположить, что r = 0, то можно считать, что функция аналитична в открытом круге за исключением центральной точки z0. Как правило, в этой точке функция бывает не определена.
Тогда точка z0 называется изолированной особой точкой функции f.
Рассмотрим следующие частные случаи:
1) Функция f(x) имеет вид: . Т.к. степенной ряд сходится во всех точках внутри круга, то его сумма f1(x) определена и непрерывно дифференцируема во всех точках круга, а, следовательно, и в центре круга z0.
В этом случае говорят, что особенность функции f в точке z0 устранима. Для устранения особой точки достаточно доопределить функцию в центре круга (f(z0) = c0) и функция будет аналитической не только в окрестности центра круга, но и в самом центре.
В этом случае для любого контура L, содержащего точку z0 и принадлежащего к кругу .
2) Функция f(x) имеет вид: .
В этом случае точка z0 называется полюсом функции f(z) порядка (кратности) m. При m = 1 точку z0 называют еще простым полюсом.
Порядок полюса может быть определен по формуле:
z0 полюс порядка т.
3) Функция f(z) имеет вид , где в ряду не равно нулю бесконечное количество коэффициентов с-k.
В этом случае говорят, что функция f(z) имеет в точке z0 существенно особую точку.
Определение. Пусть z0 изолированная особая точка функция f(z), т.е. пусть функция f(z) аналитическая в некотором круге из которого исключена точка z0. Тогда интеграл
называется вычетом функции f(z) в точке z0, где L контур в круге , ориентированный против часовой стрелки и содержащей в себе точку z0.
Вычет также обозначают иногда .
Если есть ряд Лорана функции f в точке z0, то .
Таким образом, если известно разложение функции в ряд Лорана, то вычет легко может быть найден в случае любой особой точки.
В частных случаях вычет может быть найден и без разложения в ряд Лорана.
Например, если функция , а имеет простой нуль при z = z0 , то z = z0 является простым полюсом функции f(z).
Тогда можно показать, что вычет находится по формуле
Если z = z0 полюс порядка m 1, то вычет может быть найден по формуле:
Пример. Найти вычет функции относительно точки z = 2.
Эта точка является полюсом второго порядка. Получаем:
Теорема о вычетах.
Теорема. Пусть функция f(z) аналитическая на всей плоскости z, за исключением конечного числа точек z1, z2, …, zN. Тогда верно равенство:
А интеграл от функции по контуру L, содержащему внутри себя эти точки, равен
Эти свойства применяются для вычисления интегралов. Если функция f(z) аналитическая в верхней полуплоскости, включая действительную ось, за исключением N точек, то справедлива формула
Пример. Вычислить определенный интеграл .
Подынтегральная функция является аналитической в верхней полуплоскости за исключением точки 2i. Эта точка является полюсом второго порядка.
Найдем вычет функции
Получаем
Пример. Вычислить определенный интеграл
Подынтегральная функция является аналитической в верхней полуплоскости за исключением точки i. Эта точка является полюсом второго порядка.
Найдем вычет функции
Получаем
При использовании компьютерной версии “Курса высшей математики” возможно запустить программу, которая находит вычеты задаваемой функции.
Для запуска программы дважды щелкните на значке
Примечание: Для запуска программы необходимо чтобы на компьютере была установлена программа Maple ( Waterloo Maple Inc.) любой версии, начиная с MapleV Release 4.
Операционное исчисление.
Преобразование Лапласа.
(Пьер Симон Лаплас (1749 1825) французский математик)
Рассмотрим функцию действительного переменного t, определенную при t 0. Будем также считать, что функция f(t)- кусочно - непрерывная, т.е. в любом конечном интервале она имеет конечное число точек разрыва первого рода, и определена на бесконечном интервале (-, ), но f(t) = 0 при t < 0.
Будем считать, что функция ограничена условием:
Рассмотрим функцию
где p = a + ib комплексное число.
Определение. Функция F(p) называется изображением Лапласа функции f(t).
Также функцию F(p) называют L изображением или преобразованием Лапласа.
Обозначается
При этом функция f(t) называется начальной функцией или оригиналом, а процесс нахождения оригинала по известному изображению называется операционным исчислением.
Теорема. (Теорема единственности) Если две непрерывнные функции f(x) и g(x) имеют одно и то же L изображение F(p), то они тождественно равны.
Определение. Функцией Хевисайда (Оливер Хевисайд (1850 1925) английский физик) называется функция
Свойства изображений.
Если , то справедливы следующие свойства:
1) Свойство подобия.
2) Свойство линейности.
3) Смещение изображения.
4) Дифференцирование изображения.
5) Дифференцирование оригинала.
6) Интегрирование изображения.
(Справедливо при условии, что интеграл сходится)
7) Интегрирование оригинала.
Таблица изображений некоторых функций.
Для большинства функций изображение находится непосредственным интегрированием.
Пример. Найти изображение функции f(t) = sint.
Для многих функций изображения посчитаны и приведены в соответствующих таблицах.
№ |
f(t) |
F(p) |
№ |
f(t) |
F(p) |
1 |
1 |
9 |
|||
2 |
sint |
10 |
|||
3 |
cost |
11 |
|||
4 |
e-t |
12 |
|||
5 |
sht |
13 |
|||
6 |
cht |
14 |
|||
7 |
15 |
||||
8 |
16 |
* |
* - при условии, что
Теоремы свертки и запаздывания.
Теорема. (теорема запаздывания) Если f(t) = 0 при t < 0, то справедлива формула
где t0 некоторая точка.
Определение. Выражение называется сверткой функций f1(t) и f2(t) и обозначается f1 f2.
Теорема. (теорема свертки) Преобразование Лапласа от свертки равно произведению преобразований Лапласа от функций f1(t) и f2(t) .
Теорема. (Интеграл Дюамеля (Дюамель (1797 1872) французский математик)). Если , то верно равенство
Для нахождения изображений различных функций наряду с непосредственным интегрированием применяются приведенные выще теоремы и свойства.
Пример. Найти изображение функции .
Из таблицы изображений получаем: .
По свойству интегрирования изображения получаем:
Пример. Найти изображение функции .
Из тригонометрии известна формула .
Тогда =.
Операционное исчисление используется как для нахождения значений интегралов, так и для решение дифференциальных уравнений.
Пусть дано линейное дифференциальное уравнение с постоянными коэффициентами.
Требуется найти решение этого дифференциального уравнения, удовлетворяющее начальным условиям:
Если функция x(t) является решением этого дифференциального уравнения, то оно обращает исходное уравнение в тождество, значит функция, стоящая в левой части уравнения и функция f(t) имеет (по теореме единственности) одно и то же изображение Лапласа.
Из теоремы о дифференцировании оригинала {} можно сделать вывод, что
Тогда
Обозначим
Получаем:
Это уравнение называется вспомогательным (изображающим) или операторным уравнением.
Отсюда получаем изображение , а по нему и искомую функцию x(t).
Изображение получаем в виде:
Где
Этот многочлен зависит от начальных условий. Если эти условия нулевые, то многочлен равен нулю, и формула принимает вид:
Рассмотрим применение этого метода на примерах.
Пример. Решить уравнение
Изображение искомой функции будем искать в виде:
Находим оригинал, т.е. искомую функцию:
Пример. Решить уравнение
Пример. Решить уравнение:
Изображение искомой функции
Для нахождения оригинала необходимо разложить полученную дробь на элементарные дроби. Воспользуемся делением многочленов (знаменатель делится без остатка на p 1):
p3 6p2 + 11p 6 p - 1
p3 p2 p2 5p + 6
-5p2 + 11p
-5p2 + 5p
6p - 6
6p - 6
0
В свою очередь
Получаем:
Тогда:
Определим коэффициенты А, В и С.
Тогда
Приемы операционного исчисления можно также использовать для решения систем дифференциальных уравнений.
Пример. Решить систему уравнений:
Обозначим - изображения искомых функций и решим вспомогательные уравнения:
Решим полученную систему алгебраических уравнений.
Если применить к полученным результатам формулы
то ответ можно представить в виде:
Как видно, гиперболические функции в ответе могут быть легко заменены на показательные.
Пример. Решить систему уравнений
при x(0) = y(0) = 1
Составим систему вспомогательных уравнений:
Если обозначить то из полученного частного решения системы можно записать и общее решение:
При рассмотрении нормальных систем дифференциальных уравнений этот пример был решен традиционным способом (См. Другой способ решения.). Как видно, результаты совпадают.
Отметим, что операторный способ решения систем дифференциальных уравнений применим к системам порядка выше первого, что очень важно, т.к. в этом случае применение других способов крайне затруднительно.
Криволинейные интегралы.
Определение. Кривая () называется непрерывной кусочно гладкой, если функции , и непрерывны на отрезке [a,b] и отрезок [a,b] можно разбить на конечное число частичных отрезков так, что на каждом из них функции , и имеют непрерывные производные, не равные нулю одновременно.
Если определено не только разбиение кривой на частичные отрезки точками, но порядок этих точек, то кривая называется ориентированнной кривой.
Ориетированная кривая называется замкнутой, если значения уравнения кривой в начальной и конечной точках совпадают.
Рассмотрим в пространсве XYZ кривую АВ, в каждой точке которой определена произвольная функция .
Разобьем кривую на конечное число отрезков и рассмотрим произведение значения функции в каждой точке разбиения на длину соответствующего отрезка.
Сложив все полученные таким образом произведения, получим так называемую интегральнуюсумму функции f(x, y, z).
Определение. Если при стремлении к нулю шага разбиения кривой на частичные отрезки существует предел интегральных сумм, то этот предел называется криволинейным интегралом от функции f(x, y, z) по длине дуги АВ или криволинейным интегралом первого рода.
Свойства криволинейного интеграла первого рода.
1) Значение криволинейного интеграла по длине дуги не зависит от направления кривой АВ.
2) Постоянный множитель можно выносить за знак криволинейного интеграла.
3) Криволинейный интерал от суммы функций равен сумме криволинейных интегралов от этих функций.
4) Если кривая АВ разбита на дуга АС и СВ, то
5) Если в точках кривой АВ
то
6) Справедливо неравенство:
7) Если f(x, y, z) = 1, то
S длина дуги кривой, - наибольшая из всех частичных дуг, на которые разбивается дуга АВ.
8) Теорема о среднем.
Если функция f(x, y, z) непрерывна на кривой АВ, то на этой кривой существует точка (x1, y1, z1) такая, что
Для вычисления криволинейного интеграла по длине дуги надо определить его связь с обыкновенным определенным интегралом.
Пусть кривая АВ задана параметрически уравнениями x = x(t), y = y(t), z = z(t),
t , где функции х, у, z непрерывно дифференцируемые функции параметра t, причем точке А соответствует t = , а точке В соответствует t = . Функция f(x, y, z) непрерывна на всей кривой АВ.
Для любой точки М(х, у, z) кривой длина дуги АМ вычисляется по формуле
(См. Вычисление длины дуги кривой.):
Длина всей кривой АВ равна:
Криволинейный интеграл по длине дуги АВ будет находиться по формуле:
Таким образом, для вычисления криволинейного интеграла первого рода (по длине дуги АВ) надо, используя параметрическое уравнение кривой выразить подынтегральную функцию через параметр t, заменить ds дифференциалом дуги в зависимости от параметра t и проинтегрировать полученное выражение по t.
Пример. Вычислить интеграл по одному витку винтовой линии
Если интегрирование производится по длине плоской кривой, заданной уравнением то получаем:
Криволинейные интегралы второго рода.
Пусть АВ непрерывная кривая в пространстве XYZ (или на плоскости ХОY), а точка P(x, y, z) произвольная функция, определенная на этой кривой. Разобьем кривую точками на конечное число частичных дуг. И рассмотрим сумму произведений значений функции в каждой точке на длину соответствующей частичной дуги.
;
Определение. Если при стремлении к нулю шага разбиения кривой АВ интегральные суммы имеют конечный предел, то этот предел называется криволинейным интегралом по переменной х от функции P(x, y, z) по кривой АВ в направлении от А к В.
Криволинейный интеграл второго рода, т.е. интеграл по координатам отличается от криволинейного интеграла первого рода, т.е. по длине дуги тем, что значение функции при составлении интегральной суммы умножается не на длину частичной дуги, а на ее проекцию на соответствующюю ось. (В рассмотренном выше случае на ось ОХ).
Вообще говоря, криволинейные интегралы могут считаться также и по переменным у и z.
Сумму криволинейных интегралов также называют криволинейным интегралом второго рода.
Свойства криволинейного интеграла второго рода.
1) Криволинейный интеграл при перемене направления кривой меняет знак.
2)
3)
4)
5) Криволинейный интеграл по замкнутой кривой L не зависит от выбора начальной точки, а зависит только от направления обхода кривой.
Направление обхода контура L задается дополнительно. Если L замкнутая кривая без точек самопересечения, то направление обхода контура против часовой стрелки называется положительным.
6) Если АВ кривая, лежащая в плоскости, перпендикулярной оси ОХ, то
Аналогичные соотношения справедливы при интегрировании по переменным у и z.
Теорема. Если кривая АВ кусочно- гладкая, а функции P(x, y, z), Q(x, y, z) и
R(x, y, z) непрерывны на кривой АВ, то криволинейные интегралы
существуют.
Вычисление криволинейных интегралов второго рода производится путем преобразования их к определенным интегралам по формулам:
В случае, если АВ плоская кривая, заданная уравнением y = f(x), то
Пример. Вычислить криволинейный интеграл . L контур, ограниченный параболами . Направление обхода контура положительное.
Представим замкнутый контур L как сумму двух дуг L1 = x2 и
Формула Остроградского Грина.
(Остроградский Михаил Васильевич (1861-1862) русский математик,
академик Петерб. А.Н.)
(Джордж Грин (1793 1841) английский математик)
Иногда эту формулу называют формулой Грина, однако, Дж. Грин предложил в 1828 году только частный случай формулы.
Формула Остроградского Грина устанавливает связь между криволинейным интегралом и двойным интегралом, т.е. дает выражение интеграла по замкнутому контуру через двойной интеграл по области, ограниченной этим контуром.
Будем считать, что рассматриваемая область односвязная, т.е. в ней нет исключенных участков.
y
y = y2(x)
D
A
C
B
y= y1(x)
0 x1 x2 x
Если замкнутый контур имеет вид, показанный на рисунке, то криволинейный интеграл по контуру L можно записать в виде:
Если участки АВ и CD контура принять за произвольные кривые, то, проведя аналогичные преобразования, получим формулу для контура произвольной формы:
Эта формула называется формулой Остроградского Грина.
Формула Остроградского Грина справедлива и в случае многосвязной области, т.е. области, внутри которой есть исключенные участки. В этом случае правая часть формулы будет представлять собой сумму интегралов по внешнему контуру области и интегралов по контурам всех исключенных участков, причем каждый из этих контуров интегрируется в таком направлении, чтобы область все время оставалась по левую сторону линии обхода.
Пример. Решим пример, рассмотренный выше, воспользовавшись формулой Остроградского Грина.
Формула Остроградского Грина позволяет значительно упростить вычисление криволинейного интеграла.
Криволинейный интеграл не зависит от формы пути, если он вдоль всех путей, соединяющих начальную и конечную точку, имеет одну и ту же величину.
Условием независимости криволинейного интеграла от формы пути равносильно равенству нулю этого интеграла по любому замкнутому контуру, содержащему начальную и конечную точки.
Это условие будет выполняться, если подынтегральное выражение является полным дифференциалом некоторой функции, т.е. выполняется условие тотальности.
Поверхностные интегралы первого рода.
z
Si
0
y
x
Поверхностный интеграл является таким же обобщением двойного интеграла, каким криволинейный интеграл является по отношению к определенному интегралу.
Рассмотрим поверхность в пространстве, которая произвольно разбита на n частей.
Рассмотрим произведение значения некоторой функции F в произвольной точке с координатами (, , ) на площадь частичного участка Si, содержащего эту точку.
Определение. Если при стремлении к нулю шага разбиения поверхности существует конечный предел интегральных сумм, то этот предел называется поверхностным интегралом первого рода или интегралом по площади поверхности.
Свойства поверхностного интеграла первого рода.
Поверхностные интегралы первого рода обладают следующими свойствами:
1) S площадь поверхности.
2)
3)
4) Если поверхность разделена на части S1 и S2, то
5) Если , то
6)
7) Теорема о среднем.
Если функция F(x, y, z) непрерывна в любой точке поверхности S, то существует точка (, , ) такая, что
S площадь поверхности.
Проведя рассуждения, аналогичные тем, которые использовались при нахождении криволинейного интеграла, получим формулу для вычисления поверхностного интеграла первого рода через двойной интеграл по по площади проекции поверхности на плоскость XOY.
Поверхностные интегралы второго рода.
Если на поверхности S есть хотя бы одна точка и хотя бы один не пересекающий границу поверхности контур, при обходе по которому направление нормали в точке меняется на противоположное, то такая поверхность называется односторонней.
Если при этих условиях направление нормали не меняется, то поверхность называется двухсторонней.
Будем считать положительным направлением обхода контура L, принадлежащего поверхности, такое направление, при движении по которому по выбранной стороне поверхности сама поверхность остается слева.
Двухсторонняя поверхность с установленным положительным направлением обхода называется ориентированной поверхностью.
Рассмотрим в пространстве XYZ ограниченную двухстороннюю поверхность S, состоящую из конечного числа кусков, каждый из которых задан либо уравнением вида z = f(x, y), либо является цилиндрической поверхностью с образующими, параллельными оси OZ.
Определение. Если при стремлении к нулю шага разбиения поверхности S интегральные суммы, составленные как суммы произведений значений некоторой функции на площадь частичной поверхности, имеют конечный предел, то этот предел называется поверхностным интегралом второго рода.
- поверхностный интеграл второго рода.
Свойства поверхностного интеграла второго рода аналогичны уже рассмотренным нами свойствам поверхностного интеграла первого рода.
Т.е. любой поверхностный интеграл второго рода меняет знак при перемене стороны поверхности, постоянный множитель можно выносить за знак интеграла, поверхностный интеграл от суммы двух и более функций равен сумме поверхностных интегралов от этих функций, если поверхность разбита на конечное число частичных поверхностей, интеграл по всей поверхности равен сумме интегралов по частичным поверхностям.
Если S- цилиндрическая поверхность с образующими, параллельными оси OZ, то . В случае, если образующие поверхности параллельны осям OX и OY, то равны нулю соответствующие составляющие поверхностного интеграла второго рода.
Вычисление поверхностного интеграла второго рода сводится к вычислению соответствующих двойных интегралов. Рассмотрим это на примере.
Пример. Вычислить интеграл по верхней стороне полусферы
Преобразуем уравнение поверхности к виду:
Заданная поверхность проецируется на плоскость XOY в круг, уравнение которого:
Для вычисления двойного интеграла перейдем к полярным координатам:
(См. Двойной интеграл в полярных координатах.)
,
Связь поверхностных интегралов первого и второго рода.
Поверхностные интегралы первого и второго рода связаны друг с другом соотношением:
В этой формуле cos, cos, cos - направляющие косинусы нормали к поверхности S в выбранную сторону поверхности.
Формула Гаусса Остроградского.
Формула Гаусса Остроградского является аналогом формулы Грина Остроградского. Эта формула связывает поверхностный интеграл второго рода по замкнутой поверхности с тройным интегралом по пространственной области, ограниченной этой поверхностью.
Для вывода формулы Гаусса Остроградского надо воспользоваться рассуждениями, подобными тем, которые использовались при нахождении формулы Грина Остроградского.
Рассматривается сначала поверхность, ограниченная сверху и снизу некоторыми поверхностями, заданными известными уравнениями, а сбоку ограниченную цилиндрической поверхностью. Затем рассматривается вариант когда поверхность ограничена цилиндрической поверхностью с образующими, параллельными дум доугим координатным осям.
После этого полученные результаты обобщаются, приводя к формуле Гаусса Остроградского:
Отметим, что эта формула применима для вычисления поверхностных интегралов по замкнутой поверхности.
На практике формулу Гаусса Остроградского можно применять для вычисления объема тел, если известна поверхность, ограничивающая это тело.
Тиеют место формулы:
Пример. Найти формулу вычисления объема шара.
В поперечных сечениях шара (сечения параллельны плоскости XOY) получаются окружности.
Уравнение шара имеет вид:
Найти объем шара можно по формуле:
Для решения этой же задачи можно воспользоваться преобразованием интеграла к сферическим координатам. (См. Сферическая система координат.) Это значительно упростит интегрирование.
Элементы теории поля.
Определение. Если каждой точке пространства М ставится в соответствие некоторая скалярная величина U, то таким образом задается скалярное поле U(M). Если каждой пространства М ставится в соотвтствие вектор , то задается векторное поле (М).
Пусть в пространстве М задана поверхность . Будем считать, что в каждой точке Р определяется положительное направление нормали единичным вектором .
В пространстве М зададим векторное поле, постовив в соответствие каждой точке точке пространства вектор, определенный координатами:
Если разбить каким либо образом поверхность на частичные участки i и составить сумму , где - скалярное произведение, то предел этой суммы при стремлении к нулю площадей частичных участков разбиения (если этот предел существует) будет поверхностным интегралом.
Определение. Поверхностный интеграл называется потоком векторного поля через поверхность .
Если поверхность разбита на конечное число частичных поверхностей, то поток векторного поля через всю поверхность будет равен сумме потоков через частичные поверхности.
Если преобразовать скалярное произведение в координатную форму, то получаем соотношение:
Если на области существует функция f(x, y, z), имеющая непрерывные частные производные, для которых выполняются свойства:
то такую функцию называют потенциальной функцией или потенциалом вектора .
Тогда вектор является градиентом функции f.
Потенциал может быть найден по формуле:
В этой формуле x0, y0, z0 координаты некоторой начальной точки. В качестве такой точки удобно брать начало координат.
Теорема. Для того, чтобы поле вектора , заданного в некоторой области, имело потенциал, необходимо и достаточно, чтобы выполнялось одно из двух условий:
1) Интеграл от вектора по любому кусочно гладкому контуру, принадлежащему области, равен нулю.
2) Интеграл по любому кусочно гладкому пути, соединяющему две любые точки поля не зависит, от пути интегрирования.
Формула Стокса.
(Джордж Габриель Стокс (1819 1903) английский математик)
Формула Стокса связывает криволинейные интегралы второго рода с поверхностными интегралами второго рода.
Пусть в пространстве задана некоторая поверхность S. L непрерывный кусочно гладкий контур поверхности S.
z S
L
y
l
x
Предположим, что функции P,Q и R непрерывны на поверхности S вместе со своими частными производными первого порядка. Применим формулу, выражающую криволинейный интеграл через определенный.
Введем обозначения:
Применив формулу Грина Остроградского, можно заменить криволинейный интеграл равным ему двойным интегралом. После преобразований устанавливается следуюшее соответствие между криволинейным и поверхностным интегралом:
эта формула и называется формула Стокса.
Определение. Вектор , компоненты которого равны соответственно равны
называется вихрем или ротором вектора и обозначается:
Определение. Символический вектор называется оператором Гамильтона. ( Уильям Роуан Гамильтон (1805 1865) ирландский математик) Символ - “набла”.
С учетом этого обозначения можно представить себе понятие ротора вектора как векторного произведения оператора Гамильтона на вектор .
Определение. Криволинейный интеграл, представляющий собой работу векторного поля вдоль некоторой кривой L называется линейным интегралом от вектора по ориентированной кривой L.
Если кривая L представляет собой замкнутый контур, то линейный интеграл по такому контуру называется циркуляцией вектроного поля вдоль контура L.
В векторной форме теорему Стокса можно сформулировать так:
Циркуляция вектора вдоль контура некоторой поверхности равна потоку вихря (ротора) через эту поверхность.
Отметим, что рассмотренная выше формула Грина Остроградского является частным случаем формулы Стокса.
Также при условии равенства нулю всех компонент ротора вектора, получаем, что криволинейный интеграл по любой пространственной кривой равен нулю, т.е. криволинейный интеграл не зависит от пути интегрирования.
Определение. Выражение называется дивергенцией вектора (дивергенцией векторной функции) и обозначается
Таким образом, формулу Гаусса Остроградского может быть записана в виде:
или
т.е. интеграл от дивергенции векторного поля по объему равен потоку вектора через поверхность, ограниченную этим объемом.
Определение. Векторное поле называется соленоидальным (трубчатым), если div =0 .
C помощью описанного выше оператора Гамильтона можно представить определенные нами понятия следующим образом:
Как было сказано выше (См. Уравнение Лапласа.), выражение
называется оператором Лапласа.
Справедливы следующие соотношения:
Справедливость этих равенств легко проверить непосредственной подстановкой.
Теперь рассмотрим примеры применения рассмотренных выше понятий.
Пример. Найти , если
Найдем скалярное произведение:
Найдем скалярное произведение:
Пример. Найти поток векторного поля через сторону треугольника S, вырезанного из плоскости координатными плоскостями.
z
x = 1 z z = 1 - y
x
y = 1 - x
y
Пример. Найти div(grad u), если
Пример. Определить является ли векторное поле
и найти его потенциал.
Если поле потенциально, то должны выполняться следующие условия:
Эти условия эквивалентны условию равенства нулю ротора векторного поля.справедливость этого утверждения видна из формулы ротора.
Таким образом, поле потенциальное. Потенциал находится по формуле:
Содержание КВМ Часть 1.
Содержание КВМ Часть 2.
Содержание КВМ Часть 4.
Содержание:
Обыкновенные дифференциальные уравнения.
Основные определения.
Свойства общего решения.
Теорема Коши.
Интегральные кривые.
Особое решение.
Дифференциальные уравнения первого порядка.
Уравнения вида у = f(х).
Уравнения с разделяющимися переменными.
Однородные уравнения.
Уравнения, приводящиеся к однородным.
Линейные уравнения.
Линейные однородные дифференциальные уравнения.
Линейные неоднородные дифференциальные уравнения.
Метод Бернулли.
Метод Лагранжа.
Уравнение Бернулли.
Уравнения в полных дифференциалах.
Условие тотальности.
Уравнения вида у = f(y) и x = f(y).
Уравнения Лагранжа и Клеро.
Геометрическая интерпретация решений дифференциальногоуравнения первого порядка.
Поле направлений.
Изоклины.
Численные методы решения дифференциальных уравнений.
Метод Эйлера.
Ломаная Эйлера.
Уточненный метод Эйлера.
Метод Рунге Кутта.
Дифференциальные уравнения высших порядков.
Уравнения, допускающие понижение порядка.
Уравнения вида y(n) = f(x).
Уравнения, не содержащие явно искомой функции и ее производных до порядка n-1 включительно.
Уравнения, не содержащие явно независимой переменной.
Линейные дифференциальные уравнения высших порядков.
Линейные однородные дифференциальные уравненияс произвольными коэффициентами.
Структура общего решения.
Фундаментальна система решений.
Определитель Вронского.
Общее решение линейного однородного дифференциальногоуравнения второго порядка.
Линейные однородные дифференциальные уравнения спостоянными коэффициентами.
Характеристический многочлен и характеристическое уравнение.
Линейные неоднородные дифференциальные уравнения с произвольными коэффициентами.
Метод вариации произвольных постоянных.
Линейные неоднородные дифференциальные уравнения с постоянными коэффициентами.
Уравнения с правой частью специального вида.
Нормальные системы обыкновенных дифференциальныхуравнений.
Нормальные системы линейных однородных дифференциальных уравнений с постоянными коэффициентами.
Элементы теории устойчивости.
Устойчивость по Ляпунову.
Точка покоя.
Теорема Ляпунова.
Классификация точек покоя.
Уравнения математической физики.
Уравнения в частных производных.
Линейные однородные дифференциальные уравнения в частных производных первого порядка.
Классификация основных типов уравнений математической физики.
Уравнение колебаний струны.
Граничные, начальные и краевые условия.
Решение задачи Коши методом разделения переменных. (Метод Фурье).
Решение задачи Коши методом Даламбера.
Уравнение теплопроводности.
Уравнение Лапласа.
Задача Дирихле.
Решение задачи Дирихле для круга.
Ряды.
Основные определения.
Свойства рядов.
Критерий Коши.
Ряды с неотрицательными членами.
Признак сравнения.
Признак Даламбера.
Предельный признак Даламбера.
Признак Коши.
Интегральный признак Коши.
Знакопеременные ряды.
Знакочередующиеся ряды.
Признак Лейбница.
Абсолютная и условная сходимость рядов.
Признак Даламбера и Коши для знакопеременных рядов.
Свойства абсолютно сходящихся рядов.
Функциональные последовательности.
Область сходимости.
Функциональные ряды.
Критерий Коши равномерной сходимости.
Признак Вейерштрасса.
Свойства равномерно сходящихся рядов.
Степенные ряды.
Теоремы Абеля.
Радиус сходимости.
Действия со степенными рядами.
Разложение функций в степенные ряды.
Решение дифференциальных уравнений с помощью степенных рядов.
Ряды Фурье.
Тригонометрический ряд.
Коэффициенты Фурье.
Достаточные признаки разложимости в ряд Фурье.
Разложение в ряд Фурье непериодической функции.
Ряд Фурье для четных и нечетных функций.
Ряд Фурье для функций любого периода.
Ряд Фурье по ортогональной системе функций.
Интеграл Фурье.
Преобразование Фурье.
Элементы теории функций комплексной переменной.
Свойства функций комплексной переменной.
Основные трансцендентные функции.
Производная функций комплексной переменной.
Условия Коши Римана.
Интегрирование функций комплексного переменного.
Теорема Коши.
Интегральная формула Коши.
Ряды Тейлора и Лорана.
Изолированные особые точки.
Теорема о вычетах.
Вычисление интегралов с помощью вычетов.
Операционное исчисление.
Преобразование Лапласа.
Свойства изображений.
Таблица изображений некотрых функций.
Теорема свертки и запаздывания.
Интеграл Дюамеля.
Решение дифференциальных уравнений с помощью операционногоисчисления.
Криволинейные интегралы.
Криволинейные интегралы первого рода.
Свойства криволинейных интегралов первого рода.
Криволинейные интегралы второго рода.
Свойства криволинейных интегралов второго рода.
Формула Остроградского Грина.
Поверхностные интегралы первого рода.
Свойства поверхностных интегралов первого рода.
Поверхностные интегралы второго рода.
Связь поверхностных интегралов первого и второго рода.
Формула Гаусса Остроградского.
Элементы теории поля.
Поток векторного поля.
Потенциал.
Формула Стокса.
Ротор.
Оператор Гамильтона.
Циркуляция.
Дивиргенция.
Соленоидальное поле.
Ларин Александр Александрович 2000 год. mailto: aalar@yandex.ru
EMBED Package
EMBED Package
EMBED Package
EMBED Package
EMBED Package
EMBED Package
EMBED Package