ИНФОРМАТИКА
УДК 519.713,519.718.7
АВТОМАТНАЯ ИНТЕРПРЕТАЦИЯ ЦЕЛОЧИСЛЕННЫХ ПОСЛЕДОВАТЕЛЬНОСТЕЙ
А.С. Епифанов
Саратовский государственный университет,
кафедра дискретной математики и информационных технологий
E-mail: [email protected]
Преобразование фазовых картин в геометрические образы законов функционирования автоматов, предложенное и разработанное В.А. Твердохлебовым, позволило представлять фазовые картины едиными математическими структурами -- ломаными линиями с числовыми координатами точек. В.А.Твердохлебовым показано, что последовательность элементов из конечного множества, совмещенная с линейным порядком на множестве входных слов, определяет законы функционирования дискретной детерминированной динамической системы (автомата). Это позволяет проводить анализ законов функционирования автоматов на основе исследования свойств числовых последовательностей, а также исследовать свойства числовых последовательностей на основе анализа свойств автоматов.
Ключевые слова: конечный детерминированный автомат, геометрический образ автомата, целочисленные последовательности.
The Automata Interpretation of Integer Sequences
A.S. Epifanov
Saratov State University,
Chair of Discrete Mathematic and Information Technologies E-mail: [email protected]
Transformation of phase pictures to geometrical images of laws of functioning of state machines, offered and developed by V.A.Tverdokhlebov, has allowed to represent phase pictures by uniform mathematical structures -- broken lines with numerical coordinates of points. V.A.Tverdokhlebov shows, that sequence of elements from the finite set, combined with linear order on set of input words, defines laws of functioning of the discrete determined dynamic system (state machine). It allows carry out the analysis of laws of functioning of state machines on the basis of research of properties of numerical sequences, and also it allows to research properties of numerical sequences by the analysis of properties of state machines.
Key words: finite state machine, geometrical image of state machine, integer sequences. ВВЕДЕНИЕ
Методы решения задач управления, синтеза и анализа дискретных детерминированных динамических систем используют задание законов функционирования систем. В фазовых картинах систематизированы и представлены фазовые траектории, определяющие поведение систем. Большое разнообразие фазовых картин в работе [1] сведено к единой форме — символьным и числовым графикам, точки которых размещаются на геометрических кривых линиях. Это позволяет не только сводить фазовые картины к удобной стандартной форме, но и использовать для постановок и решения задач средства непрерывной математики. Существенной оказывается возможность доопределения фазовых картин, представленных частично заданными геометрическими образами, с использованием классических ме-
тодов интерполяции. При задании автомата геометрическим образом соответствующее автоматное отображение преобразуется в точечный числовой график с точками, расположенными на геометрических кривых линиях (см. [1]). При таком задании законов функционирования автомата мощность множества состояний может быть бесконечной или, для удобства математических вычислений, полагаться бесконечной. В данной работе содержатся результаты исследований законов функционирования дискретных детерминированных динамических систем (автоматов), представленных в форме числовых последовательностей (последовательностей вторых координат точек геометрических образов) длиной до 1000000 знаков. Оценка сложности и классификация математических структур в форме последовательностей производится на основе спектра динамических параметров рекуррентного описания последовательностей.
Данная работа является продолжением работ [2-3], в которых осуществлены построение и анализ классов автоматов, законы функционирования которых заданы в форме числовых последовательностей, и содержат результаты более детального исследования 10-ти математических величин: п, е, ф (так называемое золотое сечение), лД, 32, 1п(2), 1п(10), С(3), константы Каталана и константы Эйлера. Исследуется зависимость числа состояний у минимального автомата, построенного по последовательности, от мощности входного алфавита автомата и способа доопределения функции переходов автомата.
1. ГЕОМЕТРИЧЕСКИЕ ОБРАЗЫ ЗАКОНОВ ФУНКЦИОНИРОВАНИЯ АВТОМАТОВ
Преобразование символьной формы автоматной модели в числовую структуру (геометрический образ законов функционирования автомата) включает линейное упорядочивание автоматного отображения рв = и {(р, А(з,р))} для инициального автомата Лв = (Б, X, У, 5, Л, 5), где Б, X и У — соот-
рех *
ветственно множества состояний, входных и выходных сигналов, 5 : Б хX ^ Б — функция переходов, А : Б х X ^ У — функция выходов, а 5 е Б — начальное состояние автомата. Автоматное отображение рв взаимнооднозначно преобразуется в автоматное отображение вида р'3 = и {(Р,А'(з,р))}, где
рех*
А'(з,р) — последний знак последовательности Л(5,р). После введения на множестве X* линейного порядка « множество пар р3 оказывается линейно упорядоченным множеством (рв, «1), которое при любом линейном порядке «2 на множестве Y становится графиком 7в = (р5,<1 , «2) в системе координат с осью абсцисс (X*,<1) и осью ординат (У,<2). Для преобразования множества пар рв и р'3 в графики принципиально важным является задание линейного порядка <1 на множестве X* (см., например, [1]).
Правило 1. На множестве X вводится линейный порядок <1 ): х1 х2 ... х&.
Правило 2. Порядок <1 на X распространяется до линейного порядка на множестве X*:
- для любых слов р15р2 е X* неодинаковой длины |р11 < |р21 ^ р1 р2;
- для любых слов р1 ,р2 е X*, для которых |р11 = |р21 и р1 = р2, их отношение по порядку <1 повторяет отношение ближайших слева несовпадающих букв в словах р1 и р2.
На множестве У = {у1 , у2,... } определяется некоторый линейный порядок «2. Упорядоченные множества пар (рв ,<1) и (р^,<1) дополняются линейными порядками «0 на У * и «2 на У. В результате получаем графики (р5,«1,«0) и (рв,«1 ,«2). Построенные графики размещены в системе координат с осью абсцисс (X*,«1) и осями ординат соответственно (У*,«0) и (У,«2). Замена элементов р е X*, р = Хг1 Хг2 ...Хгк , и у е У в графике = (рв ,«1 ,«2) их номерами г1(р) по порядку «1 и г2(у) по порядку «2 позволяет преобразовать символьный график 7в в числовой график в системе координат с осью абсцисс N + и осью ординат {1, 2,...,/}, где |У | = /.Возможны и другие варианты линейных порядков на X* (см., например, [1]).
В данной работе исследование законов функционирования автоматов проводится с использованием порядка «1. Из геометрического образа 7в автомата Лв выделяется последовательность вторых координат точек геометрического образа, которая взаимнооднозначно соответствует полному геометрическому образу (при выбранных множестве входных сигналов X и линейном порядке «1 на множестве входных слов X*). Произвольная последовательность элементов из конечного множества может рас-
сматриваться как последовательность вторых координат точек геометрического образа и, следовательно, как задание законов функционирования автомата. Это позволяет некоторые свойства законов функционирования автомата представлять свойствами последовательностей.
2. СПЕКТР ДИНАМИЧЕСКИХ ПАРАМЕТРОВ РЕКУРРЕНТНОГО ОПРЕДЕЛЕНИЯ ЧИСЛОВЫХ ПОСЛЕДОВАТЕЛЬНОСТЕЙ
В работе [1] предложен спектр динамических параметров, характеризующих сложность правил порождения последовательности. Пусть и = {щ, и2} — конечное множество и £ — последовательность элементов из множества и : £ =< и(1), и(2),..., и(£),... >. Спектр
0 =< О0, ,О2, О3, О4 > вводится как многоуровневая структура, в которой на каждом уровне представлены наборы характеристик использованных рекуррентных форм (г1 ,г2,..., ¿то) = ¿т+1, где т — порядок рекуррентной формы, т=1,2,... . Каждая рекуррентная форма с учетом отношения £ е и *, описывается функцией /т : и * ^ и и применяется для определения одного или нескольких отрезков в зависимости от структуры последовательности £. Исследования свойств спектров, соответствующих последовательностям, имеет общее значение для конструктивных объектов.
Определим понятие спектра.
Определение 1. Для любой последовательности £ е иу (|£| = V) наименьший порядок рекуррентной формы, определяющей последовательность £, будем обозначать т0(£).
Определение 2. Для любой последовательности £ е ии т е N +, где 1 < т < т0(£), наибольшую длину начального отрезка последовательности £, определяемого рекуррентной формой порядка т, будем обозначать йт(£).
Определение 3. Для любой последовательности £ е ии т е N +, где 1 < т < т0(£), число смен рекуррентных форм порядка т, требующихся при определении последовательности £, будем обозначать гт(£).
Определение 4. Для любой последовательности £ е иу и т е N+, где 1 < т < т0(£) и 3, где
1 < 3 < гт(£), длину з-го отрезка в определении последовательности £ будем обозначать й2(£).
Используя введенные обозначения, определим спектр параметров, характеризующих последовательность, как следующую структуру: О0(£) =< т0(£) >, О1 (£) =< й1 (£),й2(£),...,(£) >, О2(ё) =< Г1 (£),г2(£),...,га(£)_ >, Оз(£) =< О3(£),О3(£),...,О£(£) >, где а = т0(£) и От (£) =< (£ ),й22 (£),... ,42т (£) > (пт — номер последнего отрезка в определении последовательности £ как последовательности отрезков, определяемых отдельными рекуррентными формами порядка т). Формально О4(£) = ©(О3(£)), где © — оператор замены в О3(£) величин длин отрезков весами использованных рекуррентных форм для определения отрезков. В качестве примера в табл. 1 приведены числовые показатели на уровнях О0-О3 спектра О для последовательности £ =< 3,1,4,1, 5, 9, 2, 6, 5,3, 5, 8,9, 7, 9, 3, 2,3, 8,4,6, 2, 6,4,3,3,8, 3, 2, 7,9, 5,0, 2, 8, 8,4,1,9, 7,1,6,9,3, 9, 9, 3, 7, 5,1 > длины 50-ти знаков, задающей приближение числа п.
Таблица 1
Числовые показатели на уровнях П0-П3 спектра П для последовательности длины 50-ти знаков, задающей приближение числа п
Уровни спектра Числовые показатели спектра
По < то (О = 3 >
П1 <4(0 = 4,4(0 = 23,4(£) = 50 >
П2 <г1 (0 = 14, г2 (0 = 3, г3 (0 = 1 >
Оз < П1 (0 =< ¿1 (0 = 3,$(0 = 5,4(0 = 2,4(0 = 4,4(0 = 3, 4 (0 = 5,4 (0 = 3,4 (0 = 2,4 (0 = 6,4о(0 = 2,4^) = 5, 42(0 = 4,4з (0 = М14 (0=4 >, п2(0=< ¿2 (0 = 21,4(0 = 24,4 (0 = 3 >, О3(€)=< (0=47 >)
3. ПОСТРОЕНИЕ И АНАЛИЗ КЛАССОВ (Я, М, Я (Я))-АВТОМАТОВ
В работе [1] предложены и разработаны методы синтеза автомата по последовательностям и геометрическим кривым. Последовательность элементов из конечного множества, совмещенная с линейным порядком на множестве входных слов, определяет законы функционирования дискретной детерминированной динамической системы (автомата). В.А. Твердохлебовым предложен новый тип автомата — (Я, т, ¿(Я))-автомат. Законы функционирования данного типа автомата задаются числовой последовательностью Я, которая полагается последовательностью вторых координат точек геометрического образа. Рассматривается начальный отрезок длины ¿(Я) последовательности Я. Величина т — мощность входного алфавита автомата, количество выходных сигналов определяется спецификой начального отрезка последовательности Я длины ¿(Я) (число различных значений элементов в начальном отрезке длины ¿(Я)).
В данной работе содержатся результаты по построению и анализу конечных детерминированных автоматов, законы функционирования которых определены начальными отрезками геометрических образов и выбором числа входных сигналов автомата. Для этого из банка последовательностей [5] извлечены 12 последовательностей длины 1000000 знаков, задающих приближения следующих математиче-
те
ских величин: п, е, ф = (так называемое золотое сечение), л/2, 32, 1п(2), 1п(10), С(3) = ^ (хЗ),
х = 1
те (—1™) 11 1 константы Каталана С = ^ , константы Эйлера 7 = Нш (1 + 2 + 3 + ... + п — ^(п)),
п=0( ) п^те 23 п
Х^(характеристическая последовательность для чисел Фибоначчи) и хр (характеристическая последовательность для простых чисел). Каждая последовательность представлена набором начальных отрезков, имеющих длины ¿к знаков,где ¿к = 10000 ■ к, к = 1, 2,..., 100. Полученное множество из 1200 последовательностей рассматривается как множество начальных отрезков последовательностей вторых координат точек геометрических образов законов функционирования автоматов. Соответствующие последовательности первых координат точек геометрического образа определялись вариантами выбора числа входных сигналов автомата и линейным порядком «1 на множестве входных последовательностей. Рассматриваются множества входных сигналов, содержащие 2, 5, 10, 20, 50, 100, 200, 500, 1000 и 2000 элементов. Подробное описание метода синтеза автомата по последовательности, а также метода проверки эквивалентности состояний автомата по геометрическому образу содержится в монографии [1]. Отметим принципиальные моменты метода синтеза автомата по последовательности. Если в качестве задания законов функционирования автомата рассматривается последовательность £ =< и(1),и(2),..., и(£),...}, то в ее интерпретацию включается разбиение последовательности £ на отрезки длины т, где т = IX|. Каждый отрезок при таком разбиении определяет функционирование автомата для конкретного состояния памяти. Если множеством состояний автомата полагать множество Б = {5р}рех*, а функцию переходов 5 определить правилами 50 = и 5(зр, х) = 5рх, то функция 5 оказывается стандартной для всех автоматов с множеством входных сигналов X. Специфика автоматов проявляется в том, что на бесконечном (счетном) множестве состояний для каждого автомата выделяются классы эквивалентных состояний. Существенным при синтезе законов функционирования автомата является способ доопределения функции переходов 5 автомата. Возможно циклическое доопределение, доопределение в начальное состояние, генерация состояния псевдослучайным образом (из множества возможных состояний). В случае, когда ^Х^ = [^Х"^], где IX| — мощность входного алфавита автомата, а ¿(Я) — длина начального отрезка последовательности Я (по которой строятся законы функционирования автомата), доопределение требуется и для функции выходов А.
В данной работе доопределение функции переходов осуществляется всеми указанными способами, а значение мощности входного алфавита и длины начальных отрезков последовательностей выбраны таким образом, что ^Хр = [^хТ], поэтому доопределение функции А не требуется. Рассматриваемое множество из 1200 последовательностей при десяти различных значениях мощности входного алфавита и трех способах доопределения функции переходов сопоставляется с классом автоматов, состоящем из 36000 элементов. Далее осуществляется минимизация автоматов из построенного класса и разбиение на подклассы автоматов по числу состояний в автоматах после минимизации. Для построения семейства автоматов по заданной последовательности при различных значениях мощности входного алфавита, минимизации построенных автоматов по числу состояний на основе уже существующих
методов разработаны соответствующие алгоритмы, составлены и отлажены программы, реализующие данные алгоритмы. В результате получен программный комплекс, позволяющий по заданному в файле массиву последовательностей осуществлять построение семейств автоматов, минимизацию автоматов по числу состояний и классификацию по сложности на основе числа состояний в минимальном автомате, построенном по последовательности. В программном комплексе реализовано три указанных выше способа доопределения функции переходов автомата. В качестве примера в табл. 2 приведены табличные задания функций переходов и выходов автомата, построенного по начальному отрезку длины 60 числа п. В данном примере использовано циклическое доопределение функции переходов автомата, а значение мощности входного алфавита равно 5. В результате анализа класса п-автоматов, т. е. класса автоматов, построенных по последовательностям, задающим приближение фундаментальной математической величины п длины ¿к знаков, где ¿к = 10000 ■ к, к=1,2,... ,100, при рассмотрении различных вариантов мощности входного алфавита (IX| = 2, 5,10, 20, 50,100, 200, 500,1000, 2000) отмечено, что при некоторых из указанных способах доопределения функции переходов уменьшается число состояний у автоматов после минимизации.
Таблица 2
Табличное задание функций переходов 5 и выходов Л (п,5,60)-автомата
5 50 51 52 53 54 55 56 57 58 59 510 511
51 56 511 54 59 52 57 50 55 510 53 58
Х2 52 57 50 55 510 53 58 51 56 511 54 59
Х3 53 58 51 56 511 54 59 52 57 50 55 510
Х4 54 59 52 57 50 55 510 53 58 51 56 511
Х5 55 510 53 58 51 56 511 54 59 52 57 50
Л 50 51 52 53 54 55 56 57 58 59 510 511
Х1 У3 У9 У5 У3 У6 У3 У9 У8 У1 У9 У0 У9
Х2 У1 У2 У8 У2 У2 У8 У5 У4 У6 У3 У5 У7
Х3 У4 У6 У9 У3 У6 У3 У0 У1 У9 У7 У8 У4
Х4 У1 У5 У7 У8 У4 У2 У2 У9 У3 У5 У2 У9
Х5 У5 У3 У9 У4 У3 У7 У8 У7 У9 У1 У0 У4
При доопределении в начальное состояние из 3000 постренных п-автоматов менее чем у 100 автоматов после минимизации число состояний уменьшилось. Например, у (п,2,10000)-автомата, (п,2,50000)-автомата и (п,2,100000)-автомата число состояний после минимизации сокращается почти в 2 раза. У всех построенных п-автоматов при доопределении функции переходов на основе использования генератора случайных чисел (состояний выбирается из множества возможных состояний случайным образом) число классов эквивалентности совпадает с числом состояний, т. е. автоматы являются минимальными по числу состояний. Используемый аппарат геометрических образов автоматов позволяет рассматривать с автоматной интерпретацией как сами классические последовательности чисел Фибоначчи (и другие так называемые возвратные последовательности) и простых чисел, так и их характеристические последовательности. Автором статьи на основе автоматной интерпретации проведен анализ как самих классических последовательностей чисел Фибоначчи и простых чисел, так и их характеристических последовательностей. При построении автоматов по самим классическим последовательностям отмечено, что такие автоматы являются минимальными по числу состояний (при анализируемых значениях мощности входного алфавита и рассмотренных начальных отрезках таких последовательностей). Кроме того, характеристики данных последовательностей по спектру О (на рассмотренных в работе первых четырех уровнях О0-О3) также являются «тривиальными» (так как формально для определения рекуррентной формой последовательности, в которой все элементы различны (например,последовательности простых чисел), достаточно порядка рекуррентной формы, равного единице, а характеристики на первом, втором и третьем уровнях спектра «вырождаются» из наборов чисел в одно число на каждом уровне). Ввиду ограничений на объем статьи приводятся результаты исследования не самих классических последовательностей, а только характеристических последовательностей (как более интересных, по мнению автора, с точки зрения полученных результа-
тов их анализа). Проведенный анализ класса (хр,т,^(хр))-автоматов, где хр — характеристическая последовательность чисел Фибоначчи, показал, что и при циклическом доопределении функции переходов, и при доопределении в начальное состояние, а также и при доопределении функции переходов на основе использования генератора случайных чисел возможно уменьшение числа состояний у автомата после минимизации. Например, при доопределении в начальное состояние у (хр,2,10000), (хр,2,20000), (хр,2,30000), (хр,2,40000), (хр,2,50000) и (хр,2,100000)-автоматов после минимизации число состояний уменьшается более чем в 60 раз. Кроме того, проведено дополнительное исследование автоматов из класса (хр, т,^(хр))-автоматов, где хр — характеристическая последовательность простых чисел при значениях ¿(хр) < 100 и циклическом доопределении функции переходов.
Отмечены следующее свойства: (хр, 2,60)-автомат после минимизации имеет 19 состояний (число состояний до минимизации 30), (хр, 2, 70)-автомат после минимизации имеет 35 состояний (число состояний не уменьшилось после минимизации), (хр, 2,80)-автомат после минимизации имеет 31 состояние (число состояний до минимизации 40), (хр, 2,90)-автомат после минимизации имеет 45 состояний (число состояний не уменьшилось после минимизации), (хр, 2,100)-автомат после минимизации имеет 38 состояний (число состояний до минимизации 50). Аналогичное свойство отмечено при построении автоматов (при циклическом доопределении функции переходов) по характеристической последовательности чисел Фибоначчи при ¿(хр) < 100. На основе этих данных выдвигается гипотеза, что в случаях ¿(хр) < 100, когда при построении автомата используется циклическое доопределение функции переходов и величина ^ХР^ является четным числом, после минимизации число состояний у автомата уменьшается, если же ^ является нечетным числом, то число состояний после минимизации не уменьшается.
4. ОЦЕНКА СЛОЖНОСТИ И КЛАССИФИКАЦИЯ ЗАКОНОВ ФУНКЦИОНИРОВАНИЯ (Я, М, £(Я))-АВТОМАТОВ
В данной работе на основе использования аппарата геометрических образов автоматов и спектра динамических параметров, предложенных и разработанных В.А. Твердохлебовым [1], осуществляется оценка сложности и классификация законов функционирования дискретных детерминированных динамических систем (автоматов). Сравнение по сложности таких автоматов может быть сделано на основе сравнения математических структур, представляющих специфику законов функционирования автомата. В качестве таких структур могут быть использованы геометрические образы законов функционирования автоматов. Геометрический образ построен на основе, во-первых, линейного упорядочивания элементов автоматного отображения, а во-вторых, извлечения из линейно упорядоченного автоматного отображения последовательности вторых координат точек [1]. Сравнение по сложности законов функционирования автоматов производится на основе сравнения спектров последовательностей. На основании того, что на структуру характеристических последовательностей ограничения не накладываются и она может быть любой последовательностью элементов из конечного множества, метод сравнения моделей по сложности допускает интерпретацию на любых, в частности, числовых последовательностях. Уровни спектра построены иерархически по добавлению в спектр новых параметров. В связи с этим эквивалентность последовательностей представляется через ¿-эквивалентности, определяемые уровнем 0^,0 < £ < 3.
В результате вычисления и анализа спектров для рассматриваемых математических последовательностей построены классы эквивалентных по сложности последовательностей и соответствующие им классы эквивалентных по сложности законов функционирования автоматов. Например, среди классов (п, т, 1000)-автоматов, (е,т, 1000)-автоматов, (хр, т, 1000)-автоматов и (хр, т, 1000)-автоматов с точки зрения используемого спектра О самым сложным оказался класс (хр,т, 1000)-автоматов, а самым простым класс (п,т, 1000)-автоматов. При этом на нулевом уровне 00 спектра О отмечена следующая специфика: в классе (п,т, 1000)-автоматов значение 00 =< 6 >, в классе (е,т, 1000)-автоматов 00 =< 7 >, в классе (хР,т, 1000)-автоматов 00 =< 60 >, а в классе (хР,т, 1000)-автоматов 00 =< 377 >. Также отмечено, что при увеличении параметра ¿(Я) в классах (п,т,^(п))-автоматов и (е,т,^(е))-автоматов от 1000 до 10000 сложность с точки зрения спектра динамических параметров существенно не увеличивается, кроме того, и в классе (п,т,10000)-автоматов О0 =< 8 >, и в классе (е,т,10000)-автоматов 00 =< 8 >, т.е. при таком значении величины ¿(Я) данные клас-
сы автоматов имеют одинаковую сложность. Подробное описание проведенного исследования всех 150000 целочисленных последовательностей из банка [5] содержится в монографии [2].
ВЫВОДЫ
Использование геометрических образов позволяет сопоставлять последовательности (а соответственно и той величине, приближение которой она задает) дискретную детерминированную динамическую систему (автомат) и изучать свойства последовательностей, а также объектов, процессов и явлений, представленных в виде последовательностей, на основе анализа свойств автоматов.
В работе построены и проанализированы классы автоматов, законы функционирования которых представлены в форме последовательностей из банка целочисленных последовательностей [5]. Проведена оценка сложности последовательностей на основе числа состояний в минимальном автомате, построенном по последовательности. Для рассматриваемых последовательностей построены и проанализированы спектры числовых показателей. На основе совпадения показателей на соответствующих уровнях спектра проведена оценка сложности и классификация последовательностей (и соответствующих им законов функционирования автоматов).
Библиографический список
1. Твердохлебов, В.А. Геометрические образы законов функционирования автоматов / В.А. Твердохлебов. -Саратов, 2008. - 183 с.
2. Епифанов, А.С. Анализ фазовых картин дискретных динамических систем / А.С. Епифанов. - Саратов, 2008. - 156 с.
3. Епифанов, А.С. Анализ геометрических образов законов функционирования автоматов / А.С. Епифанов //
Управление большими системами. - 2009. - Вып. 24. - С. 81-98.
4. Епифанов, А.С. Построение и анализ классов (H, m,d(H))-автоматов / А.С. Епифанов // Управление большими системами: V Всерос. школа-семинар молодых ученых. - Липецк, 2008. - С. 23-30.
5. The On-Line Encyclopedia of Integer Sequences [Электронный ресурс]. URL: www.research.att.com (дата обращения: 11.09.2010).