CC BY
59
УДК 621.396
АЛГОРИТМ ПОСТРОЕНИЯ СЛУЧАЙНЫХ РЕАЛИЗАЦИЙ НОРМАЛЬНОГО МАРКОВСКОГО ПОЛЯ НА ПОВЕРХНОСТИ ТОРА А.С. Мазманишвили
Сумской государственный университет,
Сумы, Украина, e-mail: [email protected]
Рассмотрена задача построения алгоритма численной генерации нормальных двумерных марковских полей на цилиндрической поверхности и на поверхности тора. Эти поля обладают свойством однородности по отношению к сдвигам вдоль поверхностей, соответственно, цилиндра и тора. Описаны возможные перспективы применения нормальных марковских полей, реализованных на поверхности тора.
Ключевые слова: алгоритм, нормальное марковское поле, случайная реализация, тор, цилиндр.
1. Введение
Задачи численного стохастического моделирования начинаются с синтеза алгоритмов генерации случайных величин с необходимыми для данной задачи свойствами. При моделировании случайных полей, принимается допущение о евклидовой структуре пространства, в котором разыгрываются случайные события. Это допущение приводит к тому, что эти поля определяются на отрезке, в прямоугольнике, или, наконец, в многомерном гиперкубе [1, 2]. По этой причине, как правило, разрабатываемые алгоритмы опираются на прямоугольную декартову систему даже в тех случаях, когда исследуемый объект мало удобен при описании в этой системе. Именно такая ситуация имеет место при построении случайных полей на торе или цилиндре [2]. В то же время, алгоритмы генерации двумерных марковских полей, представленные в [3, 4], таковы, что они приспособлены только для генерации случайных полей, задаваемых на прямоугольнике, стороны которого параллельны декартовым осям на плоскости. Однако, в задачах машиностроения [5], технической электродинамики [6], навигации [7] и других в большинстве объекты не являются плоскими или прямоугольными, и их описание в прямоугольной декартовой системе координат неестественно и, вообще говоря, приводит к неудовлетворительным результатам.
Настоящая работа посвящена попытке построения алгоритма численной генерации нормального марковского случайного поля на поверхности тора. Поперечными сечениями такого поля, реализующимися вдоль большого радиуса Я тора, являются марковские процессы Орнштейна-Уленбека, определённые на отрезке, равном длине соответствующей окружности, а затем продолженные периодически на всю числовую ось изменения угловой координаты (т.н. "броуновский мост"). Сечениями же случайного поля, которые реализуются на окружности заданного малого радиуса тора г, являются случайные процессы типа броуновского моста, построенные по такой же схеме но с периодом, равным соответствующей длине окружности. Тогда, внутри замкнутого прямоугольника, накрывающего тор, реализуется стационарное марковское гауссовское поле. В нашей работе построение алгоритма численной реализации указанного случайного поля основывается на иерархическом принципе. При этом на нижнем уровне иерархии используется стандартный генератор нормального "белого" шума п(Ь) с интенсивностью (среднем квадратичным отклонением) а. На его основе строятся случайные реализации процессов вдоль образующих тора. Последовательное заполнение зацепленными броуновскими мостами и позволяет сгенерировать численно реализации марковского поля на поверхности тора.
Рис. 1: Идеальный недеформированный тор с параметрами: Я = 6, г = 3, N = 40, М = 40.
2. Алгоритм генерации нормальных двумерных марковских полей
на поверхностях цилиндра и тора
Обозначим Дх - величину шага вдоль оси х на плоскости при численной реализации случайного поля, а Ду - величину шага вдоль оси у. Соответственно, N - число узлов по х, а М - число узлов вдоль оси у. Обозначим Н -
случайное поле, которое является проекцией рассматриваемого нами случайного поля на совокупность узлов {кщт}, п = 0 ^ N, т = 0 ^ М. Обозначим, далее, ух, иу - декременты затухания нормального марковского поля, соответственно, вдоль продольной и поперечной координатных осей. Корреляторы этого поля в каждых двух соседних узлах равны
р = ехр(-МхДх) , Я = ехр(-^уДу) . (1)
Рассмотрим сначала сечение случайного поля при у = 0, когда т = 0. Совместная плотность распределения вероятностей случайных значений поля в узлах с номерами п = 0 ^ N имеет, согласно [3, 5], вид
/„ ({/Чо}) = [27г;72(1 _ р2)]ЛУ2 ехр |“2<гЦ1 - Vі) ^(/Ч° — Р*«-1,о)21 ■
(2)
При этом справедливо следующее правило нумерации, обеспечивающее периодичность случайного поля = ^о,о, +1 = Н\ и т.д. Пусть ип,о - генератор
белого шума в рассматриваемых целочисленных точках, обладающий нулевым средним и (иП0) = ^2, п = 1,..., N. Тогда из вида квадратичной формы в показателе экспоненты в плотности (2) следует система связей
Ь,п,о = рйп-1,о + аип,о , п = 1,N (3)
с постоянной ___________________
а = д/ (1 -Р2)(1 - РШ)~1, (4)
которая необходима для нормировки вторых моментов поля. Набор значений поля в нулевом по х сечении (т = 0) может быть выражен явно через величины ип,0 следующим образом
^і,о = а (иі,о + ри2,о + р2из,о + ... + рм-1иж,о) ,
^2,о = а (рм-1иі,о + и2,о + риз,о + ... + рМ-2иж,о) ,
^з,о = а (рм-2иі,о + рн-1и2,о + из,о + ... + рн-зи^,о) , (5)
^ж,о = а (риі,о + р2и2,о + р3из,о + ... + и^,о) .
Таким образом, на основе (5), построен броуновский мост при п = 0. При
выбранном значении а из (2) и (5) следуют выражения для безусловных средних
(^п,0^п+к,о) = рк, к = 0, (6)
Из формулы (6) вытекает свойство стационарности броуновского моста в нулевом сечении.
Рассмотрим сечение поля при у = Ду, т.е. при т = 1. При фиксированном значении в нулевом сечении вдоль оси у, совместная плотность распределения вероятностей случайных значений поля имеет вид [3, 5]
/„({/41}) = [2тп72(1 - р2)( 1 - д2)] ~м/2ехр |-2(г2(1 _ ^)(1 _ д2) | • (7)
где ^1 - квадратичная форма относительно {^пд}
22 ^1 = (^1,1 - Лд - д^,о + рдЦо) +(^2,1 - Лд - ^2,0 + рд^,о) + ...
... + (^N,1 - -1,1 - ^^N,0 + -1,о)2 . (7а)
При этом также действует циклическое правило нумерации ^жд = ^о,1 и т.д. Из вида квадратичной формы (7а) вытекает, что значение поля ^пд связано с тремя "предыдущими" значениями. Тогда на основе плотности (7) можно прийти к следующей системе связей для броуновского моста в первом сечении
К,1 = Р^п-1,1 + Л,о - р?^п-1,о + а^п,1, п = 1,N. (8)
Из (8) видно, что система величин йпд = ^пд - р^п,о, где п = 1,...,Ж,
подчиняется тем же соотношениям, что и величины (3). Следовательно, два
шага вычислений
51.1 = а (щд + РМ2,1 + Р2М3,1 + ... + -1Мжд) ,
52.1 = а (рж-1М1,2 + РМ2,2 + Р2М3,2 + ... + -2Мж,^ ,
53.1 = а (рж-2М1,3 + -1М2,3 + риз,з + ... + -3иж,^ , (9)
5ж,1 = а (рМ1,ж + Р2^2,ж + Р3^3,ж + ... + иж,ж) .
и
К,1 = Л,о + 5п,1 , п = 1, ..., N (10)
дадут искомый набор значений броуновского моста в первом сечении.
Рассмотрим сечение за номером m. Для него плотность распределения значений поля в узлах (n,m) следующая:
/я({Л»,»})= [27Г(72(1-р2)(1-д2)]~‘¥/2ехр|-2^2(1 _д2)| , (11)
где Qm - квадратичная форма относительно {hn,m},
Qm = (hl,m — ph0,m — qhl,m-l + pqh0,m-l) +
2
+ (h2,m — phl,m — qh2,m-1 + pqhl,m-l) + ...
2
... + (hN,m — phN —l,m — qhN,m-l + pqhN-l,m-l) . (11a)
Теперь, если в (9) и (lO) везде заменить в y-ом индексе 0 на m — 1 и l на m, то получим алгоритм нахождения набора значений броуновского моста в m-ом сечении. А именно, для произвольного m-го сечения заполнение его узлов осуществляется следующими двумя шагами:
sl,m+l = а (ul,m + pu2,m + p2u3,m + ... + pN luN,m) ’
52,m+l = а (pN—lUl,m + №,m + p2U3,m + ... + pN—2UN,m) ,
53,m+l = а (pN—2Ul,m + PN —lU2,m + №,m + ... + pN—3UN,m) , (12)
5Ж,т+1 = а (ри1,ш + р2^2,ш + р3и3,т + ... + иЖ,ш) .
и
^п,т — р^п,т-1 + 5п,т > п — 1? ...; N . (13)
Используя векторно-матричные обозначения, выражения (9) и (11) можно записать в компактной форме, введя циркулярные матрицы [9].
Объединяя формулу (5) и М раз повторенные формулы (11) и (12), получим итоговый алгоритм генерации значений поля на прямоугольнике, накрывающем поверхность цилиндра и образованном из N х (М + 1) узлов с учетом их периодичности вдоль оси х. При выбранном значении а, из этих формул, плотностей (2), (7) и правила нумерации, учитывающего периодичность, получаем для математических ожиданий
(^п,т^п+к,т+^') = /^^2 , к,; = 1, ..., N , (^П,т) = ^ . (14)
Таким образом, во всех сечениях, сгенерированное на цилиндрической поверхности поле характеризуется свойствами однородности. Свойство марковости непосредственно вытекает из вида плотности условных вероятностей перехода (7). Построенное на цилиндре поле обладает свойствами нормальности и марковости в том смысле, что плотности условных вероятностей перехода имеют гауссовскую форму вдоль выбранных направлений.
Рис. 2: Вид сбоку сгенерированного случайного поля на торе с параметрами: р = 0.5, д = 0, 5, а = 0.2 в виде случайной деформации тора, пропорциональной значению поля в каждой из точек. Сеткой показана связь случайных процессов в соседних узлах.
Теперь, с целью построения случайного поля на поверхности тора, осуществим замыкание построенного нами поля на поверхности цилиндра вдоль оси у. Заметим, что если в выражениях, описывающих замыкание вдоль оси х, везде выполнить замены по правилу {Ж ^ М,р ^ д}, то будет получено замыкание вдоль оси у. А именно, для произвольного п-го сечения заполнение его узлов осуществляется следующими двумя шагами:
£п+і,і = в (ип,і + дип,2 + д2мп,з + ... + 3м-1мп,м) ,
£п+2,2 = в -1ип,1 + дип,2 + 32Мп,з + ... + 3м-2Мп,м) ,
£п+і,з = в (зм-2ип,1 + 3м-1ип,2 + дмп,з + ... + 3м-3мп,м) , (15)
£п+1,м = в (^ип,1 + ^2^п,2 + ^3^п,3 + ... + ип,м)
и
При этом введена новая нормировочная постоянная, необходимая для сохранения вторых моментов амплитуды вещественного поля,
/3= \/(1-<г'2)(1-<г2М)-1- (17)
В результате указанного замыкания будет построено случайное поле на поверхности цилиндра, ориентированного вдоль оси х.
Совместное действие двух замыканий приведет к построению стохастического поля на поверхности тора. При этом для сохранения вторых моментов значений поля должна быть использована новая нормировочная постоянная
7 = а/3 = ^(1 -г>2)(1 -^)-!(1 _92)(1 - ^2М)—1, (18)
Рис. 3: Сгенерированная случайная реализация поля на торе (вид сбоку).
3. Применение алгоритмов генерации нормальных двумерных марковских полей на поверхностях цилиндра и тора
Рассмотрим геометрическую интерпретацию нормального двумерного марковского поля на цилиндрической и торической поверхностях. Тор характеризуется двумя радиусами - азимутальным Я и радиусом г кругового сечения плоскостью, проходящей через центр тора перпендикулярно к азимутальной плоскости. Соответственно при заданных Я и г уравнения, определяющие тор, можно параметризовать с помощью двух азимутальных углов: ^ - угол в плоскости ХОУ, проходящей перпендикулярно оси OZ, и ф - фазовый угол кругового сечения. Задав сетку из N х М узлов, множество узлов на торе
Рис. 4: Сгенерированная случайная реализация поля на торе (вид сверху).
можно определить формулами
Xn,m = (R + Г Sin фп) COS m ,
Yn,m = (R + Г Sin фп) Sin , (19)
Zn,m — r COS фт •
Наложение случайного поля H со случайными значениями {hn,m} в узлах выбранной сетки осуществляется на основе выражений
Xn,m = [R + (Г + hn,m) Sin фп] COS ^m ,
Ynm — [R + (r + hn,m) Sin Фп] Sin ^m , (20)
Zn,m — (r + hn,m) COS фm •
Случайное поле H, синтезированное в узлах сетки, стационарно, поскольку для него, с учетом нормировки (15), справедливо соотношение (13), гарантирующее инвариантность средних относительно сдвигов вдоль координатных осей. При программировании оказываются более удобными аналоги формул (2), (9), (10), (12), выполненные в терминах векторов и циркулярных матриц [9]. Из рисунков 2-4 можно заметить самосогласованное поведение значений
поля, изменение сечений - случайных процессов, а также их связь в соседних
узлах.
5. Выводы
Таким образом, иерархический подход синтеза реализаций случайного поля на двумерной поверхности дает возможность численно генерировать марковские поля на поверхности тора. По существу развиваемый метод отличается от метода генерации полей на плоскости [9] только действием двойного правила периодичности, а в роли стационарного распределения вероятностей выступают броуновские мосты. Область применения подобных случайных полей достаточно широка. Они возникают всякий раз, когда необходим учет случайного изменения некоторой "идеальной" поверхности. Недостатком реализованного подхода является неаналитическое введение азимутальных периодичностей броуновских мостов. Возможность применения подобного подхода для пространственных тел более сложной структуры, например, таких как сфера, остается открытой.
Литература
1. Бусленко Н.П. Метод статистических испытаний / Н.П.Бусленко, Ю.А.Шрейдер. - М.: Физматгиз,1961. - 280с.
2. Харин Ю.С. Практикум на ЭВМ по математической статистике / Ю.С.Харин, М.Д.Степанова. - Минск: "Университетское", 1987. - 304с.
3. Habibi A. Two-Dimensional Bayesian Estimate of Image // Proc. IEEE. -1972. - 60;7. - P.878-883.
4. Хусу А.П. Шероховатость поверхностей / А.П.Хусу, Ю.Р.Витенберг, В.А.Пальмов. - М.: Наука,1975. - 344с.
5. Шифрин Я.С. Вопросы статистической теории антенн / Я.С.Шифрин.
- М.: Советское радио,1970. - 384с.
6. Ярлыков М.С. Статистическая теория навигации / М.С.Ярлыков. - М.: Радио и связь,1985. - 344с.
7. Рытов С.М. Введение в статистическую радиофизику / С.М.Рытов. -М.: Наука,1966. - 404с.
8. Мазманишвили А.С., Щербань В.Е. Моделирование марковских случайных последовательностей и алгоритм генерации однородного двумерного марковского поля // Электронное моделирование. - 1996. - 18;2. - С.93-95.
9. Воеводин В.В. Матрицы и вычисления / В.В.Воеводин, Ю.А.Кузнецов.
- М.: Наука. - 1984. - 320с.
10. Мазманишвили А.С. Алгоритм генерации нормального марковского поля на поверхности идеального цилиндра // Электронное моделирование.
- 1998. - 20;6. - С.65-69.
THE GENERATION ALGORITHM FOR RANDOM REALIZATIONS OF THE NORMAL MARKOVIAN FIELD ON TORUS SURFACE A. S.Mazmanishvili
Sumy State University, Sumy, Ukraine, e-mail: [email protected]
The paper deals with the generation problem of normal two-dimensional Markovian fields on cylindrical and torus surfaces. The generation algorithm for such fields is developed. Generated fields are homogeneous. The numerical example of generated field is given.
Key words: algorithm, normal markov field, random realization, torus, cylinder.
