УДК 519.23
К СТАТИСТИКЕ МОРФОЛОГИИ СЛУЧАЙНЫХ ПРОСТРАНСТВЕННО-ВРЕМЕННЫХ СИГНАЛОВ, СФОРМИРОВАННЫХ ДВУМЕРНЫМ ТОЧЕЧНЫМ ПОТОКОМ ВОССТАНОВЛЕНИЯ
Д.В. Дубинин, А.И. Кочегуров*, В.Е. Лаевский (V. Geringer)**
Томский государственный университет систем управления и радиоэлектроники E-mail: dima@info.tusur.ru *Томский политехнический университет E-mail: kai@cc.tpu.edu.ru ** Baden-W?rttemberg Cooperative State University, Faculty of Engineering,
Friedrichshafen, Germany E-mail: geringer@DHBW-Ravensburg.de
Исследуются характеристики морфологии пространственно-временных сигналов, полученных при помощи двумерного рекуррентного точечного потока восстановления [1] обладающего рядом свойств: стационарностью, отсутствием последствий, ординарностью. Исследования проводились с помощью программного комплекса стохастического моделирования «КИМ СП» [2]. Выбор групп составных элементов контурного рисунка поля, определяющих его морфологию (форму контурной структуры), значения финальных вероятностей появления составляющих элементов, задается оператором.
Ключевые слова:
Исследования на моделях, стохастическое моделирование, моделирование пространственно-временных сигналов, двумерный поток восстановления, поток Пальма, одноуровневая цепь Маркова.
Key words:
Research on models, stochastic simulation, modelling the space-time signals, two-dimensional renewal process, ordinary point process of Palms, one-step Markov chains.
Актуальность проблемы
В настоящее время широкое применение находят системы получения информации, включающие обработку пространственно-временных сигналов. При решении задач синтеза, доводки существующих комплексов, их оптимизация немыслима без поддержки компьютерных систем, в частности прикладного моделирования сложных процессов и статистического анализа многомерных данных.
Большое значение имеет выбор адекватной математической модели при построении наблюдаемых данных. На данный момент не существует универсального способа формирования случайных полей с произвольно заданными характеристиками, а полное решение проблемы описания реальных изображений в общем случае отсутствует. Поэтому известные модели случайных полей соответствуют реальным пространственно-временным сигналам лишь по ограниченному числу характеристик [3] (форма корреляционной функции, радиус корреляции, контрастность, размеры и т. п.). В работах [2, 4, 5] приведена методика построения случайных яркостных полей, аппроксимированных двумерным рекуррентным точечным потоком восстановления.
В данной работе приводятся результаты статистического анализа полей, полученных по методике [2]: исследование закона распределения числа сегментов изображения, зависимость его от морфологии поля, определение среднего значения и дисперсии числа сегментов.
Результаты обработки статистических данных
В ходе компьютерного моделирования были построены яркостные поля на различных типах
морфологий изображений: мозаики типа В, ВВ, ВС, ВБ, ВВ. Каждое изображение строилось построчно. Оно содержало пхп простейших элементов, где п -число отсчетов. Совокупность простейших элементов для каждого типа мозаики приведена в работе [2]. На рис. 1 показаны моменты восстановления в структуре получаемого изображения.
Мозаики каждого типа были получены для размерностей изображений 30, 40, 50, 60, 70, 80, 90, 100, 150, 200 отсчетов. На рис. 2 приведены примеры некоторых мозаик типа В, ВВ, ВС, ВО, ВВ при числе отсчетов 200.
Для каждого типа мозаики и каждой размерности было построено тысяча изображений, для каждого из которых было получено распределение количества сегментов. На рис. 3 приведены примеры распределения числа сегментов для мозаик типа ВС, ВБ при количестве отсчетов 60.
Для каждого типа мозаик при различном числе отсчетов была проведена проверка гипотезы о нормальном распределении числа сегментов. Проверка проводилась по критерию согласия Пирсона с коррекцией Йетса [6]:
где к - количество групп, оцениваемое по формуле Стерджеса [7]; м] - экспериментальное количество значений в у-й группе; р - теоретическая вероятность попадания значений в у-ю группу при нормальном распределении; п - размер выборки.
Гипотеза о нормальном распределении Н0 считалась истинной, если выполнялось неравенство
ё < х2(а, к - 2),
Количество отсчётов
"П--ГГ-ГГ--Г“
12 3 4
П
О
1
Рис. 1. Моменты восстановления рекуррентного точечного потока в растре изображения
Рис. 2. Примеры сегментов мозаик типов F, FB, FC, FD, FE
где а - уровень значимости, (к-2) - число степеней свободы.
При а=0,05, к=11 квантиль распределения хи-квадрат:
Х2(а, к - 2) = х2(0,05, 9) = 16,91.
Гипотеза о нормальности распределения Н0 нашла подтверждение в 73 % случаях для мозаик типа В, в 84 % - для мозаик типа ВВ, в 95 % -для мозаик типа ВС, в 78 % - для мозаик типа ВБ, в 89 % - для мозаик типа ВВ.
300.
250.
200.
150,
100.
50
Мозаика типа РБ
теоретическая экспериментальн ая
300.
250.
200.
150.
100.
50
Мозаика типа РО
1338 3 5 8 3 78
300
250
200
150.
100,
50
1418 438 458 478 498 1518 538 559
Количество сегментов
666 675 685 694 704 713 723 732 742 751 761 771
Количество сегментов
Мозаика типа РО
теоретическая
экспериментальная
250 258 266 274 282 290 298 306 314 322 330 339
Количество сегментов
Рис. 3. Распределение количества сегментов мозаик F, FB, FC, FD, FE
550.
500.
450.
400.
350.
300.
250,
200.
150.
100.
50
'(п)
. п
550
500
450
400
350
300
250
200
150.
100,
50
7 9 100
Количество отсчётов
10 30
• - Тип Р
• - Тип РБ (О(п) /10)
• - Тип РО (О(п) /10)
Рис. 4. Зависимость среднего числа сегментов от количества отсчетов
'(п)
п
10 30
• - Тип РО ■ - Тип РЕ
7 9 100
Количество отсчётов
5
На основании полученных данных было проведено исследование зависимости среднего числа сегментов от размерности изображения. Данные показывают, что с увеличением размерности изображения среднее число сегментов возрастает по квадратичному закону. На рис. 4 приведены зависимости для мозаик всех типов.
Для записи этой зависимости в аналитическом виде была применена аппроксимация результатов по методу наименьших квадратов. Коэффициенты квадратного полинома вычислялись по данным восьми размерностей изображений: 30, 40, 50, 60, 70, 80, 90, 100. Данные размерностей 150 и 200 использовались для вычисления ошибок аппроксимации. В частности, для мозаики типа ВС зависимость среднего числа сегментов от числа отсчетов изображения СВС(п) определяется соотношением
СРС (п) = 0,3646 • п2 + 3,4536 • п - 48,5357.
Погрешность экстраполяции для числа отсчетов 150 и 200 составила менее 1 %.
Подобные зависимости среднего _числа сегментов от числа отсчетов изображения С(п) получены для мозаик других типов. Квадратичный характер зависимости сохраняется во всех случаях. Погрешности аппроксимации не превышают 5 %. В табл. 1 приведены аналитические зависимости среднего числа сегментов С(п) от количества отсчетов п для всех типов мозаик, а также относительные погрешности аппроксимации этих зависимостей при числе отсчетов 150 и 200.
Таблица 1. Зависимость среднего числа сегментов от количества отсчетов
Тип мозаи- ки Зависимость среднего значения числа сегментов С(п) от количества отсчетов п Относительная погрешность, %
при 150 отсчетах при 200 отсчетах
F 0,069п2+0,8589п-0,1488 1,6 1,89
FB 0,1695п2+1,9685п-2,7798 1,0 1,36
FC 0,3646п2+3,4536п-48,5357 0,31 0,89
FD 0,0667п2+0,6476п-7,9881 2,6 3,0
FE 0,1366п2+0,5065п-9,0655 3,15 3,91
Анализ полученных данных позволил построить зависимости среднеквадратических отклонений количества сегментов от числа отсчётов т(п). Данные показывают, что с увеличением размерности изображения СКО числа сегментов возрастает по линейному закону. На рис. 5 приведены зависимости т(п) для всех типов мозаик и максимальные погрешности аппроксимации на всём промежутке количества отсчётов.
Для записи зависимости в аналитическом виде был применен аналогичный подход, описанный выше. Отличие заключалось в выборе аппроксимирующей функции. Результаты, полученные в ходе моделирования, показывают, что исследуемая зависимость носит линейный характер. В частности, для мозаики FC зависимость СКО числа сегментов (TFC(n) от количества отсчетов изображения n описывается соотношением
tfc (n) = 0,5661 • n + 1,069.
Рис. 5. Зависимость СКО числа сегментов от количества отсчетов
Погрешность аппроксимации для числа отсчетов 150 и 200 составила менее 0,6 %.
Зависимости СКО числа сегментов <г(п) от количества отсчетов изображения п получены для мозаик других типов и представлены в табл. 2. Линейный характер зависимости сохраняется во всех случаях. Погрешности аппроксимации для всех типов мозаик не превышают 5 %.
Таблица 2. Зависимость СКО числа сегментов от количества отсчетов
Зависимость СКО Относительная
Тип мозаики числа сегментов погрешность, %
и(п) от количества при 150 при 200
отсчетов n отсчетах отсчетах
F 0,2822n+1,8142 1,87 1,15
FB 0,2893n-0,4766 2,66 2,07
FC 0,5661n+1,0690 0,14 0,57
FD 0,2206n+0,3159 0,5 3,74
FE 0,4005n+1,5842 1,9 0,14
Заключение
По критерию согласия Пирсона проверена гипотеза о нормальном законе распределения числа
сегментов изображений, построенных надвумер-ном рекуррентном точечном потоке восстановления. Показано, что в 73 % случаев гипотеза является истинной для мозаики типа В, и в более 78 % -для мозаик других типов. Гипотеза не получила подтверждение для изображений небольших размеров (в 30 и 40 отсчетов). При возрастании размеров изображения распределение числа сегментов стремится к нормальному закону для мозаик всех типов.
Получены зависимости среднего числа сегментов от количества отсчетов изображений для мозаик типов В, ВВ, ВС, ВБ, ВВ и установлено, что зависимость имеет квадратичный характер. Проведена аппроксимация данных зависимостей по методу наименьших квадратов. Результаты аппроксимации дают удовлетворительные результаты (погрешность менее 5 %) в случае изображений больших размерностей.
Проведена оценка поведения СКО числа сегментов от количества отсчетов изображений для мозаик типов В, ВВ, ВС, ВБ, ВВ. Анализ показал, что зависимость имеет линейных характер, а погрешность аппроксимации не превышает 5 %.
СПИСОК ЛИТЕРАТУРЫ
1. Гнеденко Б.В., Коваленко И.Н. Введение в теорию массового обслуживания. - М.: Наука, 1966. - 431 с.
2. Дубинин Д.В., Лаевский В.Е., Кочегуров А.И. Методика моделирования случайных яркостных полей, аппроксимированных однородными, одноуровневыми марковскими цепями // Проблемы информатики. - 2011. - № 4 (12). - С. 128-136.
3. Фамилия И.О. Название статьи // Техническое зрение в системах управления мобильными объектами-2010: Труды научнотехнической конференции. Вып. 4 / под ред. Р.Р. Назирова. -М.: КДУ, 2011. - 328 с.
4. Лаевский В.Е. Алгоритм построения одноуровневых марковских полей // Известия Томского политехнического университета. - 2006. - Т. 309. - № 8. - С. 32-36.
5. Dubinin D., Geringer V. (Laevski V.), Kochegurov A. Eine Methode zur Erzeugung stochastischer Helligkeitsfelder durch homogene, einstufige Markoff-Ketten // Oldenbourg Wissenschaftsverlag, Technisches Messen. - 2012. - № 5. - S. 271-276.
6. Yates J. Contingency table involving small numbers and the £ test // Royal Statistical Society. - 1934. - V. 1. - P. 217-235.
7. Sturges H.A. The choice of a class interval // Journal of the American Statistical Association. - 1926. - № 21. - P. 65-66.
Поступила 25.09.2012 г.