CC BY
26
УДК 519.2
МАНИПУЛИРОВАНИЕ УЧАСТКАМИ ФУНКЦИИ ПЛОТНОСТИ ПРИ МОДЕЛИРОВАНИИ СЛУЧАЙНЫХ ПОСЛЕДОВАТЕЛЬНОСТЕЙ С РАВНОМЕРНЫМ РАСПРЕДЕЛЕНИЕМ
С.В.Гарбарь
Институт электронных и информационных систем НовГУ, sergey.garbar@gmail.com
Предложен способ моделирования стандартной равномерно распределенной случайной последовательности с заданным коэффициентом корреляции между соседними членами.
Ключевые слова: моделирование случайных последовательностей, равномерная случайная последовательность
The method of modeling sequences of standard uniformly distributed random variables with specified correlation coefficient between its adjacent members is proposed.
Keywords: random sequences modeling, uniformly distributed random sequence
Введение
Возможность моделирования случайных процессов может быть использована для имитации различных реальных ситуаций, явлений, объектов. Стационарные случайные последовательности можно использовать для описания временных рядов, течение которых стабилизировалось и происходит в неизменных условиях. Они применяются в радиотехнике, теории связи, механике жидкости и газа, океанологии, метеорологии и т.д.
Примером задачи, где применяется моделирование случайных последовательностей, может служить проведение виброиспытаний с заданными спектрами случайных процессов. Так как подготовка вибростенда к испытаниям требует усилий, проведение предварительной имитации на ЭВМ может сэкономить время и материальные ресурсы [1]. Теорема Винера—Хинчина утверждает, что автокорреляционная функция и энергетический спектр стационарного в широком смысле случайного процесса связаны между собой преобразованиями Фурье, поэтому возможность моделирования случайных последовательностей заданного вида с заданными корреляционными свойствами облегчает эту имитацию.
Постановка задачи
Рассмотрим задачу моделирования стационарной последовательности стандартных равномерно распределенных случайных величин {Т(/)} с заданным значением коэффициента корреляции между соседними членами (т. е. с заданным значением автокорреляционной функции в точке 1).
Под стандартной равномерно распределенной случайной величиной понимается случайная величина, равномерно распределенная на отрезке [0;1], имеющая функцию плотности
0, если у < 0, g (У) = 11, если у е[0;1],
0, если у > 1.
Обычно задача моделирования случайных последовательностей решается на основе алгорит-
мов авторегрессии [2], однако поиск вида независимой случайной составляющей в общем случае — довольно сложная задача. В работе [3] найдено решение данной задачи для последовательностей равномерно распределенных случайных величин с заданным значением коэффициента корреляции между соседними членами в случае, если требуе-
мый коэффициент корреляции равен
1
n +1
где n
— натуральное. Моделирование последовательности с иным значением коэффициента корреляции требует другого подхода.
Манипулирование участками функции плотности
Линейная комбинация Y ^) = a ■Т ^ -1) + Ь ■ X двух стандартных равномерно распределенных случайных величин Т ^ -1) и X имеет плотность распределения
0 при у < 0,
1
g ( У ) =
—y при y е[0; a], ab
1 при y e[a;b], b
a + b 1 r, , гл
—;------------т у при y e[b; a + b],
ab ab
0 при у > Ь,
при условии, что а < Ь . Зафиксируем а = 1 и будем считать, что Ь > 1 (иначе будет нарушено условие а < Ь). График функции плотности g(у) будет иметь вид, представленный на рис.1.
Рис.1. Плотность распределения линейной комбинации равномерно распределенных случайных величин
У ^) = а ■ У ^ -1) + Ь ■ X
Проведем следующее преобразование полученной случайной величины: если значение случайной величины у будет принадлежать отрезку
1
0;-
2 J
Если же у є
то следует вместо него взять значение 1 - у.
Ь + -; Ь +1 2
то нужно поменять его на
2 • Ь - у + 1. Таким образом поменяются вероятности появления значений случайной величины в отрезках [0;1] и [Ь;Ь +1], и, как видно по графику на рис.1, функция плотности ее распределения примет вид
ё (У) =
-, если у є
2'
2 J
0, если у > Ь +—,
что есть плотность равномерно распределенной на
отрезке
1; Ь +1
случайной величины. Взяв вместо
. 2 2 J
полученного значения у линейно преобразован-
1
2
у-
ное
Ь
получим стандартное равномерное рас-
пределение.
Если же вместо полученного на последнем шаге алгоритма значения у взять в качестве значения
случайной величины 1 - у, то знак коэффициента корреляции поменяется на противоположный (так можно поступить, так как функция плотности симметрична относительно точки 1 ). метрична относительно точки —).
Определение значения коэффициента Ь
Для получения наперед заданного требуемого значения коэффициента корреляции следует подобрать соответствующее значение коэффициента Ь. Для вычисления коэффициента корреляции между случайными величинами Т ^ -1) и Т () при различных значениях коэффициента Ь использовались численные методы (многократное моделирование соответствующих последовательностей на ЭВМ). График зависимости коэффициента корреляции от значения коэффициента Ь представлен на рис.2.
Как видно по графику, описанным способом можно получить коэффициенты корреляции, лежащие в пределах от 0 (при значениях коэффициента Ь, стремящихся к бесконечности) до 0,63 (при Ь = 1).
Для получения наперед заданного необходимого коэффициента корреляции следует подобрать соответствующее значение коэффициента Ь . Для значений, принадлежащих промежутку (0,1;0,63), можно воспользоваться приближением
г(Т^ - 1),Т^)) = (0,62 + 0,02 ■ Ь) - 0,02,
Рис.2. Зависимость значения коэффициента корреляции между соседними членами последовательности от значения коэффициента Ь
где г(Т ^ - 1),Т^)) — необходимый коэффициент корреляции. Абсолютная погрешность при этом составляет не более 0,01 (см. рис.2).
Выводы
Таким образом, алгоритм моделирования требуемой последовательности состоит из следующих шагов:
1. Определить значение коэффициента Ь, обеспечивающее требуемый коэффициент корреляции.
2. Получить промежуточное значение
у* = у^ -1) + Ьх.
3. Если у* є ние у* = 1-у * . Если у* є
0;
то присвоить новое значе-
1
Ь +—; Ь +1 2
то у* = 2Ь - у * +1.
4. Если требуемый коэффициент корреляции положительный, то линейно преобразовать у* , что-
бы получить У (ґ): у(ґ) =
у
* - 2
Ь
Если требуемый коэффициент корреляции от-
рицательный, то у(ґ) = 1 -
Богданов А.В., Кирьянов Б.Ф. Формирование случайных процессов для задания механических вибраций в вибростендах с помощью микро-ЭВМ // Проектирование радиоэлектронных систем: Сб. М., 1983. С.77-81.
Кирьянов Б. Ф. Микро-ЭВМ как средства имитации и обработки случайных процессов в радиоэлектронных системах: Монография. Новгород: НПИ, 1986. 217 с.
Гарбарь С.В. // Вестник НовГУ. Сер: Техн. науки. 2010. №55. С.53
2
Ь
