Формула Симпсона для приближенного вычисления определенного интеграла

Метод парабол

Формула Симпсона для приближенного вычисления определенного интеграла

где: - длина каждого из маленьких отрезков или шаг;

f(x_i) - значения подынтегральной функции в точках x₀,x₁,x₂,x₃,…,x_2n-2,x_2n-1,x_2n.

Детализируя это нагромождение, разберу формулу подробнее:

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

- сумма членов, с чётными индексами умножаемая на 2.

- сумма членов с нечётными индексами умножается на 4.

На основании полученных данных строим график (рисунок 2), который показывает погрешность:

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

Решение:

Метод левых прямоугольников

Метод правых прямоугольников

Метод трапеции

Метод Симпсона

4. Приближенное решение обыкновенных дифференциальных уравнений первого порядка. Задача Коши

Многие задачи науки и техники сводятся к решению обыкновенных дифференциальных уравнений (ОДУ). ОДУ называются такие уравнения, которые содержат одну или несколько производных от искомой функции. В общем виде ОДУ можно записать следующим образом: , где x - независимая переменная, yⁱ - i-ая производная от искомой функции. n - порядок уравнения. Общее решение ОДУ n-го порядка содержит n произвольных постоянных c₁,.., c_n,т.е. общее решение имеет вид y=ц(x, c₁, …, c_n).

Для выделения единственного решения необходимо задать n дополнительных условий. В зависимости от способа задания дополнительных условий существуют два различных типа задач: задача Коши и краевая задача. Если дополнительные условия задаются в одной точке, то такая задача называется задачей Коши. Дополнительные условия в задаче Коши называются начальными условиями. Если же дополнительные условия задаются в более чем одной точке, т.е. при различных значениях независимой переменной, то такая задача называется краевой. Сами дополнительные условия называются краевыми или граничными.

Ясно, что при n=1 можно говорить только о задачи Коши.

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

4.1 Численные методы решения задачи Коши для ОДУ первого порядка

Постановка задачи. Найти решение ОДУ первого порядка на отрезке [x₀, x_n] при условии y(x₀)=y₀.

При нахождении приближенного решения будем считать, что вычисления проводятся с расчетным шагом , расчетными узлами служат точки xi=x0+ih, (i=0,1,…,n) промежутка [x₀, x_n].

Целью является построение таблицы.

Таблица 2

Т.е. ищутся приближенные значения y в узлах сетки.

Интегрируя уравнение на отрезке [x_i,x_i+1]получим

(4.5)

Вполне естественным (но не единственным) путем получения численного решения является замена в нем интеграла какой-либо квадратурной формулой численного интегрирования. Если воспользоваться простейшей формулой левых прямоугольников первого порядка

(4.6)

то получим явную формулу Эйлера:

(4.7)

Порядок расчетов:

Зная , находим , затем т.д..

4.2 Геометрическая интерпретация метода Эйлера

Пользуясь тем, что в точке x₀ известно решение y(x₀) = y₀ и значение его производной , можно записать уравнение касательной к графику искомой функции y=y(x)

в точке (x₀,y₀):

При достаточно малом шаге h ордината , этой касательной, полученная подстановкой в правую часть значения , должна мало отличаться от ординаты y(x₁) решения y(x) задачи Коши. Следовательно, точка (x1,y1) пересечения касательной с прямой x = x₁ может быть приближенно принята за новую начальную точку. Через эту точку снова проведем прямую , которая приближенно отражает поведение касательной к y=y(x) в точке (x1, y(x1)). Подставляя сюда x2=x1+h (т.е. пересечение с прямой x = x₂), получим приближенное значение y(x) в точке x₂: , и т.д. В итоге для i-й точки получим формулу Эйлера.

Рисунок 7. Метод Эйлера

Явный метод Эйлера имеет первый порядок точности или аппроксимации. Если использовать формулу правых прямоугольников:

то придем к методу

(4.8)

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

Неявный метод Эйлера имеет первый порядок точности или аппроксимации.

Модифицированный метод Эйлера: в данном методе вычисление y_i+1 состоит из двух этапов:

(4.9)

(4.10)

Данная схема называется еще методом предиктор - корректор (предсказывающее - исправляющее). На первом этапе приближенное значение предсказывается с невысокой точностью (h), а на втором этапе это предсказание исправляется, так что результирующее значение имеет второй порядок точности.

Решение

Метод Эйлера