Лосев-КлячинМиклюков МА в КИ
.pdf§10. Приближенные методы вычисления определенного интеграла |
211 |
В качестве отмеченных точек выберем точки ξi = (xi+xi+1)/2, т.е. середины отрезков [xi, xi+i]. Тогда справедливо
|
n−1 f(ξ |
)∆x |
= |
b − a |
n−1 f |
|
xi + xi+1 |
, |
|
|
||||||
|
i=0 |
i |
i |
|
|
n |
i=0 |
2 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и приходим к формуле прямоугольников |
|
|
|
|
||||||||||||
|
b |
|
|
b − a |
n−1 f |
xi + xi+1 |
|
|
|
|
||||||
f(x) dx = |
+ R |
. |
|
|||||||||||||
|
|
|
||||||||||||||
|
|
|
|
n |
|
i=0 |
|
|
|
2 |
|
n |
|
(1) |
||
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь Rn – остаточный член и для практических применений формулы (1) важно знать его оценки.
Покажем, что если f имеет на [a, b] непрерывную произ-
водную второго порядка, то |
|
|
|
|
||||
| |
R |
n| ≤ |
(b − a)3 |
M, |
M = max |
f (x) |
. |
(2) |
|
24n2 |
[a,b] | |
| |
|
Оценим сначала интеграл
h
f(x) dx,
−h
считая, что f дважды непрерывно дифференцируема на отрезке [−h, h].
Пусть F – первообразная по отношению к f. Тогда
h
f(x) dx = F (h) − F (−h).
−h
Так как F трижды непрерывно дифференцируема на [−h, h] (что следует из того, что f – дважды непрерывно дифференцируема, и теоремы о дифференцируемости интеграла по переменному верхнему пределу), то, пользуясь формулой Тейлора с остаточным членом в форме Лагранжа, получаем
F (x) = F (0) + F (0)x + 12F (0)x2 + 16F (ξx)x3,
где ξx – некоторая точка между точками 0 и x. Заметим, что
F (0) = f(0), F (0) = f (0), F (ξx) = f (ξx).
214 Глава 10. Свойства интегрируемых функций
Тем самым, мы получаем приближенную формулу |
n−1 |
|||||
a |
≈ n |
|
2 |
1 · · · |
|
|
b |
b − a |
|
f(a) + f(b) |
|
|
|
f(x)dx |
|
|
|
+ f(x ) + |
+ f(x ) , |
называемую формулой трапеций.
10.3.Формула парабол (формула Симпсона)
Пусть f – непрерывная на [−h, h] функция. Заменим при-
ближенно f квадратичной параболой y = a x2 + b x + c так, чтобы y(−h) = f(−h), y(0) = f(0) и y(h) = f(h).
|
Найдем коэффициенты a, b и c. Имеем |
|
|
|
|
|
|
|
||||||||||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
b = f(h) − f(−h) |
|
|
|||||||||||
|
ah − bh + c = f(−h) |
|
|
|
|
|
|
|
|
|
|
|
|
2h |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
2 |
|
c |
= |
f(0) |
|
|
|
c = f(0) |
− |
. |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
a = |
|
|
|
|
|
− |
2h2 |
|||||||
|
ah |
+ bh + c |
= |
f(h), |
|
|
|
|
|
|
|
f( |
|
h) + f(h) |
|
2f(0) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Выполним некоторые |
вычисления: |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
h |
|
|
|
|
|
|
|
3 |
|
h |
|
|
2 |
|
|
h |
|
|
|
|
||||
|
|
h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
ax |
|
− |
|
bx |
|
|
|
− |
|
|
|
|
||||||
|
− |
|
|
|
|
|
|
|
|
|
|
h + 2 |
|
|
|
|
|
|
+ cx|−h h = |
|
||||||
|
|
(ax2 + bx + c) dx = 3 |
|
|
|
|
|
|
|
|
h |
|
||||||||||||||
|
|
a |
|
3 |
|
3 |
b |
2 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
= |
|
|
(h |
|
+ h ) + |
|
(h |
− h ) + c(h + h) = |
|
|
||||||||||||||
|
|
3 |
|
|
2 |
|
|
§10. Приближенные методы вычисления определенного интеграла |
215 |
|||||||||||||||
|
= |
f(−h) + f(h) − 2f(0) |
|
|
2h3 |
|
+ f(0)2h = |
|
||||||||
|
|
|
|
|
|
|
||||||||||
|
3 |
3 |
2h2 |
|
3 |
|
− |
|
|
|||||||
3 |
3 |
|
|
|
|
|
||||||||||
= h |
f(h) |
+ |
4f(0) |
+ |
f(−h) |
|
= |
h |
(f( |
|
h) + 4f(0) + f(h)) . |
|||||
|
|
|
|
|
|
Заменяя теперь интеграл от f эквивалентным интегралом от квадратного трехчлена y = ax2 + bx + c, получаем
h
f(x)dx ≈ h3 [f(−h) + 4f(0) + f(h)].
−h
В случае, если f непрерывна на [a, b], данное соотношение влечет
b |
|
b − a |
|
|
a + b |
|
|
|
f(x)dx |
|
f(a) + 4f |
+ f(b) . |
|
||||
|
|
|
|
|
||||
a |
≈ |
6 |
|
2 |
|
(3) |
Разобьем [a, b] на n равных частей точками
a = x0 < x2 < · · · < x2n = b и обозначим через x2i+1 середину [x2i, x2i+2]. Тогда
|
b |
n−1 |
x2i+2 |
|
|
|
|
|
|
|
f(x) dx = i=0 |
f(x) dx |
||
a |
|
x2i |
|
и, применяя к каждому из интегралов в сумме приближенную формулу (3), находим
|
|
b |
|
|
n−1 |
|
b − a |
|
|
|
|
|
|
|
|
|
|
|
|
f(x) dx |
|
(f(x ) + 4f(x |
|
|
) + f(x |
|
)) , |
|
|||||||
|
|
|
|
|
|
|
|
||||||||||
|
a |
|
|
|
|
|
2i |
|
2i+1 |
|
|
2i+2 |
|
|
|
||
|
|
|
|
|
≈ i=0 6n |
|
|
|
|
|
|
||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
n−1 |
|
|
|
n−1 |
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
≈ |
b − a |
|
|
|
|
|
|
|
|
|
|
|
|||
|
f(x) dx |
|
f(a) + f(b) + 2 f(x |
|
|
) + 4 |
f(x |
|
) . |
||||||||
|
|
|
|
6n |
|
|
|
i=1 |
|
2i |
i=0 |
|
2i+1 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Это есть известная формула Симпсона.
218 |
Глава 11. Приложения определенного интеграла |
через каждую точку квадрата (кривые Пеано).
Схематический рисунок кривой Пеано
УПРАЖНЕНИЕ 1. Сравнить приведенный рисунок со знаменитой картиной Малевича.2
УПРАЖНЕНИЕ 2. Выяснить, не могла ли являться картина Малевича "Черный квадрат" попыткой изобразить кривую Пеано?
1.1.Уравнения касательной и нормали к кривой
Пусть задана дифференцируемая кривая (или дуга) γ посредством вектор-функции
r(t) = (ϕ(t), ψ(t)),
определенной на промежутке a, b . Фиксируем некоторое значение параметра t0 a, b . Придадим t0 приращение ∆t так, чтобы параметр t0 + ∆t a, b . Тогда вектор-функция получит приращение
r(t0 + ∆t) − r(t0) = (ϕ(t0 + ∆t) − ϕ(t0), ψ(t0 + ∆t) − ψ(t0)).
Предельный вектор
lim r(t0 + ∆t) − r(t0) = r (t0) = (ϕ (t0), ψ (t0))
∆t→0 ∆t
называется производной вектор-функции r(t) в точке t0.
2Малевич Казимир Северинович (1878–1935 гг.) – художник. Учился в Училище живописи, ваяния и зодчества (1904–1905 гг.). Участник выставок художественных группировок "Бубновый валет", "Ослиный хвост" и др. В 1919-1922 гг. преподавал в Народной художественной школе "нового революционного образца" в Витебске. В 1923-1927 гг. – директор Ленинградского гос. института художественной культуры. Один из основателей абстрактного искусства, отошёл от отражения реальных вещей и явлений, отказался от конкретной сюжетной содержательности произведений, трактовал предметную форму как комбинации контрастных по цвету геометрических элементов ("Чёрный квадрат", 1913 г.; "Полёт аэроплана", 1915 г.; "Красный квадрат", 1917 г.) Автор книг: "Искусство, церковь, фабрика", "Супрематизм", "От кубизма к супрематизму".
220 Глава 11. Приложения определенного интеграла
ЛЕММА 1.1. Предположим, что кривая (или дуга) γ задана посредством пары непрерывных функций x = ϕ(t), y = ψ(t), t [a, b].
Если кривая (или дуга) γ является простой, то мелкость µ(P ) → 0 тогда и только тогда, когда мелкость
µ(T ) → 0.
Доказательство. Утверждение о том, что µ(T ) → 0 влечет µ(P ) → 0, следует из теоремы Кантора о равномерной непрерывности для ϕ и ψ.
Для доказательства обратного утверждения предположим противное. В таком случае найдется последовательность пар
точек {pkpk}, лежащих на γ и таких, что diam pkpk |
→ 0, одна- |
||||||||||||||||||||||||||||||
ко |
для последовательности их |
|
прообразов |
|
|
|
|
выполнено |
|||||||||||||||||||||||
| |
t |
− |
t |
|
ε |
> 0 (k = 1, 2, . . .) |
|
|
|
ε |
– |
некоторое фикси- |
|||||||||||||||||||
|
k |
k| ≥ |
|
|
|
|
|
|
, где |
|
|
|
|
|
|
" |
|
|
|
|
|||||||||||
рованное число. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
Так как последовательность {tk} лежит на отрезке [a, b], то |
|||||||||||||||||||||||||||||
она ограничена. По теореме Больцано-Вейерштрасса из нее |
|||||||||||||||||||||||||||||||
можно извлечь подпоследовательность {tkm } → t |
|
[a, b]. |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
} в свою |
|||
Из соответствующей ей подпоследовательности { km |
|||||||||||||||||||||||||||||||
очередь также извлечем подпоследовательность |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
{ |
t |
|
|
t |
|
[a, b] |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
kmn } → |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
Мы имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
p |
|
ϕ t |
|
, ψ |
t |
|
p |
|
|
ϕ |
|
t |
) |
, ψ |
( |
t |
|
, |
|
||||||||
и |
|
|
kmn = ( ( |
kmn ) |
|
( kmn )) → |
|
= ( ( |
|
|
|
|
|
)) |
|
|
|||||||||||||||
|
p |
|
ϕ |
t |
, ψ t |
p |
|
|
ϕ |
t |
) |
, ψ |
( |
t |
)) |
. |
|||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
kmn |
= ( ( |
|
kmn ) ( kmn )) → |
|
|
= ( ( |
|
|
|
|
|
|
|
||||||||||||||
Но diam p p |
|
|
0, а потому p |
совпадает с p . |
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
kmn |
kmn → |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
С другой стороны, |t − t | ≥ ε и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(ϕ(t ), ψ(t )) = p , (ϕ(t ), ψ(t )) = p .
Тем самым, мы нашли два различных значения параметра t = t , переходящих в одну и ту же точку p = p на γ,
причем в случае замкнутости кривой γ найденная точка не может совпадать с концевой точкой A = B. Это противоречит предположению, что γ – простая кривая (дуга).
УПРАЖНЕНИЕ 3. Показать, что лемма 1.1 будет неверна, если мелкость определить формулой
µ˜(P ) = max |pi − pi+1|.
0≤i≤n−1