6.4. Выбор шага интегрирования для равномерной сетки

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

Известно два подхода к решению данной задачи:

1) выбор шага hпо теоретическим оценкам (6.22);

2) по косвенным схемам (эмпирическим оценкам).

6.4.1. Выбор шага интегрирования по теоретическим оценкам погрешностей

Пусть требуется вычислить интеграл с точностью . Тогда, используя формулу дляR, выбирают шаг так, чтобы

| R | </2.

Учитывается также число знаков после запятой, чтобы погрешность округления не превышала /2.

Пример 6.2. С помощью формулы Симпсона вычислить значение интегралас точностью= 10^–3.

Решение. Выберем шагh.

;[a,b], т. е.[/4,/2].

Согласно соотношениям (6.22) получим

< 0,510^–3.

Вычислим f ^(IV)(x):

. (6.23)

Оценим | f ^(IV)| на отрезке [/4,/2]. Воспользуемся величинами из (6.23)и. Они положительные и убывают, следовательно, их максимальное значение в точкеx=/4.

При этом +< 81. Таким образом,< 0,510^–3;h⁴< 1410^–4;h0,19.

Однако для данного метода hвыбирается с учетом того, чтобы [/4,/2] делился на четное число отрезков. Этим двум требованиям отвечаетh =/24 = = 0,13 < 0,19, при которомn = = 6. Тогда, чтобы погрешность округления не превысила 0,510^–3, достаточно вычисления выполнить с четырьмя знаками после запятой.

Составим таблицу y=sinx/xсh=/24 = 730´ = 0,1309:

i	xi0	xi	sin x	y0, y6	y2m	y2m–1
0	45 00´	0,7854	0,7071	0,9003
1	52 30´	0,9163	0,7934			0,8659
2	60 00´	1,0472	0,8660		0,8270
3	67 30´	1,1781	0,9239			0,7843
4	75 00´	1,3090	0,9659		0,7379
5	82 30´	1,4399	0,9914			0,6885
6	90 00´	1,5708	1,0000	0,6366
Сумма				1,5369	1,5649	2,3386

Для n= 6 по формуле Симпсона

.

6.4.2. Выбор шага интегрирования по эмпирическим схемам

Двойной пересчет. В связи с тем что вычисления максимального значения по абсолютной величинеk-й производной приводят к громоздкости расчетов, на практике прибегают к искусственным приемам достижения заданной точности. А именно, определенный интеграл вычисляют по какой-либо квадратурной формуле дважды с шагомhиh/2, что удваивает числоn.

Определяют:

 если | I_n–I_2n| <, тоI=I_2n;

 если | I_n–I_2n| >, то берут шагh/4; (6.24)

 если | I_2n – I_4n | <  , то I = I_4n.

В качестве начального шага hможно рекомендоватьh=, гдеm = 2 для формул среднего и трапеций,m = 4 – для формулы Симпсона.

Схема Эйткина. На практике для повышения точности численного интегрирования широко используетсясхема Эйткина.

Расчеты интегралов проводятся с шагами h₁,h₂иh₃, при этом соотношение между ними. Получают три значенияI₁,I₂иI₃.

Далее производится уточнение по эмпирической формуле

. (6.25)

Порядок точности =.

Правило Рунге. Это наиболее популярное практическое правило, разработанное в предположении, чтоf(x)C⁴[a,b] для квадратурных формул прямоугольников и трапеций,f(x)C⁶[a,b] – для формулы Симпсона. В этом случае можно показать, что погрешностиR(h,f) имеют следующие представления приh0:

(6.26)

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

; (6.27)

где k= 2,m= 2 – для формул прямоугольников и трапеций;k= 4,m= 2 – для формулы Симпсона.

Другие оценки погрешности

1. Приближенной оценкой погрешности могут быть:

– для формул трапеций и прямоугольников;

– для формулы Симпсона.

2. Следует заметить, что эмпирические формулы (6.24), (6.25), (6.27) предполагают и автоматическое изменение шага интегрирования h. Для этой цели имеется другая схема расчета, заключающаяся в следующем.

Анализ составных формул (6.19), (6.20), (6.21) для вычисления интегралов I_П,I_Т,I_Споказывает, что точное значение интеграла находится междуI_ПиI_Т, при этом имеет место соотношение

I_С= (2I_П +I_Т) / 3. (6.28)

Соотношение (6.28) используется и для контроля погрешности вычисления. Если │I_С–I_Т│ ≥ ε, то шаг уменьшают вдвое и расчет повторяют. Если точность достигнута, то окончательное значение интеграла получают по формуле (6.28).