Решение уравнений способом эйлера

Дифференциальное уравнение Эйлера и методы его решения

Более общее уравнение Эйлера имеет вид:
.
Это уравнение подстановкой t = ax+b приводится к более простому виду, которое мы и будем рассматривать.

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

Рассмотрим уравнение Эйлера:
(1) .
Оно сводится к линейному уравнению с постоянными коэффициентами подстановкой:
x = e t .
Действительно, тогда
;
;
;

;
;
.

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

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

Рассмотрим однородное уравнение Эйлера:
(2) .
Ищем решение уравнения (2) в виде
.
;
;
.
.
Подставляем в (2) и сокращаем на x k . Получаем характеристическое уравнение:
.
Решаем его и получаем n корней, которые могут быть комплексными.

Рассмотрим действительные корни. Пусть ki – кратный корень кратности m . Этим m корням соответствуют m линейно независимых решений:
.

Рассмотрим комплексные корни. Они появляются парами вместе с комплексно сопряженными. Пусть ki – кратный корень кратности m . Выразим комплексный корень ki через действительную и мнимую части:
.
Этим m корням и m комплексно сопряженным корням соответствуют 2 m линейно независимых решений:
;
;
.
.

После того как получены n линейно независимых решений, получаем общее решение уравнения (2):
(3) .

Примеры

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

Рассмотрим неоднородное уравнение Эйлера:
.
Метод вариации постоянных (метод Лагранжа) также применим и к уравнениям Эйлера.

Сначала мы решаем однородное уравнение (2) и получаем его общее решение (3). Затем считаем постоянные функциями от переменной x . Дифференцируем (3) n – 1 раз. Получаем выражения для n – 1 производных y по x . При каждом дифференцировании члены, содержащие производные приравниваем к нулю. Так получаем n – 1 уравнений, связывающих производные . Далее находим n -ю производную y . Подставляем полученные производные в (1) и получаем n -е уравнение, связывающее производные . Из этих уравнений определяем . После чего интегрируя, получаем общее решение уравнения (1).

Пример

Неоднородное уравнение Эйлера со специальной неоднородной частью

Рассмотрим уравнение Эйлера со специальной неоднородной частью:
(4)
,
где – многочлены от степеней и , соответственно.

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

Автор: Олег Одинцов . Опубликовано: 14-08-2013 Изменено: 24-10-2020

Источник

Решение уравнений способом эйлера

Сначала находим общее решение однородного уравнения Эйлера;

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

Общее решение неоднородного уравнения будет представлять собой сумму общего решения однородного уравнения (шаг \(1\)) и частного решения неоднородного уравнения (шаг \(2\)).

Сделаем подстановку \(x = .\) Поскольку \[y» = >\left( <\frac<<y>><>> — \frac<><

>> \right),\] то дифференциальное уравнение принимает вид: \[ <4\cancel>\cancel>\left( <\frac<<y>><>> — \frac<><
>> \right) + y = 0,>\;\; <\Rightarrow 4\frac<<y>><>> — 4\frac<><
> + y = 0.> \] Вычислим корни соответствующего характеристического уравнения: \[ <4— 4k + 1 = 0,>\;\; <\Rightarrow D = 16 - 16 = 0,>\;\; <\Rightarrow k = \frac<4><<2 \cdot 4>> = \frac<1><2>.> \] Уравнение имеет один корень второго порядка. Тогда общее решение для функции \(y\left( t \right)\) будет определяться выражением \[y\left( t \right) = \left( <+ t> \right)<2>\normalsize>>.\] Решение для исходной функции \(y\left( x \right)\) можно записать в виде: \[ + \ln x> \right)><2>\normalsize>> > = <\left( <+ \ln x> \right)<2>\normalsize>> > = <\left( <+ \ln x> \right)\sqrt x ,> \] где \(,\) \(\) − произвольные действительные числа.

Для решения данного уравнения воспользуемся вторым способом, т.е. будем искать решение в форме \(y = .\) Тогда \[y’ = k>,\;\;y» = k\left( \right)>.\] Подставляя в исходное дифференциальное уравнение, получаем \[ <k\left( \right)> — xk> — 8 = 0,>\;\; <\Rightarrow k\left( \right) — k — 8 = 0,>\;\; <\Rightarrow \left[ \right) — k — 8> \right] = 0.> \] Соответствующее характеристическое уравнение имеет корни: \[ \right) — k — 8 = 0,>\;\; <\Rightarrow — 2k — 9 = 0,>\;\; <\Rightarrow D = 4 + 32 = 36,>\;\; <\Rightarrow > = \frac<<2 \pm 6>> <2>= 4, — 2.> \] Следовательно, общее решение для промежуточной функции \(y\left( t \right)\) выражается формулой \[y\left( t \right) = > + >.\] Возвращаясь к переменной \(x,\) получаем окончательный ответ: \[ > + > > = < + > > = < + \frac<<>><<>>.> \] Здесь \(,\) \(\) − произвольные действительные постоянные.

Читайте также:  Revive удобрение способ применения

Сначала построим общее решение однородного уравнения: \[y» + xy’ + y = 0.\] Сделаем подстановку: \[ ,\;\;y’ = >\frac<><

>,>\;\; >\left( <\frac<<y>><>> — \frac<><
>> \right).> \] В результате однородное уравнение примет вид: \[ <\cancel>\cancel>\left( <\frac<<y>><>> — \frac<><
>> \right) + \cancel\cancel>\frac<><
> + y = 0,>\;\; <\Rightarrow \frac<<y>><>> — \cancel<\frac<><
>> + \cancel<\frac<><
>> + y = 0,>\;\; <\Rightarrow \frac<<y>><>> + y = 0.> \] Решим характеристическое уравнение: \[ + 1 = 0,\;\; \Rightarrow > = \pm i.\] Как видно, корни характеристического уравнения мнимые. Поэтому общее решение однородного уравнения записывается в виде \[\left( t \right) = \cos t + \sin t,\] где \(\) и \(\) − действительные постоянные.

Теперь определим частное решение неоднородного уравнения \[\frac<<y>><>> + y = 5>.\] Принимая во внимание структуру правой части, будеи искать частное решение в виде \(\left( t \right) = a>,\) где \(a\) − некоторый постоянный коэффициент. Тогда \[\frac<>><

> = 2a>,\;\;\frac<<>><>> = 4a>.\] Подставим данную функцию вместе с ее производными в уравнение и найдем коэффициент \(a:\) \[ <4a> + a> = 5>,>\;\; <\Rightarrow 5a> = 5>,>\;\; <\Rightarrow a = 1.>\] Итак, частное решение неоднородного уравнения определяется выражением \[\left( t \right) = >.\] Теперь можно записать общее решение исходного неоднородного уравнения: \[ \left( t \right) + \left( t \right) > = <\cos t + \sin t + >.> \] Возвращаясь обратно к переменной \(x,\) получаем \[y\left( x \right) = \cos \left( <\ln x>\right) + \sin \left( <\ln x>\right) + >.\] Поскольку \(> = >> = ,\) окончательный ответ записывается в виде \[y\left( x \right) = \cos \left( <\ln x>\right) + \sin \left( <\ln x>\right) + .\]

Сначала найдем решение однородного уравнения: \[y» — 2xy’ + 2y = 0.\] Используя подстановку \(x = ,\) можно преобразовать последнее уравнение в уравнение с постоянными коэффициентами: \[\frac<<y>><>> — 3\frac<><

> + 2y = 0.\] Вычислим корни характеристического уравнения и запишем общее решение \(\left( t \right)\) однородного дифференциального уравнения: \[ <— 3k + 2 = 0,>\;\; <\Rightarrow D = 9 - 4 \cdot 2 = 1,>\;\; <\Rightarrow > = \frac<<3 \pm 1>> <2>= 2,1.> \] Следовательно, \[\left( t \right) = > + .\] Теперь рассмотрим неоднородное уравнение, которое можно записать через переменную \(t\) в виде \[ <\frac<<y>><>> — 3\frac<><
> + 2y = 6> + 4\ln \left( <> \right),>\;\; <\Rightarrow \frac<<y>><>> — 3\frac<><
> + 2y = 6> + 4t.> \] В показателе экспоненциальной функции в правой части содержится коэффициент \(2,\) который совпадает с одним из корней характеристического уравнения. Поэтому будем искать частное решение в форме \[\left( t \right) = at> + bt + c,\] где \(a, b\) и \(c\) − некоторые (пока еще) неизвестные числа.

Первая и вторая производные данной функции будут равны \[\frac<>><

> = a> + 2at> + b,\] \[ <\frac<<>><>> = 2a> + 2a> + 4at> > = <4a> + 4at>.> \] Подставим это в неоднородное уравнение: \[ <4a> + 4at> > — <3\left( > + 2at> + b> \right) > + <2\left( > + bt + c> \right) > = <6> + 4t,> \] \[ <\Rightarrow 4a> + \cancel<4at>> — 3a> > — <\cancel<6at>> — 3b > + <\cancel<2at>> + 2bt + 2c > = <6> + 4t,> \] Последнее равенство является тождественным, то есть оно справедливо для любых значений \(t.\) Приравнивая коэффициенты при членах с одинаковыми степенями в левой и правой части, получаем: \[ <\left\< \begin a = 6\\ 2b = 4\\ — 3b + 2c = 0 \end \right.,>\;\; <\Rightarrow \left\< \begin a = 6\\ b = 2\\ c = 3 \end \right..> \] Таким образом, частное решение неоднородного уравнения описывается выражением \[\left( t \right) = 6t> + 2t + 3.\] Теперь мы можем записать общее решение неоднородного уравнения Эйлера: \[ \left( t \right) + \left( t \right) > = <> + + 6t> + 2t + 3.> \] Поскольку \(t = \ln x,\) то окончательный ответ записывается как \[y\left( x \right) = + x + 6\ln x + 2\ln x + 3,\] где \(,\) \(\) − произвольные действительные числа.

Читайте также:  Засолить перцы холодным способом

Источник

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

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

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

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

Функция , стоящая под знаком интеграла, зависит от трех аргументов . Будем считать ее определенной и дважды непрерывно дифференцируемой по аргументу для всех значений, по аргументам же и — в некоторой области плоскости . Ниже предполагается, что мы всегда будем находиться внутри этой области.

Под понимается некоторая функция от

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

Геометрически функцию можно изобразить в плоскости некоторой линией , лежащей над отрезком (рис. 3).

Интеграл (9) является обобщением интегралов (3) и (6), с которыми мы встретились в задачах о линии наискорейшего ската и поверхности вращения наименьшей площади. Значение его зависит от выбора функции или от линии , и задача о его минимуме имеет следующий смысл.

Дано некоторое множество функций (10) (линий ). Среди них нужно найти ту функцию (линию ), для которой интеграл имеет наименьшее значение.

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

1) функция непрерывно дифференцируема на отрезке ;

2) на концах отрезка функция принимает заданные наперед значения

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

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

Рассмотрим семейство функций, зависящее от численного параметра

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

Интеграл (9), вычисленный для , будет некоторой функцией параметра

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

Так как функция дает минимальное значение интегралу, то функция должна иметь минимум при , и производная от нее в этой точке обязана обращаться в нуль

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

и придать условию (14) другую форму

Читайте также:  Народные способы лечения онемения пальцев ног

Может быть доказана следующая простая лемма.

Пусть выполняются условия:

1) функция непрерывна на отрезке ;

2) функция непрерывно дифференцируема на отрезке и на концах отрезка обращается в нуль.

Если при любой такой функции интеграл равен нулю, то отсюда следует, что .

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

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

Всегда можно построить функцию , непрерывно дифференцируемую на , положительную на и равную нулю всюду вне (рис. 4).

Такой будет, например, , определенная равенствами

Но для такой функции

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

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

или после вычисления производной по переменной

Равенство это является дифференциальным уравнением 2-го порядка относительно функции . Оно называется уравнением Эйлера .

Мы можем сделать следующее заключение.

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

Сами решения уравнения (17) обладают тем свойством, что производная

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

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

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

Их нужно определить так, чтобы интегральная кривая проходила через точки и , что доставляет два уравнения для нахождения постоянных и

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

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

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

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

Пример . Ранее было установлено, что [url]задача о линии наискорейшего ската[/url] может быть приведена к нахождению минимума интеграла

на множестве функций, удовлетворяющих граничным условиям .

Источник

Оцените статью
Разные способы