Задание 3

Приближенное интегрирование с заданным шагом

1. Цель работы

Изучение способов приближенного интегрирования

2. Основные теоретические положения

1. Постановка задачи

Пусть необходимо вычислить определенный интеграл

I = . (13)

Методы приближенного интегрирования основаны на использовании геометрической интерпретации значения определенного интеграла, как площади криволинейной трапеции, ограниченной осью абсцисс, прямыми x = a, x = b и кривой f (x) (рис.2).

f(x) f(x)

x 0 a b

Рис.2

Для вычисления интересующей нас площади (см. рис.3) разобьем область интегрирования на n равных частей точками:

x = a, x₁, x₂, ... , x_i, x_i+1, ... , x _n = b. (14)

f(x) f(x) Рис.3

I_i

x

0 х= a х_i x_i+1 х= b

Тогда I = ,

где I_i = . (15)

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

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

Как правило, в результате эксперимента получают дискретные данные, т.е. в узлах х_i производят измерение значений некоторой функции y_i, (см.работу 1).

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

Рассмотрим три простейших метода приближенного интегрирования.

3. Типы формул интегрирования

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

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

Согласно методу прямоугольников, область между точками разбиения интервала интегрирования [a,b] заменяется прямоугольником, высота которого соответствует координате Y одной из точек, а ширина равна расстоянию между точками. Значение интеграла определяется по следующей формуле:

I= . (16)

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