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

Об одной формуле Христиана Гюйгенса

Вавилов В. В.

9

Об одной формуле Христиана Гюйгена

Кандидат физико-математических наук
В. В. ВАВИЛОВ

Длина окружности радиуса 1/2 может быть определена как предел последовательности длини перimetров правильных многоугольников, вписанных в данную окружность при неограниченном увеличении числа его стороны. Определение таким образом число обозначется символом л. Аналогично можно определить число л и как предел перimetров правильных описанных около этой окружности многоугольников (рис. 1). Точнее, если через pnp_n обозначить перimetр вписанного, а через qnq_n перimetр описанного правильного nn-угольника, то имеют место неравенства
pn<π<qn,n⩾3.p_n < \pi < q_n, \quad n \geqslant 3.
Когда nn возрастает, обе последовательности pnp_n и qnq_n монотно приближаются к числу л (рис. 2); при перехождение от nn к n+1n+1 мы получаем меньше ошибку в приближениях, которое числа pnp_n и qnq_n дают для числа л.
В своей значении соинения «Измерение кругга» Архимед (287—212 гг. до н. э.) дал следующее приближение для числа л:
31071<π<31070,3\frac{10}{71} < \pi < 3\frac{10}{70},
то есть
3,14084507<π<3,14285714.3,14084507 < \pi < 3,14285714.
Этот результат Архимед получил при помощи вычисления перimetра правильного вписанного 96-угольника.
Многоековая история числа л связана с многими именами известных математиков. Однако из последовательей

10

Именно это основное неравенство и позволило Гюйгенсу пойти дальше Архимеда.

Ниже мы докажем приведенные результаты Архимеда и Х. Гюйгена и рассмотрим некоторые смежные вопросы.

Доказательство основного неравенства

Доказательство неравенства (1) мы разобем на несколько шагов. Левую часть этого неравенства будем считать известной и поэтому докажем только справедливость второго неравенства.

1°. Если в некоторых сегмент круга ACB приscan равноведенный треугольник ABC и в точках A, B и C проведены касательные к окружности так, как показано на рисунке 4, то

SABC<2SKLM.S_{ABC} < 2S_{KLM}.

(2) Действенный, из соответствий ∣KL∣>∣LC∣|KL| > |LC| и ∣LC∣=∣AL∣|LC| = |AL| следует ∣AK∣<2∣KL∣|AK| < 2|KL|. Так как треугольники KLM и KAB подобны, их площади относятся как квадраты соответствующих сторонах, справедливоно,

SKAB/SKLM=∣KA∣2/∣KL∣2<4,S_{KAB} / S_{KLM} = |KA|^2 / |KL|^2 < 4,

откуда SKAB<4SKLMS_{KAB} < 4S_{KLM}.

С другой стороны, равноведенные треугольники KAB и ABC имеют общее основание AB, и поэтому их площади относятся так же, как их высоты, или (по теореме Фалеса) как длины соответствующих сторонах. Следовательно,

SABCSKAB=∣AL∣∣AK∣,\frac{S_{ABC}}{S_{KAB}} = \frac{|AL|}{|AK|},

и, так как ∣AL∣<(1/2)∣AK∣|AL| < (1/2)|AK|, получаем SABC<(1/2)SKABS_{ABC} < (1/2)S_{KAB}.

Отсюда и из предыдущего неравенства для площадей получаем (2).

2°. Если обозначить через ScermS_{cerm} площадь сегмента ACB (см. рис. 4), то

Scerm≤23SABK.S_{cerm} \leq \frac{2}{3} S_{ABK}.

(3)

Для доказательства этого неравенства воспользуемся следующим сообщением. Обозначим (см. рис. 5) через S0S_0 и A0A_0 площади треугольников ABC и KLM соответственно. Расмотрим теперьдугу AC, разделим ее точкой C1C_1 пополам и построим точно таким же образом треугольники площадей S2S_2 и A2A_2. Точка C3C_3 — серединадуги AC2^2 — порождает аналогичным образом треугольники площадей S3S_3 и A3A_3. Продолжая этот процесс, мы получим две последовательности {Sn}\{S_n\} и {An}\{A_n\} положительных чисел — площадей соответствующих треугольников. По доказанному в 1° — см. (2), — имеем следующее:

S0<2A0,S1<2A1,S2<2A2,S3<2A3,…,Sn<2An.S_0 < 2A_0, \quad S_1 < 2A_1, \quad S_2 < 2A_2, \quad S_3 < 2A_3, \dots, \quad S_n < 2A_n.

Ясно, что (рис. 6)

S0+2S1+22S2+23S3+⋯+2nSn+⋯=ScermS_0 + 2S_1 + 2^2S_2 + 2^3S_3 + \dots + 2^nS_n + \dots = S_{cerm}

11

и
A0+2A1+22A2+23A3+⋯+2nAn+⋯=SABK−Sсегм.A_0 + 2A_1 + 2^2A_2 + 2^3A_3 + \dots + 2^nA_n + \dots = S_{ABK} - S_{\text{сегм}}.
Из приведенной цепочки неравенства и двух последних соответствий следует, что *)
Sсегм≤2(SABK−Sсегм),S_{\text{сегм}} \leq 2(S_{ABK} - S_{\text{сегм}}),
то есть утверждение (3).
3°. Следующая геометрическая лемма X. Гойгена представляет и самостоятельный интерес.
Лемма. На рисунке 7: О — центр кругга, прямые AE и CD — касательные к окружности, AB — произвольная секция, пересекающая прямую CD, l_AC — длина дуги AC. Тогда справедливо следующее неравенство:
lAC≤23∣CD∣+13∣AF∣.l_{\text{AC}} \leq \frac{2}{3} |CD| + \frac{1}{3} |AF|.

По свойству касательных, проведенных из одной точки к окружности, |AE| = |EC|; следовательно, треугольник ACE равноведенный и углы EAC и ECA равны. Отсюда следует, что треугольник AED также равноведенный (так как CAD=90∘CAD = 90^\circ, DAE=90∘−EACDAE = 90^\circ - EAC, ADE=90∘−ECAADE = 90^\circ - ECA) и тем самым |AE| = |ED| = |CE|. Из (3) следует, что Sсегм≤(2/3)SAECS_{\text{сегм}} \leq (2/3)S_{AEC}, где SсегмS_{\text{сегм}} — площадь сегмента кругга, закрашен-

*) Нетрудно доказать, что если последовательности {an}\{a_n\} и {bn}\{b_n\} имеют предель и an<bna_n < b_n, n=1,2,…n=1, 2, \dots, то lim⁡n→∞an≤lim⁡n→∞bn\lim_{n \to \infty} a_n \leq \lim_{n \to \infty} b_n.

ного на рисунке 7. Прибавив к обеим частиям этого неравенства площадь SAOCS_{AOC} треугольника OAC, мы получим
Sсегт≤23SAEC+SOAC,S_{\text{сегт}} \leq \frac{2}{3} S_{AEC} + S_{OAC},
где SсегтS_{\text{сегт}} — площадь сектора OAC, откуда
Sсегт≤23SOAC+13SOAC,S_{\text{сегт}} \leq \frac{2}{3} S_{OAC} + \frac{1}{3} S_{OAC},
где SOACS_{OAC} — площадь четырехугольника OAEC. Так как
SOAC=2⋅SOEC=2⋅12∣OC∣⋅12∣CD∣=12∣OC∣⋅∣DC∣,SOAC=12∣OC∣⋅∣AF∣,Sсегт=12∣OC∣⋅lAC,S_{OAC} = 2 \cdot S_{OEC} = 2 \cdot \frac{1}{2} |OC| \cdot \frac{1}{2} |CD| = \frac{1}{2} |OC| \cdot |DC|, S_{OAC} = \frac{1}{2} |OC| \cdot |AF|, S_{\text{сегт}} = \frac{1}{2} |OC| \cdot l_{AC},
из последнего неравенства следует нужное соответствие (4).
4°. Для завершения доказательства неравенства (1) сделаем следующие построение (рис. 8): пусть AA1AA_1 — сторона правильно вписанного nn-угольника, а LL1LL_1 — сторона правильно описанного многоугольника с таким же числом сторон; проведем прямую AB и из точки A опустим перпендикуляр AK на касательную LL_1.

Так как AOC=2ABCAOC = 2ABC (одни из углов центральный, а второй опирается на эту же дуг), мы видим, что AD — bicсектриса в треугольнике LAK, и по теореме о биссектрисе имеем ∣AL∣/∣LD∣=∣AK∣/∣KD∣|AL| / |LD| = |AK| / |KD|.

Но (AK)⊥(KL)(AK) \perp (KL), поэтому ∣AK∣<∣AL∣|AK| < |AL| и из предыдущей пропорции получаем ∣KD∣<∣LD∣|KD| < |LD|. Так как ∣AF∣=∣KC∣|AF| = |KC|, получаем ∣KD∣=∣CD∣|KD| = |CD|

12

| nn | pnp_n | qnq_n |
| :— | :— | :— |
| 3 | 2,59807621 | 5,19615242 |
| 6 | 3,00000000 | 3,46410161 |
| 12 | 3,10582854 | 3,21539030 |
| 24 | 3,13262861 | 3,15965994 |
| 48 | 3,13935020 | 3,14608621 |
| 96 | 3,14103195 | 3,14271459 |
| 192 | 3,14145247 | 3,14187304 |
| 384 | 3,14155760 | 3,14166274 |
| 768 | 3,14158389 | 3,14161017 |
| 1536 | 3,14159046 | 3,14159703 |
| 3072 | 3,14159210 | 3,14159374 |

−∣AF∣-|AF|. Кроме того, ∣LD∣=∣CL∣−|LD|=|CL|-
−∣CD∣-|CD| и, следовательно,
∣CD∣−∣AF∣<∣CL∣−∣CD∣,|CD|-|AF|<|CL|-|CD|,
2∣CD∣<∣AF∣+∣CL∣,2|CD|<|AF|+|CL|,
2∣CD∣+∣AF∣<2∣AF∣+∣CL∣2|CD|+|AF|<2|AF|+|CL|.

Используя теперь (4), получаем
lAC<23∣AF∣+13∣CL∣.l_{AC} < \frac{2}{3} |AF| + \frac{1}{3} |CL|.
Аналогично можно написать, что
lA1C<23∣A1F∣+13∣CL1∣l_{A1C} < \frac{2}{3} |A1F| + \frac{1}{3} |CL1|
или, наконец,
lAC+lA1C<23∣AA1∣+13∣LL1∣.l_{AC} + l_{A1C} < \frac{2}{3} |AA1| + \frac{1}{3} |LL1|.

Умножив полученное неравенство на nn, мы и получим неравенство (1).

Дальнейшее исследование основного неравенства

Нами установлено, что число лежит в первой трети интеграла pnp_n, qnq_n при всех n⩾3n \geqslant 3. Для того чтобы уточнить расположение числа л в этом новом интеграле, рассмотрим отношение длины интегралов pnp_n, π\pi при pp, qnq_n [. Вычисления показывают *] (см. таблицы 1, 2), что это отношение длины, то есть значения дробей
(qn−π)/(π−pn),n=3,6,12,24,(q_n-\pi)/(\pi-p_n), n=3, 6, 12, 24,
достаточно близки к 2. На основании этих вычислений мы с большой степенью уверенности можем предположить, что в действительности имеет место соответствие

*) Для желающих проверки вычисления на калькуляторе приведем приближение равенство: π≈3,14159265\pi \approx 3,14159265.

lim⁡n→∞qn−ππ−pn=2.\lim_{n \to \infty} \frac{q_n-\pi}{\pi-p_n} = 2. (5)

Для доказательства соответствия (5) заметим, что (рис. 9))

pn=n⋅sin⁡πn,qn=n⋅tgπn,n⩾3,p_n=n \cdot \sin \frac{\pi}{n}, q_n=n \cdot \tg \frac{\pi}{n}, n \geqslant 3,
и, следовательно,
qn−ππ−pn=1cos⁡πn⋅nsin⁡πn−πcos⁡πnπ−nsin⁡πn.\frac{q_n-\pi}{\pi-p_n} = \frac{1}{\cos \frac{\pi}{n}} \cdot \frac{n \sin \frac{\pi}{n}-\pi \cos \frac{\pi}{n}}{\pi-n \sin \frac{\pi}{n}}.
Для анализа полученного выражения нам понадобятся неравенства *)
x−x36<sin⁡x<x−x36+x5120,x>0.x-\frac{x^3}{6} < \sin x < x-\frac{x^3}{6} + \frac{x^5}{120}, x > 0. (6)

Чтобы доказать левое неравенство в (6), положим
f(x)=sin⁡x−x+x36.f(x)=\sin x-x+\frac{x^3}{6}.

Тогда имеем
g1(x)=f′(x)=cos⁡x−1+x22,g_1(x)=f'(x)=\cos x-1+\frac{x^2}{2},
g2(x)=g1′(x)=−sin⁡x+x.g_2(x)=g'_1(x)=-\sin x+x.

*) Неравенство (6) значительное улучшает известное неравенство sin x<xx<x при x>0x>0; см. по этому поводу: «Квант», № 5, 1979, c. 5.

13

Так как sin x < x при x > 0, получим g2(x)>0g_2(x) > 0 при x > 0. Тем самым функция g1(x)g_1(x) возвращает при x > 0. Но g1(0)=0g_1(0) = 0 и, следовательно, g1(x)>0g_1(x) > 0 при x > 0. Функция g1(x)g_1(x) является произвольной для функции f(x)f(x), для которого также f(0)=0f(0) = 0. Но при x > 0 имеем g1(x)>0g_1(x) > 0, поэтому функция f(x)f(x) также возвращает и, следовательно, принимает только положительные значения, то есть f(x)>0f(x) > 0 при x > 0, что и утверждалось.

Аналогично установливается правила част неравенства (6), а также неравенства (доказите их самостояльно!)

1−x22<cos⁡x<1−x22+x424,x>0.(7)1 - \frac{x^2}{2} < \cos x < 1 - \frac{x^2}{2} + \frac{x^4}{24}, \quad x > 0. \quad (7)

Из неравенств (6) и (7) вытекают следующие приближенные формулы:

sin⁡πn≈πn−π36n3,\sin \frac{\pi}{n} \approx \frac{\pi}{n} - \frac{\pi^3}{6n^3},

cos⁡πn≈1−π22n2.\cos \frac{\pi}{n} \approx 1 - \frac{\pi^2}{2n^2}.

Следовательно,

qn−ππ−pn≈2⋅(1−π22n2)−1,\frac{q_n - \pi}{\pi - p_n} \approx 2 \cdot \left(1 - \frac{\pi^2}{2n^2}\right)^{-1},

что и завершает доказательство соответствия (5), так как π22n2\frac{\pi^2}{2n^2} стремится к нулю при n→∞n \to \infty.

Равенство (5) позволяет сделать следующий качественный вывод: число π\pi, находясь при любом n⩾3n \geqslant 3 в интервале ]pn,23pn+13qn[]p_n, \frac{2}{3}p_n + \frac{1}{3}q_n[, при всех достаточно больших значениях nn близке к правому концу этого интервала, чем к левому.

Формула Гойгенаса и ее эффективность

Архимед использовал для вычисления числа π\pi приближенную формулу

π≈pn,n⩾3.\pi \approx p_n, \quad n \geqslant 3.

Гойгенс в своей работе, в частности, получил другую приближенную формулу

π≈23pn+13qn,n⩾3,\pi \approx \frac{2}{3}p_n + \frac{1}{3}q_n, \quad n \geqslant 3,

то есть взял в качестве приближения для числа π\pi правую часть неравенства (1).

Бóльшую эффективность формулы Гойгенаса по сравнению с формулой Архимеда можно обнаружить неосредственными вычислениями на микрокалькуляторе, сравив результаты таких вычислений, сведенных в таблицах 1, 2 и 3 (цветом на этих таблицах выделены точные цифры числа π\pi). Отметим, что проведится также вычисления — увлекательная и непростая задача.

Можно сравнить эффективность формул Архимеда и Гойгенаса другим методом, не произвольной конкретных вычислений для pnp_n и qnq_n. Можно использовать так называемые аппирорные оценки для точности этих формул, то есть такие неравенства, которые позволяют оценить разности π−pn\pi - p_n и 23pn+13qn−π\frac{2}{3}p_n + \frac{1}{3}q_n - \pi для всех n⩾3n \geqslant 3 заранее, без вычислений периметров. Ценность таких аппирорных оценок состоит так же в том, что они позволяют в точности этих оценки числа nn, для которых нужно найти точность будет достигуть. Используя неравенства (6) и (7) и произвольные простые аппирорные вычисления, нетрудно получить (мы оставляем это читателю в качестве упражнения!) следующие аппирорные оценки:

5n2<π−pn<6n2,n⩾3,(8)\frac{5}{n^2} < \pi - p_n < \frac{6}{n^2}, \quad n \geqslant 3, \quad (8)

и

23pn+13qn−π<64n4,n⩾3.(9)\frac{2}{3}p_n + \frac{1}{3}q_n - \pi < \frac{64}{n^4}, \quad n \geqslant 3. \quad (9)

Непосредственно из этих неравенств

14

замечаем, что «скорость», с которой разность 23pn+13qn−π\frac{2}{3} p_n + \frac{1}{3} q_n - \pi стремится к нулю, больше «скорости», с которой стремится к нулю разность π−pn\pi - p_n при n→∞n \to \infty.

Для сравнения этих «скоростей» найдем те значения nn, которые нам наверняка обеспечат точность вычисления числа π\pi, например, равную 10−210^{-2}. Так как π−pn<6n2\pi - p_n < \frac{6}{n^2}, для вычисления по формуле Архимеда с этой точностью достаточно взять nn таким, чтобы 6/n2<10−26/n^2 < 10^{-2}, то есть n≥25n \geq 25. Используя же формулу Гойгенаса, нам достаточно выбрать nn таким, чтобы 64/n4<10−264/n^4 < 10^{-2}, то есть n≥9n \geq 9.

Как видно из разобранного примера, а также из неравенств (8) и (9), формула Гойгенаса быстрое приводит к заданной точности.

Две задачи

Мы помещаем здесь материал, который может стать началом небольшого математического исследования.

Задача 1. Пусть ACAC — хорда, ∣BB1∣=2R|BB_1| = 2R — диаметр окружности. OO — центр круга, (BB1)⊥(AC)(BB_1) \perp (AC), BB — середина дуги ACAC, EE и FF — середины дуг ABAB и BCBC соответственно (рис. 10). Доказите, что

a) ∣AB∣2=∣BD∣⋅∣BB1∣|AB|^2 = |BD| \cdot |BB_1|, ∣BE∣2=∣BG∣×∣BB1∣|BE|^2 = |BG| \times |BB_1|;

b) ∣BD∣<4⋅∣BG∣|BD| < 4 \cdot |BG|;

c) ∣AC∣<2⋅∣EF∣|AC| < 2 \cdot |EF|;

d) SABC<4(SABE+SBFC)S_{ABC} < 4(S_{ABE} + S_{BFC});

e) πR2>43S2n−13Sn\pi R^2 > \frac{4}{3} S_{2n} - \frac{1}{3} S_n;

где SnS_n — площадь правильного в окружность nn-угольника (S2nS_{2n} — соответствующее 2n2n-угольника);

ж) 2πR>43p2n−13pn2\pi R > \frac{4}{3} p_{2n} - \frac{1}{3} p_n,

где pnp_n — перimetр правильного в окружность nn-угольника (p2np_{2n} — соответствующее 2n2n-угольника);

задача 2. Докажите, что все нижелейные предыдущего существуют и найдете их:

a) n2(π−pn)n^2(\pi - p_n);

b) n4(23pn+13qn−π)n^4\left(\frac{2}{3} p_n + \frac{1}{3} q_n - \pi\right);

b) lim⁡n→∞pn2−pn−1pn+12pn−pn−1−pn+1\lim_{n \to \infty} \frac{p_n^2 - p_{n-1}p_{n+1}}{2p_n - p_{n-1} - p_{n+1}};

r) n4(π−43p2n+13pn)n^4\left(\pi - \frac{4}{3} p_{2n} + \frac{1}{3} p_n\right);

d) lim⁡n→∞πR2−Sn(n)Sq(n)−πR2\lim_{n \to \infty} \frac{\pi R^2 - S_n(n)}{S_q(n) - \pi R^2},

где RR — радиус круга; Sp(n),Sq(n)S_p(n), S_q(n) — площади правильных виссанных и описанных многоугольников.

В заключение отмечим, что наблюдение Х. Гойгена о том, что число nn, равное длине окружности радиуса 1/21/2, находится в два раза больше к перimetр виссанного в эту окружность правильного nn-угольника, чем к перimetр описанного вокруг этой окружности правильного nn-угольника, было не так давно обобщено на случай произвольных выпуклых фигур, ограниченных гладками (без изломов) кривыми. Именно, имеет место следующие утверждения.

Пусть LL обозначает длину выпуклой замкнутой кривой без изломов, PnP_n — перimetр nn-угольника с максимальным перimetром, виссанного в данной кривой и QnQ_n — перimetр описанного вокруг данной кривой с минимальным перimetром (такие существуют!). Тогда справедливо соответствие утверждений

lim⁡n→∞Qn−LL−Pn=2.\lim_{n \to \infty} \frac{Q_n - L}{L - P_n} = 2.

Доказательство этого утверждения и многое другое интересные результаты, связанные с ним, можно найти в книге [4].

Литература

  1. Рудно Ф. О кадаратуре круза (Архимед, Гойгенс, Лежандр, Ламберт), Матезис, 1911.
  2. Веселовский И. Н. ХристIAN Гойгенс. — М.: Учредгиз, 1959.
  3. Кыпман Ф. История числа π\pi, — М.: Наука, 1971.
  4. Тот Л. Ф. Расположения на плоскости, на сфере и в пространстве, — М.: Физматгиз, 1958.