Среднее арифметико-геометрическое
(
арифметико-геометрическое среднее
,
АГС
) — величина, определяющаяся для двух величин
a
{\displaystyle a}
и
b
{\displaystyle b}
как предел последовательности
{
a
N
}
{\displaystyle \{a_{N}\}}
,
{
b
N
}
{\displaystyle \{b_{N}\}}
, где:
a
0
=
a
b
0
=
b
{\displaystyle a_{0}=a\quad \quad \quad \quad \quad \quad \quad \quad b_{0}=b}
a
1
=
a
0
+
b
0
2
b
1
=
a
0
b
0
{\displaystyle a_{1}={\frac {a_{0}+b_{0}}{2}}\quad \quad \quad \quad b_{1}={\sqrt {a_{0}b_{0}}}}
a
2
=
a
1
+
b
1
2
b
2
=
a
1
b
1
{\displaystyle a_{2}={\frac {a_{1}+b_{1}}{2}}\quad \quad \quad \quad b_{2}={\sqrt {a_{1}b_{1}}}}
…
a
N
=
a
N
−
1
+
b
N
−
1
2
b
N
=
a
N
−
1
b
N
−
1
{\displaystyle a_{N}={\frac {a_{N-1}+b_{N-1}}{2}}\quad \quad b_{N}={\sqrt {a_{N-1}b_{N-1}}}}
имеют при
N
→
+
∞
{\displaystyle N\to +\infty }
один и тот же предел:
lim
N
→
∞
a
N
=
lim
N
→
∞
b
N
=
M
(
a
,
b
)
{\displaystyle \lim _{N\to \infty }a_{N}=\lim _{N\to \infty }b_{N}=M(a,b)}
.
АГС может быть применено для быстрого вычисления точного периода
математического маятника
.
Часто используется сокращение
M
(
x
)
=
M
(
x
,
1
)
{\displaystyle M(x)=M(x,1)}
. В частности
M
(
x
,
y
)
=
y
M
(
x
/
y
,
1
)
=
y
M
(
x
/
y
)
{\displaystyle M(x,y)=yM(x/y,1)=yM(x/y)}
.
Модифицированное арифметико-геометрическое среднее
(
МАГС
) двух величин
x
{\displaystyle x}
и
y
{\displaystyle y}
— (общий) предел (убывающей) последовательности
{
x
n
}
n
=
1
∞
{\displaystyle \{x_{n}\}_{n=1}^{\infty }}
и (возрастающей) последовательности
{
y
n
}
n
=
1
∞
{\displaystyle \{y_{n}\}_{n=1}^{\infty }}
, где
x
0
=
x
{\displaystyle x_{0}=x}
,
y
0
=
y
{\displaystyle y_{0}=y}
и
z
0
=
0
{\displaystyle z_{0}=0}
.
x
n
+
1
=
x
n
+
y
n
2
{\displaystyle x_{n+1}={\frac {x_{n}+y_{n}}{2}}}
y
n
+
1
=
z
n
+
(
x
n
−
z
n
)
(
y
n
−
z
n
)
{\displaystyle y_{n+1}=z_{n}+{\sqrt {(x_{n}-z_{n})(y_{n}-z_{n})}}}
z
n
+
1
=
z
n
−
(
x
n
−
z
n
)
(
y
n
−
z
n
)
{\displaystyle z_{n+1}=z_{n}-{\sqrt {(x_{n}-z_{n})(y_{n}-z_{n})}}}
МАГС может быть применено для быстрого вычисления длины нити в линейном параллельном поле сил отталкивания.
МАГС выразимо посредством АГС
[
как?
]
, такое опосредованное вычисление МАГС предпочтительно при вычислении длины периметра
эллипса
L
{\displaystyle L}
с полуосями
a
{\displaystyle a}
и
b
{\displaystyle b}
:
L
=
2
π
N
(
a
2
;
b
2
)
M
(
a
;
b
)
,
{\displaystyle L={\frac {2\pi N(a^{2};b^{2})}{M(a;b)}},}
где
M
(
x
;
y
)
{\displaystyle M(x;y)}
— АГС чисел
x
{\displaystyle x}
и
y
{\displaystyle y}
, а
N
(
x
;
y
)
{\displaystyle N(x;y)}
— МАГС чисел
x
{\displaystyle x}
и
y
{\displaystyle y}
. Тем самым, такая формула выражает метод Гаусса, с
квадратичной сходимостью
, для вычисления полного эллиптического интеграла второго рода.
Приложения
С использованием АГС и МАГС можно вычислять значения некоторых
трансцендентных функций
и
числа
π
{\displaystyle \pi }
.
Например, по формуле Гаусса — Саламина
:
π
=
2
M
(
1
;
2
)
2
1
−
∑
j
=
1
∞
2
j
c
j
2
,
{\displaystyle \pi ={\frac {2M\left(1;{\sqrt {2}}\right)^{2}}{1-\sum _{j=1}^{\infty }2^{j}c_{j}^{2}}},}
где
c
j
=
1
2
(
a
j
−
1
−
b
j
−
1
)
{\displaystyle c_{j}={\frac {1}{2}}\left(a_{j-1}-b_{j-1}\right)}
,
a
0
=
1
{\displaystyle a_{0}=1}
,
b
0
=
2
{\displaystyle b_{0}={\sqrt {2}}}
.
В то же время, если взять:
a
0
=
1
,
b
0
=
cos
α
{\displaystyle a_{0}=1,\quad \quad \quad b_{0}=\cos \alpha }
,
то
lim
N
→
∞
a
N
=
π
2
K
(
sin
α
)
{\displaystyle \lim _{N\to \infty }a_{N}={\frac {\pi }{2K(\sin \alpha )}}}
,
где
K
(
α
)
{\displaystyle K(\alpha )}
есть полный
эллиптический интеграл
K
(
α
)
=
∫
0
π
2
(
1
−
α
2
sin
2
θ
)
−
1
2
d
θ
{\displaystyle K(\alpha )=\int _{0}^{\frac {\pi }{2}}(1-\alpha ^{2}\sin ^{2}\theta )^{-{\frac {1}{2}}}d\theta }
.
То есть
π
{\displaystyle \pi }
выражается формулой:
π
=
M
(
2
)
2
N
(
2
)
−
1
{\displaystyle \pi ={\frac {M({\sqrt {2}})^{2}}{N(2)-1}}}
,
где
M
(
x
)
{\displaystyle M(x)}
— АГС 1 и
x
{\displaystyle x}
, а
N
(
x
)
{\displaystyle N(x)}
— МАГС 1 и
x
{\displaystyle x}
.
Пользуясь этим свойством, а также
преобразованиями Ландена
,
Брент
предложил
первые АГС-алгоритмы для быстрого вычисления простейших трансцендентных функций (
e
x
,
cos
x
,
sin
x
{\displaystyle e^{x},\cos x,\sin x}
). В дальнейшем исследование и использование АГС-алгоритмов было продолжено многими авторами
Примечания
B. C. Carlson.
(англ.)
//
Amer. Math. Monthly
: journal. — 1971. —
Vol. 78
. —
P. 496—505
. —
doi
:
.
An algorithm for computing logarithms and arctangents // Math.Comp.. —
Т. 26
,
№ 118
. —
С. 543—549
. —
doi
:
.
|язык=en |тип=journal |автор=B. C. Carlson |год=1972}}
↑
Adlaj, Semjon (September 2012),
(PDF)
,
Notices of the AMS
,
76
(8): 1094—1099,
doi
:
,
ISSN
от 6 мая 2016 на
Wayback Machine
.
Computation of
π
{\displaystyle \pi }
using arithmetic-geometric mean
(англ.)
//
: journal. — 1976. —
Vol. 30
,
no. 135
. —
P. 565—570
. —
doi
:
.
Landen J.
(англ.)
// Philosophical Transactions of the Royal Society of London. — 1775. —
Vol. 65
. —
P. 283—289
. —
ISSN
. —
doi
:
.
[
]
R.P. Brent
.
Fast Multiple-Precision Evaluation of Elementary Functions
(англ.)
// J. Assoc. Comput. Mach. : journal. — 1976. —
Vol. 23
,
no. 2
. —
P. 242—251
. —
doi
:
.
and
.
Pi and the AGM
(англ.)
. — New York: Wiley, 1987. —
ISBN 0-471-83138-7
.