Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лосев-КлячинМиклюков МА в КИ

.pdf
Скачиваний:
864
Добавлен:
27.05.2015
Размер:
3.66 Mб
Скачать

§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).

212

 

Глава 10. Свойства интегрируемых функций

Теперь имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

F (h) − F (−h) = F (0) + f(0)h +

 

f (0)h2

+

 

 

f (ξh)h3

2

6

 

1

 

 

 

 

1

 

 

 

 

 

 

− F (0) − f(0)h +

 

 

f (0)h2

 

 

f (ξ˜h)h3

=

 

2

6

 

 

= 2f(0)h + f (ξh) + f (ξ˜h)

h3

 

 

 

 

,

 

 

 

˜

6

 

 

 

где ξh расположена между точками −h и 0.

 

 

 

 

 

 

Таким образом, для всякого отрезка [xi, xi+1] выполнено

xi+1

 

 

 

 

 

 

 

 

 

 

f(x) dx = 2f

xi + xi+1

 

xi+1 xi

+

 

2

2

 

xi

 

 

+[f (ξi) + f (ξ˜i)]

(xi+1 − xi)3

.

48

 

 

˜

(Здесь ξi и ξi — некоторые точки, расположенные между серединой [xi, xi+1] и точками xi+1 и xi соответственно.)

Тем самым, мы приходим к следующей формуле

 

b

n−1 f

 

 

f(x) dx =

b − a

xi + xi+1

 

 

 

a

 

n

 

 

2

 

 

i=0

(b − a)3 + 24n3

n−1

˜

i

f (ξi) + f (ξi)

.

=0

2

 

 

 

Остаточный член Rn здесь имеет вид

 

 

3 n−1

˜

Rn =

(b − a)

 

f (ξi) + f (ξi)

.

24n3

=0

2

 

 

 

 

i

 

 

Поэтому, предполагая, что |f (x)| ≤ M для всех x [a, b], находим

R

(b − a)3

n−1

M =

(b − a)3

M,

 

| n| ≤

24n3

i

=0

 

24n2

 

 

 

 

 

что и требовалось доказать.

10.2.Формула трапеций

Предположим, что f – непрерывная на [a, b] функция и T

разбиение отрезка

[a, b]

на n равных частей точками

a = x0 < x1 < · · ·

< xn

= b. На каждом из элементарных

§10. Приближенные методы вычисления определенного интеграла

213

отрезков [xi, xi+1] рассмотрим линейную функцию, совпада-

ющую с f в xi и xi+1.

 

трапеции,

 

расположенной над

Площадь элементарной

 

[xi, xi+1], дается формулой

 

 

 

 

 

 

 

f(xi) + f(xi+1)

(x

i+1

x

) =

b − a

 

f(xi) + f(xi+1)

,

2

n

2

 

i

 

 

а площадь Str(T ) фигуры, состоящей из элементарных трапеций, равна

 

 

Str(T ) =

b − a

n−1

f(xi) + f(xi+1)

=

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

n

2

 

 

 

 

 

 

 

=0

 

 

 

 

 

=

b − a

 

f(a) + f(b)

+ f(x1) +

· · ·

+ f(xn

1) .

 

 

 

n

 

2

 

 

 

 

 

Заметим, что s(T ) ≤ Str(T ) ≤ S(T ), где s(T ) и S(T ) – нижняя и верхняя суммы Дарбу, соответствующие разбиению T .

Так как функция f непрерывна, то, как было доказано в теореме об интегрируемости непрерывной функции, величины

s(T ) и S(T ) стремятся к интегралу

b

 

 

 

f(x) dx,

откуда следует

 

a

 

 

 

 

 

 

 

 

 

 

 

 

b

lim

S

tr

(T ) =

f(x) dx.

µ(T )0

 

 

 

a

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это есть известная формула Симпсона.

Глава 11

Приложения определенного интеграла

§1. Кривые и дуги

Определим некоторые геометрические понятия, связанные с теорией кривых на плоскости. Более строгое изложение данных вопросов будет дано позже и после того, как мы приобретем некоторый опыт работы с этими понятиями.

Предположим, что в плоскости R2 задана декартова прямоугольная система координат xOy. Пусть x = ϕ(t), y = ψ(t)

– пара функций, заданных на промежутке a, b . В этом случае говорят, что на a, b задана вектор-функция

r(t) = (ϕ(t), ψ(t)).

Вектор-функция непрерывна (дифференцируема), если каждая из функций ϕ, ψ непрерывна (дифференцируема).

При изменении параметра t от a до b точка (ϕ(t), ψ(t)) описывает некоторую кривую γ в плоскости (x, y). Кривая γ называется замкнутой, если ее начало и конец совпвдают, т.е.

(ϕ(a), ψ(a)) = (ϕ(b), ψ(b)).

Кривая γ называется разомкнутой кривой (или дугой), если ее начало (ϕ(a), ψ(a)) не совпадает с концом (ϕ(b), ψ(b)).

§1. Кривые и дуги

217

Может случиться, что через одну и ту же точку на плоскости кривая γ проходит несколько раз, т.е. нескольким различным значениям параметра t1, t2, . . . соответствует одна единственная точка (x0, y0) плоскости. Эта точка называется кратной точкой кривой γ.

Если кривая (или дуга) γ не имеет кратных точек, то она называется простой кривой (простой дугой).

Аналогично, простой замкнутой кривой называется замкнутая кривая, не имеющая точек пересечения кроме начала и конца.

Если простая дуга γ является графиком некоторой1 функции f, то говорят, что γ задана в непараметрическом виде.

Существуют примеры непрерывных кривых, проходящих

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 г.) Автор книг: "Искусство, церковь, фабрика", "Супрематизм", "От кубизма к супрематизму".

§1. Кривые и дуги

219

 

Уравнение прямой,

проходящей через точки

 

(t0) и

r

 

(t0 + ∆t), имеет вид

 

 

 

r

 

 

 

x ϕ(t0) = ϕ(t0 + ∆t) ϕ(t0), y − ψ(t0) ψ(t0 + ∆t) − ψ(t0)

или

(x−ϕ(t0)) (ψ(t0 +∆t)−ψ(t0)) = (y−ψ(t0)) (ϕ(t0 +∆t)−ϕ(t0)).

Разделим обе части равенства на ∆t и устремим ∆t → 0. Приходим к уравнению касательной к кривой γ в точке (ϕ(t0), ψ(t0)), имеющему вид

(x − ϕ(t0)) ψ (t0) = (y − ψ(t0)) ϕ (t0).

Уравнение нормали к γ в точке (ϕ(t0), ψ(t0)) записывается по обычным правилам

(x − ϕ(t0)) ϕ (t0) = (y − ψ(t0)) ψ (t0).

1.2.Длина дуги

Пусть r = r(t), t a, b , – непрерывная вектор-функция, определяющая некоторую простую дугу или кривую γ, замкнутую или разомкнутую. При изменении параметра t от a до b точка r(t) пробегает вдоль γ от одного ее конца A до другого B, т.е. на γ определено направление.

Если

T= {a = t0 ≤ t1 ≤ . . . ≤ tn = b}

разбиение a, b , то точки

A = p0 = r(a), p1 = r(t1), . . . , pn = r(b) = B

образуют разбиение P дуги γ. Соединив их последовательно одну за другой прямолинейными отрезками pipi+1, мы получим ломаную линию p0p1 . . . pn.

Каждый из отрезков pipi+1 является хордой для некото-

рой, стягиваемой им дуги pipi+1. Обозначим через diam pipi+1 диаметр этой дуги, т.е.

diam pipi+1 = sup{p , p pipi+1 : |p − p |}.

Здесь символом |p − p | обозначено обычное евклидово расстояние между точками p , p на плоскости.

Пусть

µ(P ) = max diam pipi+1

0≤i≤n−1

– мелкость разбиения P дуги (или кривой) γ.

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