О решении линейных m-разностных уравнений порядков, меньших m Текст научной статьи по специальности «Математика»

Sorokin S.P. POSITIONAL NECESSARY OPTIMALITY CONDITIONS AND NUMERICAL METHOD FOR DISCRETE OPTIMAL CONTROL PROBLEMS IN LINEAR SYSTEM

The paper is devoted to new numerical method for discrete optimal control problems in system which is linear wrt the state variable. The method is based on the positional necessary optimality conditions in Pontryagin form. The designed method proved itself quite effective in numerical analysis of optimal control problems with bilinear dynamical systems.

Key words: optimal control; discrete linear systems with controlled coefficients; necessary optimality conditions; feedback controls; numerical method.

УДК 517.929.2

О РЕШЕНИИ ЛИНЕЙНЫХ т -РАЗНОСТНЫХ УРАВНЕНИЙ ПОРЯДКОВ,

МЕНЬШИХ т

© Д.Н. Спичкин

Ключевые слова: т -разностные уравнения; функции Виленкина - Крестенсона.

Получена структура решения линейных т -разностных уравнений с постоянными коэффициентами порядка д<т при комплексных решениях характеристического уравнения.

Обозначим через [а, Ь] = {а, а + 1,... ,Ь} упорядоченное множество неотрицательных целых чисел. Пусть х = (х\... хп)т, р = (р1... рп)т — п -разрядные т -ичные представления этих чисел.

Определение! Операцией т -сдвига двух чисел х, р € [а, Ь], обозначаемой х 0 р,

т

назовем их поразрядную разность по модулю т.

О п р е д е л е н и е 2. Линейное разностное уравнение, где в качестве сдвига аргумента берется т -сдвиг, назовем разностным уравнением с модулярной арифметикой или линейным т -разностным уравнением.

Линейное т -разностное уравнение порядка д с постоянными коэффициентами имеет вид

у(х 0 д) + кгу(х 0 (д — 1)) + ... + кд-1у(х 0 1) + кду(х) = 0, (1)

т т т

где к^ € С, ] = 1,д.

Решениями уравнений (1) являются функции Виленкина - Крестенсона (ВКФ) [1], записываемые в форме Пэли

п

У. Рп+1-^ х^

Ра1 (р, х) = Wj=1 , (2)

- 2П

где W = вгт, г — мнимая единица, х — аргумент, р — параметр, причем х,р € [0, тп — 1].

Известно [1], что при решении линейных т -разностных уравнений часть разрядов пара-

метра р фиксируется, а часть остается произвольной.

Пусть, далее, порядок д уравнения (1) удовлетворяет неравенствам 1 ^ д<т. Тогда число д в т -ичной п -разрядной форме представления запишется как д = (0 0 . . . д)т.

Подставим выражение (2) в уравнение (1) и проведем ряд упрощений. Так как выражения (хп — д)т и (хп — (д — Ь))т в модулярной арифметике, соответственно, равняются выражениям хп + т(5о Хп + Хп + ... + 5(д-1) хп ) — д и хп + т(§0 Хп + Хп + ... + 5(д-г-1) хп ) —

2687

—(д — Ь) (где ¿ав — символ Кронекера) в обычной арифметике и, учтя, что Vа € К

п-1

Рп+1-^ ху

Wта = 1, а также разделив обе части получаемого уравнения на W^=1 =0, полу-

чаем

г-р 1д

г=1

Пусть Л = W-р1, тогда выражение (3) перепишется в виде

W-piq + ktW-pi(q-t) = 0. (3)

Aq + ^2 kt\q-t = 0. (4)

t=i

Это уравнение является характеристическим уравнением уравнения (1) в случае 1 ^ д<т .

Решением исходного уравнения (1) будет ВКФ (2), в которой разряд р1 параметра р фиксируется из характеристического уравнения (4) (если учесть, что Л = W-р1 ),аосталь-ные разряды — произвольны.

Заметим также, что в силу равенства \Ш'-р1| = 1, необходимо имеем |Л| = 1 для всех значений Л (как вещественных, так и комплексных). Поэтому справедлива

Л е м м а 1. Решения Лг характеристического уравнения (4) для т -разностного уравнения (1) порядка д<т с постоянными коэффициентами ^, ] = 1,д, должны удовлетворять условию |ЛГ | = 1.

Для вещественных корней характеристического многочлена справедлива [2]

Т е о р е м а 1. При вещественном однократном решении Лг характеристического уравнения (4) необходимыми условиями существования решения уравнения (1) в базисе ВКФ (2) при д<т являются условия:

1. при Лг > 0 р1 = 0 , параметр р имеет представление: р = (0 р2 р3 ... рп).

2. при Лг < 0 р1 = т , параметр р имеет представление: р = ( т р2 р3 ... рп), т — четное.

Обратимся к случаю простых комплексных корней характеристического многочлена.

Пусть Лг = аг + гЬг € С, |ЛГ| = 1, является корнем уравнения (4). Тогда

2п

ЬпЛг = г(агг(аг + гЬг) + 2пк) = —г—р1, к € 2,

т

откуда

т

/ m \

p1 = I —— arg(ar + ibr) — mk I € [0, m — 1].

V 2n )

Т е о р е м а 2. Пусть уравнение (1) имеет однократное комплексное решение \r характеристического уравнения (4). Тогда для существования решения в базисе ВКФ необходимо, чтобы в параметре p = (pi p2 ■ ■ ■ Рп)т разряд p1 = (—m arg(ReAr + iImAr) — mk) и p1 € [0, m — 1], k € Z (остальные разряды могут быть произвольными).

ЛИТЕРАТУРА

1. Трахтман А.М., Трахтман В.А. Основы теории дискретных сигналов на конечных интервалах. М.: Сов. радио, 1975. 239 с.

2. Спичкин Д.Н. О разностных уравнениях на конечных интервалах произвольного порядка с постоянными коэффициентами // Материалы 9 студенческой научной конференции кафедры ПМИ ИжГТУ. Ижевск: Изд-во ИжГТУ, 2007. С. 44-45.

Spichkin D.N. SOLUTION OF LINEAR m -DIFFERENCES EQUATIONS ORDER LESS THAN m

The structure of decision of order q in q<m linear m -differences equations with constant coefficients for complex solutions of the characteristic equation is received.

Key words: m -differences equations; Vilenkin-Chrestenson functions.

2688

УДК 519.254

БЫСТРОЕ ВЫЧИСЛЕНИЕ ЦИКЛИЧЕСКОЙ СВЕРТКИ, ОСНОВАННОЕ

НА ФАКТОРИЗАЦИИ ТЕНЗОРОВ

© Т.М. Спичкина

Ключевые слова: циклическая свертка; дискретные сигналы; факторизация тензоров. Приведены результаты о вычислении циклической свертки дискретных сигналов, заданных тензорами первого ранга в контравариантных координатах, разложенных в произведение тензоров другой структуры, обладающих рядом эффективных вычислительных свойств.

Основные определения

Пусть V — N -мерное линейное пространство с базисом (ві), и пусть V* — сопряженное линейное пространство к V с базисом (ві). Будем рассматривать тензоры над V типа

(р, д) [1], т. е. элементы пространства (V)хр х (V*)хд . Тензор типа (р, д) называется также

р раз контравариантным и д раз ковариантным.

Тензору А1 поставим в соответствие тензор

В = АЛ, (1)

где ак, вк Є Ъ, индексы аі* — ві, а2] — в2 при г,] =0, N — 1 берутся по модулю N и принимают все значения из множества {0,1,... N — 1}. Если коэффициенты ак, вк в индексах

а1г — ві, а2] — в2 такие, что ак, вк = 0, N — 1, а г = ] = N, то в результате получим тензор вида

в2 _ л а2 N—в2 (2)

Ваі ві = АагМ—@1. (2)

Можно показать, что определение тензоров (1) и (2) — корректно [1].

Для тензора (2) определим операцию, которую назовем обобщенной операцией свертки, следующим образом:

N—1

__ Вв в2 _ \ ' ДвN—в2 или Ва2

Вві = ві =2AвN—ві, или, Ваі

в=0

где все индексы берутся по модулю N2.

По аналогии с транспонированием матриц введем операцию транспонирования тензора. Тензор (1) не допускает перестановки индексов, поэтому под транспонированием тензора по смешенным координатам будем понимать описанное далее преобразование тензора.

Пусть дан тензор А3І. Умножим исходный тензор на метрический тензор О в ковари-антных координатах, заданных в базисе (ві), получим А^ = д^кАк и умножим результат на метрический тензор в контравариантных координатах в базисе (ві), тогда А^ = дікА^. В дальнейшем под транспонированием по ковариантной координате г и контравариантной координате ] тензора А1 будем понимать введенное преобразование, обозначаемое А^ =

=А.

Вычисление свертки с применением тензорных операций

N—1

в _ \ ' дa2N—в

Ваі в =2.^ Аа^—в, в=0

2689