Решение линейных уравнений методом эйлера

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

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

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

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

;
;
.

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

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

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

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

Рассмотрим комплексные корни. Они появляются парами вместе с комплексно сопряженными. Пусть k_i – кратный корень кратности m . Выразим комплексный корень k_i через действительную и мнимую части:
.
Этим 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

Метод Эйлера: для чего он нужен, порядок действий и упражнения

Метод Эйлера: для чего он нужен, порядок действий и упражнения — Наука

Содержание:

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

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

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

Самый общий способ написать уравнение первой степени:

Что такое метод Эйлера?

Идея метода Эйлера заключается в нахождении численного решения дифференциального уравнения в интервале между X₀ и X_F .

Сначала интервал дискретизируется на n + 1 балл:

Которые получаются так:
Икс_я= х₀+ я

Где h — ширина или шаг подынтервалов:

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

Эта производная представляет собой наклон касательной к кривой функции y (x) точно в точке:

Затем делается приблизительный прогноз значения функции y (x) в следующей точке:

Затем была получена следующая приближенная точка решения, которая будет соответствовать:

Процедура повторяется для получения последовательных баллов.

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

Решенные упражнения

Упражнение 1

я) Пусть дифференциальное уравнение имеет вид:

При начальном условии x = a = 0; Y_к= 1

Используя метод Эйлера, получить приближенное решение Y в координате X = b = 0,5, разбивая интервал [a, b] на n = 5 частей.

Решение

Численные результаты резюмируются следующим образом:

Из чего делается вывод, что решение Y для значения 0,5 составляет 1,4851.

Примечание: для проведения расчетов, Smath studio, бесплатная программа для бесплатного использования.

Упражнение 2.

II) Продолжая работу с дифференциальным уравнением из упражнения I), найдите точное решение и сравните его с результатом, полученным методом Эйлера. Найдите ошибку или разницу между точным и приблизительным результатом.

Решение

Точное решение найти не очень сложно. Производная функции sin (x) известна как функция cos (x). Следовательно, решение y (x) будет:

Для выполнения начального условия и (0) = 1 константа C должна быть равна 1. Затем точный результат сравнивается с приблизительным:

Сделан вывод, что в расчетном интервале аппроксимация имеет три значащих цифры точности.

Упражнение 3.

III) Рассмотрим дифференциальное уравнение и его начальные условия, указанные ниже:

При начальном условии x₀ = 0; Y₀ = 1

Используйте метод Эйлера, чтобы найти приблизительные значения решения у (х) в интервале х = [0, 1,5]. Используйте шаг h = 0,1.

Решение

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

В электронной таблице на рисунке показаны три столбца (A, B, C), первый — это переменная Икс , второй столбец представляет переменную Y, а третий столбец — производная Y ‘.

Строка 2 содержит начальные значения Икс, Y, Y ‘ .

Шаг значения 0,1 помещен в ячейку абсолютного положения ($ D $ 4).

Начальное значение y0 находится в ячейке B2, а y1 — в ячейке B3. Чтобы вычислить y₁ используется формула:

Эта формула электронной таблицы будет иметь вид B3: = B2 + $ D $ 4 * C3.

Точно так же y2 будет в ячейке B4, и его формула показана на следующем рисунке:

На рисунке также показан график точного решения и точки A, B,…, P приближенного решения по методу Эйлера.

Ньютоновская динамика и метод Эйлера

Классическая динамика была разработана Исааком Ньютоном (1643 — 1727). Первоначальная мотивация Леонарда Эйлера (1707 — 1783) к разработке своего метода заключалась именно в решении уравнения второго закона Ньютона в различных физических ситуациях.

Второй закон Ньютона обычно выражается в виде дифференциального уравнения второй степени:

куда Икс представляет положение объекта в момент т. Указанный объект имеет массу м и подвергается силе F. Функция F связана с силой и массой следующим образом:

Для применения метода Эйлера требуются начальные значения времени. т, скорость v и положение Икс.

В следующей таблице объясняется, как, начиная с начальных значений t1, v1, x1, можно получить приближение скорости v2 и положения x2 в момент t2 = t1 + Δt, где Δt представляет небольшое увеличение и соответствует шагу в методе Эйлер.

Упражнение 4.

IV) Одна из фундаментальных проблем механики — это проблема блока массы M, привязанного к пружине (или пружине) с постоянной упругостью K.

Второй закон Ньютона для этой проблемы будет выглядеть так:

В этом примере для простоты мы возьмем M = 1 и K = 1. Найдите приблизительные решения для положения Икс и скорость v методом Эйлера на временном интервале [0, π / 2], разбивая интервал на 12 частей.

Возьмите 0 в качестве начального момента, начальную скорость 0 и начальное положение 1.

Решение

Численные результаты представлены в следующей таблице:

Также отображаются графики положения и скорости между временами от 0 до 1,44.

Предлагаемые упражнения для дома

Упражнение 1

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

y ’= — Exp (-y) с начальными условиями x = 0, y = -1 в интервале x = [0, 1]

Начните с шага 0,1. Постройте результат.

Упражнение 2.

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

y ’’ = -1 / y² с начальным условием t = 0; и (0) = 0,5; у ‘(0) = 0

Найдите решение на интервале [0,5; 1.0] с шагом 0,05.

Постройте результат: y vs t; y ‘vs t

Ссылки

1. Метод Эрлера Взято с wikipedia.org
2. Решатель Эйлера. Взято с en.smath.com

3 различия между пандемией и эпидемией (и примеры)

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

Решение линейных уравнений методом эйлера

Системой дифференциальных уравнений называется система вида

где x — независимый аргумент,

y_i — зависимая функция, ,

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

Численные методы решения систем дифференциальных уравнений.

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

Метод Рунге-Кутта четвертого порядка.

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

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

Численно ищется частное решение уравнения (2), которое удовлетворяет заданным начальным условиям, то есть решается задача Коши.

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

.

(3)

Функция f₂(x, y₁, y) в систему (3) введена формально для того, чтобы методы, которые будут показаны ниже, могли быть использованы для решения произвольной системы дифференциальных уравнений первого порядка. Рассмотрим несколько численных методов решения системы (3). Расчетные зависимости для i+1 шага интегрирования имеют следующий вид. Для решения системы из n уравнений расчетные формулы приведены выше. Для решения системы из двух уравнений расчетные формулы удобно записать без двойных индексов в следующем виде:

Метод Рунге-Кутта четвертого порядка.

где h — шаг интегрирования. Начальные условия при численном интегрировании учитываются на нулевом шаге: i=0, x=x₀, y₁=y₁₀, y=y₀.

Контрольное задание по зачетной работе.

Колебания с одной степенью свободы

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

Задание. Численно и аналитически найти:

1. закон движения материальной точки на пружинке х(t),
2. закон изменения силы тока I(t) в колебательном контуре (RLC — цепи) для заданных в табл.1,2 режимов. Построить графики искомых функций.

Свободные незатухающие колебания

Затухающее колебательное движение

Предельное апериодическое движение

Вынужденное колебание без сопротивления

Вынужденное колебание без сопротивления, явление резонанса

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

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