25Формула Ньютона — Лейбница

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

25.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

fig, (ax1, ax2) = plt.subplots(2, 1, figsize=(6, 6), sharex=True)
x = np.linspace(-2, 2, 200)
ax1.plot(x, x ** 2, label='$y=x^2$')
ob.center_spines(grid=False, minor_ticks=False, ax=ax1)
ob.settle_axes(xmin=-1.7, xmax=1.7, ymin=-1.2, ymax=2.8, 
               xlabel="x", ylabel="y", ax=ax1) 
ax1.legend()
ax1.set_xticks([])
ax1.set_yticks([])

ax2.plot(x, x ** 3 / 3, label='$y=x^3/3$', color='C1')
ax2.plot(x, x ** 3 / 3 + 0.5, label='$y=x^3/3+1/2$', color='C2')
ax2.plot(x, x ** 3 / 3 + 1, label='$y=x^3/3+1$', color='C3')
ax2.plot(x, x ** 3 / 3 - 0.5, label='$y=x^3/3-1/2$', color='C4')
ax2.set_yticks([])
ax2.legend()
ob.center_spines(grid=False, minor_ticks=False, ax=ax2)
ob.settle_axes(xmin=-1.7, xmax=1.7, ymin=-2.2, ymax=2.8, 
               xlabel="x", ylabel="y", ax=ax2) 
Рис. 25.1: Различные первообразные функции отличаются друг от друга добавлением константы, что соответствует вертикальному сдвигу графика

Определение 1. Пусть функция непрерывна на отрезке . Её первообразной (англ. antiderivative) называется такая функция , что для всякого ,
где производные на концах отрезка понимаются в смысле односторонних производных.

Пример 1. Любая функция вида является первообразной для функции . Любая функция вида является первообразной к функции . Любая функция вида является первообразной к функции .

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

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

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

25.2Связь интегралов и первообразных

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

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

Пример 2. Пусть и . Тогда при интеграл
равен площади прямоугольника шириной и высотой . Значит, .

Пример 3. Пусть и . Тогда при интеграл
равен площади трапеции высотой и основаниями и . Из геометрии мы знаем, что площадь такой трапеции равна произведению высоты на полусумму оснований, то есть

25.2.2Формулировка и доказательство формулы Ньютона — Лейбница

Теорема 1. (Формула Ньютона — Лейбница.) Пусть непрерывна на отрезке и
Тогда справедливы два утверждения:
  1. является первообразной , то есть для всех ,
  2. Для любой первообразной функции справедливо утверждение:

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

Приращение функции. Запишем определение производной для функции в некоторой точке (на концах нужно взять соответствующую одностороннюю производную):

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

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

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

По определению минимума и максимума, для всех ,
Неравенства можно интегрировать (см. утверждение 3 из предыдуещй лекции). Проинтегрируем (25.10) по отрезку . Получим:
В самой левой и самой правой частях неравенства подынтегральные функции — константы (они зависят от , но интегрирование происходит по переменной , и с точки зрения этого интегрирования, они константы). Соответствующие интегралы — это просто площади прямоугольников, они равны и соответственно. Значит
Геометрически это соответствует тому, что мы ограничили нашу узкую полоску двумя прямоугольниками: она содержится целиком внутри прямоугольника, высота которого равна , а прямоугольник с высотой содержится внутри неё. Значит, её площадь находится между площадями этих двух прямоугольников.

Предельный переход. Разделим неравенство (25.11) на (мы сейчас предполагаем, что , обратный случай рассматривается аналогично):

Посередине написано выражение, стоящее под знаком предела в определении производной функции (см. (25.5) с учётом (25.6)). Согласно (25.8),
и . По теореме о двух милиционерах, из этого следует, что при . В силу непрерывности функции по теореме о пределе сложной функции, отсюда следует, что
при . Аналогичное утверждение справедливо и для :
Применяя теорему о двух милиционерах к неравенствам (25.12) при , имеем:
Выражение посередине является производной , и в силу полученных неравенство, , как и ожидалось.

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

Осталось найти эту константу. Подставим в это равенство . В левой части будет — то есть интеграл от по отрезку от до . Он равен нулю (см. (24.6)). Справа будет . Таким образом,
то есть . Значит для любого :
и значит
Формула Ньютона — Лейбница доказана.

Замечание 2. У нас есть некоторый произвол в выборе функции — она может быть любой первообразной для . Однако, этот произвол не влияет на результат: все первообразные отличаются друг от друга на константу, и если к добавить эту константу, то в правой части формулы (25.4) она появится дважды с разными знаками и значение не поменяется.

25.2.3Нахождение интегралов с помощью первообразных

Формулу Ньютона — Лейбница также называют фундаментальной теоремой анализа (англ. fundamental theorem of calculus), поскольку она связывает два главных раздела этой науки — дифференциальное и интегральное исчисление. Наиболее популярным её применением — по крайней мере, в учебных курсах типа нашего — является нахождение интегралов с помощью первообразных. Приведём несколько примеров.

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

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

Пример 6. Найдём
Мы знаем, что первообразной для является . Значит

Замечание 3. При вычислении значений интегралов с помощью формулы Ньютона — Лейбница часто приходится вычислять разность значений одной и той же функции в двух разных точках. Чтобы сократить запись, используют такое обозначение:
Например, решение последнего примера обычно записывается так:

25.3Заключение

Мы доказали формулу Ньютона — Лейбница, и обнаружили, что задача отыскания площади сводится к задаче отыскания первообразной — то есть такой функции, чья производная равна подыинтегральной функции. Как находить первообразные? До сих пор мы это делали методом угадывания, и с совсем простыми функциями (степенными, тригонометрическими, экспонентой) этот фокус проходил. Но что делать, если функция, которую нужно проинтегрировать, устроена сложнее? Плохая новость: универсальных методов отыскания первообразных не существует, и, более того, первообразная вообще может не записываться в виде формулы, содержащей функции, которые мы проходили до сих пор. Хорошая новость: иногда эта задача всё-таки решается, и в следующей лекции мы обсудим некоторые методы интегрирования.