считаем, что эта масса сосредоточена в центре масс прямоугольника (совпадает с его центром) с координатами
x, 2y
dK y
m
. Тогда dKx
y
f(x)
x
0
a x
y xdx K x
b
f (x)dx; X c
a
1 |
y ydx |
1 |
|
y 2dx и |
||
2 |
2 |
|||||
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y/2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||||
|
|
x+dx |
b x |
|
x |
|
||||||
Рис. 49. |
|
|
|
|
|
|||||||
|
|
|
|
|
||||||||
|
|
|
b |
|
|
|
|
b |
||||
|
|
|
|
|
|
f (x) 2 dx;K y |
|
x f (x) dx; |
||||
2 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
a |
|
|
|
|
a |
||||
K y |
|
;Yc |
K x |
|
|
|
||||||
|
m |
|
m |
|
|
|||||||
|
|
|
|
|
|
|
|
|||||
111
5. ПРИБЛИЖЕННОЕ ВЫЧИСЛЕНИЕ ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ
Если функция f(x) непрерывна на отрезке a, b
и
известна еѐ первообразная F(х), то определѐнный интеграл можно вычислить по формуле Ньютона-Лейбница:
b |
|
f (x)dx F (b) F (a) . |
(1) |
a
Однако во многих случаях первообразная функция F(x) не выражается через элементарные или находится слишком сложно. Поэтому возникает задача о приближѐнном вычислении интеграла.
Численные методы приближѐнного вычисления интеграла основаны на приближѐнном представлении определѐнного интеграла интегральной суммой I n или
величиной, близкой к I n : b
f (x)dx I n .
a
5.I. Формулы прямоугольников b
Пусть надо вычислить f (x)dx . Разобьѐм отрезок a, b
|
a |
|
|
на n равных частей длиной h |
|
b a |
точками |
|
n |
||
|
|
|
|
x0 a, x1, x2 ,... xn b . Вычислим значения функции f(x) в этих точках:
y0
f (a); y1
f ( x1 );... yk
f ( xk );... yn 1
f ( xn 1 ); yn
f (b)
.
Составим две суммы:
112
y
у0 |
|
|
|
|
М |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
уn |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
уn-1 |
|
|||
0 |
у0 |
|
|
|
|
|
у1 |
|
|
|
у2 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
х1 |
|
|
х2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
||||
|
х0=a |
|
|
|
... |
|
|
|
|
|
хn-1 |
|
xn=b |
||||||||||||
|
|
|
|
|
|
|
|
|
|
Рис. 50. |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
1 |
|
|
I1 |
h y0 |
y1 |
... |
|
yk ... |
yn 1 |
|
|
f (xk )h |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
1 |
|
|
|
I1 |
h y1 |
y2 |
... |
|
yk ... |
yn |
|
f (xk |
1)h . |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
0 |
|
|
|
Каждая из этих сумм является частным случаем интегральной суммы для функции f(x) и отрезка a, b : k xk 1. Поэтому
при n
эти суммы дают значение определѐнного интеграла (1). Если n достаточно велико (h мало), получаются суммы, мало отличающиеся от определѐнного интеграла. (площади ступенчатых фигур мало отличаются от площади криволинейной трапеции). Поэтому:
b
f (x)dx |
b a |
y |
|
y ... y |
|
... y |
|
(2) |
|
0 |
k |
n 1 |
|||||
|
n |
1 |
|
|
||||
|
|
|
|
|
|
|
||
a
113
b |
b a |
|
|
|
|
|
|
|
|
|
f (x)dx |
y |
y |
|
... y |
|
... |
y |
|
(3) |
|
|
2 |
k |
n |
|||||||
|
n |
1 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
Полученные |
формулы |
|
|
носят |
|
название |
формул |
|||
прямоугольников. |
|
|
|
|
|
|
|
|
|
|
Если функция f(x) возрастает на отрезке |
a, b , то формула (2) |
|||||||||
даѐт значение интеграла с недостатком, а (3) – с избытком.
Можно доказать, что если |
f (x) |
существует и ограничена на |
|||||||||||||
a, b , то погрешность |
формул |
(16) |
и |
(17) можно |
оценить |
||||||||||
|
|
|
|
|
|
b |
|
|
|
n 1 |
|
|
M1 |
(b a)2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
R |
|
|
f (x)dx |
h |
|
yk |
|
|||||
|
|
|
|
|
|
||||||||||
следующим образом: |
|
|
|
|
|
|
, |
||||||||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
a |
|
|
|
k 0 |
|
|
|
2h |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е. |
R |
порядка h. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(x) |
|
|
||||||||||||
Здесь M1 - наибольшее значение |
f |
на отрезке |
a, b . В |
||||||||||||
более точных способах вместо интегральных сумм рассматривают суммы другого вида.
5.2. Формула трапеций
Как и в предыдущем случае, делим отрезок a, b |
на n |
|||
равных частей: h |
b a |
. Каждую элементарную |
дугу |
|
n |
||||
|
|
|
||
заменим соответствующей хордой, в результате чего криволинейная трапеция заменяется суммой площадей элементарных прямоугольных трапеций. Геометрически очевидно, что такие трапеции более точно выражает искомую площадь, чем прямоугольники. Площадь каждой элементарной трапеции равна:
114
y
у0 |
|
|
М |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
у2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
у1 |
|
|
|
|
|
|
|
|
|
|
|
уk |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
0 |
|
|
у0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
х0=a |
х1 |
|
|
|
х2 |
|
... |
|
|
|
|
|
|
хk |
|
... |
|
хn-1 |
|
|
xn=b |
|
|
x |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 51. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
y0 y1 |
h; |
|
y1 |
y2 |
|
h;... |
yn 1 |
|
yn |
|
h . |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Складывая элементарные площади, получим : |
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
I3 |
h |
|
y0 |
|
|
|
y1 |
|
y2 ... |
|
yn |
1 |
|
|
yn |
|
|
|
|
|
|
|
|||||||||||||||
|
|
2 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
n |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
||
|
|
|
h |
|
|
|
|
yk |
|
|
h y0 |
|
|
yn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
k |
0 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
b |
|
|
|
b |
a y0 |
|
|
yn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
f (x)dx |
|
|
|
|
y1 |
|
y2 ... |
yn 1 |
. |
(4) |
||||||||||||||||||||||||||||
|
|
|
|
|
n |
|
|
|
|
|
|
|
2 |
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
a
Полученная формула называется формулой трапеций. Можно доказать, что если на a, b существует ограниченная производная f (x) , то ошибка оценивается формулой:
|
R |
|
M |
2 b |
a 3 |
|
R |
|
h |
2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
, т.е. |
|
|
|
. |
|
|
|||
|
|
|
12n2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
f (x) |
|
на a, b . |
|
Здесь М2 – наибольшее значение |
|
|
||||||||||||
|
|
|
|
|
115 |
|
|
|
||||||