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

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

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

1. Формуле левых прямоугольников:

2. Формуле правых прямоугольников:

Формуле прямоугольников(средних):

Составные квадратурные формулы

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

1. Для левых прямоугольников:

2. Для правых прямоугольников:

3. Для средних прямоугольников:

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

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

Сравнение применения различных формул прямоугольников
Формула левых прямоугольников   Формула средних прямоугольников   Формула правых прямоугольников

Составные формулы для равномерных сеток

Равномерную сетку можно описать следующим набором формул:

где — шаг сетки.

Для равномерных сеток формулы прямоугольников можно записать в виде следующих формул Котеса:

1. Составная формула левых прямоугольников:

2. Составная формула правых прямоугольников:

Погрешность метода

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

Для формулы прямоугольников (средних)

Для составных формул правых и левых прямоугольников на равномерной сетке:

Для составной формулы прямоугольников:

Пример реализации

Формула средних прямоугольников для аналитически заданной функции, написанная на С

#include <stdio.h>

#include <math.h>

double f(double x){ //Подынтегральная функция

return sin(x); //Например, sin(x)

}

double rectangle_integrate(double a, double b, int n, double (*f)(double) ){

double result, h;

int i;

h = (b-a)/n; //Шаг сетки

result = 0.0;

for(i=1; i <= n; i++){

result += f( a + h * (i - 0.5) ); //Вычисляем в средней точке и добавляем в сумму

}

result *= h;

return result;

}

int main(void){

double integral;

integral=rectangle_integrate(0,2,100,f);

printf("The value of the integral is: %lf \n", integral);

return 0;

}