27Интегрирование функций с разрывами и несобственные интегралы

До сих пор мы интересовались в первую очередь интегралами непрерывных функций. Именно для такого случая мы доказали формулу Ньютона — Лейбница: непрерывность подынтегральной функции является важным её условием. Что бывает, когда подынтегральная функция терпит разрывы? Мы не сможем дать полный ответ на этот вопрос, но обсудим некоторые важные частные случаи.

27.1Интегрирование функций с конечным числом скачков

27.1.1Скачки интегралу не помеха

Пусть функция имеет на отрезке конечное число ( штук) разрывов в точках , во всех остальных точках отрезка она непрерывна, и все разрывы являются скачками, то есть в них определены пределы справа и слева. Можно ли посчитать интеграл по от такой функции?

Если воспринимать задачу интегрирования геометрически, это не выглядит проблемой. Интеграл — это площадь, а площадь обладает свойством аддитивности — если разрезать фигуру, площадь от которой мы находим, на несколько непересекающихся кусочков, площадь всей фигуры должна равняться сумме площадей этих кусочков. Соответственно, можно разбить интеграл в сумму интегралов по промежуткам, на которых функция непрерывна (а значит интегрируема), и всё сложить. Именно так мы и сделаем! Но сперва докажем полезную лемму.

Лемма 1. Пусть функции и отличаются друг от друга лишь в конечном числе точек . Тогда их интегралы совпадают (и одновременно существуют или не существуют).

Доказательство. Изменения конечного числа значений не влияет на ограниченность, следовательно, функции и либо одновременно неограниченны, либо одновременно ограничены. В первом случае они обе неинтегрируемы и доказывать нечего. Рассмотрим второй случай, пусть обе функции ограничены некоторым числом (по модулю).

Возьмём произвольное размеченное разбиение , у которого границами отрезков разбиения являются точки и отмечены точки . Поскольку функции и различаются лишь в конечном числе точек, эти точки могут попасть лишь в конечное число отрезков разбиения. (В худшем случае точка может оказаться общим концом двух отрезков разбиения, значит, число затронутых отрезков не более чем в два раза больше числа точек различия.) Пусть затронуты отрезки с номерами , . Модуль разности между интегральными суммами для функций и равен

Модуль разности между значениями функций и не превосходит , поскольку сами значения ограничены по модулю числом . (Мы используем неравенство .) Длина отрезка разбиения не превосходит диаметра разбиения . Всего слагаемых штук, где . Значит разница между двумя интегральными суммами не превосходит и стремится к нулю при стремящемся к нулю. Значит, пределы интегральных сумм совпадают.

Теорема 1. Пусть функция имеет конечное число точек разрывов на отрезке и все разрывы являются скачками, то есть в каждой точке разрыва есть конечные односторонние пределы функции справа и слева. В остальных точках функция непрерывна. Тогда она интегрируема на отрезке .

Доказательство. Воспользовавшись свойством аддитивности (24.5), разобьем интеграл в сумму интегралов:
Рассмотрим любой из получившихся отрезков интегрирования, обозначим его концы через и . Функция непрерывна на , за исключением, быть может, концов отрезка. Но в концах существуют односторонние пределы (поскольку все точки разрывов — скачки), и следовательно для вычисления интеграла по функцию можно переопределить в и , таким образом, чтобы сделать её непрерывной на всём отрезке (как обычно, в концах требуется односторонняя непрерывность). По лемме, значение интегралов от этого не поменяется, равно как и факт интегрируемости. Но функция, непрерывная на отрезке, интегрируема на этом отрезке, и значит все интегралы в правой части равенства (27.1) существует и функция интегрируема на всём отрезке .

Пример 1. Пусть
Найдём
Возможно три варианта: , или . В первом случае наш интеграл — это площадь прямоугольника высотой и шириной , взятая со знаком «минус» (потому что значение функции отрицательно), то есть (см. анимацию на рис. 27.1). При можно переопределить функцию таким образом, чтобы в точке она принимала значение и тем самым сделать её непрерывной на всём отрезке интегрирования. Если же , наш интеграл разбивается в сумму двух интегралов: по отрезку и по отрезку . Первый равен , второй является площадью прямоугольника с шириной и высотой , то есть .
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

def nansign(x):
    sign = np.sign(x)
    return np.where(sign != 0, sign, np.nan)
fig, (ax1, ax2) = plt.subplots(2, 1, sharex=True)
camera = Camera(fig)

x = np.linspace(-3, 3, 201)

for t in np.linspace(-2, 3, 11):
    ax1.plot(x, nansign(x), label=r'$y=\mathrm{sign}(x)$', 
             color='gray')
    ax1.plot([0], [0], 'o', color='gray')
    ax1.plot([0, 0], [1, -1], 'o', color='gray',
        markerfacecolor='white')
    t_span_neg = np.linspace(-2, min(t, 0), 200)
    ax1.fill_between(t_span_neg, -np.ones_like(t_span_neg), color='red', 
                     alpha=0.5)
    if (t > 0):
        t_span_pos = np.linspace(0, t, 200)
        ax1.fill_between(t_span_pos, np.ones_like(t_span_pos), color='green', 
                         alpha=0.5)
    ax1.annotate('$t$', (t, 0), xytext=(t, 0.1))

    ob.center_spines(grid=False, minor_ticks=False, ax=ax1)
    ob.settle_axes(xmin=-3, xmax=3, ymin=-1.3, ymax=1.3, 
                   xlabel="x", ylabel="y", ax=ax1) 
    ax1.set_xticks([-2, -1, 0, 1, 2])
    ax1.set_yticks([-1, 1])


    t_span = np.linspace(-2, t)
    ax2.plot(t_span, -2+np.abs(t_span), color='C0')
    ob.center_spines(grid=False, minor_ticks=False, ax=ax2)
    ob.settle_axes(xmin=-3, xmax=3, ymin=-3.2, ymax=3.2, 
                   xlabel="t", ylabel="y", ax=ax2) 
    camera.snap()
animation = camera.animate()
Рис. 27.1: Интеграл от как функция верхнего предела. Сверху: подынтегральная функция и площадь; снизу: значение интеграла.
Итак,
Это выражение можно записать короче: .

Заметим, что во всех точках , производная совпадает с подынтегральной функцией (проверьте!):

Это неудивительно: первая часть формулы Ньютона — Лейбница выполняется для всех точек, в которых подынтегральная функция непрерывна, а разрыв у нас только в нуле. В нуле же функция не является дифференцируемой, но является непрерывной. Оказывается, так будет всегда.

27.1.2Интеграл с переменным верхним пределом непрерывен

Утверждение 1. Пусть функция интегрируема на отрезке . Тогда функция
непрерывна на отрезке .

Доказательство. Действительно, для любых точек
Иными словами, разность значений в точках и равна интегралу от по отрезку . Поскольку функция под интегралом интегрируема, она ограничена, и значит интеграл по маленькому отрезку от неё будет маленьким. Следовательно, значения функции в близких точках близки.

Аккуратное доказательство выглядит так. Поскольку функция ограничена, найдётся такое , что для всех

Интегрируя эти неравенства по отрезку (будем считать, что , противоположный случай разбирается аналогично), имеем:
Считая интегралы в левой и правой части, получаем:
В пределе при , левая и правая части стремятся к нулю, значит и интеграл посередине стремится к нулю. Но этот интеграл равен разности , следовательно эта разность стремится к нулю, и значит
то есть непрерывна в точке .

Теперь мы можем заключить, что интеграл от кусочно-непрерывной функции с конечным числом скачков, если рассматривать его как функцию от верхнего предела, является непрерывной функцией, дифференцируемой во всех точках, кроме точек разрывов ; в точках разрывов у интеграла точки излома и производная не существует (если только эти разрывы не являются устранимыми, устранимые разрывы интеграл «не замечает» в силу леммы).

27.2Несобственные интегралы

До сих пор мы рассматривали интегралы, которые отвечали площадям ограниченнух фигур — то есть таких фигур, которые можно поместить в достаточно большой прямоугольник. Действительно, в горизонтальном направлении мы были ограничены отрезком интегрирования, а в вертикальном — ограниченностью подынтегральной функции — мы всегда говорили, что раз функция не является ограниченной, она и не интегрируемая, и даже однажды это доказали.

Однако, можно представить себе фигуру, которая не будет помещаться ни в какой прямоугольник, но тем не менее будет иметь конечную площадь. Как бы нам её найти?

27.2.1Несобственные интегралы с бесконечными пределами интегрирования

Начнём с примера.

Пример 2. Пусть мы хотим найти такой интеграл:
Область интегрирования теперь — не отрезок, а бесконечный луч, см. рис. 27.2.
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

fig, ax1 = plt.subplots()

x = np.linspace(-1.2, 3, 201)
ax1.plot(x, 1 / x ** 2)

x_clip = np.linspace(1, 3, 201)
ax1.fill_between(x_clip, 1 / x_clip ** 2, alpha=0.5)

ob.center_spines(grid=False, minor_ticks=False, ax=ax1)
ob.settle_axes(xmin=-1.2, xmax=3, ymin=-0.7, ymax=3, 
               xlabel="x", ylabel="y", ax=ax1) 

ax1.set_xticks([1, 2])
ax1.set_yticks([])
Рис. 27.2: Интеграл по лучу
Использовать обычное определение нельзя — луч нельзя разбить в конечное число отрезков. Теоретически, можно было бы расширить понятие разбиения, разрешив разбиения на бесконечное число отрезков, но тогда в интегральных суммах было бы бесконечное число слагаемых, и надо было думать, как эти суммы определять…

В общем, мы сделаем иначе. Положим по определению, что интеграл (27.2) равен такому пределу:

То есть мы считаем интеграл до какой-то конечной границы , для разных получаем разные значения интеграла, и дальше устремляем к бесконечности, см. анимацию на рис. 27.3.

Это выглядит логичным: если мы будем всё сильнее увеличивать верхний предел интегрирования , мы будем захватывать всё большую и большую часть нашей неограниченной фигуры, так что в пределе ничего не останется не захваченным, и значит площадь захваченного куска должна стремиться к площади всей фигуры.

import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

fig, ax1 = plt.subplots()
camera = Camera(fig)
x = np.linspace(-1.2, 3, 201)

for t in np.linspace(1, 3, 11):
    ax1.plot(x, 1 / x ** 2, color='C0')

    x_clip = np.linspace(1, 3, 201)
    ax1.fill_between(x_clip, 1 / x_clip ** 2, alpha=0.5, color='C0')

    x_clip_t = np.linspace(1, t, 201)
    ax1.fill_between(x_clip_t, 1 / x_clip_t ** 2, alpha=0.5,
        color='C1')

    ob.center_spines(grid=False, minor_ticks=False, ax=ax1)
    ob.settle_axes(xmin=-1.2, xmax=3, ymin=-0.7, ymax=3, 
                   xlabel="x", ylabel=None, ax=ax1) 
    ax1.annotate('$y$', (0.03, 3 - 0.15))

    ax1.annotate('$t$', (t, 0), xytext=(t + 0.03, -0.25))
    ax1.set_xticks([1, 2])
    ax1.set_yticks([])
    camera.snap()
animation = camera.animate()
Рис. 27.3: Несобственный интеграл как предел
Попробуем найти, что получится. Заметим, что первообразной является . Воспользуемся формулой Ньютона — Лейбница:
Получилось! Искомая площадь равна .

Определение 1. Пусть функция интегрируема на любом отрезке , . Интеграл
называется несобственным интеграл по бесконечному промежутку. По определению он равен

Всегда ли предел будет существовать? Конечно, нет — легко придумать пример, когда не будет. Скажем,

Но тут история понятная — подынтегральная функция является константой, наша фигура — «бесконечный прямоугольник», неудивительно, что у него бесконечная площадь. Но что если подынтегральная функция стремится к нулю. Может ли быть так, чтобы интеграл по лучу не существовал? Давайте рассмотрим ещё один пример.

Пример 3. Найдём
Действуя так же, как и раньше, запишем этот интеграл в виде предела. Напомним, что первообразная при положительных — это . Имеем:

Определение 2. Если предел в определении несобственного интеграла существует, говорят, что интеграл сходится, в противном случае он расходится.

Итак, интеграл по одному и тому же лучу от функции сходится, а от функции расходится. Почему так? Казалось бы, они обе стремятся к нулю при , и на глаз их графики очень похожи. Однако, мало стремиться к нулю, нужно делать это ещё достаточно быстро. Функция стремится к нулю медленнее, чем , и этого стремления не хватает, чтобы площадь оказалась конечной.

Является ли требование стремления к нулю необходимым для сходимости интеграла? Оказывается, нет.

Упражнение 1. Придумайте пример функции, которая не стремится к нулю при , но интеграл от которой по лучу сходится.

Определение 3. Легко определить интеграл по лучу :
Интеграл по всей прямой определяется так:
где промежуточную точку можно выбирать произвольным образом, значение интеграла от этого не изменится (почему?).

27.2.2Признак сравнения

Часто бывает важно даже не найти значение несобственного интеграла, а просто показать, что он сходится. В этом случае полезной оказывается следующая теорема.

Теорема 2. Пусть функции и интегрируемы на любом отрезке , . Пусть также для всякого
и интеграл
сходится. Тогда
тоже сходится.

Доказательство. Утверждение выглядит естественным. Фигура, соответствующая интегралу (27.5), содержится внутри фигуры, соответствующей интегралу (27.4). Логично предположить, что если площадь первой конечна, то и второй тоже.

Формальное доказательство выглядит так. Введём обозначения:

Заметим, что обе функции и неубывают на : это следует из того факта, что их производные, равные подынтегральным функциям, неотрицательны. Также интегрируя неравенства (27.3) по отрезку , имеем:
По предположению, имеет предел при . Обозначим его через . Все значения неубывающей функции не превосходят предела: действительно, если для какого-то , , то найдётся такое , что (можно в качестве взять половину разности , тогда — середина отрезка ), и значит для всех выполняется такое же неравенство (в силу неубывания ), а значит предел не меньше , то есть строго больше — противоречие. Значит для всех :
Итак, функция неубывает и ограничена. Можно легко адаптировать теорему Вейерштрасса для функций (сделайте это!) и получить, что существует предел
то есть интеграл (27.5) сходится.

Пример 4. Докажем, что
сходится.

Действительно, разобьём этот интеграл в сумму двух:

С первым слагаемым никаких проблем нет: это обычный (собственный) интеграл Римана, подынтегральная функция непрерывна, значит, интеграл существует. Рассмотрим второе слагаемое. Заметим, что при , и
Таким образом, если интеграл
сходится, то и наш интеграл сходится. Но легко проинтегрировать явно: первообразная равна . Таким образом,
Мы нашли этот интеграл явно, значит, он сходится, и значит наш исходный интеграл тоже сходится.

В том виде, в котором мы его сформулировали, признак сравнения работает только с неотрицательными подынтегральными функциями.

Вопрос 1. Покажите, что условие неотрицательности здесь существенно. Приведите пример функций и , для которых условие выполняется для всех , сходится, а интеграл расходится.

Следующее утверждение позволяет расширить признак сравнения на случай знакопеременных подынтегральных выражений:

Утверждение 2. Если сходится интеграл
то сходится и интеграл

Мы не будем сейчас доказывать это утверждение. Однако, из него мгновенно следует, что если существует такая функция , что для всех , и сходится, то и тоже сходится.

27.2.3Несобственные интегралы и вертикальные асимптоты

Итак, нам удалось расширить понятие интеграла на случай бесконечных промежутков интегрирования. Может быть, и с неограниченными функциями удастся что-то сделать?

До сих пор мы много раз повторяли мантру «раз функция неограничена, значит, она и не интегрируема» — и даже доказали утверждение на эту тему. Однако, значит ли это, что соответствующие площади на самом деле не определены? Нам уже удалось найти площади некоторых неограниченных фигур, как в примере 2 выше. Может быть, и с площадями под графиками неограниченных функций ситуация может быть аналогичной — и наш реузльтат о неинтегрируемости вызван просто неудачным определением интеграла? Можно ли модифицировать это определение каким-то образом, чтобы хотя бы некоторые неограниченные функции стали интегрируемыми? Давайте попробуем.

Пример 5. Попробуем найти интеграл
Подынтегральная функция стремится к бесконечности при , и значит не является ограниченной — стало быть интеграл в обычном смысле не определён. Однако, давайте воспользуемся тем же приёмом, которым мы воспользовались для нахождения интеграла по бесконечному промежутку. Положим по определению
Для всякого интеграл под знаком предела определён. Более того: его можно найти явно — подынтегральная функция может быть записана в виде , а все степенные функции мы имеем интегрировать. Первообразная имеет вид (проверьте дифференцированием!) и значит
Интеграл сошёлся, ура!

Определение 4. Пусть функция имеет вертикальную асимптоту в точке и для всякого функция интегрируема по Риману на отрезке . Несобственный интеграл
определяется следующим образом:

Как обычно, предел может существовать (в этом случае говорят, что интеграл сходится), или не существовать (значит, интеграл расходится).

Пример 6. Сходится ли интеграл
Мы уже считали интеграл от этой функции, но по другой области — по бесконечному промежутку — и он сошёлся. Может быть, и в этот раз нам повезёт?
Увы, этот интеграл расходится.

Почему такая разница с функцией ? Дело в том, что стремится к нулю гораздо медленнее, чем , и значит стремится к бесконечности гораздо медленнее, чем . Поэтому площадь под графиком оказывается конечной, а под графиком — бесконечной.

В определении 4 вертикальная асимптота находится на левой границе отрезка интегрирования. Конечно, можно с тем же успехом определить несобственный интеграл для функции с вертикальной асимптотой на правой границе области отрезка интегрирования.

Определение 5. Пусть функция имеет вертикальную асимптоту в точке и для всякого функция интегрируема по Риману на отрезке . Несобственный интеграл
определяется следующим образом:

А что если точка бесконечного разрыва попала куда-то внутрь нашего отрезка? Ничего страшного — просто разобьем интеграл на два.

Определение 6. Пусть функция имеет вертикальную асимптоту в точке и интегрируема на любой отрезке , и на любом отрезке , . Тогда
где интегралы в правой части — несобственные. Если хотя бы один из этих двух интегралов расходится, интеграл в левой части тоже считается расходящимся.

Пример 7. Что вы можете сказать про интеграл
Очень хочется сказать, что он равен нулю: ведь это интеграл от нечётной функции по симметричному отрезку, и значит положительная и отрицательные площади должны взаимно уничтожиться. Однако, это не так: по определению,
и оба интеграла в правой части расходятся, значит и исходный интеграл расходится. (Подробнее о том, почему это правильный подход, будет обсуждаться на семинаре.)

27.2.4Несобственные интегралы с несколькими особенностями

Каждое из определений несобственных интегралов, которые мы обсуждали, рассматривает только одну «особенность» — бесконечный разрыв или бесконечную границу интегрирования. Если у интересующего нас интеграла несколько особенностей, его следует разбить в сумму интегралов, у каждого из которых лишь одна особенность, и воспользоваться нашими определениями.

Пример 8. Пусть у функции бесконечные разрывы в точках , и , и мы хотим посчитать от неё интеграл по всей числовой прямой. Его следует представить в таком виде:
Каждый из получившихся интегралов — один из несобственных интегралов, которые мы обсуждали. Промежуточные точки , , , выбраны произвольным образом — можно выбрать какие-то другие точки, лежаще на тех же промежутках, от этого ничего не изменилось бы. Если хотя бы один из интегралов в правой части расходится, значит, весь исходный интеграл тоже расходится.

Замечание 1. Теперь если вы видите запись
вам нужно посмотреть, как ведёт себя функция на отрезке . Если она ограничена на этом отрезке, перед нами обычный интеграл Римана (собственный интеграл). Если же у неё есть точки бесконечных разрывов, значит, это несобственный интеграл, и его нужно искать, представляя интеграл как сумму несобственных интегралов разобранных нами типов.

27.3Заключение

Мы вышли за пределы уютного и безопасного мира, в котором все подыинтегральные функции непрерывны. И тут много всего интересного! Во-первых, если функция имеет разрывы типа «скачок», и их конечное число, никаких проблем для интегрирования это не представляет. Правда, интеграл, если думать о нём как о функции от его верхнего предела, будет иметь изломы в точках разрывов подынтегральной функции — но по крайней мере сам останется непрерывным. Случай бесконечных разрывов подынтегральной функции представляет большую опасность: несобственные интегралы, которые мы придумали специально для этого случая, существуют не всегда, иногда они сходится, а иногда расходится. Но тут уж ничего не поделать — площадь под графиком неограниченной функции может быть бесконечной. Мы также обсудили интегралы по бесконечным промежуткам — особенно важный для нас случай, поскольку такие интегралы сплошь и рядом встречаются в теории вероятностей, которую вы будете проходить совсем скоро.

На этом наш краткий экскурс в теорию интегрирования заканчивается. Конечно, мы успели обсудить не всё, но фундамент заложен. Время идти дальше!