У этого термина существуют и другие значения, см.
АГС .
Среднее арифметико-геометрическое (арифметико-геометрическое среднее , АГС ) — величина, определяющаяся для двух величин
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 }
один и тот же предел[ 1] [ 2] :
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)}
.
АГС может быть применено для быстрого вычисления точного периода математического маятника [ 3] .
Часто используется сокращение
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}
. Тем самым, такая формула выражает метод Гаусса, с квадратичной сходимостью , для вычисления полного эллиптического интеграла второго рода[ 3] .
Приложения
С использованием АГС и МАГС можно вычислять значения некоторых трансцендентных функций и числа
π π -->
{\displaystyle \pi }
.
Например, по формуле Гаусса — Саламина[ 4] :
π π -->
=
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}
[ 3] .
Пользуясь этим свойством, а также преобразованиями Ландена [ 5] , Брент предложил[ 6] первые АГС-алгоритмы для быстрого вычисления простейших трансцендентных функций (
e
x
,
cos
-->
x
,
sin
-->
x
{\displaystyle e^{x},\cos x,\sin x}
). В дальнейшем исследование и использование АГС-алгоритмов было продолжено многими авторами[ 7] .
Примечания
↑ B. C. Carlson. Algorithms involving arithmetic and geometric means (англ.) // Amer. Math. Monthly : journal. — 1971. — Vol. 78 . — P. 496—505 . — doi :10.2307/2317754 .
↑ An algorithm for computing logarithms and arctangents // Math.Comp.. — Т. 26 , № 118 . — С. 543—549 . — doi :10.2307/2005182 . |язык=en |тип=journal |автор=B. C. Carlson |год=1972}}
↑ 1 2 3 Adlaj, Semjon (September 2012), "An eloquent formula for the perimeter of an ellipse" (PDF) , Notices of the AMS , 76 (8): 1094– 1099, doi :10.1090/noti879 , ISSN 1088-9477 , Архивировано из оригинала (PDF) 6 мая 2016 , Дата обращения: 14 июля 2016
↑ E. Salamin [англ.] . Computation of
π π -->
{\displaystyle \pi }
using arithmetic-geometric mean (англ.) // Math. Comp. [англ.] : journal. — 1976. — Vol. 30 , no. 135 . — P. 565—570 . — doi :10.2307/2005327 .
↑ Landen J. XXVI. An investigation of a general theorem for finding the length of any arc of any conic hyperbola, by means of two elliptic arcs with some other new and useful theorems deduced therefrom (англ.) // Philosophical Transactions of the Royal Society of London. — 1775. — Vol. 65 . — P. 283—289 . — ISSN 0261-0523 . — doi :10.1098/rstl.1775.0028 . [исправить ]
↑ R.P. Brent . Fast Multiple-Precision Evaluation of Elementary Functions (англ.) // J. Assoc. Comput. Mach. : journal. — 1976. — Vol. 23 , no. 2 . — P. 242—251 . — doi :10.1145/321941.321944 .
↑ J. M. Borwein [англ.] and P. B. Borwein [англ.] . Pi and the AGM (англ.) . — New York: Wiley, 1987. — ISBN 0-471-83138-7 .