Методы решения дифференциальных уравнений обзор

Методы решения дифференциальных уравнений

Здесь мы рассмотрим методы решения обыкновенных дифференциальных уравнений. Это уравнения, зависящие от одной независимой переменной, зависимой переменной, и ее производных:
.
Основные определения, относящиеся к дифференциальным уравнениям, изложены на странице Основные понятия и определения дифференциальных уравнений.

Мы считаем, что уравнения имеют решения в области задания переменных; Функции, заданные неявно можно разрешить относительно одной из переменной. Мы не проводим исследования этих и подобных вопросов. Здесь мы рассматриваем только методы решения.

Мы часто будем делить, и умножать уравнения на какие-то функции. В таких операциях нужно соблюдать осторожность. От этого могут появляться дополнительные решения, или исчезать имеющиеся. Например, если мы умножим все части уравнения на , то может появиться новое решение . Если мы разделим все части уравнения на , то может исчезнуть решение , если оно имелось в исходном уравнении. То есть, если мы умножаем или делим уравнение на некоторую функцию f , то всегда нужно особо рассматривать случай f = 0 . Здесь мы не будем заострять на этом внимание.

Дифференциальные уравнения первого порядка

Дифференциальное уравнение и его интеграл

Далее, если это особо не оговорено, мы считаем, что x – это независимая переменная, а y – зависимая. То есть y есть функция от x : . Однако, в уравнениях первого порядка, мы можем легко менять роли переменных. То есть можно считать y независимой переменной, а x – зависимой. Но по умолчанию, x – это независимая переменная, а y – зависимая.

Пусть у нас есть дифференциальное уравнение первого порядка, разрешенное относительно производной. Запишем его в следующем виде:
(1) .
Здесь p и q – заданные функции двух переменных.

Далее рассмотрим уравнение:
(2) ,
где φ – некоторая функция двух переменных; C – постоянная, то есть число. Положим, что y есть функция от x : . Тогда будет уже сложной функцией от одной переменной x . Обозначим ее буквой : .
Перепишем уравнение (2), выразив левую часть через переменную x :
(3) .
Дифференцируем это уравнение по x , применяя правило дифференцирования сложной функции:
;
.
Мы получили дифференциальное уравнение первого порядка, имеющее тот же вид, что и уравнение (1). Отсюда следует, что если , то функция , определяемая из уравнения , является решением исходного уравнения (1).

Заметим, что левая часть уравнения является производной от функции :
.
Тогда сама функция является интегралом по отношению к уравнению (1), точнее – к его левой части. По этой причине решение уравнения, записанного в виде , называется интегралом уравнения, а сам процесс решения называется интегрированием дифференциального уравнения.

Уравнения в дифференциалах

Воспользуемся свойством дифференциалов, согласно которому
.
Перепишем уравнение (1) и умножим его на dx :
(1) ;
(4) .
Мы получили уравнение, связывающее дифференциалы переменных x и y . По этой причине такие уравнения называются дифференциальными уравнениями. Такая форма записи называется уравнением в дифференциалах, или дифференциальной формой уравнения. Уравнения (1) и (4) эквивалентны. Можно использовать любую из этих форм.

Пусть
(5) ,
где – некоторая функция двух переменных. Подставим в (4):
(6) .
Отсюда видно, что левая часть уравнения (6) является дифференциалом функции : . Тогда уравнение (6) можно переписать в виде равенства нулю дифференциала:
.
Отсюда следует, что функция равняется постоянной, которую обозначим буквой C . Тогда общий интеграл уравнения (4), при условии (5), имеет вид:
(7) .

Уравнения в полных дифференциалах

Итак, мы нашли, что если в уравнении
(4) ,
функции p и q являются частными производными
(5)
от некоторой функции φ , то уравнение (4) имеет интеграл
(7) .
Такие уравнения называются дифференциальными уравнениями в полных дифференциалах.

Как правило интеграл уравнения (7) нам не известен, а известно лишь само уравнение, то есть известны функции и . Возникает вопрос, как по известным функциям p и q определить, что левая часть уравнения является полным дифференциалом? Оказывается, что сделать это достаточно просто. Для того, чтобы уравнение было в полных дифференциалах, необходимо и достаточно, чтобы выполнялось условие:
(8) .
Доказательство

Зная, что уравнение относится к классу уравнения в полных дифференциалах, мы можем найти функцию , применяя несколько методов. Рассмотрим метод последовательного выделения дифференциала. В этом методе, мы применяем формулы дифференцирования, записанные в дифференциальной форме:
;
;
;
.
Здесь и могут быть любыми функциями от и . Рассмотрим применение этого метода на конкретном примере.

Пример

Дано уравнение:
(П1) .
Требуется проверить, является ли это уравнение в полных дифференциалах. И если является, то решить его.

В нашем случае . Проверим, является ли это уравнение в полных дифференциалах. Находим частные производные.
;
.
Видно, что . То есть это уравнение в полных дифференциалах. Решаем его, последовательно выделяя дифференциал.

.
Итак, мы нашли эквивалентное (П1) уравнение
.
Отсюда получаем его общий интеграл:
.

Решать подобные уравнения можно также и методом последовательного интегрирования. Решение этим методом можно найти на странице Дифференциальные уравнения в полных дифференциалах

Интегрирующий множитель

Итак, мы научились решать дифференциальные уравнения первого порядка
(4)
при условии
(8) .
Но если существует единственное решение уравнения (4), и условие (8) не выполняется, то оказывается, что существует такая функция , умножив на которую уравнение (4), оно становится уравнением в полных дифференциалах. Причем существует бесконечное множество таких функций. Доказательство

В качестве примера рассмотрим уравнение:
(П2) .
Перепишем его, сгруппировав члены:
.
Заметим, что . Поэтому разделим уравнение на , чтобы выделить полный дифференциал . При имеем:
.
Выделяем полный дифференциал:
;
.
Отсюда получаем общий интеграл исходного уравнения:
.
В уравнении (П2), интегрирующий множитель равен . Когда мы умножили на него уравнение, то оно стало уравнением в полных дифференциалах, которое мы и решили.

Заметим, что при умножении уравнения на множитель , мы получили другое уравнение. Оно эквивалентно исходному за исключением точек, в которых и . Уравнение корней не имеет. Поэтому этот случай отпадает. А уравнение имеет корень :
.
Поэтому умножение уравнения на множитель дает эквивалентное уравнение, за исключением точек . Другими словами, поскольку мы разделили уравнение на , то нужно проверить случай . Подстановкой в (П2) убеждаемся, что также является решением исходного уравнения. Поэтому общее решение имеет вид:
; .

В этом примере мы угадали, что если уравнение умножить на , то можно выделить полный дифференциал. Не смотря на то, что для любого уравнения, при условии существования его решения, интегрирующий множитель существует, у нас нет общего метода, который позволяет найти его для любого дифференциального уравнения. Можно попытаться это сделать, но для произвольного уравнения нет гарантии, что мы найдем интегрирующий множитель, и решим уравнение. К счастью есть несколько классов уравнений, для которых это сделать можно. Эти типы уравнений мы и рассмотрим.

Уравнения с разделяющимися переменными

Рассмотрим уравнение
(9) ,
где – некоторые заданные функции. Перепишем это уравнение в дифференциалах:
.
Разделим его на . При имеем:
(11) .
Уравнение имеет вид суммы, каждое слагаемое которой зависит только от одной переменной. Говорят, что переменные разделились, а уравнение (9), по этой причине, называют дифференциальным уравнением с разделяющимися переменными.

Нетрудно видеть, что уравнение (11) в полных дифференциалах. Действительно, поскольку множитель при dx не зависит от y , то . Поскольку множитель при dy не зависит от x , то .
Видно, что необходимое и достаточное условие для полных дифференциалов выполняется:
.

Таким образом мы нашли интегрирующий множитель: . Это позволяет нам выделить дифференциал и получить решение в квадратурах:
;
;
.
Отсюда получаем общий интеграл:
.

Пример

Решить уравнение:
(П3) .

Перепишем (П3) в дифференциалах:
.
Разделим на . При имеем:
.
Переменные разделились. Общий интеграл имеет вид:
.
Далее, см. Дифференциальные уравнения с разделяющимися переменными

Уравнения, приводящиеся к разделяющимся переменным

К уравнению с разделяющимися переменными приводятся уравнения вида
,
где f – функция; a, b, c – постоянные. Для решения подобного уравнения нужно от переменной y перейти к новой переменной u , сделав подстановку .

Пример

Решить уравнение:
(П4.1)

От переменной y перейдем к переменной u . Делаем подстановку:
(П4.2) .
Здесь и – функции от x . Дифференцируем (П4.2) по x , и подставляем (П4.1):
;
.
Тем самым мы получили дифференциальное уравнение с разделяющимися переменными.
См. далее Дифференциальные уравнения, приводящиеся к уравнениям с разделяющимися переменными

Однородные уравнения

Однородные дифференциальные уравнения первого порядка имеют вид
.
Чтобы определить, является ли уравнение однородным, нужно сделать замену . Здесь t – постоянная. Если t сократится, то это однородное уравнение. Для его решения нужно от переменной y перейти к переменной u , сделав подстановку . После этого, уравнение сводится к уравнению с разделяющимися переменными.

Пример

Проверим, является ли это уравнение однородным. Сделаем замену . Считаем, что постоянная :
;
;
.
Постоянная t сократилась. Она также сократится, если считать . Это однородное уравнение. Переходим от переменной y к переменной u . Для этого делаем подстановку , где u – функция от x . Дифференцируем по x :
.
Подставляем в(П5):
;
;
.
При , берем знак ′+′ . При – знак ′–′ . Мы получили уравнение с разделяющимися переменными, решать которое мы уже умеем.
Далее см. Однородные дифференциальные уравнения первого порядка

Уравнения, приводящиеся к однородным

Уравнение вида
.
приводится к однородному подстановками
,
где – новые переменные; – постоянные, которые выбираются из условий
.

Пример

От переменных x и y , переходим к переменным t и u . Делаем подстановку . Тогда ;
;
;
.
Решаем систему из двух линейных уравнений

Определив и , получаем однородное уравнение:
.
Метод решения такого уравнения мы только что рассмотрели. См. Дифференциальные уравнения первого порядка, приводящиеся к однородным

Обобщенные однородные уравнения

К однородным уравнениям приводятся уравнения вида
.
Чтобы определить, является ли дифференциальное уравнение обобщенным однородным, нужно ввести постоянную t и сделать замену: . Если удастся выбрать такое значение α , при котором постоянная t сократится, то это – обобщенное однородное дифференциальное уравнение. Для решения этого уравнения, нужно от переменной y перейти к переменной u , сделав подстановку . При этом уравнение сводится к разделяющимся переменным.

Пример

Проверим, является ли уравнение (П7) обобщенным однородным. Делаем замену: .
.
Подставляем в (П7):
.
Делим на :
.
Отсюда видно, что t сокращается, если положить .

Итак, мы нашли, что это обобщенное однородное уравнение с . Решаем его. От переменной y переходим к переменной u , выполняя подстановку .
;
;
.
Подставляем в (П7):
(П7) ;
;
;
;
.
Мы получили уравнение с разделяющимися переменными.
См. далее Обобщенные однородные дифференциальные уравнения первого порядка

Линейные уравнения

Дифференциальные уравнения, вида
(11)
называются линейными дифференциальными уравнениями первого порядка.

Решение с помощью интегрирующего множителя

Уравнение (11) имеет интегрирующий множитель .
См. Решение линейных дифференциальных уравнений первого порядка
Продемонстрируем это на примере.

Пример

Это линейное уравнение первого порядка. Решаем его с помощью интегрирующего множителя. Разделим (П8.1) на x :
(П8.2) .
Тогда . Находим интегрирующий множитель :
; .
Пусть . Тогда . Умножаем (П8.2) на и выделяем полный дифференциал:
;
;
;
.
Отсюда , или .

Мы нашли интегрирующий множитель полагая, что . После умножения на него, мы получили уравнение в полных дифференциалах как при , так и при . При решении мы нигде не полагали, что . Это предположение нам потребовалось, только чтобы выбрать интегрирующий множитель. На самом деле, любое уравнение первого порядка имеет бесконечное число интегрирующих множителей. Поэтому, если бы мы в самом начале взяли , то получили бы множитель . И с его помощью, получили то же самое решение.

Решение методом Бернулли

Линейное уравнение первого порядка можно решить красивым приемом, введя две функции и , зависящие от переменной x . Сделаем подстановку . Тогда . Подставим в исходное уравнение (11):
(11) ;
;
(12) .
Наложим условие
(13) .
Уравнение (13) с разделяющимися переменными. Решаем его, и возьмем любое, отличное от нуля частное решение. Так мы определим функцию . Учитывая (13), уравнение (12) примет вид:
.
Теперь здесь уже известная функция, и это уравнение с разделяющимися переменными. Решая его, найдем общее решение . Вместе с этим получаем общее решение исходного уравнения (11): .
Подробнее, см. Решение линейного ДУ первого порядка методом Бернулли

Решение методом Лагранжа

Метод Лагранжа интересен тем, что указывает путь поиска решения от простого к сложному. Рассмотрим линейное уравнение:
(11) .
Давайте его упростим. Сначала рассмотрим однородное уравнение – то есть уравнение с :
(14) .
Это уравнение с разделяющимися переменными, и мы можем его решить:
;
;
;
;
.
Заменим постоянную на C . Тогда общее решение примет вид:
(15) , где .

Теперь вернемся к исходному неоднородному уравнению (11). Попытаемся найти его решение, используя решение более простого, однородного уравнения (14). Для этого в (15) заменим постоянную C на функцию, зависящую от переменной x : . То есть будем искать решение в виде
.
Подставляя в (11), получим для дифференциальное уравнение с разделяющимися переменными, которое решается в квадратурах. Решив его, получаем решение исходного уравнения. Такой метод решения называется методом вариации постоянных, или методом Лагранжа.
См. Решение линейных ДУ первого порядка методом Лагранжа

Дифференциальное уравнение Бернулли

Дифференциальное уравнение Бернулли имеет вид:
,
где и – заданные функции от x . Можно убедиться, что оно сводится к линейному уравнению подстановкой .
См. Дифференциальное уравнение Бернулли и методы его решения

Однако его легче решать методом двух функций Бернулли. Для этого вводим две функции и . Ищем решение в виде . Одну из этих функций выбираем так, чтобы уравнение для другой функции превратилось в уравнение с разделяющимися переменными.

Пример

Это уравнение Бернулли. Решаем методом Бернулли. Ищем решение в виде произведения двух функций: . Тогда
. Подставляем в (П9.1):
;
(П9.2) .
Одну из этих функций мы можем выбрать произвольным образом. Выберем v так, чтобы выражение в круглых скобках равнялось нулю:
(П9.3) .
Тогда уравнение (П9.2) превратится в уравнение с разделяющимися переменными.

Решаем уравнение (П9.3). Разделяем переменные.
;
;
;
;
.
Возьмем решение , или .

Подставим в (П9.2), учитывая (П9.3), и разделяем переменные:
(П9.2) ;
;
;
.
При имеем:
;
;
;
;
;
.
Заменим постоянную интегрирования: . Тогда решение уравнения (П9.1) примет вид:
.

Теперь рассмотрим случай . Нетрудно увидеть, что это также решение уравнения (П9.2). Тогда является решением исходного уравнения. Получаем общее решение исходного уравнения:
.

Уравнения, не разрешенные относительно производной

Существует несколько типов уравнений, не разрешенных относительно производной, которые допускают решение. При этом они должны быть разрешены относительно одной из переменной. Далее перечислены типы этих уравнений, и даны ссылки на страницы с методами их решений.

Дифференциальные уравнения второго и высших порядков

Уравнения, допускающие понижение порядка

Уравнения, не содержащие y в явном виде

Рассмотрим уравнения вида
.
Если сделать подстановку , то . То есть мы понизили на единицу порядок такого уравнения.
См. Дифференциальные уравнения, не содержащие функцию в явном виде

Уравнения, не содержащие x в явном виде

Рассмотрим уравнения, которые не содержат независимую переменную x в явном виде:
.
Мы можем понизить порядок таких уравнений, если от переменных x и y перейдем к независимой переменной y и зависимой переменной y′ . То есть, считаем, что все производные являются функциями от y .

Пример

Это уравнение не содержит независимую переменную x в явном виде. Переходим к новым переменным. Пусть независимой переменной является y , а зависимой y′ . Введем для нее обозначение:
. Тогда по правилу дифференцирования сложной функции, имеем:
.

Подставляем в (П10.1):
.
Мы получили дифференциальное уравнение первого порядка с разделяющимися переменными.
Далее, см. Дифференциальные уравнения высших порядков, не содержащие переменную в явном виде

Уравнения, однородные относительно функции и ее производных

Это уравнения вида
.
Чтобы распознать такое уравнение, нужно сделать замены , и т.д. Если постоянная t сократится, то это уравнение однородное относительно функции и ее производных.

Для решения, мы от зависимой переменной y переходим к новой зависимой переменной u с помощью подстановки
,
где – функция от x .

Пример

Проверим, является ли это уравнение однородным относительно функции и ее производных. Заменим в исходном уравнении y на ty , y′ на ty′ , y′′ на ty′′ :
;
.
Постоянная t сокращается. Значит это уравнение однородное относительно функции и ее производных.

Делаем подстановку , где – функция от x .
.
(П11.1) ;
.
Делим на . При имеем:
;
;
.
Мы получили линейное дифференциальное уравнение первого порядка. См. далее ДУ высших порядков, однородные относительно функции и ее производных

Линейные уравнения с постоянными коэффициентами

В линейных уравнениях с постоянными коэффициентами, допускающими решение в аналитическом виде, можно сделать линейную подстановку, и понизить порядок уравнения. Однако проще воспользоваться свойствами линейных уравнений и решать их более простым методом.

Общие свойства линейных уравнений

Рассмотрим линейное однородное дифференциальное уравнение n-го порядка с постоянными коэффициентами:
(Л1) ,
где – постоянные, то есть не зависящие от переменной x коэффициенты (числа). При этом . Это уравнение имеет n линейно независимых решений:
(Л2) .
Они называются фундаментальной системой решений. Когда n линейно независимых решений найдены, то общее решение однородного уравнения (Л1) имеет вид:
.

Теперь рассмотрим более общее – линейное неоднородное дифференциальное уравнение n-го порядка с постоянными коэффициентами:
(Л3) ,
где – непрерывная функция на некотором отрезке . Тогда, на этом отрезке, уравнение (Л3) имеет решение, удовлетворяющее заданным начальным условиям:
, где – любые действительные числа; .

Пусть есть частное (любое) решение уравнения (Л3). Тогда общее решение неоднородного уравнения (Л3) равно сумме частного решения неоднородного уравнения, и общего решения однородного:
,
где – общее решение однородного уравнения (Л1).

Если, в уравнении (Л3), неоднородную часть можно представить в виде суммы p слагаемых:
,
то частное решение равно сумме отдельных частных решений: . Здесь – частное решение уравнения
.

Решение однородного уравнения

Рассмотрим линейное однородное ДУ n-го порядка с постоянными коэффициентами:
(Л1) .
Чтобы найти его общее решение, нам нужно найти n линейно независимых решений. Или, как говорят, найти фундаментальную систему решений. Ищем решение в виде . Подставляя в (Л1), получаем уравнение степени n, которое называют характеристическим уравнением:
(Л4) .
Оно имеет n корней , и может быть записано в виде:
.
Каждому корню соответствует частное решение, входящее в состав фундаментальной системы. При этом корни могут быть кратными и комплексными. Рассмотрим правила составления линейно независимых решений.

Действительному единственному корню соответствует решение .
Действительному корню кратности p , соответствуют p линейно независимых решений:
.
Если есть единственный комплексный корень , то имеется и комплексно сопряженный корень . Им соответствуют два линейно независимых решения
.
Если есть кратный комплексный корень кратности p , то имеется и комплексно сопряженный корень, кратности p : . Им соответствуют 2 p линейно независимых решений
;
;
;
.
.

Пример

Ищем решение в виде . Составляем характеристическое уравнение и преобразуем его:
;
(П12.2) .
Решаем квадратное уравнение :
.
Перепишем характеристическое уравнение (П12.2) в эквивалентном виде:
.
Корням кратности 2 соответствуют два линейно независимых решения:
;
.
Комплексно сопряженным корням , соответствуют решения
.
Общее решение:
.

Решение уравнений со специальной неоднородностью

Рассмотрим, часто встречающееся в приложениях, линейное неоднородное уравнение с постоянными коэффициентами со специальной неоднородностью:
(Л5.1) ,
где правая часть представлена в виде произведений степенной функции, экспоненты, косинусов и синусов:
(Л5.2) .
Здесь – многочлены степеней и , соответственно.

Общее решение (Л5.1) – (Л5.2) имеет вид:
.
Здесь – общее решение однородного уравнения (с ); – частное решение неоднородного уравнения (Л5.1)–(Л5.2). Как найти общее решение , мы рассмотрели в предыдущем пункте. Изложим метод нахождения частного решения.

Ищем методом неопределенных коэффициентов. Известно, что для уравнения (Л5.1) – (Л5.2), частное решение имеет следующий вид:
(Л6) .
Здесь ; и – многочлены степени s . Если среди корней характеристического уравнения (Л4) ⇑ нет корня , то . Если такой корень есть, то m – его кратность.

Метод нахождения частного решения заключается в том, что мы ищем решение в виде (Л6). Для этого записываем многочлены в общем виде:
;
.
Здесь коэффициенты (числа), которые нужно определить. Далее мы выписываем (Л6) в общем виде:

.
Находим n производных , и подставляем их выражения в исходное уравнение (Л5.1) – (Л5.2). В левой части мы получим сумму из членов с множителями , и членов с множителями . Здесь . В правой, неоднородной части, также имеется членов с множителями , и членов с множителями . При этом часть из них может равняться нулю. Приравнивая коэффициенты при этих множителях, получим систему из уравнений, решая которую, мы определяем неизвестные коэффициенты .

Если , то правая часть имеет более простой вид:
(Л5.3) .
Тогда частное решение содержит неопределенных коэффициентов:
.
Здесь , если характеристическое уравнение (Л4) ⇑ не имеет действительного корня . Если характеристическое уравнение имеет действительный корень , то m – его кратность.
Далее находим выражения для n производных , и подставляем их в исходное уравнение (Л4.1) – (Л4.2). В левой части мы получим сумму из членов с множителями . В правой, неоднородной части, также имеется членов с множителями . Часть из них может равняться нулю. Приравнивая коэффициенты при этих множителях, получим систему из уравнений, решая которую, мы определяем неизвестные коэффициенты .

Наконец, если и , и , то
(Л5.4) .
Частное решение, как и в предыдущем случае, имеет неопределенных коэффициентов:
.
Если характеристическое уравнение не имеет действительного корня , то . Если характеристическое уравнение имеет такой корень, то m – его кратность. Находим выражения для n производных ; подставляем их в исходное уравнение (Л4.1) – (Л4.2). В левой части мы получим сумму из членов с множителями . В правой, неоднородной части, также имеется членов с множителями . Часть из них может равняться нулю. Приравнивая коэффициенты при этих множителях, получим систему из уравнений, решая которую, определяем неизвестные коэффициенты .

Решение неоднородных уравнений общего вида

Теперь рассмотрим методы решения линейного неоднородного ДУ n-го порядка с постоянными коэффициентами с неоднородностью общего вида:
(Л3) .
В отличие от предыдущего случая со специальной неоднородностью, в этом разделе мы считаем, что неоднородность имеет произвольный вид.

Решение методом Бернулли

Метод Бернулли заключается в том, что мы ищем решение уравнения
(Л3)
в виде произведения двух функций и , зависящих от переменной x :
.
Если в качестве v взять частное решение однородного уравнения
,
то такая подстановка приводит к понижению порядка исходного уравнения (Л3).

Пример

Ищем решение в виде произведения двух функций; подставляем в уравнение (П13.1) и группируем члены:
(П13.2) ;
;
.
(П13.1) ;
;
(П13.3) ;

Решаем однородное уравнение
(П13.4) .
Ищем решение в виде . Составляем и решаем характеристическое уравнение:
;
.
Получаем два кратных корня . Общее решение уравнения (П13.4):
.
В качестве v мы можем взять любое, отличное от нуля решение. Поэтому положим
. Тогда ; .

Понижение порядка линейной подстановкой

Порядок линейного уравнения с постоянными коэффициентами можно понизить с помощью подстановки . Более подробно этот материал изложен на странице «Понижение порядка в линейных неоднородных ДУ с постоянными коэффициентами». Здесь мы рассмотрим пример применения этого метода.

Пример

Решить уравнение, применяя линейную подстановку
(П14.1)

Перепишем левую часть уравнения (П14.1), введя оператор дифференцирования :
.
Подставим в (П14.1). Исходное уравнение принимает вид
.
Сделаем подстановку . В результате для переменной u получаем уравнение первого порядка:
;
.

Итак, подстановкой
(П14.2) ,
мы получили уравнение первого порядка:
(П14.3) .

Решаем уравнение (П14.3), умножая его на интегрирующий множитель :
;
;
;
;
.

Подставляем в уравнение (П14.2) и решаем его с помощью интегрирующего множителя .
;
;
;
;
;
.

Метод вариации постоянных Лагранжа

Выпишем еще раз линейное неоднородного ДУ n-го порядка с постоянными коэффициентами:
(Л3) .
Метод вариации постоянных, который мы применили для уравнения первого поряддка, также применим и для уравнений произвольного порядка.

Для решения уравнения (Л3), мы вначале решаем однородное уравнение
.
Получаем его общее решение, которое имеет вид:
(Л7) .
Далее мы считаем, что постоянные являются функциями от x . То есть заменяем постоянные на некоторые, пока не известные, функции . Подставляем в (Л7), и ищем решение исходного уравнения (Л3) в следующем виде:
(Л8) .
Подставляем (Л8) в (Л3). При этом на функции накладываем дополнительные ограничения:
.
В результате получаем систему n линейных уравнений относительно неизвестных . Решая эту систему, получаем значения производных , как функций от x . Интегрируя, получаем выражения для самих функций . Подставляя в (Л8), получаем общее решение исходного уравнения (Л3).

Уравнение Эйлера

Уравнение
(Л9)
называется дифференциальным уравнением Эйлера. Подстановкой
(Л10)
оно приводится к уравнению с постоянными коэффициентами.

В некоторых случаях, уравнение Эйлера проще решать напрямую, не прибегая к подстановке (Л10). При решении неоднородного уравнения, к уравнению Эйлера применимы методы двух функций Бернулли, и метод вариации постоянных Лагранжа.

Рассмотрим однородное уравнение Эйлера:
(Л11) .
Его тоже проще решить без подстановки (Л10). Для этого мы ищем решение в виде . Находим производные, и подставляем в уравнение (Л11). В результате получаем характеристическое уравнение степени n . Оно имеет n корней.

Действительному корню , кратности p , соответствуют p линейно независимых решений
;
.
Если есть комплексный корень кратности p , то есть и комплексно сопряженный корень кратности p . Им соответствуют линейно независимых решений
;
;
.
.
Определив фундаментальную систему решений , получаем общее решение однородного уравнения (Л11):
.

Использованная литература:
В.В. Степанов, Курс дифференциальных уравнений, «ЛКИ», 2015.
Л.Э. Эльсгольц, Дифференциальные уравнения и вариационное исчисление, М., 1969.

Автор: Олег Одинцов . Опубликовано: 31-05-2020

Курсовая работа: Решение дифференциальных уравнений. Обзор

Нижегородский государственный технический университет

Кафедра «Общеобразовательные и общепрофессиональные дисциплины»

«Решение дифференциальных уравнений. Обзор»

Выполнила: Аверина Л.А

Группа. ТМв 151001-09

Проверила: Ловыгина М.Б

1 Обзор методов решения в Excel

1.1 Метод Рунге-Кутта четвертого порядка для решения уравнения первого порядка

1.3 Метод Эйлера

1.4 Модифицированный метод Эйлера

1.5 Практическая часть

2 Решение дифференциальных уравнений с помощью Mathcad

2.1 Метод Эйлера

2.2 Метод Эйлера с шагом h/2

2.3 Метод Рунге – Кутты

Введение

Уравнение называется обыкновенным дифференциальным n-го порядка, если F определена и непрерывна в некоторой области и, во всяком случае, зависит от . Его решением является любая функция u(x), которая этому уравнению удовлетворяет при всех x в определённом конечном или бесконечном интервале. Дифференциальное уравнение, разрешенное относительно старшей производной имеет вид

Решением этого уравнения на интервале I=[a,b] называется функция u(x).

Решить дифференциальное уравнение у / =f(x,y) численным методом — это значит для заданной последовательности аргументов х₀ , х₁ …, х_n и числа у₀ , не определяя функцию у=F(x), найти такие значения у₁ , у₂ ,…, у_n , что у_i =F(x_i )(i=1,2,…, n) и F(x₀ )=y₀ .

Таким образом, численные методы позволяют вместо нахождения функции y=F(x) (3) получить таблицу значений этой функции для заданной последовательности аргументов. Величина h=x_k -x_k _-1 называется шагом интегрирования.

Метод Эйлера относиться к численным методам, дающим решение в виде таблицы приближенных значений искомой функции у(х). Он является сравнительно грубым и применяется в основном для ориентировочных расчетов. Однако идеи, положенные в основу метода Эйлера, являются исходными для ряда других методов.

Метод Эйлера для обыкновенных дифференциальных уравнений используется для решений многих задач естествознания в качестве математической модели. Например задачи электродинамики системы взаимодействующих тел (в модели материальных точек), задачи химической кинетики, электрических цепей. Ряд важных уравнений в частных производных в случаях, допускающих разделение переменных, приводит к задачам для обыкновенных дифференциальных уравнений – это, как правило, краевые задачи (задачи о собственных колебаниях упругих балок и пластин, определение спектра собственных значений энергии частицы в сферически симметричных полях и многое другое)

1 Обзор методов решения в Excel

1.1 Метод Рунге-Кутта четвертого порядка для решения уравнения первого порядка

Идея Рунге-Кута состоит в том, чтобы использовать метод неопределённых коэффициентов. Наиболее употребительным методом Рунге-Кутта решения уравнения первого порядка y’ = F(x,y) (1) является метод четвертого порядка, в котором вычисления производятся по формуле:

yk+1 = yk +(k1 +2k2 +2k3 +k4 )/6, (2)

k₁ = Fk h = F(xk , yk )h

Рассмотрим задачу Коши для уравнений первого порядка на отрезке [a,b]:

, (4)

Разобьём промежуток [a,b] на N частей . Обозначим , где u(x) –точное решение задачи Коши, и через значения приближенного решения в точках . Существует 2 типа численных схем :

1. явные: ) (5)

2. неявные: (6)

Здесь F некоторая функция, связывающая приближения. В явных схемах приближенное значение в точке определяется через некоторое число k уже определённых приближенных значений. В неявных схемах определяется не рекурентным способом, как в явных схемах, а для его определения возникает уравнение, поскольку равенство (6) представляет из себя именно уравнение на . Явные схемы проще, однако зачастую неявные схемы предпочтительнее

1.3 Метод Эйлера

Таким образом, численные методы позволяют вместо нахождения функции У=F(x) получить таблицу значений этой функции для заданной последовательности аргументов. Величина h=x_k -x_k _-1 называется шагом интегрирования.

Рассмотрим дифференциальное уравнение первого порядка (7) с начальным условием

Требуется найти решение уравнения (7) на отрезке [а,b].

Разобьем отрезок [a, b] на n равных частей и получим последовательность х₀ , х₁ , х₂ ,…, х_n , где x_i =x₀ +ih (i=0,1,…, n), а h=(b-a)/n-шаг интегрирования.

В методе Эйлера приближенные значения у(х_i )»y_i вычисляются последовательно по формулам у_i +hf(x_i , y_i ) (i=0,1,2…).

При этом искомая интегральная кривая у=у(х), проходящая через точку М₀ (х₀ , у₀ ), заменяется ломаной М₀ М₁ М₂ … с вершинами М_i (x_i , y_i ) (i=0,1,2,…); каждое звено М_i M_i ₊₁ этой ломаной, называемой ломаной Эйлера, имеет направление, совпадающее с направлением той интегральной кривой уравнения (7), которая проходит через точку М_i . Если правая часть уравнения (7) в некотором прямоугольнике R<|x-x₀ |£a, |y-y₀ |£b>удовлетворяет условиям:

|df/dx|=|df/dx+f(df/dy)| £ M (M=const),

то имеет место следующая оценка погрешности:

где у(х_n )-значение точного решения уравнения (7) при х=х_n , а у_n — приближенное значение, полученное на n-ом шаге.

Формула (13) имеет в основном теоретическое применение. На практике иногда оказывается более удобным двойной просчет: сначала расчет ведется с шагом h, затем шаг дробят и повторный расчет ведется с шагом h/2. Погрешность более точного значения у_n * оценивается формулой

Метод Эйлера легко распространяется на системы дифференциальных уравнений и на дифференциальные уравнения высших порядков. Последние должны быть предварительно приведены к системе дифференциальных уравнений первого порядка.

1.4 Модифицированный метод Эйлера

Рассмотрим дифференциальное уравнение (7) y / =f(x,y) с начальным условием y(x₀ )=y₀ . Разобьем наш участок интегрирования на n равных частей. На малом участ интегральную кривую заменим прямой линией.

Рисунок 1 Метод Эйлера в графическом виде

Получаем точку М_к (х_к ,у_к ). Через М_к проводим касательную:

Делим отрезок (х_к ,х_к1 ) пополам

Получаем точку N_k / . В этой точке строим следующую касательную:

Из точки М_к проводим прямую с угловым коэффициентом α_к и определяем точку пересечения этой прямой с прямой Х_к1 . Получаем точку М_к / . В качестве у_к+1 принимаем ординату точки М_к / . Тогда:

(14)-рекурентные формулы метода Эйлера.

Сначала вычисляют вспомогательные значения искомой функции у_к+1/2 в точках х_к+1/2 , затем находят значение правой части уравнения (11) в средней точке y / _k _+1/2 =f(x_k _+1/2 , y_k _+1/2 ) и определяют у_к+1 .

Для оценки погрешности в точке х_к проводят вычисления у_к с шагом h, затем с шагом 2h и берут 1/3 разницы этих значений:

где у(х)-точное решение дифференциального уравнения.

Таким образом, методом Эйлера можно решать уравнения любых порядков. Например, чтобы решить уравнение второго порядка y // =f(y / ,y,x) c начальными условиями y / (x₀ )=y / ₀ , y(x₀ )=y₀ , выполняется замена

Тем самым преобразуются начальные условия

Здесь решается уравнение dy/dx = 2x-y+x 2 на интервале [0,2], начальное значение y(0)=0, для оценки точности задано также точное решение в виде функции u(x)=x 2 . Оценка погрешности делается в нормеL₁ , как и принято в данном случае

2 Решение дифференциальных уравнений с помощью Mathcad

Mathcad имеет ряд встроенных функций, предназначенных для решения обыкновенных дифференциальных уравнений (ОДУ). При решении ОДУ искомой величиной является функция. При использовании любых методов численного интегрирования необходимо, чтобы были заданы по крайней мере следующие величины:

набор точек в которых нужно найти решение;

само дифференциальное уравнение, записанное в некотором специальном виде, который будет описан ниже.

Один из наиболее эффективных алгоритмов интегрирования ОДУ основан на численном методе Рунге-Кутты четвертого порядка. Функция, реализующая этот метод, имеет вид rkfixed (y,x₁ ,x₂ , npoints,D)

y-вектор начальных условий размерности n, где n- порядок дифференциального уравнения или число уравнений в системе (если решается система уравнений);

x₁ , x₂ – граничные точки интервала, на котором ищется решение дифференциального уравнения. Начальные условия ,заданные в векторе y,- это значение решения в точке x₁ ;

npoints- число точек (не считая начальной точки), в которых ищется приближенное решение. При помощи этого аргумента определяется число строк (1+npoints) в матрице, возвращаемой функцией rkfixed;

D(x,y) – функция,возвращающая значение в виде вектора n элементов, содержащих первые производные неизвестных функций.

Если задачу об отыскании всех решений дифференциального уравнения удается свести к конечному числу алгебраических операций, операций интегрирования и дифференцирования известных функций, то говорят, что уравнениеинтегрируется в квадратурах. В приложениях крайне редко встречаются уравнения, интегрируемые в квадратурах. Поэтому для исследования дифференциальных уравнений широко используются приближенные, численные методы их решения.

Численное решение на отрезке [a, b] задачи Коши

состоит в построении таблицы приближенных значений

решенияy(x)в узлах сетки

a=x₀ , щелкните по соответствующей позиции в панели Matrix или введите с клавиатуры символ («точка с запятой»)

Определим шаг формулы Эйлера — шаг интегрирования

Для того чтобы ввести нижний индекс переменной, щелкните по соответствующей позиции в панели Matrix или в панели Calculator

Определим по формулам Эйлера значения приближенного решения в узлах сетки

Выведем в рабочий документ вычисленные значения решения

Построим график найденного решения y(x)

Для того чтобы вывести значение переменной в рабочий документ, введите имя переменной, знак равенства и щелкните по рабочему документу вне выделяющей рамки

Для того чтобы построить график приближенного решения, щелкните в панели Graph по пиктограмме декартова графика, введите в помеченной позиции возле оси абсцисс обозначение компонент вектора, содержащего значения узлов сетки, а в позиции возле оси ординат — обозначение компонент вектора, содержащего значения приближенного решения в узлах сетки; затем щелкните по свободному месту в рабочем документе вне поля графиков.

Метод Эйлера допускает простуюгеометрическую интерпретацию. Пусть известна точка (x_i ,y_i ) интегральной кривой уравненияy’=f(x, y).

Касательная к интегральной кривой уравнения, проходящая через эту точку, определяется уравнением

Следовательно, вычисленная методом Эйлера точка (x_i+ 1 ,y_i+ 1 ),

Найдем методом Эйлера на [0, 1] с шагом h=0.2 и с шагом h=0.1 приближенное решение задачи Коши

y’ = sin x – cosy,y(0)=1.

Изобразим приближенные решения графически.

Расчетные формулы метода Эйлера для решения этой задачи имеют вид

x0=0, y0= 1, xi+1 = xi + 0.2, yi+1 = yi + 0.2(sinxi — cosyi), i =0, 1, . 4

xi+1 = xi + 0.2, yi+1 = yi + 0.2(sinxi — cosyi), i =0, 1, . 9

Определим правую часть уравнения

Знак присваивания можно ввести щелчком по соответствующей позиции в панели Evaluation.

Определим диапазон изменения номера точки i=0,1, . 4 для вычислений с шагом h=0.2

Для того чтобы ввести символ диапазона изменения индекса , щелкните по соответствующей позиции в панели Matrix или введите с клавиатуры символ («точка с запятой»)

При решении задачи с шагом h=0.2 назовем шаг h1, аргумент — x1, а решение — y1.

Определим начальное условие

Определим шаг формулы Эйлера — шаг интегрирования

Определим по формулам Эйлера значения приближенного решения в узлах сетки

Выведем в рабочий документ вычисленные значения решения

Построим график найденного решения y1(x1)

Определим диапазон изменения номера точки i=0,1, . 9 для вычислений с шагом h=0.1

При решении задачи с шагом h=0.1 назовем шаг h2, аргумент — x2, а решение — y2.

Определим начальное условие

Определим шаг формулы Эйлера — шаг интегрирования

Определим по формулам Эйлера значения приближенного решения в узлах сетки

Выведем в рабочий документ вычисленные значения решения. Для сравнения рядом выведены значения решения, вычисленные с большим шагом

Построим график решения y2(x2)

Построим на одном графике оба приближенные решения

Для того чтобы одновременно построить графики нескольких функций от разных аргументов, щелкните в панели Graph по пиктограмме декартова графика, введите в помеченной позиции у оси абсцисс имя первого аргумента, запятую, имя второго аргумента, и т.д., разделяя имена аргументов запятой.

Аналогично, в позиции возле оси ординат введите имя функции первого аргумента, запятую, имя функции второго аргумента и т.д.разделяя имена функций запятой.

Когда функции определены, щелкните по рабочему документу вне поля графиков.

Методом Рунге-Кутты четвертого порядкаточности называют одношаговый метод, относящийся к широкому классу методов Рунге-Кутты. В этом методе величиныy_i+ 1 вычисляются по следующим формулам:

Найдем на [0, 1]приближенноерешение задачи Кошиy’ = sinx– cosy,y(0)=1методом Рунге-Кутты 4-го порядка с шагом h=0.2 и методом Эйлера с тем же шагом.Изобразим оба приближенные решения графически

Для решения задачиметодом Рунге-Кутты воспользуемся функциейrkfixed

Определим начальное условие — решение в начальной точке

Определим правую часть уравнения

Знак присваивания можно ввести щелчком по соответствующей позиции в панели Evaluation.

Вычислим приближенное решение на отрезке [0,1], выполнив n=1/h=5 одинаковых шагов, методом Рунге-Кутты 4-го порядка; обозначим приближенное решение Y

Выведем в рабочий документ вычисленное приближенное решение

Для того чтобы вывести значение переменной в рабочий документ, введите имя переменной, знак равенства и щелкните по рабочему документу вне выделяюшей рамки

В первом столбце приведены значения x, во втором столбце — соответствующие значения приближенного решения

Решим ту же задачу методом Эйлера

Выведем в рабочий документ вычисленное приближенное решение, и, для сравнения, решение, вычисленное методом Рунге-Кутты

Построим графики приближенных решений

Для того чтобы построить график приближенного решения, щелкните в панели Graph по пиктограмме декартова графика, введите в помеченной позиции возле оси абсцисс имя первого столбца матрицы Y, содержащего значения x в узлах сетки, а в позиции возле оси ординат — имя второго столбца, содержащего значения приближенного решения в узлах сетки; затем щелкните по свободному месту в рабочем документе вне поля графиков.

Для того чтобы одновременно построитьграфики нескольких функций от разных аргументов, щелкните в панели Graph по пиктограмме декартова графика, введите в помеченной позиции у оси абсцисс имя первого аргумента, запятую, имя второго аргумента, и т.д., разделяя имена аргументов запятой.

Аналогично, в позиции возле оси ординат введите имя функции первого аргумента, запятую, имя функции второго аргумента и т.д. разделяя имена функций запятой.

Когда функции определены, щелкните по рабочему документу вне поля графиков.

Для того чтобы ввести номер столбца, щелкните по соответствующему символу в панели Matrix

Для того чтобы изменить стиль изображения, щелкните дважды по полю графиков и установите в окне соответствующие параметры

Миллионы людей занимаются математическими расчетами, иногда в силу влечения к таинствам математики и ее внутренней красоте, а чаще в силу профессиональной или иной необходимости, не говоря уже об учебе. Ни одна серьезная разработка в любой отрасли науки и производства не обходится без трудоемких математических расчетов. Система Mathcadпользуется огромной популярностью во всем мире, позволяя готовить вполне профессиональные документы, имеющие вид статей и книг по математике.

Программа MicrosoftExcel входит в офисный пакет MicrosoftOfficeи предназначена для подготовки и обработки электронных таблиц под управлением операционной системой Windows. MicrosoftExcel – это многофункциональный, мощный редактор электронных таблиц. Он представляет возможность производить различные расчеты, составлять списки, сметы и что немаловажно, строить наглядные графики и диаграммы.

MathCAD – это мощная и в то же время простая универсальная среда для решения задач в различных отраслях науки и техники, финансов и экономики, физики и астрономии, строительства и архитектуры, математики и статистики, организации производства и управления… Она располагает широким набором инструментальных, информационных и графических средств. Сегодня MathCAD – одна из самых популярных математических систем. Она пользуется большим спросом у студентов, инженеров, экономистов, менеджеров, научных работников и всех тех, чья деятельность связана с количественными методами расчета.

Microsoft Excel ‑ средство для работы с электронными таблицами, намного превышающее по своим возможностям существующие редакторы таблиц, первая версия данного продукта была разработана фирмой Microsoft в 1985 году. Microsoft Excel ‑ это простое и удобное средство, позволяющее проанализировать данные и, при необходимости, проинформировать о результате заинтересованную аудиторию, используя Internet. Microsoft ® Excel разработан фирмой Microsoft, и является на сегодняшний день самым популярным табличным редактором в мире. Кроме стандартных возможностей его отличает следующие возможности, он выводит на поверхность центральные функции электронных таблиц и делает их более доступными для всех пользователей. Для облегчения работы пользователя упрощены основные функции, создание формул, форматирование, печать и построение графиков.

Данная курсовая работа позволила мне более близко познакомится с пакетом прикладных программ MathCAD и MicrosoftExcel. Мной было рассмотрено несколько способов решения дифференциальных уравнений.

Всё это позволило в полном объеме усвоить лекционный материал и понять перспективы использования вычислительной техники при решении различных задач практического характера.

1. Индейкин В. В. Табличный редактор Microsoft Excel. Учебное пособие. – Казань, 1999. – 75с.

2. Кудрявцев Е. М. MathCAD 2000 Pro. – М.: ДМК Пресс, 2001. – 571с.

Методические рекомендации для преподавателей математики и студентов средних специальных учебных заведений по теме «Дифференциальные уравнения»

Разделы: Математика

I. Обыкновенные дифференциальные уравнения

1.1. Основные понятия и определения

Дифференциальным уравнением называется уравнение, связывающее между собой независимую переменную x, искомую функцию y и её производные или дифференциалы.

Символически дифференциальное уравнение записывается так:

Дифференциальное уравнение называется обыкновенным, если искомая функция зависит от одного независимого переменного.

Решением дифференциального уравнения называется такая функция , которая обращает это уравнение в тождество.

Порядком дифференциального уравнения называется порядок старшей производной, входящей в это уравнение

1. Рассмотрим дифференциальное уравнение первого порядка

Решением этого уравнения является функция y = 5 ln x. Действительно, , подставляя y’ в уравнение, получим – тождество.

А это и значит, что функция y = 5 ln x– есть решение этого дифференциального уравнения.

2. Рассмотрим дифференциальное уравнение второго порядка y» — 5y’ +6y = 0. Функция – решение этого уравнения.

Действительно, .

Подставляя эти выражения в уравнение, получим: , – тождество.

А это и значит, что функция – есть решение этого дифференциального уравнения.

Интегрированием дифференциальных уравнений называется процесс нахождения решений дифференциальных уравнений.

Общим решением дифференциального уравнения называется функция вида ,в которую входит столько независимых произвольных постоянных, каков порядок уравнения.

Частным решением дифференциального уравнения называется решение, полученное из общего решения при различных числовых значениях произвольных постоянных. Значения произвольных постоянных находится при определённых начальных значениях аргумента и функции.

График частного решения дифференциального уравнения называется интегральной кривой.

1.Найти частное решение дифференциального уравнения первого порядка

xdx + ydy = 0, если y = 4 при x = 3.

Решение. Интегрируя обе части уравнения, получим

Замечание. Произвольную постоянную С, полученную в результате интегрирования, можно представлять в любой форме, удобной для дальнейших преобразований. В данном случае, с учётом канонического уравнения окружности произвольную постоянную С удобно представить в виде .

— общее решение дифференциального уравнения.

Частное решение уравнения, удовлетворяющее начальным условиям y = 4 при x = 3 находится из общего подстановкой начальных условий в общее решение: 3 2 + 4 2 = C 2 ; C=5.

Подставляя С=5 в общее решение, получим x 2 +y 2 = 5 2 .

Это есть частное решение дифференциального уравнения, полученное из общего решения при заданных начальных условиях.

2. Найти общее решение дифференциального уравнения

Решением этого уравнения является всякая функция вида , где С – произвольная постоянная. Действительно, подставляя в уравнения , получим: , .

Следовательно, данное дифференциальное уравнение имеет бесконечное множество решений, так как при различных значениях постоянной С равенство определяет различные решения уравнения .

Например, непосредственной подстановкой можно убедиться, что функции являются решениями уравнения .

Задача, в которой требуется найти частное решение уравнения y’ = f(x,y) удовлетворяющее начальному условию y(x₀) = y₀, называется задачей Коши.

Решение уравнения y’ = f(x,y), удовлетворяющее начальному условию, y(x₀) = y₀, называется решением задачи Коши.

Решение задачи Коши имеет простой геометрический смысл. Действительно, согласно данным определениям, решить задачу Коши y’ = f(x,y) при условии y(x₀) = y₀,, означает найти интегральную кривую уравнения y’ = f(x,y) которая проходит через заданную точку M₀(x₀,y₀).

II. Дифференциальные уравнения первого порядка

2.1. Основные понятия

Дифференциальным уравнением первого порядка называется уравнение вида F(x,y,y’) = 0.

В дифференциальное уравнение первого порядка входит первая производная и не входят производные более высокого порядка.

Уравнение y’ = f(x,y) называется уравнением первого порядка, разрешённым относительно производной.

Общим решением дифференциального уравнения первого порядка называется функция вида , которая содержит одну произвольную постоянную.

Пример. Рассмотрим дифференциальное уравнение первого порядка .

Решением этого уравнения является функция .

Действительно, заменив в данном уравнении, его значением, получим

то есть 3x=3x

Следовательно, функция является общим решением уравнения при любом постоянном С.

Найти частное решение данного уравнения, удовлетворяющее начальному условию y(1)=1 Подставляя начальные условия x = 1, y =1 в общее решение уравнения , получим откуда C = 0.

Таким образом, частное решение получим из общего подставив в это уравнение, полученное значение C = 0 – частное решение.

2.2. Дифференциальные уравнения с разделяющимися переменными

Дифференциальным уравнением с разделяющимися переменными называется уравнение вида: y’=f(x)g(y) или через дифференциалы , где f(x) и g(y)– заданные функции.

Для тех y, для которых , уравнение y’=f(x)g(y) равносильно уравнению, в котором переменная y присутствует лишь в левой части, а переменная x- лишь в правой части. Говорят, «в уравнении y’=f(x)g(y разделим переменные».

Уравнение вида называется уравнением с разделёнными переменными.

Проинтегрировав обе части уравнения по x, получим G(y) = F(x) + C– общее решение уравнения, где G(y) и F(x) – некоторые первообразные соответственно функций и f(x), C произвольная постоянная.

Алгоритм решения дифференциального уравнения первого порядка с разделяющимися переменными

Производную функции переписать через её дифференциалы
Разделить переменные.
Проинтегрировать обе части равенства, найти общее решение.
Если заданы начальные условия, найти частное решение.

Решить уравнение y’ = xy

Решение. Производную функции y’ заменим на

разделим переменные

проинтегрируем обе части равенства:

Ответ:

Найти частное решение уравнения

Это—уравнение с разделенными переменными. Представим его в дифференциалах. Для этого перепишем данное уравнение в виде Отсюда

Интегрируя обе части последнего равенства, найдем

Подставив начальные значения x₀ = 1, y₀ = 3 найдем С 9=1-1+C, т.е. С = 9.

Следовательно, искомый частный интеграл будет или

Составить уравнение кривой, проходящей через точку M(2;-3) и имеющей касательную с угловым коэффициентом

Решение. Согласно условию

Это уравнение с разделяющимися переменными. Разделив переменные, получим:

Проинтегрировав обе части уравнения, получим:

Используя начальные условия, x = 2 и y = — 3 найдем C:

Следовательно, искомое уравнение имеет вид

2.3. Линейные дифференциальные уравнения первого порядка

Линейным дифференциальным уравнением первого порядка называется уравнение вида y’ = f(x)y + g(x)

где f(x) и g(x) — некоторые заданные функции.

Если g(x)=0 то линейное дифференциальное уравнение называется однородным и имеет вид: y’ = f(x)y

Если то уравнение y’ = f(x)y + g(x) называется неоднородным.

Общее решение линейного однородного дифференциального уравнения y’ = f(x)y задается формулой: где С – произвольная постоянная.

В частности, если С =0, то решением является y = 0 Если линейное однородное уравнение имеет вид y’ = ky где k — некоторая постоянная, то его общее решение имеет вид: .

Общее решение линейного неоднородного дифференциального уравнения y’ = f(x)y + g(x) задается формулой ,

т.е. равно сумме общего решения соответствующего линейного однородного уравнения и частного решения данного уравнения.

Для линейного неоднородного уравнения вида y’ = kx + b,

где k и b— некоторые числа и частным решением будет являться постоянная функция . Поэтому общее решение имеет вид .

Пример. Решить уравнение y’ + 2y +3 = 0

Решение. Представим уравнение в виде y’ = -2y — 3 где k = -2, b= -3 Общее решение задается формулой .

Следовательно, где С – произвольная постоянная.

Ответ:

2.4. Решение линейных дифференциальных уравнений первого порядка методом Бернулли

Нахождение общего решения линейного дифференциального уравнения первого порядка y’ = f(x)y + g(x) сводится к решению двух дифференциальных уравнений с разделенными переменными с помощью подстановки y=uv, где u и v — неизвестные функции от x. Этот метод решения называется методом Бернулли.

Алгоритм решения линейного дифференциального уравнения первого порядка

1. Ввести подстановку y=uv.

2. Продифференцировать это равенство y’ = u’v + uv’

3. Подставить y и y’ в данное уравнение: u’v + uv’ = f(x)uv + g(x) или u’v + uv’ + f(x)uv = g(x).

4. Сгруппировать члены уравнения так, чтобы u вынести за скобки:

5. Из скобки, приравняв ее к нулю, найти функцию

Это уравнение с разделяющимися переменными:

Разделим переменные и получим:

Откуда . .

6. Подставить полученное значение v в уравнение (из п.4):

и найти функцию Это уравнение с разделяющимися переменными:

7. Записать общее решение в виде: , т.е. .

Найти частное решение уравнения y’ = -2y +3 = 0 если y =1 при x = 0

Решение. Решим его с помощью подстановки y=uv, .y’ = u’v + uv’

Подставляя y и y’ в данное уравнение, получим

Сгруппировав второе и третье слагаемое левой части уравнения, вынесем общий множитель u за скобки

Выражение в скобках приравниваем к нулю и, решив полученное уравнение, найдем функцию v = v(x)

Получили уравнение с разделенными переменными. Проинтегрируем обе части этого уравнения: Найдем функцию v:

Подставим полученное значение v в уравнение Получим:

Это уравнение с разделенными переменными. Проинтегрируем обе части уравнения: Найдем функцию u = u(x,c) Найдем общее решение: Найдем частное решение уравнения, удовлетворяющее начальным условиям y = 1 при x = 0:

Ответ:

III. Дифференциальные уравнения высших порядков

3.1. Основные понятия и определения

Дифференциальным уравнением второго порядка называется уравнение, содержащее производные не выше второго порядка. В общем случае дифференциальное уравнение второго порядка записывается в виде: F(x,y,y’,y») = 0

Общим решением дифференциального уравнения второго порядка называется функция вида , в которую входят две произвольные постоянные C₁ и C₂.

Частным решением дифференциального уравнения второго порядка называется решение, полученное из общего при некоторых значениях произвольных постоянных C₁ и C₂.

3.2. Линейные однородные дифференциальные уравнения второго порядка с постоянными коэффициентами.

Линейным однородным дифференциальным уравнением второго порядка с постоянными коэффициентами называется уравнение вида y» + py’ +qy = 0, где pи q— постоянные величины.

Алгоритм решения однородных дифференциальных уравнений второго порядка с постоянными коэффициентами

1. Записать дифференциальное уравнение в виде: y» + py’ +qy = 0.

2. Составить его характеристическое уравнение, обозначив y» через r 2 , y’ через r, yчерез 1: r 2 + pr +q = 0

3.Вычислить дискриминант D = p 2 -4q и найти корни характеристического уравнения; при этом если:

а) D > 0; следовательно, характеристическое уравнение имеет два различных действительных корня . Общее решение дифференциального уравнения выражается в виде , где C₁ и C₂ — произвольные постоянные.

б) D = 0; следовательно, характеристическое уравнение имеет равные действительные корни . Общее решение дифференциального уравнения выражается в виде

Общее решение

Дифференцируя общее решение, получим

Составим систему из двух уравнений

Подставим вместо ,и заданные начальные условия:

Таким образом, искомым частным решением является функция

2. Найти частное решение уравнения

Название: Решение дифференциальных уравнений. Обзор
Раздел: Рефераты по информатике, программированию
Тип: курсовая работа Добавлен 01:47:37 18 января 2011 Похожие работы
Просмотров: 4852 Комментариев: 21 Оценило: 2 человек Средний балл: 5 Оценка: неизвестно Скачать

2. а)

б)

в)

г)

источники:

http://www.bestreferat.ru/referat-197208.html

http://urok.1sept.ru/articles/527195