ИСПОЛЬЗОВАНИЕ СВОЙСТВ МАТРИЦ НАД КОНЕЧНЫМИ ПОЛЯМИ ДЛЯ ОБНАРУЖЕНИЯ

НЕУСТОЙЧИВЫХ ЦИКЛОВ И НЕПОДВИЖНЫХ СОСТОЯНИЙ ДВОИЧНЫХ ДИНАМИЧЕСКИХ СИСТЕМ

А. А. Мельников, Е. В. Рукуйжа, А. В. Ушаков

Введение. Постановка задачи.

Рассматривается линейная двоичная динамическая система [1–3], векторноматричное описание которой задается в формах рекуррентной

последовательность, Y – m мерная выходная последовательность, A – (nxn) матрица состояния, B – (nxr) матрица входа, C – (mxn) матрица выхода, H – (mxr) матрица "вход-выход", V– (k) вектор стратегии управления состоянием ЛДДС, задаваемый на первых к – тактах в форме

В дальнейшем используется для оценки свойств пары матриц (A,B) матрица управляемости (5) для случая k = n, в силу чего она представляется выражением

Ставится задача обнаружения неустойчивых циклов и неподвижных состояний двоичных динамичных систем вида (1) - (3) на основе анализа свойств пары матриц (AB). Анализ опирается на то обстоятельство, что часть свойств матриц над конечными полями идентичны свойствам матриц над бесконечными [4], а другая часть свойств порождена модулярной природой конечных полей.

Свойства матриц состояния и входа над конечными полями

Анализ свойств матрицы состояния проведем опираясь на понятие матричной функции от матрицы с учетом модулярной специфики конечного поля, а анализ свойств матрицы B, в котором в основном ограничимся случаем r = 1, будем проводить на предмет установления факта принадлежности матрицы В геометрическому спектру собственных векторов матрицы А, а также установления свойства управляемости пары

(A, B).

Основные результаты изложены в виде системы определений, утверждений и их доказательств.

243

Определение 1 (О.1) Матричной функцией f ( A) , представляющей собой (nxn)

матрицу от (nxn) матрицы А, заданных над конечным полем GF(2) называется матричный ряд, порождаемый рядом от скалярной переменной α вида

Доказательство утверждения опирается на положение теоремы Гамильтона–Кэли [4], сформулированное для случая характеристического полинома матриц над бесконечным полем [1]. ■

Определение 2 (О.2) Квадратная (nxn) матрица A обладает свойством принадлежности показателю μ , где μ - целое положительное, если выполняется

условие

Определение 3 (О.3) Матрица A размерности (nxn) называется нильпотентной с индексом нильпотентности ν , где ν - целое положительное, если выполняется условие

Доказательство утверждения строится на свойстве (10) матрицы обнулять свой характеристический полином. ■

Утверждение 4 (У.4) Индекс нильпотентности ν (nxn)–нильпотентной матрицы

244

Доказательство утверждения строится на канонической форме нильпотентной матрицы. Если каноническая форма матрицы A представляет собой (nxn) клетку Жордана [4] с нулевыми диагональными элементами, то индекс нильпотентности ν ее равен значению n и является максимальным, что доказывает справедливость правой части неравенства (17). Если отсечением согласованных по размеру левых нулевых столбцов и нижних нулевых строк получается Жорданова клетка размера ν <n, при этом индекс нильпотентности ν такой матрицы устанавливается прямым возведением в степень ν матрицы A, то очевидно, что минимальный размер Жордановой клетки с нулевыми диагональными элементами будет равен двум, что доказывает справедливость левой части неравенства (17).

■

Утверждение 5 (У.5) Если матрица входа B имеет размерность (nx1) и при этом она является собственным вектором матрицы состояния A ЛДДС (1), то пространство управляемости ЛДДС имеет размерность равную единице. □

Доказательство утверждения строится на том, что размерность пространства управляемости системы (1) определяется рангом ее матрицы управляемости

Воспользуемся свойством собственного вектора матрицы A. Если некоторый вектор ξ является собственным вектором матрицы A, то над конечным полем

Подставим в (19) вместо ξ матрицу-столбец B , тогда на основе (19) получим

матричные равенства

AB = B;

Матрица управляемости вида (21) обладает тем свойством, что ранг ее равен единице. ■

Утверждение 6(У.6) Пусть ЛДДС (1) характеризуется передаточной функцией

245

Доказательство утверждения строится на применении к векторно–матричному описанию (1) ЛДДС прямогоD - преобразования [1,2] при нулевом начальном ее состоянии, тогда с использованием свойств D - преобразования получим

где X (d) – Д образы состояния x(k) входной последовательности U (k) и выходной

Нетрудно видеть, что (27) представляет собой передаточную функцию ЛДДС, записанную в мультипликативной форме так, что

Основной результат

Для решения задачи, вынесенной в название статьи, введем и сформулируем дополнительные определения и утверждения.

Определение 4(О.4) ЛДДС (1) называется асимптотически устойчивой при

Утверждение 7 (У.7) Если матрица A состояния ЛДДС (1) обладает индексом нильпотентности ν , то эта двоичная динамическая система будет асимптотически

По условию утверждения матрица A является нильпотентной с индексом нильпотентности ν , в силу чего для k ≥ν , выполняется равенство

тому, что если матрица состояния не является нильпотентной, то такая ЛДДС не обладает свойством асимптотической устойчивости, в ней возможны неподвижные состояния, замкнутые неустойчивые циклы и особые режимы поведения в случае, если матрица B оказывается собственным вектором матрицы A.

Определение 5 (О.5) Пусть ЛДДС (1) характеризуется (nxn) матрицей состояния

будет называться неустойчивым, если условие

246

Условие (35) может быть реализовано следующими двумя способами, которые определим утверждениями.

Утверждение 8(У.8) ЛДДС (1), с ( nxn) матрицей состояния A, при любой ее реализации всегда имеет в качестве неподвижного состояния нулевое состояние

реализации имеет в качестве неподвижного состояния состояние, которое является ее

Необходимо сделать следующие примечания.

Примечание 1 (П.1)

Нетрудно видеть, что условию утверждения У.9 удовлетворяет любой вектор ξ ,

принадлежит матрица A.

Если характеристический полином D(λ) =det(λI + A) является неприводимым, то, как известно [5], показатель μ в этом случае связан с размерностью n соотношением

Таким образом, в случае удовлетворяющем (40) все множество мощностью 2n состояний ЛДДС (1), (36) разбивается на два подмножества, где одно представлено нулевым неподвижным состоянием, а второе образовано μ состояниями, где μ

удовлетворяет (40), образующими замкнутый неустойчивый цикл.

Если характеристический полином D(λ) , является приводимым, то в зависимости от конкретной его реализации μ принимает одно из значений удовлетворяющее неравенству

n ≤μ < 2n −1

247

В этом случае 2n возможных состояний ЛДДС (1), (36) образует сумму подмножеств. Первое из них образуется нулевым неподвижным состоянием. Второе–

В ключе поставленной проблемы обнаружения неподвижных состояний и замкнутых неустойчивых циклов, рассмотрим теперь поведение ЛДДС (1) для случая ненулевой входной последовательности. При этом выделим ситуацию, когда матрица B размерности (nx1) является собственным вектором матрицы A. Для этой ситуации

сформулируем следующее утверждение.

Утверждение 10 (У.10) Пусть входная последовательность представляет собой унитарный код, так что

Тогда, под действием такой входной последовательности в пространстве

Примеры

В качестве примеров, иллюстрирующих основные положения работы, рассмотрим три версии ЛДДС третьего порядка (n =3) со следующими реализациями матриц

Анализ информации, представленный (47) показывает, что матрица A (47) обладает неприводимым характеристическим полиномом D(λ) поэтому она

принадлежит показателю μ = 2n −1=7 . Как следствие, структура состояний ЛДДС с такой матрицей при U (k) ≡ 0 представлена двумя подмножествами. Первое

подмножество состоит из нулевого состояния, являющегося неподвижным, второе подмножество образует неустойчивый цикл из оставшихся семи состояний. Матрица A

248

(47) с матрицей B (50) образует полностью управляемую пару, так что матрица управляемости

W =[B | AB | A2 B]= row {[111]T , [110]T , [101]T }

обладает рангом равным 3, что соответствует выполнению условия управляемости rangW = n . Установленное свойство полной управляемости ЛДДС гарантирует

соответствующим выбором входной последовательности перевод системы из любого начального состояния в любое конечное за фиксированное число тактов.

Анализ информации представленный (48) обнаруживает, что матрица A принадлежит показателю μ = 4 . Структура пространства ЛДДС для этого случая при

состояниями образующими замкнутый неустойчивый цикл с периодом равным 4. Третье подмножество образовано неподвижным состоянием x =[111]T , которое

является собственным вектором матрицы A. Четвертое подмножество представлено состояниями, образующими замкнутый неустойчивый цикл с периодом Т = 2, каждое

из которых x(k) =[010]T и x(k) =[101]T таково, что для них выполняются соотношения x(k + 2) = A2 x(k) = x(k) , является собственным вектором матрицы A2 .

A2 = col{[001], [111], [100]} ξ1 =[010]T ,ξ2 =[101]T

Матрица A (48) и матрица B (50) образуют не полностью управляемую пару, т.к. матрица B (50) является собственным вектором матрицы A, как следствие при входной последовательности в виде унитарного входа U (k) =1(k) в структуре состояний ЛДДС

возникает неустойчивый вынужденный цикл периода 2, содержащий нулевое состояние x =[000]T и состояние x = B = [111]T .

Анализ информации представленный соотношением (49) показывает, что матрица A ЛДДС принадлежит показателю μ =3. Как следствие, при U (k) ≡ 0 структура

пространства ЛДДС представлена четырьмя подмножествами. Первое подмножество образовано нулевым неподвижным состоянием. Второе и третье подмножества представлены неустойчивыми неподвижными циклами периода T = μ =3 . Четвертое

подмножество, представленное неподвижным состоянием x =[111]T , которое является

собственным вектором матрицы A. Пара матриц (49), (50) является не полностью управляемой, т.к. B-собственный вектор A. Как следствие, при входной последовательности в виде унитарного кода в поведении ЛДДС обнаруживается

неустойчивый вынужденный цикл периода Т = 2 на состояниях x =[000]T и x =[111]T .

В заключении необходимо отметить, что полученная на основе анализа свойств матриц A и B в рассматриваемых примерах структура неподвижных состояний и неустойчивых циклов также обнаруживается таблично-графическим способом с использованием таблиц и графов переходов ЛДДС [6].

Литература

1.Гилл А. Линейные последовательностные машины. М.: Наука, 1974.

2.Фараджев Р. Г. Линейные последовательностные машины. М.: Сов. радио, 1975.

3.Бохман Д., Постхофф Х. Двоичные динамические системы. М.: Энергоатомиздат, 1986.

4.Гантмахер Ф. Д. Теория матриц. М.: Наука, 1988.

5.Питерсон У., Уэлдон Э. Коды, исправляющие ошибки. / Пер. с англ.М.: Мир, 1976.

6.Алгебраические методы в теории устройств дискретной автоматики и телемеханики // Труды лаборатории телемеханики СПб ГИТМО (ТУ) / под ред. А. В. Ушакова. СПб: СПбГИТМО (ТУ), 2001.

249