CC BY
51
МОДЕЛЮВАННЯ І УПРАВЛІННЯ
УДК 519.2: 530.1: 600.1 И.И. ГОРБАНЬ
КРИТЕРИИ И ПАРАМЕТРЫ СТАТИСТИЧЕСКОЙ НЕУСТОЙЧИВОСТИ
Анотація. Розширено поняття статистичної нестійкості процесу (послідовності). Введено визначення статистично стійкого процесу у вузькому та широкому сенсі. Розвинута методика оцінки статистичної нестійкості на обмеженому інтервалі спостереження. Запропоновано додаткові параметри статистичної нестійкості, що характеризують порушення статистичної стійкості по відношенню до вибіркового середньоквадратичного відхилення.
Ключові слова: статистична стійкість, параметр статистичної нестійкості, теорія гіпервипа-дкових явищ.
Аннотация. Расширено понятие статистической неустойчивости процесса (последовательности). Введены определения статистически устойчивого процесса в узком и широком смыслах. Развита методика оценки статистической неустойчивости на ограниченном интервале наблюдения. Предложены дополнительные параметры статистической неустойчивости, характеризующие нарушение устойчивости по отношению к выборочному среднеквадратическому отклонению.
Ключевые слова: статистическая устойчивость, параметр статистической неустойчивости, теория гиперслучайных явлений.
Abstract. Statistical instability conception of the process (of the sequence) has been extended. Definitions of process’s statistical stability in the narrow and wide sense have been introduced. The methodology of statistical instability estimation for limited observation interval has been developed. Additional parameters of statistical instability characterizing the disturbance of statistical stability to standard deviation have been proposed.
Keywords: statistical stability, parameter of statistical instability, theory of hyper-random phenomena.
1. Введение
Одним из удивительных физических явлений окружающего мира является феномен статистической устойчивости частоты массовых событий. Долгое время полагали, что при неограниченном увеличении объема выборки частота любого массового события стремится к пределу - определенной вероятности.
Однако экспериментальные исследования на больших интервалах наблюдения процессов различной физической природы показали [1], что феномен статистической устойчивости частоты носит не идеальный характер: вначале с увеличением объема выборки дисперсия флуктуации выборочного среднего, как правило, уменьшается, но затем, достигнув определенного значения, практически перестает меняться, а в некоторых случаях начинает расти.
Таким образом, вероятность, оказывается, не более как математическая абстракция, удобная для описания реальных массовых явлений окружающего мира при относительно небольших интервалах наблюдения, однако не имеющая смысла при больших интервалах.
Для описания реальных физических событий, величин, процессов и полей с учетом ограниченного характера их статистической устойчивости была разработана физикоматематическая теория гиперслучайных явлений [1, 2].
В рамках этой теории формализовано понятие статистической устойчивости процесса. К статистически устойчивым отнесены процессы, у которых выборочные дисперсии
© Горбань И.И., 2012
ISSN 1028-9763. Математичні машини і системи, 2012, № 4
среднего стремятся к нулю при неограниченном увеличении объема выборки, а к неустойчивым - процессы, не удовлетворяющие этому требованию.
Для оценки статистической устойчивости реальных процессов была разработана специальная методика, основанная на расчете параметров статистической неустойчивости, характеризующих степень нарушения устойчивости на ограниченном интервале наблюдения.
С ее помощью на больших интервалах наблюдения были зафиксированы нарушения статистической устойчивости целого ряда физических процессов: колебаний напряжения городской сети, магнитного поля Земли, курса валют, высоты и периода морских волн, температуры воздуха, температуры воды в океане, скорости ветра и других процессов [1, 3-10].
На основе результатов проведенных экспериментальных исследований была выдвинута гипотеза [1, 2], что нарушение статистической устойчивости присуще всем реальным физическим явлениям. Исключение могут составлять, возможно, лишь фундаментальные физические постоянные, такие как скорость света, гравитационная постоянная и пр., рассматриваемые современной наукой как мировые константы.
Дальнейшие исследования показали, что в природе, хотя и редко, но встречаются физические явления, параметры статистической неустойчивости которых практически неотличимы от параметров идеально статистически устойчивых процессов. Таковыми оказались, например, колебания количества осадков на протяжении столетия [5] и колебания интенсивности излучения некоторых пульсаров на протяжении полутора десятков лет [10].
Естественно, возникает вопрос: можно ли считать эти физические процессы статистически устойчивыми? В соответствии с принятым в [1] критерием устойчивости, они формально устойчивы, однако не исключено, что при другой формализации понятия статистической устойчивости ответ на поставленный вопрос может оказаться иным.
Целью настоящей статьи является расширение понятия статистической неустойчивости процесса и развитие методики оценки статистической неустойчивости на ограниченном интервале наблюдения.
2. Известные критерии и параметры статистической неустойчивости
Определение 1а. По определению [1] последовательность Х1,Х2,... случайных величин (случайная выборка) считается статистически устойчивой (статистически стабильной), если при устремлении объема выборки N к бесконечности математическое ожидание выборочной дисперсии
- 1 N
°г, = лц£ (У -шм )2 О)
флуктуации выборочного среднего
1 п ------------
У = - £ X. (п = 1, М) (2)
П г=1
_ 1 N
стремится к нулю, где Шу =— £ Уп - выборочное среднее флуктуации среднего. После-
М М п=1
довательности, не удовлетворяющие этому условию, статистически неустойчивы.
Уменьшение выборочного среднего по мере увеличения объема данных может быть вызвано не только стабилизацией среднего, но также уменьшением дисперсии исходного процесса. Для подавления влияния изменения дисперсии процесса в работе [10] предложено следующее определение понятия статистической устойчивости процесса.
Определение 1б. Последовательность случайных величин Х1, Х2,... названа статистически устойчивой, если при устремлении объема выборки N к бесконечности параметр статистической неустойчивости
М Г Бг ] М Г Бг ]
.. = 1_ Ум J = Ь (3)
ND Бх ( }
Ум xм
стремится к нулю, где Dy^ - дисперсия выборочного среднего:
1 N
D =А- £ D
Ум д т2 х
¿У п=1
М[ ] - оператор математического ожидания, Dx = М[(Хп — шх ) ] - дисперсия случайной величины Хп, а шх = М[Хп ] - ее математическое ожидание, D^ - средняя дисперсия
____ _ 1 N
случайных величин Хп (п = 1, ^ : Dx = — £ Dnn
N п=1
При выполнении условий Определения 1а или Определения 1б последовательность естественно называть статистически устойчивой по отношению к среднему, а при невыполнении этих условий - статистически неустойчивой по отношению к среднему.
Процесс обычно описывается последовательностью значений в фиксированные моменты времени, поэтому, говоря о процессах, можно отождествлять их с соответствующими последовательностями.
Заметим, что в частном случае постоянной дисперсии понятия статистически устойчивого процесса, согласно Определениям 1 а и 1 б, совпадают.
По результатам наблюдения процесса на конечном интервале невозможно точно установить факт нарушения статистической устойчивости. Однако можно количественно оценить степень флуктуации выборочного среднего в фиксированные моменты времени и, анализируя динамику происходящих изменений, выявить некоторые тенденции, ведущие к нарушению статистической устойчивости.
Такие тенденции характеризуют параметр статистической неустойчивости уN и параметр цN , связанный с параметром уN, соотношением
Ц N = >/ У N1(1 + У N ) . (4)
В отличие от параметра уN, ограниченного лишь снизу нулевым значением, параметр ^ ограничен как снизу, так и сверху: минимально возможное его значение равно нулю, а максимально возможное - единице.
В качестве единицы измерения статистической неустойчивости параметра yМ при
фиксированном N может выступать величина у0М [8], представляющая собой параметр уN, рассчитанный для идеальной статистически устойчивой последовательности N некоррелированных отсчетов с постоянной дисперсией Dx = Dx и нулевым математическим ожиданием. Поэтому параметром статистической неустойчивости является также параметр
^ = У N IУ 0 N
[8], диапазон изменения которого - [0, ¥) .
■ | * * 1 *
Для практических расчетов можно использовать оценки уN, цN, ^ параметров
У N , Ц N , ^ .
Оценка yМ может быть вычислена по множеству реализаций или как y*N = , где
DXм
- 1 N
—хк = — £ —*Х - среднее оценок дисперсии —X , сформированное по отдельным фраг-
М N п=1 " "
ментам реализации последовательности Хп (п = 1, N). Оценки же и ИХ могут быть
^ < * / * ! s л * \ 1 * * /
рассчитаны при наличии оценки уN по формулам ^N =^уN/ (1 + уN ) , hN = yN/ у0N .
Факт нарушения устойчивости процесса устанавливается по тенденции изменения параметров уN, цN, hN (или y*N, ^N, hN ). Если при больших значениях N наблюдается рост этих параметров или стабилизация значения на ненулевом уровне, процесс считается неустойчивым, в противном случае - устойчивым.
Известная методика оценки статистической неустойчивости [1] основана на расчете этих параметров.
3. Альтернативные критерии и параметры статистической неустойчивости
Рассмотренные критерии и параметры статистической неустойчивости процесса отслеживают динамику изменения важнейшей величины процесса - выборочного среднего. Однако при этом остаются без внимания функция распределения и другие выборочные моменты.
Стремление к нулю математического ожидания выборочной дисперсии среднего процесса (или какого-либо связанного с ним параметра) не гарантирует сходимость выборочной функции распределения к какой-либо функции распределении. Поэтому Определения 1 а и 1 б нельзя признать безукоризненными.
Принимая во внимание сказанное, можно предложить несколько альтернативных вариантов определения понятия статистически устойчивой последовательности (процесса).
Рассмотрим последовательность случайных величин X1,X2,... (с различными в общем случае законами распределения). Пусть N (х) - количество членов последовательности {XN} = X1,X2,...,XN объема N, меньших х, а F*(х) = - эмпирическая (выбо-
рочная) функция распределения, представляющая собой неубывающую ступенчатую функцию.
При неограниченном увеличении числа членов N эмпирическая функция распределения FN* (х) не обязательно сходится к определенной функции.
Определение 1в. Последовательность случайных величин X1, X2,... назовем статистически устойчивой в узком смысле почти наверное (с вероятностью единица), если существует случайная величина X, описываемая функцией распределения F (х), к которой сходится почти наверное (с вероятностью единица) эмпирическая функция распределения F*(х) при неограниченном увеличении N : lim F*(х) = F(х) , т.е.
P{lim sup FN(х)-F(х)| = o} = 1. (5)
[N ®¥-¥< х<¥ J
Если такая случайная величина X не существует, последовательность будем называть статистически неустойчивой в узком смысле.
Рис. 1. Статистически устойчивые и неустойчивые случайные последовательности (процессы) в узком и широком смыслах
Поскольку сходимость с вероятностью единица более сильная, чем сходимость по вероятности, то в случае, когда FN*( х) сходится к ^ (х) почти наверное, имеет место сходимость Г*( х) к ^ (х) и по вероятности (рис. 1).
Определение 1в может быть полезным при проведении теоретических исследований, однако оно мало пригодно для оценки нарушений статистической устойчивости реальных процессов, поскольку оценка сходимости эмпирической функции распределения представляет собой практически непреодолимую проблему.
Задачу можно существенно упростить, если ограничиться оценкой математических ожиданий выборочных дисперсий двух первых выборочных моментов: среднего и дисперсии (или среднеквадратического отклонения (СКО)).
Определение 1г. Последовательность Х1, Х2,... случайных величин назовем статистически устойчивой в широком смысле, если при устремлении объема выборки N к бесконечности
1) математическое ожидание выборочной дисперсии (1) флуктуации выборочного среднего (2) и
2) математическое ожидание выборочной дисперсии
- 1 N
—г, = ^ £ (г, - )2 (6)
выборочного СКО
г =
_ 1 N
стремятся к нулю, где тг =------£ 2п - среднее флуктуации выборочного СКО.
М N — 1 п=2
Последовательности, не удовлетворяющие этому условию, будем называть статистически неустойчивыми в широком смысле.
Таким образом, статистически устойчивый процесс в широком смысле - процесс, статистически устойчивый по отношению к среднему и к СКО.
п
1 п
Е (X
1 ,=1
■Уя )2 (п = 2, М)
(7)
Словосочетания «в узком смысле» и «широком смысле» используются здесь и далее по аналогии с тем, как они применяются применительно к понятию стационарного процесса [11].
Обратим внимание, что области статистической устойчивости (как в широком, так и в узком смысле) располагаются внутри области существования функции распределения Е (х), а области неустойчивости охватывают не только части области, в которой существует функция распределения Е (х), но и область, в которой она не существует (рис. 1).
Не всякий процесс, устойчивый в узком смысле, является устойчивым в широком смысле, и, наоборот, не всякий процесс, устойчивый в широком смысле, является устойчивым в узком смысле.
Для подавления влияния изменения дисперсии исследуемого процесса можно предложить еще один, более конструктивный, вариант определения понятия статистической устойчивости в широком смысле.
Определение 1д. Последовательность случайных величин Х1, Х2,... назовем статистически устойчивой в широком смысле, если при устремлении объема выборки N к бесконечности стремятся к нулю
1) параметр статистической неустойчивости по отношению к среднему (3) и
2) параметр статистической неустойчивости по отношению к СКО:
М Г Ё>7 1 М Г Ё>7 1
Г" =~Ж~=~^Г. (8)
хм
Последовательности, не удовлетворяющие этому условию, будем называть статистически неустойчивыми в широком смысле (рис. 1).
Наряду с параметром ГМ можно использовать параметр ММ = ^jгN/(1+ ГN ) , аналогичный параметру цN, и параметр HN = ГN /у0N , аналогичный параметру ^ .
Параметры ^, ММ и HN - безразмерные величины. Параметры ^ и HN ограничены снизу нулевым значением, а параметр МN - нулевым значением снизу и единичным значением сверху.
Для практических расчетов на конечном интервале наблюдения вместо параметров
^, ММ и HN целесообразно использовать их оценки: ГN = , MN =^'N1 (1 + ^) и
-х„
К =Г'м/ То N .
Таким образом, методика оценки статистической неустойчивости физических процессов, основанная на расчете параметров статистической неустойчивости по отношению
'''
к выборочному среднему уМ , цN , ^, дополняется расчетами параметров статистической неустойчивости по отношению к выборочному СКО ^, MN, К .
Обратим внимание на одно важное обстоятельство. При расчете оценок параметров статистической неустойчивости yМ, цN, Кг и ^, MМ, ^ на основе выборочных дисперсий не используется тот факт, что исследуемые величины Х1, Х2,... являются случайными (т.е. характеризуются вероятностной мерой). Поэтому понятие статистической устойчивости в широком смысле можно обобщить на любые последовательности, в том числе и неслучайные.
Определение 1е. Последовательность величин Х1, Х2,... назовем статистически устойчивой в широком смысле, если оценки параметров статистической неустойчивости по
* 4 *
отношению к среднему YN = и по отношению к СКО Гд, = при неограниченном
XМ XМ
увеличении объема выборки N стремятся к нулю: lim yn = 0 , lim TN = 0.
N
Последовательности, не обладающие этим свойством, будем считать статистически неустойчивыми в широком смысле.
4. Исследование статистической устойчивости реальных процессов
Для иллюстрации информативности дополнительных параметров оценки статистической устойчивости TN, MN и HN на рис. 2 и 3 приведены результаты расчетов оценок этих параметров для двух астрофизических аккрецирующих источников рентгеновского излучения: GRS 1915+105 (рис. 2) и PSRJ 1012+5307 (рис. 3).
Рис. 2. Выборочное среднее Yn и выборочное среднеквадратическое отклонение Zn колебания интенсивности излучения источника GRS 1915+105 (а) и соответствующие параметры статистической неустойчивости yN , TN (б), ^N, MN (в) и h*, H* (г)
Как и в статье [10], исходные данные взяты с сайта [12]. Измерения интенсивности излучения осуществлялись в период с 1 января 1996 г. по 31 декабря 2011 г. Средняя периодичность измерений составляла 2,7 ч для GRS 1915+105 и 2,8 ч для PSRJ 1012+5307.
На рис. 2а и 3 а сплошными линиями изображены колебания во времени выборочных средних Yn и выборочных СКО ZN, на рис. 2б, 3б - зависимости от времени параметров yN , TN, на рис. 2в, 3в - зависимости от времени параметров ^N , MN, на рис. 2г, 3г -зависимости от времени параметров hNN , HNN .
Пунктирными линиями на рис. 2б, 3б, а также 2в, 3в и 2г, 3г представлены соответственно параметры yn , ЦN и hN для идеального устойчивого процесса, а точечными линиями - СКО от них.
Для источника GRS 1915+105 (рис. 2) сильные нарушения статистической устойчивости по отношению к среднему сопровождаются значительными нарушениями устойчи-
вости по отношению к СКО. Степень нарушения устойчивости по отношению к среднему выше, чем по отношению к СКО.
10~2 101 10° 101 t, лет ю‘2 ю1 10° 101 t, лет
Рис. 3. Выборочное среднее Yn и выборочное среднеквадратическое отклонение Zn колебания интенсивности излучения пульсара PSR J1012+5307 (а) и соответствующие параметры статистической неустойчивости yN , С (б), ^N, MN (в) и h*, H* (г)
Для источника PSRJ 1012+5307 ситуация иная (рис. 3): нарушения статистической устойчивости по отношению к среднему не наблюдаются, однако имеют место нарушения устойчивости по отношению к СКО.
Отсюда следует, что колебания обоих рассматриваемых источников излучения статистически неустойчивы в широком смысле.
Аналогичные исследования флуктуации количества осадков в нескольких регионах выявили нарушения статистической устойчивости по отношению к СКО.
Таким образом, претендентами на роль статистически устойчивых процессов в широком смысле остаются пока лишь фундаментальные физические постоянные (мировые константы).
5. Выводы
1. Расширено понятие статистической неустойчивости процесса (последовательности). Введены определения статистически устойчивого процесса в узком и широком смыслах.
2. Развита методика оценки статистической неустойчивости на ограниченном интервале наблюдения. Предложены дополнительные параметры статистической неустойчивости, характеризующие нарушение устойчивости по отношению к выборочному среднеквадратическому отклонению.
3. На примере астрофизических источников излучения продемонстрирована информативность предложенных дополнительных параметров статистической неустойчивости.
4. Установлено, что претендентами на роль статистически устойчивых процессов в широком смысле остаются пока лишь фундаментальные физические постоянные (мировые константы).
СПИСОК ЛИТЕРАТУРЫ
1. Горбань И.И. Теория гиперслучайных явлений: физические и математические основы [Электронный ресурс] / Горбань И.И. - К.: Наукова думка, 2011. - 318 с. - Режим доступа: http://www.immsp.kiev.ua/perspages/ gorban i i/index.html.
2. Горбань И.И. Теория гиперслучайных явлений [Электронный ресурс] / Горбань И.И. - К.: ИПММС, 2007. - 181 с. - Режим доступа: http://www.immsp.kiev.ua/perspages/ gorban i i/index.html.
3. Горбань И.И. Нарушение статистической устойчивости физических процессов / И.И. Горбань // Математичні машини і системи. - 2010. - № 1. - С. 171 - 184.
4. Gorban I.I. Disturbance of statistical stability / I.I. Gorban // Information Models of Knowledge. -Sofia: ITHEA, 2010. - P. 398 - 410.
5. Г орбань И.И. Статистическая устойчивость колебаний температуры воздуха и осадков в районе Москвы / И.И. Горбань // Математичні машини і системи. - 2011. - № 3. - С. 97 - 104.
6. Исследование статистической устойчивости колебаний температуры шельфовой зоны окраинных морей / И.И. Горбань, Н.И. Горбань, В.В. Новотрясов [и др.] // Седьмой Всероссийский симпозиум «Физика геосфер». - Владивосток, 2011. - С. 542 - 547.
7. Gorban I.I. Disturbance of statistical stability (part II) / I.I. Gorban // International Journal “Information Theories and Applications”. - 2010. - Vol. 18, N 4. - P. 321 - 334.
8. Горбань И.И. Статистическая неустойчивость физических процессов / И.И. Горбань // Известия вузов. Радиоэлектроника. - 2011. - Т. 54, № 9. - С. 40 - 52.
9. Горбань И.И. Исследование нарушений статистической устойчивости колебаний скорости ветра в Чернобыле / И.И. Горбань, А.Д. Скорбун // Системи підтримки прийняття рішень. Теорія і практика: зб. доп. наук.-прак. конф. з міжнар. участю. «СППР’2012». - К., 2012. - С. 39 - 42.
10. Горбань И.И. Статистическая устойчивость излучения астрофизических объектов / И.И. Горбань // Математичні машини і системи. - 2012. - № 2. - С. 155 - 160.
11. Горбань 1.1. Теорія ймовірностей і математична статистика для наукових працівників та інженерів [Електронний ресурс] / І.І. Горбань. - К.: ІПММС НАН України, 2003. - 245 с. - Режим доступу: http://www.immsp.kiev.ua/perspages/gorban i i/index.html.
12. All-Sky Monitor (ASM) team at the Kavli Institute for Astrophysics and Space Research at the Massachusetts Institute of Technology [Електронний ресурс]. - Режим доступу:
http://xte.mit.edu/ASM lc.html.
Стаття надійшла до редакції 27.07.2012
