Среднее значение круговых величин

редактировать

В математике среднее значение круговых величин - это среднее значение, который иногда лучше подходит для таких величин, как углы, дневное время и дробные части действительных чисел. Это необходимо, поскольку большинство обычных средств может не подходить для круглых количеств. Например, среднее арифметическое 0 ° и 360 ° равно 180 °, что вводит в заблуждение, поскольку для большинства целей 360 ° - то же самое, что 0 °. В качестве другого примера, «среднее время» между 23:00 и 1:00 - это либо полночь, либо полдень, в зависимости от того, являются ли эти два времени частью одной ночи или частью одного календарного дня. Это один из простейших примеров статистики неевклидовых пространств.

Содержание

1 Среднее значение углов
2 Свойства
3 Пример
4 См. Также
5 Ссылки
6 Внешние ссылки

Среднее значение углов

Поскольку среднее арифметическое не всегда подходит для углов, следующий метод может использоваться для получения как среднего значения, так и меры для отклонения углов:

Преобразование всех углов в соответствующие точки на единичной окружности, например, $α {\ displaystyle \ alpha}$ $\ alpha$ в $(соз ⁡ α, грех ⁡ α) {\ Displaystyle (\ соз \ альфа, \ грех \ альфа)}$ $(\ cos \ alpha, \ sin \ alpha)$ . То есть преобразовать полярные координаты в декартовы координаты. Затем вычислите среднее арифметическое этих точек. Полученная точка будет находиться внутри единичного диска. Преобразуйте эту точку обратно в полярные координаты. Угол - это разумное среднее значение входных углов. Результирующий радиус будет равен 1, если все углы равны. Если углы равномерно распределены по окружности, то результирующий радиус будет 0, а среднего кругового нет. (Фактически, невозможно определить непрерывную операцию среднего на окружности.) Другими словами, радиус измеряет концентрацию углов.

Учитывая углы $α 1,…, α n {\ displaystyle \ alpha _ {1}, \ dots, \ alpha _ {n}}$ $\ alpha _ {1}, \ dots, \ alpha _ {n}$ общая формула среднее значение

α ¯ = atan2 ⁡ (1 n ∑ j = 1 n sin ⁡ α j, 1 n ∑ j = 1 n cos ⁡ α j) = atan2 ⁡ (∑ j = 1 n sin ⁡ α j, ∑ j = 1 N соз ⁡ α j) {\ displaystyle {\ bar {\ alpha}} = \ operatorname {atan2} \ left ({\ frac {1} {n}} \ sum _ {j = 1} ^ {n } \ sin \ alpha _ {j}, {\ frac {1} {n}} \ sum _ {j = 1} ^ {n} \ cos \ alpha _ {j} \ right) = \ operatorname {atan2} \ left (\ sum _ {j = 1} ^ {n} \ sin \ alpha _ {j}, \ sum _ {j = 1} ^ {n} \ cos \ alpha _ {j} \ right)}

{\ displaystyle {\ bar {\ alpha}} = \ operatorname {atan2} \ left ({\ frac {1} {n}} \ sum _ {j = 1} ^ {n} \ sin \ alpha _ {j}, {\ frac {1} {n}} \ sum _ {j = 1} ^ {n} \ cos \ alpha _ {j} \ right) = \ operatorname {atan2} \ left (\ сумма _ {j = 1} ^ {n} \ sin \ alpha _ {j}, \ sum _ {j = 1} ^ {n} \ cos \ alpha _ {j} \ right)}

используя вариант atan2 функции арктангенс, или

α ¯ = arg ⁡ (∑ j = 1 n exp ⁡ (i ⋅ α j)) {\ displaystyle {\ bar {\ alpha}} = \ arg \ left (\ sum _ {j = 1} ^ {n} \ exp (i \ cdot \ alpha _ {j}) \ right)}

{\ displaystyle {\ bar {\ alpha}} = \ arg \ left (\ sum _ {j = 1} ^ {n} \ exp (i \ cdot \ alpha _ {j}) \ right)}

с использованием комплексных чисел. Чтобы соответствовать приведенному выше выводу с использованием среднего арифметического баллов, суммы должны быть разделены на $n {\ displaystyle n}$ $n$ . Однако масштабирование не имеет значения для $atan2 {\ displaystyle \ operatorname {atan2}}$ ${\ displaystyle \ operatorname {atan2}}$ и $arg {\ displaystyle \ arg}$ ${\ displaystyle \ arg}$ , поэтому его можно опустить..

Это вычисление дает результат, отличный от среднего арифметического, причем разница больше, когда углы широко распределены. Например, среднее арифметическое трех углов 0 °, 0 ° и 90 ° составляет (0 + 0 + 90) / 3 = 30 °, а среднее арифметическое составляет 26,565 °. Более того, с помощью среднего арифметического круговая дисперсия определяется только ± 180 °.

Свойства

Круговое среднее $α ¯ {\ displaystyle {\ bar {\ alpha}}}$ ${\ bar {\ alpha}}$

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

α ¯ = argmin β ∑ j = 1 nd (α j, β), где d (φ, β) = 1 - cos ⁡ (φ - β). {\ displaystyle {\ bar {\ alpha}} = {\ underset {\ beta} {\ operatorname {argmin}}} \ sum _ {j = 1} ^ {n} d (\ alpha _ {j}, \ beta), {\ текст {где}} d (\ varphi, \ beta) = 1- \ cos (\ varphi - \ beta).}

{\ displaystyle {\ bar {\ alpha}} = {\ underset {\ beta} {\ operatorname {argmin}}} \ sum _ {j = 1} ^ {n} d (\ alpha _ {j}, \ beta), {\ text {where}} d (\ varphi, \ beta) = 1- \ cos (\ varphi - \ beta).}

расстояние

d (φ, β) {\ displaystyle d ( \ varphi, \ beta)}

d (\ varphi, \ beta)

равно половине квадрата евклидова расстояния между двумя точками на единичной окружности, связанной с

φ {\ displaystyle \ varphi}

\ varphi

β {\ displaystyle \ beta}

\ beta

Пример

Простой способ вычислить среднее значение серии углов (в интервале [0 °, 360 °)) - это вычислить среднее значение косинусов и синусов каждого угла и получить угол, вычислив арктангенс. В качестве примера рассмотрим следующие три угла: 10, 20 и 30 градусов. Интуитивно, вычисление среднего значения потребовало бы сложения этих трех углов вместе и деления на 3, что в данном случае действительно привело бы к правильному среднему углу в 20 градусов. При повороте этой системы против часовой стрелки на 15 градусов три угла становятся 355 градусов, 5 градусов и 15 градусов. Наивное среднее значение теперь составляет 125 градусов, что является неправильным ответом, поскольку должно быть 5 градусов. Среднее векторное $θ ¯ {\ textstyle {\ bar {\ theta}}}$ ${\ textstyle {\ bar {\ theta}}}$ можно вычислить следующим образом, используя средний синус $s ¯ {\ textstyle {\ bar { s}}}$ ${\ textstyle {\ bar {s}}}$ и средний косинус $c ¯ ≠ 0 {\ textstyle {\ bar {c}} \ not = 0}$ ${\ textstyle {\ bar {c}} \ not = 0}$ :

s ¯ = 1 3 (sin ⁡ (355 ∘) + грех ⁡ (5 ∘) + грех ⁡ (15 ∘)) = 1 3 (- 0,087 + 0,087 + 0,259) ≈ 0,086 {\ displaystyle {\ bar {s}} = {\ frac {1} {3}} \ left (\ sin (355 ^ {\ circ}) + \ sin (5 ^ {\ circ}) + \ sin (15 ^ {\ circ}) \ right) = {\ frac {1} {3}} \ влево (-0,087 + 0,087 + 0,259 \ вправо) \ приблизительно 0,086}

\ bar s = \ frac { 1} {3} \ left (\ sin (355 ^ \ circ) + \ sin (5 ^ \ circ) + \ sin (15 ^ \ circ) \ right) = \ frac {1} {3} \ left (- 0,087 + 0,087 + 0,259 \ справа) \ приблизительно 0,086

c ¯ = 1 3 (cos ⁡ (355 ∘) + cos ⁡ (5) + cos ⁡ (15)) = 1 3 (0,996 + 0,996 + 0,966) ≈ 0,986 {\ displaystyle {\ bar {c}} = {\ frac {1} {3}} \ left (\ cos (355 ^ {\ circ}) + \ cos (5 ^ {\ circ }) + \ cos (15 ^ {\ circ}) \ right) = {\ frac {1} {3}} \ left (0,996 + 0,996 + 0,966 \ right) \ приблизительно 0,986}

\ bar c = \ frac {1} {3} \ left (\ cos (355 ^ \ circ) + \ cos (5 ^ \ circ) + \ cos (15 ^ \ circ) \ right) = \ frac {1} {3} \ left (0,996 + 0,996 + 0,966 \ right) \ приблизительно 0,986

θ ¯ = {arctan ⁡ (s ¯ c ¯) s ¯>0, c ¯>0 arctan ⁡ (s ¯ c ¯) + 180 ∘ c ¯ < 0 arctan ⁡ ( s ¯ c ¯) + 360 ∘ s ¯ < 0, c ¯>0} = arctan ⁡ (0,086 0,986) = arctan ⁡ (0,087) = 5 ∘. {\ displaystyle {\ bar {\ theta}} = \ left. {\ begin {case} \ arctan \ left ({\ frac {\ bar {s}} {\ bar {c}}} \ right) {\ bar {s}}>0, \ {\ bar {c}}>0 \\\ arctan \ left ({\ frac {\ bar {s}} {\ bar {c}}} \ right) + 180 ^ { \ circ} {\ bar {c}} <0\\\arctan \left({\frac {\bar {s}}{\bar {c}}}\right)+360^{\circ }{\bar {s}}<0,\ {\bar {c}}>0 \ end {case}} \ right \} = \ arctan \ left ({\ frac {0.086} {0.986}} \ right) = \ arctan (0,087) = 5 ^ {\ circ}.}

\bar \theta = \left. \begin{cases} \arctan \left( \frac{\bar s}{ \bar c} \right) \bar s>0, \ \ bar c>0 \\ \ arctan \ left (\ frac {\ bar s} {\ bar c} \ right) + 180 ^ \ circ \ bar c < 0 \\ \arctan \left (\frac{\bar s}{\bar c} \right)+360^\circ \bar s <0,\ \bar c>0 \ end {case} \ right \} = \ arctan \ left (\ frac {0.086} {0.986} \ right) = \ arctan (0,087) = 5 ^ \ circ.

Это может быть более лаконично понимая, что направленные данные на самом деле являются векторами единичной длины. В случае одномерных данных эти точки данных могут быть удобно представлены как комплексные числа единичной величины $z = cos ⁡ (θ) + i sin ⁡ ( θ) знак равно ei θ {\ displaystyle z = \ cos (\ theta) + i \, \ sin (\ theta) = e ^ {i \ theta}}$ $z = \ cos (\ theta) + i \, \ sin (\ theta) = е ^ {я \ тета}$ , где $θ {\ displaystyle \ theta}$ $\ theta$ - измеренный угол. Среднее значение результирующего вектора для выборки будет:

ρ ¯ = 1 N ∑ n = 1 N zn. {\ displaystyle {\ overline {\ mathbf {\ rho}}} = {\ frac {1} {N}} \ sum _ {n = 1} ^ {N} z_ {n}.}

\ overline {\ mathbf {\ rho}} = \ frac {1} {N} \ sum_ {n = 1} ^ N z_n.

Тогда выборочный средний угол равен аргумент среднего результирующего:

θ ¯ = A rg (ρ ¯). {\ Displaystyle {\ overline {\ theta}} = \ mathrm {Arg} ({\ overline {\ mathbf {\ rho}}}).}

\ overline {\ theta} = \ mathrm {Arg} (\ overline {\ mathbf {\ rho}}).

Длина результирующего вектора выборочного среднего составляет:

R ¯ = | ρ ¯ | {\ displaystyle {\ overline {R}} = | {\ overline {\ mathbf {\ rho}}} |}

\ overline {R} = | \ overline {\ mathbf {\ rho}} |

и будет иметь значение от 0 до 1. Таким образом, результирующий вектор выборочного среднего может быть представлен как:

ρ ¯ = R ¯ ei θ ¯. {\ displaystyle {\ overline {\ mathbf {\ rho}}} = {\ overline {R}} \, e ^ {i {\ overline {\ theta}}}.}

\ overline {\ mathbf {\ rho}} = \ overline {R} \, e ^ {i \ overline {\ theta}}.

Подобные вычисления также используются для определения круговая дисперсия.

См. также

Ссылки

Джаммаламадака, С. Рао и СенГупта, А.. (2001). Темы циркулярной статистики, раздел 1.3, World Scientific Press, Сингапур. ISBN 981-02-3778-2

^Кристофер М. Бишоп: Распознавание образов и машинное обучение (информатика и статистика), ISBN 0-387-31073-8

Внешние ссылки

Математика и статистика круговых значений с C ++ 11, инфраструктура C ++ 11 для круговых значений (углы, время суток и т. Д.) математика и статистика