CC BY
38
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Meurant G. Computer solution of large linear systems. ELSEVIER, Amsterdam, 1999. - 753 p.
2. Виноградова C.A., Крукиер Л.А. Решение трехмерной стационарной задачи конвекции-диффузии в анизотропной среде методом неполного LU-рюложения // Труды XIV Молодежной конференции-школы с международным участием «Современные проблемы математического моделирования», пос. Дюрсо, 12-17 сентября 2011 г. - С. 74-90.
3. Самарский А.А. Теория разностных схем. - М.: Наука, 1977. - 656 с.
4. Самарский А.А., Николаев КС. Методы решения сеточных уравнений. - М.: Наука, 1978. - 592 с.
5. . ., . ., . . //
техника полупроводников. - 1997. - Т. 31, № 11. - С. 1291-1298.
6. Роуч П. Вычислительная гидродинамика. - М.: Мир, 1980. - 618 с.
7. . . // -
. - 1996. - 11. - . 37-46.
Статью рекомендовал к опубликованию д.ф.-м.н. АЛ. Чикин.
Виноградова Светлана Александровна - Южно-Российский региональный центр информатизации (ЮГИНФО) Южного федерального университета; email: svetlavi@sfedu.ru; 344090, г. Ростов-на-Дону, пр. Стачки, 200/1, корп. 2; тел.: +78632199734; лаборатория вычислительного эксперимента на супер-ЭВМ; младший научный сотрудник.
Vinogradova Svetlana Alexandrovna - Computer Center SFU; email: svetlavi@sfedu.ru; 200/1, Stachki av., build. 2, Rostov-on-Don, 344090, Russia; phone: +78632199734; laboratory of computational experiments on super-computers; junior research fellow.
УДК 519.254.1
Д.В. Иванов, O.B. Усков
РЕКУРРЕНТНОЕ ОЦЕНИВАНИЕ БИЛИНЕЙНЫХ ДИНАМИЧЕСКИХ СИСТЕМ С ПОМЕХАМИ НАБЛЮДЕНИЯ ВО ВХОДНОМ И ВЫХОДНОМ
СИГНАЛАХ
, -метров для динамических систем с помехами во входных и выходных сигналах. В работе предложен рекуррентный алгоритм на основе стохастической аппроксимации для оценивания параметров билинейных динамических систем с помехами наблюдения во входном и выходном. Доказана сильная состоятельность, получаемых оценок параметров, при ограничениях, не требующих знания закона распределения помех. Предложенный алгоритм был реали-Matlab
и рекуррентным методом расширенных инструментальных переменных. Результаты моделирования подтвердили высокую эффективность предложенного алгоритма.
; ; -ременных; билинейные системы; метод наименьших квадратов.
D.V. Ivanov, O.V. Uskov
RECURSIVE ESTIMATION BILINEAR DYNAMICAL SYSTEMS WITH ERRORS-IN-VARIABLES
It is well known that the method of least squares gives biased estimates of parameters for dynamic systems when the observed input-output data are corrupted with noise. Recursive algorithm based on stochastic approximation is proposed for identification single-input-single-output (SISO) bilinear dynamic systems with errors-in-variables. The estimates are proved to be convergent to the true values with probability one, the resulting estimates of parameters under con-
straints that do not require knowledge of density. The proposed algorithm was realized in Matlab and compared with the recursive least-squares method and the recursive method extended instrumental variables. The results of a simulated example indicate that the proposed algorithm provides good estimates.
Recursive estimation; stochastic approximation; errors-in-variables model; bilinear systems; least square method.
.
. -
нейных систем находит применение во многих областях науки, таких как ядерная физика, электрические сети, химическая кинетика, гидродинамика и т.д. [1].
Проблема идентификации динамических систем с помехами во входных и , , , -становка задачи в регрессионном анализе, когда зашумленным является только . -ных методов идентификации систем с помехами во входных и выходных сигналах [2-4]. Подавляющее большинство работ посвящено линейным динамическим сис-.
В работе [5] предложен метод идентификации билинейных систем с помехой наблюдения в выходном сигнале, являющийся обобщением метода наименьших квадратов [6]. В данной работе дано обобщение результатов работы [5] на случай билинейной системы с помехами во входных и выходных сигналах. Разработан
,
параметров билинейных систем, и не требует априорной информации о распреде-.
.
стохастическими уравнениями с дискретным временем i = ... —1,0,1... r r r2 r3(m)
Zi — Ёb0m)Zi—m = Ё a0m)xi—m + Ё Ё c0mk)xi —mZi—k , (1)
m=1 m=0 m=0 k =1
yi = Zi +^1(І), wi = xi +^2(i).
где Zi, yi - ненаблюдаемая и наблюдаемая выходные переменные; xi , w- — ненаблюдаемая и наблюдаемая входные переменные; B,x(l), ^О')- помехи наблюдения в выходном и входном сигналах.
Пусть выполняются следующие предположения:
10. Множество В, которому априорно принадлежат истинные значения па, ,
.
20. Помехи {£(i')} и {^2(i)} статистически не зависят между собой c
E{£(i)/F(1)} = 0, E{^(i)/ F-(2)} = 0,
E{£ (i) / F-(1)} < Wi°) < ~, E{(£22 (i) / F-(2)} = Wi(2) < ~,
где F-(1), F/2) — <г-адгебры, индуцированные семействами случайных величин {1(t),t eT } И {2(t),t e T}, } = {t, t < i, t e Zc } Zc - множество целых чисел, щw, W(2)- случайные величины eЩ-(1) , E(w/2"*)<^£ , где E - оператор
.
30. {£(')} ,{^2 (')} статистически не зависит от {х. }.
40. Последовательности |х. }-стационарные в узком смысле с дробно-
рых л> 0: |хг-| <ЛХ п.н.
Требуется определять оценки неизвестных коэффициентов динамической системы, описываемой уравнением (1) по наблюдаемым последовательностям у., , при известных порядках г, г1, г2 и г3 определить оценки истинных значе-
.
Рекуррентный алгоритм оценивания параметров. Система может быть записана как линейная регрессия
(2) , -нее значение, а из предположения (3) - что ее локальная дисперсия с вероятностью
1 будет равна:
Тогда оценки неизвестного вектора д можно получить с помощью стохастически градиентного алгоритма минимизации суммы взвешенных квадратов
2 _____2
рациональной плотностью случайные сигналы с Е|(хг-) } = &Х > 0. Для некото-
У1 =ф1д + £; ,
где ф. = (фту (') | Ф1 (') | Ф1у ('))У , фу (0 = (у._1,_у._г У,
(2)
г г
£г = £(') _ ЁЪ0т)£(а _т) _ Ё«0т)^2(._т)_
т=1
т=0
(3)
т=0 к=1
2 2 Л Л Л
обобщённых ошибок (Ь0, а0, С0,1) с весом <7£ (Ъ('), й(1), С(')) , т.е.:
д(1 +1) = д(1) _ д (у.+1 _ Ф+дО'))2/&е (КО,а(г), С(0) , (4)
где аа - последовательность, удовлетворяющая условиям:
5°. = ~, а( > а.+1, £аа <-, I >1
а=0 а=0
^ «>
б0. Ха#х(1)<~, Ха&О’)<те. п^.
а=1 а=1
Теорема. Пусть динамическая система описывается уравнениями (1) и выполняются предположения 1°-6°, тогда оценки, определяемые алгоритмом (4),
УЧ л
либо д(0---------------:-— д0 и.н., либо 0(1)-:-->^.
а —— ^ а —— ^
— 2 — 2
В формуле (4) используются дисперсии помех наблюдения <7^ ,<72 И дис-
-о- 2 2
персии входного и выходного ненаблюдаемых сигналов <7Х , <7г . Согласно теоре-
2 ^ 2 —2 —2
ме Манна-Вольда [б]: если случайные величины <71 , <72 , <7Х , <7 сходятся
__2 ____2 ___2 _____2
почти к постоянным <1 , <2, <х, , т0 любая непрерывная функция
т/-72 -72 -72 -72ч „
J ( <7х , <7 , 0^ , <72 ) сходится почти наверное к постоянной
3 ( , <г2, <12, <22):
-72 пм —2 -72 пм —2 -72 пм —2 -72 пм —2
х — <х , — , °1 — <1 , 2 — <2 ,
т/—2 —2 —2 —2ч п.н. тг—2 —2 —2 —2ч
3( ^х , ^ , ^1 , ^2)----------— 3(<х, < , <1 , <2). (5)
Следовательно, если заменить в (3) <1, <, <, < оценками <712, ,
, <7г2, оценки параметров 0 останутся сильно состоятельными. Состоятель-
- ;72 -72 _
ные и несмещенные оценки дисперсии , < может быть получена как
щ+1 = щ + (+1 _ щ )/у +1),
(а+1)=оЩ(а+1)+((щ _ щ)2 _ <7'Щ(а+1)У/ а,
Уг +1 = У +(У/+1 _ У /У + l),
<7У2(а +1) = <7У2(а +1) +((Уа _ уа)2 _ <7У2(а + о)/а.
Дисперсии ненаблюдаемых и наблюдаемых входных сигналов связаны от:
—2 —2 , —2 —2 —2 , —2
+ <2 , <у = + <1 . (б)
Таким образом, для получения сильно состоятельных оценок параметров 0
- —2 —2 достаточно знать не все значения дисперсии, а только два, например: <1 , <2 или
__2 ___2
< х , , ли®0 знание отношений сигнал-шум на входе и выходе. Дальнейшим
направлением исследований должна стать разработка метода, не требующего знания дисперсий помех и сигналов.
. (4)
в МаНаЪ и сравнен с рекуррентным алгоритмом наименьших квадратов (РМНК) и
( ).
Динамическая система описывается уравнениями:
со
).
(7)
га _ 0,7 г;_1 + 0,4г;_2 =
= 0,3ха + 0,7 ха_1 + 0,2 х{ _ 2 + 0,2 х{г{ _1 На вход подавался входной сигнал:
х« + 0,5 • х«_1 = £. + 0,8 • £_ + 0,6 • £_2,
где ^ - белый шум.
Отношение “помеха-сигнал”: <1 /< ~0,5, <2/<х = 0,5. Начальные значения 0. . 1 -метров, определяемые по формуле: §д. = 0 _дJ/||(д° /•100%.
.
параметров билинейных систем с помехами наблюдения во входном и выходном сигналах. В среде МаЙаЪ создано программное обеспечение, результаты моделирования подтверждают эффективность работы предложенного алгоритма. Полученные результаты могут послужить основой для создания высокоэффективных
( ).
Дальнейшие исследования могут быть направлены на построение алгоритмов идентификации при автокоррелированных помехах.
Рис. 1. График погрешности оценок: 1 - РМНК; 2- РМИП; 3 - Алгоритм (4)
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Mohler R.R. Bilinear Control Processes: With Applications to Engineering, Ecology, and Medicine. New York: Academic Press, 1973.
2. Soderstrom T., Soverini U., Mahata K. Perspectives on errors-in-variables estimation for dynamic systems // Signal Processing. - 2002. - Vol. 82, № 8. - P. 1139-1154.
3. Soderstrom T. Errors-in-variables methods in system identification. 14th IFAC Symposium on System Identification // Invited plenary paper. Newcastle, Australia, March, 29-31, 2006.
4. Linden J. G., Vinsonneau B., Burnham K.J. Review and comparison of some identification methods in the errors-in-variables framework. In Int. Conf. of System Engineering 2006. - P. 243-254. Coventry, UK, 2006.
5. . .
// : IV
конференции в 2-х т. Т. 1. - Челябинск: Изд-во Челяб. гос. ун-та, 2011. - С. 56-57.
6. Mann, H.B., Wald A. On stochastic limit and order relationship // The Annals of Mathematical Statistics. - 1943. - Vol. 14 (3). - P. 217-226.
Статью рекомендовал к опубликованию д.т.н., профессор А.Н. Митрофанов.
Иванов Дмитрий Владимирович - ФГБОУ ВПО «Самарский государственный университет путей сообщения»; e-mail: dvi85@list.ru; 443066 г. Самара 1-й Безымянный пер., 18; тел.: 884629995460; кафедра мехатроники в автоматизированных производствах; старший преподаватель; к.ф.-м.н.
Уеков Олег Владимирович - e-mail: quentyn@bk.ru; кафедра мехатроники в автоматизированных производствах; аспирант.
Ivanov Dmitriy Vladimirovich - Samara State University of Transport; e-mail:dvi85@list.ru; 18, 1-st Bezymyanniy lane, Samara, 443066, Russia; phone: +784629995460; the department of mechatronics in automated production; senior lecturer; cand of phis.-math. sc.
Uskov Oleg Vladimirovich - e-mail: quentyn@bk.ru; the department of mechatronics in automated production; postgraduate student.
УДК 519.713.4
Н.И. Крайнюков
ПАРАЛЛЕЛЬНЫЕ АЛГОРИТМЫ МИНИМИЗАЦИИ КОНЕЧНЫХ АВТОМАТОВ ДЛЯ СУБРЕГУЛЯРНЫХ ЯЗЫКОВ
Субрегулярные языки - это подклассы регулярных языков, обладающие определенными свойствами, например к наиболее важным, с точки зрения приложений относятся ко, , -префиксные суффиксные, факториальные регулярные языки. Конечные автоматы, распознающие эти языки, представляют более узкий класс автоматов, состояния, таких конечных автоматов обладают определенными свойствами. Рассматриваются параллельные алгоритмы минимизации конечных автоматов для таких языков. Распараллеливание достигается применением функции перехода для каждого класса эквивалентных состояний. Для исследования возможностей параллельных алгоритмов используется технология MPI -параллельные процессы с передачей сообщений - наиболее распространенная технология
.
Субрегулярные языки; параллельные алгоритмы минимизации; конечные автоматы; технологии параллельного программирования.
N.I. Krayinyukov
A PARALLEL ALGORITHMS FOR MINIMIZATION OF FINITE AUTOMATA FOR SUBREGULARY LANGUAGES
Subregular languages are a subclass of regular languages with certain properties, the most important from the point of view of applications are finite languages, the languages are closed under taking prefixes, suffixes and subwords - prefix suffix, factorial regular languages. Finite automata that recognize these languages, representing a subclass of automata, the states of finite automata have certain properties. We consider parallel algorithms for the minimization of finite automata for these languages. Parallelization is achieved by using the transition functions for each class of equivalent states. To investigate the possibilities of parallel algorithms uses a technology called MP! - Parallel processes with message passing - the most common parallel programming technology.
Subregular languages; parallel algorithms for minimization; finite state machines; parallel programming technologies.
. -
зации конечных автоматов, представляющих субрегулярные языки. К наиболее важным, с точки зрения приложений, относятся конечные языки, языки, замкнутые относительно взятия префиксов, суффиксов и подслов - префиксные суф-, . -
І92
