Курс лекций по основам алгоритмизации и программирования задач машиностроения. Кравченко Д.В. - 81 стр.

UptoLike

Составители: 

Рубрика: 

79
Главный
член
погрешности
формулы
средних
прямоугольников
(22)
на
каждом
отрезке
[x
i-1
, x
i
]
равен
fh
24
1
3
i
(x
i-1/2
).
Для формулы трапеций он равен fh
12
1
3
i
(x
i
), т. е. примерно вдвое
больше и имеет другой знак.
На основании этого можно записать уточненную формулу для вычисле-
ния определенного интеграла с использованием значений I1 и I2, вычислен-
ных по методам прямоугольников и трапеций:
I (2*I
прямоугольников
+ I
трапеций
) / 3 . (24)
Так как погрешность численного интегрирования определяется шагом раз-
биения, то уменьшая его, можно добиться большей точности вычисления инте-
грала.
Однако увеличение числа точек разбиения интервала не всегда возможно.
Особенно если функция задана в табличном виде.
Поэтому в такой ситуации повысить точность численного интегрирования
можно за счет повышения степени используемых интерполяционных много-
членов, как в методе Симпсона.
В методе Симпсона отрезок интегрирования [a, b]
разбивают на четное число n равных частей с шагом h.
На каждом отрезке [x
0
, x
2
], [x
2
, x
4
], …, [x
i-1
, x
i+1
], …, [x
i-2
,
x
n
] подынтегральную функцию f(x) заменяют интерполя-
ционным многочленом второй степени (рис. 5.5):
f(x) φ
i
(x) = a
i
x
2
+ b
i
x + c
i
,
x
i-1
x x
i+1
Коэффициенты этих квадратных трехчленов могут
быть найдены из условий равенства многочлена в точках
x
i
соответствующим табличным значениям функции y
i
.
В качестве φ
i
(x) можно принять интерполяционный многочлен Лагранжа
второй степени, проходящий через точки М
i-1
(x
i-1
, y
i-1
), М
i
(x
i
, y
i
), М
i+1
(x
i+1
, y
i+1
):
( )
(
)
(
)
( )( )
( )( )
( )( )
( )( )
( )( )
.y
xxxx
xxxx
y
xxxx
xxxx
y
xxxx
xxxx
x
1i
i1i1i1i
i1i
i
1ii1ii
1i1i
1i
1i1ii1i
1ii
i
+
++
+
+
+
+
×
+×
+
+×
=ϕ
Элементарная площадь S
i
(рис. 5.5) может быть вычислена с помощью оп-
ределенного интеграла.
Рис. 5.5. Вид
криволинейной
трапеции