Квант
Основные статьи

Замечательные пределы, порождаемые классическими средними

Крейн М. Г., Нудельман А. А.

13

Замечательные предыды, порождаемые классическими средними

Для двух положительных чисел aa и bb, как известно, средним арифметическим называется число a+b2\frac{a+b}{2}, средним геометрическим — число ab\sqrt{ab}. Режке, чем эти средние, встречается среднее гармоничское 2aba+b\frac{2ab}{a+b} (ему равна, например, средняя скорость автоматчика, проехавшего первую половину пути со скоростью aa, вторую — со скоростью bb). Каждое из этих средних расположено между числами aa и bb (убедитесь в этом).

Легко проверяется справедливоность herавенств

2aba+b<ab<a+b2\frac{2ab}{a+b} < \sqrt{ab} < \frac{a+b}{2}

для любых положительных a≠b∗a \neq b^*.

Мы заимемся здесь следующим вопросом. Пусть данные два положительных числа aa и bb, a<ba < b. Вычислив какую-нибудь пару средних этих чисел, получим числа a1a_1 и b1b_1. Затем для чисел a1a_1 и b1b_1 вычислим те же средние — получим числа a2a_2 и b2b_2. С ними повторим ту же процедуру и так далее. В результате получим последовательности (an)(a_n) и (bn)(b_n).

Например, взяв среднее геометрическое и среднее арифметическое и отправлясь от чисел 1 и 3, получаем

a1=3≈1.732050808,b1=2;a_1 = \sqrt{3} \approx 1.732050808, \quad b_1 = 2;
a2≈1.861209718,b2≈1.866025404;a_2 \approx 1.861209718, \quad b_2 \approx 1.866025404;
a3≈1.863616006,b3≈1.863617561;a_3 \approx 1.863616006, \quad b_3 \approx 1.863617561;
a4≈1.863616784,b4≈1.863616784;a_4 \approx 1.863616784, \quad b_4 \approx 1.863616784;
и т. д.

Мы видим, что в нашем примере последовательности (an)(a_n) и (bn)(b_n) очень быстро сбижаются. Всетда ли так будет? Оказывается, подобные последовательности всегда имеют общий предел — это доказывается сравнительно легко. Но как найти значение этого предела?

Арифметико-гармоничское среднее

Мы начнем со случай, когда выбранная пара средних — это среднее гармоничское и среднее арифметическое; таким образом, члены последовательностей (an)(a_n) и (bn)(b_n) определяются формулами

an+1=2anbnan+bn,bn+1=an+bn2a_{n+1} = \frac{2a_nb_n}{a_n+b_n}, \quad b_{n+1} = \frac{a_n+b_n}{2}

(n=0,1,2,… ;a0=a,b0=bn=0, 1, 2, \dots; a_0=a, b_0=b).

Задача о нахождения их общего предела представляет на II командиной олимпиде школа города Омска («Квант», 1979, № 12).

Из неравенств (1) следует

a<an<an+1<bn+1<bna < a_n < a_{n+1} < b_{n+1} < b_n

то есть последовательность (an)(a_n) возратает «навстречу» убивающей последовательности (bn)(b_n).

Таким образом, обе последовательности монотонны и ограничены; следовательно, по теореме Вейерштраса («Алгебра и начала анализа 9», п. 32) они имеют предель α=lim⁡n→∞an,β=lim⁡n→∞bn\alpha = \lim_{n \to \infty} a_n, \beta = \lim_{n \to \infty} b_n. Перейдя к пределу в одном из равенств (2), получим

β=lim⁡n→∞bn+1=lim⁡n→∞an+bn2=\beta = \lim_{n \to \infty} b_{n+1} = \lim_{n \to \infty} \frac{a_n+b_n}{2} =
=12(lim⁡n→∞an+lim⁡n→∞bn)=12(α+β),= \frac{1}{2} \left( \lim_{n \to \infty} a_n + \lim_{n \to \infty} b_n \right) = \frac{1}{2} (\alpha + \beta),

откуда α=β\alpha = \beta: последовательности (an)(a_n), (bn)(b_n) имеют общий предел. Этот предел называется арифметико-гармоничским средним чисел aa и bb. Найдем его. Из (2)

14

an+1⋅bn+1=an⋅bn=…=a1b1=ab,a_{n+1} \cdot b_{n+1} = a_n \cdot b_n = \ldots = a_1 b_1 = ab,
позтому
a2=(lim⁡n→∞an)2=lim⁡n→∞an×lim⁡n→∞an=lim⁡n→∞an×bn=ab,a^2 = \left( \lim_{n \to \infty} a_n \right)^2 = \lim_{n \to \infty} a_n \times \lim_{n \to \infty} a_n = \lim_{n \to \infty} a_n \times b_n = ab,
откуда
α=ab=β,\alpha = \sqrt{ab} = \beta,
(арифметико-гариончское среднее совпадает со средним геометрическим).

Упражнение 1. Докажите, что bn−an<b−a2nb_n - a_n < \frac{b-a}{2^n}.

Таким образом, последовательности (an)(a_n) и (bn)(b_n) достаточно быстро сходятся к ab\sqrt{ab}. Поэтому они могут оказаться полезными для приближенного вычисления квадратных корней. Для вычисления c\sqrt{c} нужно начиная последовательности (an)(a_n), (bn)(b_n) с таких чисел aa и bb, что c=abc = ab (например, a=1,b=ca = 1, b = c), причем процесс будет сходиться тем быстрее, чем меньше разность между этим множительями (например, для вычисления 56\sqrt{56} лучше быть не a=1a = 1, b=56b = 56, а a=7a = 7, b=8b = 8). Легко убедиться, что последовательности (an)(a_n) и (bn)(b_n) получаются по формулам
bn+1=12(bn+cbn),an+1=cbn.b_{n+1} = \frac{1}{2} \left( b_n + \frac{c}{b_n} \right), a_{n+1} = \frac{c}{b_n}.
Для иллюстрации вычисления 12\sqrt{12}, положив b=4b = 4. Получаем b1=12(4+124)=3,5b_1 = \frac{1}{2} \left( 4 + \frac{12}{4} \right) = 3,5, b2≈3,464285715b_2 \approx 3,464285715, b3≈3,464101620b_3 \approx 3,464101620, b4≈3,464101615b_4 \approx 3,464101615 и далее все знаки стабилизируются:
12≈3,464101615.\sqrt{12} \approx 3,464101615.

Арифметико-геометрическое среднее Четырнадцатилетный Карл Фридрих Гаусс*) обнаружил на числовых примерах, что при вычислении последовательности (an)(a_n) и (bn)(b_n) с помощью геометрических иарифметических средних:
an+1=anbn,bn+1=an+bn2a_{n+1} = \sqrt{a_n b_n}, b_{n+1} = \frac{a_n + b_n}{2}

(n=0,1,2,…;a0=a,b0=bn = 0, 1, 2, \ldots; a_0 = a, b_0 = b) эти последовательности очень быстро сходятся с ближайшими.

*) Биографию «Короля математиков» и интересные сведения о его научной деятельности читатель найдет в статье С. Гиндикна «Карл Фридрих Гаусс» («Квант», 1977, № 8).

Упражнение 2. Докажите существование и совпадение пределов последовательности (3).

Этот общий предел называется арифметико-геометрическим средним чисел aa и bb и обозначается через μ(a,b)\mu(a, b). Найти явное выражение для μ(a,b)\mu(a, b) через aa и bb очень не просто. впервые оно было получено Гауссом в результате необсячайно оструемых и виртуозных расчудений и преобразований, использующих свойства так называемых эллиптических интерголов. Приведем это выражение без доказательства*):

μ(a,b)=π2∫0π/2dxa2sin⁡2x+b2cos⁡2x.\mu(a, b) = \frac{\pi}{2} \int_0^{\pi/2} \frac{dx}{\sqrt{a^2 \sin^2 x + b^2 \cos^2 x}}.

Геометрическо-гариончское среднее

Если строить последовательности (an)(a_n) и (bn)(b_n) с помощью средних грамонических и средних геометрических:
an+1=2anbnan+bn,bn+1=anbna_{n+1} = \frac{2a_n b_n}{a_n + b_n}, b_{n+1} = \sqrt{a_n b_n}

(n=0,1,2,…;a0=a,b0=bn = 0, 1, 2, \ldots; a_0 = a, b_0 = b),

то в этом случае легко доказать, что они сходятся к общему пределу. Назовем его геометрическо-гариончческим чисел aa и bb и обозначим через ν(a,b)\nu(a, b). Однако ничего существенного по сравнению с последовательностями (3) здесь нет, так как из (4) усматриваем, что
1an+1=1an+1bn,1bn+1=1an⋅1bn.\frac{1}{a_{n+1}} = \frac{1}{a_n} + \frac{1}{b_n}, \frac{1}{b_{n+1}} = \sqrt{\frac{1}{a_n} \cdot \frac{1}{b_n}}.
Позтому
1ν(a,b)=μ(1b,1a),\frac{1}{\nu(a, b)} = \mu\left(\frac{1}{b}, \frac{1}{a}\right),
или же, в силу формулы Гауса.

*) Если это выражение использовать «в обратном направлении», получится быстро способ приближения вычисления интеграла, стоящего в знаменателе: он равен π2μ(a,b)\frac{\pi}{2\mu(a, b)}, а среднее μ(a,b)\mu(a, b) вычисляется приближенно с помощью последовательностей (an)(a_n) и (bn)(b_n). В наш век ЭВМ эти вычисления можно проводится с огромной скоростью и точностью.