Численное решение дифференциального уравнения методом Эйлера
Поможем в ✍️ написании учебной работы
Поможем с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой

Явный метод Эйлера:

На малом отрезке  представим производную в виде конечно-разностного отношения. Тогда

; Данное соотношение позволяет выразить значение y в конце отрезка  через значение в начале отрезка. В частности,  и т.д. (Рис.16)

Основной вопрос, который возникает при построении приближенного решения следующий: будет ли приближенная зависимость стремиться к точной, если . Исследуем данный вопрос на частном примере – расчёте переходной характеристики RC-цепи. (Рис.17)  разряжен

Явный метод Эйлера: (Рис. 18)

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

Из проведенного анализа следует, что от выбранного шага разбиения времени на интервалы  будет зависеть сходимость приближённого решения к асимптотическому значению : если , построенное по методу Эйлера решение не только не стремится к значению , а вообще не стремится к какому-то конечному значению (верхняя зависимость на рисунке). Асимптотическое приближение к  будет только при выборе . В этом случае говорят, что метод неустойчив при  и устойчив при . Накапливающаяся ошибка при большом шаге разбиения связана с тем, что производная на этом шаге существенно меняется и замена её конечно-разностным отношением недопустима.

Неявный метод Эйлера

В данном случае производная по-прежнему аппроксимируется конечно-разностным отношением, однако  берется не в начале, а в конце интервала разбиения, т.е. . Поскольку  содержится в обеих частях уравнения, то чтобы его найти, нужно решить данное уравнение. Иными словами  связан с  неявно. Проведём исследование устойчивости данного метода также на примере RC-цепи.

 – метод устойчив при любом

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

 

Рис. 16 Рис. 17 Рис.18

 

 

Дата: 2019-03-05, просмотров: 192.