УДК 539.3
Вестник СибГАУ Том 17, № 3. С. 595-601
ОБ УПРАВЛЕНИИ ^-ПРОЦЕССАМИ
А. В. Медведев
Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнева Российская Федерация, 660037, г. Красноярск, просп. им. газ. «Красноярский рабочий», 31
Е-mail: saor_medvedev@sibsau.ru
Исследуется проблема управления дискретно-непрерывными процессами, которые имеют «трубчатую» структуру в пространстве входных-выходных переменных. Подобные процессы протекают в пространствах дробной размерности. Моделирование и управление процессами этого класса существенно отличаются от общепринятых параметрических моделей, представляющих собой поверхности в том же пространстве. При построении обучающихся параметрических моделей «трубчатых» процессов необходимо использование соответствующих непараметрических индикаторов. Многомерный Н-процесс включает в себя как управляемые, так и неуправляемые векторные переменные, при этом неуправляемые входные переменные контролируются в процессе функционирования многомерной системы. Естественно задающие воздействия также представляют собой вектор. При управлении Н-процессами необходимо компоненты вектора задающих воздействий согласовывать между собой, а не задавать их произвольно, как это принято в теории автоматического управления. В результате возникает необходимость синтеза Am -регулятора, а не A -регулятора. На первом этапе необходимо выделить некоторую область пересечения всех Н-процессов по каждой компоненте выходной переменной, для этого используется соответствующий индикатор. В результате получаем некоторую подобласть согласования в пространстве входных-выходных переменных. Только из этой подобласти можно определить задающее воздействие для каждой компоненты выхода многомерной системы. Следующий этап - это вычисление управляющих воздействий при фиксированном значении векторов входных неконтролируемых переменных. Анализируется случай, когда задающие воздействия для компонент выходных переменных не удается в конкретной ситуации согласовать или согласовать частично. Приводятся алгоритмы определения компонент вектора задающих воздействий. Даны также непараметрические алгоритмы управления многомерными Н-процессами. В результате такое управляющие устройство можно отнести к матричным системам управления. Таким образом, последний реализует как алгоритм определения задающих воздействий, так и алгоритм управления объектом. Он-то и является Am -регулятором (управляющей системой) для многомерного Н-процесса. Естественно, могут возникать случаи, когда Н-процессы по всем компонентам выхода оказываются непересекающимися. В подобной ситуации необходим поиск приемлемых компромиссов с целью поэтапного приведения состояния управляемого процесса к желаемому.
Ключевые слова: Н-процесс, H-модели, безынерционный объект, пространство дробной размерности, индикатор, подпространство согласования, Am -регулятор.
Sibirskii Gosudarstvennyi Aerokosmicheskii Universitet imeni Akademika M. F. Reshetneva. Vestnik Vol. 17, No. 3, P. 595-601
ON THE CONTROL BY ^-PROCESSES
A. V. Medvedev
Reshetnev Siberian State Aerospace University 31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation Е-mail: saor_medvedev@sibsau.ru
The paper considers the control problem of discrete-continuous processes with a "tubular " structure in the space of input-output variables. The similar processes occur in spaces of the fractional dimension. The modeling and control by this class of processes differ very much from the conventional parametric models that are surfaces in the same space. It is necessary to apply the appropriate non-parametric indicators for the development of the learning parametric models of the "tubular" processes. The multidimensional H-process includes both controlled and uncontrolled vector variables; the uncontrollable input variables are supervised while the operation of the multidimensional system. The defining impact is also a vector. It is necessary to coordinate the vector components of the defined controls but not to define them randomly as it is applied in practice of the automated control theory. As a result the necessity for synthesis a Am -controller but not a A -controller arises. Firstly, it is necessary to allocate
a certain intersection area of all H-processes for each component of the output variable the corresponding indicator which is applied for. As a result, we get some subarea of coordination in the space of input-output variables. Only from this subarea one can determine the defining impact for each component of the output of the multidimensional system. The following stage is the calculation of the defining impacts with a fixed vector value of the input uncontrolled variables. We analyze the case when the defining impacts for the component of the output variables are not possible in a particular situation to agree absolutely or partially. Algorithms of component definitions of a vector of the setting influences are given. Also nonparametric control algorithms for the multidimensional H-processes are given. As a result such a control unit can be referred to a matrix control systems. Thus, the latter is implemented as the algorithm for determining the defining impacts and a control algorithm by the object. It is a Am -controller (a control system) for the multivariate H-process. Of course there might be cases when H-processes according to all components of the output are disjoint. In such a situation it is necessary to search for the acceptable compromise to bring the state of the controlled process to the desired one.
Keywords: H-process, H-models, instantaneous object, space of the fractal dimension, indicator, subspace matching, Am -controller.
Введение. Идентификация многих стохастических объектов часто сводится к идентификации статических систем с запаздыванием. Обусловлено это тем, что некоторые выходные переменные объекта контролируются через значительно большие интервалы времени, чем входные, и существенно превышают постоянную времени объекта. Например, ряд переменных измеряется электрическим способом (в этом случае дискретность контроля At может быть достаточно мала), а другие переменные контролируются в результате химического анализа или физико-механических испытаний (в этом случае дискретность контроля АТ велика, т. е. АТ >> At). В результате, переходный процесс, протекающий в объекте, устанавливается за значительное меньшее время, чем измерение выходной переменной.
Идентификация статической системы. Пусть и = (м1, ..., ик) еП (и) с Як , ^ = (ц1, ..., цп) ей с Я" , а х = (х1, ..., хт) ей (х) с Ят. Вообще говоря, каждая компонента вектора и е [а{; Ь ], i = 1, к, ^ е С;Л] ], } = !,п , а хч е ;е?], д = 1,т . При исследовании реальных процессов значения коэффициентов , Ь, с], Л], ^, ед |, i = 1, к, ] = 1, п, д = 1, т ,
всегда известны. В технологических процессах значения этих коэффициентов регламентируются технологическим регламентом (картой). В дальнейшем, без нарушения общности, эти интервалы примем единичными [1], тогда й(и) - единичный гиперкуб,
йк (и) = [0;1], йп (ц) = [0; 1], т. е. и е [0;1], це [0;1],
йк+п+т (и' ^ х) = (и' ^ х) е йк+п+т (и' ^ х).
Задачу идентификации часто сводят к параметрической, состоящей из двух основных этапов. Первый этап - выбор (определение) параметрической модели в виде х(и, ц) = / (и, а), где а - вектор параметров,
и на втором этапе - последующая оценка параметров а на основании поступающих элементов выборки (и, ц, х!),(и2,ц2, х2), ...,(и, ,ц,, х,), т. е. получение оценки а,. Адаптивная модель в этом случае будет выглядеть следующим образом:
Х(и, ц) = / (и, а). (1)
Оценка параметров а при наличии обучающей выборки может быть осуществлена методом наименьших квадратов или методом стохастических аппроксимаций. Такова общая схема решения задач параметрической идентификации. Отметим только, что наиболее слабым местом здесь является выбор параметрической структуры модели. Если на первом этапе допущена достаточно грубая ошибка, то, в итоге, полученная модель вряд ли будет удовлетворительной. Эта проблема достаточно подробно обсуждалась в [2; 3]. Там же предложен новый класс ^-моделей, учитывающий в комплексе знание фундаментальных законов, другую априорную информацию об объекте, в том числе разнотипную. Обратим внимание на то, что модели класса (1) представляют собой гиперповерхности в пространстве входных-выходных переменных объекта, т. е. (и, ц, х) е й(и, ц, х) с Як+"+т . Процессы трубчатой структуры, или Я-процессы [4], представлены на рис. 1. Поскольку значение входных неуправляемых переменных влияет на расположение реального Я-процесса при отличающихся то это обстоятельство отражено на рис. 1. На рис. 1, а представлен Я-процесс при ц = ц', на рис. 1, б - при ц = ц".
Если исследуемый процесс имеет «трубчатую» структуру, то модель (1) необходимо подкорректировать следующим образом [2]:
хх (и, ц) = I, (и, ц) / (и, ц, а ,) (2)
либо:
N
х, (и, ц) = I, (и, ц)^а у ф; (и, ц), (3)
]=1
где ф] (и,^) - система линейно независимых функций; индикатор I, (и, имеет вид
1 ( ) {1, если и ейн (и), цей1^ (ц);
I, (и, ц) = <! (4)
|0, если и ейН(и),цейН (ц).
Заметим лишь, что, вообще говоря, область й^ (и, ц) нам не известна, а известна лишь выборка
|х1, и{, ц, / = 1, . Это обстоятельство особенно хорошо видно при анализе рис 1.
П( и х )
X н Ш и: Л"
"у
0.( и. х)
Рис. 1. Трубчатый процесс
Если индикатор равен нулю, то оценка х(и, ц), х!! (и, ц) не может быть вычислена, т. е. при таких значениях компонент вектора и е О(и), ц е ЗДц) процесс протекать не может. Если индикатор 1!! (и, ц) при любом значении и е ЗДи), ц е О(ц) равен единице, то модель (2) совпадает с (1). В качестве оценки индикатора 1!! (и, ц) можно принять следующее приближение:
I* (и) = sgn(scs )-1 X Ф(с- (X* (и) - X )) X
г=1
к т
хП Ф(с-1 (и1 - и/ ))П Ф(с- (ц? - ц?)), (5)
1=1 ?=1
где
£ к т
X* (и) = X X П Ф(с-1(и1 - и/ ))П Ф(с-1(ц? -ц?))/
1=1 1=1 ?=1
* к т
/ ХП Ф(с-1(и1 - и/ ))П Ф(с-1(ц? -ц? )), (6)
г=1 1=1 ?=1
а параметр размытости с* и колоколообразная функция Ф(-) удовлетворяют условиям сходимости [5; 6]. Подобные непараметрические оценки достаточно подробно исследовались в [7; 8].
Таким образом, при известных значениях и = и' е О(и), ц = ц ' е О(ц) сначала строится оценка х!, (и = и' , ц = ц') по формуле (6), затем вычисляется индикатор 1ц (и), и только на следующем этапе используются модели (2) или (3), если индикатор оказался равным единице. Если же индикатор равен нулю, то это означает, что хотя и' е 0.(и), ц ' е О(ц),
но и' е 0.н (и), ц ' е 0.н (ц), т. е. компоненты вектора и = и' = (и{, ...,и'к), ц = ц' = (ц1, ...,цк) определены неверно, иными словами, реально протекающий «трубчатый» процесс не соответствует совокупности
заданных значений компоненты векторов и = и , ц = ц '. Причины этого могут состоять в том, что компоненты векторов и = и' = (и{, ..., и'к), ц = ц ' = (ц1, ..., цк) выбраны неверно либо измерены с значительной погрешностью типа «выброс». Конечно же, это справедливо только при условии, что мы располагаем
представительной выборкой |хг-, ui, ц, г = 1, . Следует
заметить, что использование традиционных моделей типа (1) позволит получить оценку х(и = и', ц = ц'),
которая, естественно, будет далека от реальности.
Естественно считать, что процесс идентификации объекта в параметрической постановке также следует осуществлять с учетом «трубчатой» структуры объекта. Примем модель «трубчатого» процесса в виде (из соображений простоты ^ отпущено)
X (и) = I (и)Ха 1 ф 1 (и) ;
1=1
(7)
где ф 1 (и), 1 = 1,N - система выбранных линейно независимых функций; I* (и) - индикаторная функция (5).
Сформируем критерий оптимальности:
Л (а) = М <!(х(и) -1(и)^а^ф^ (и))
(8)
1=1
Наша цель состоит в отыскании таких а* = (а*,
Д что
Я(а*) = шшЯ(а).
(9)
Решение задачи (9) дается системой рекуррентных соотношений:
( N \
1 1-1 . 1
1=1
Ф1 (и*К и), 1 = 1, ...,N.
б
а
X
В качестве оценки I (м5) примем приближения:
I, м) = )-1 ХП Ф
¿=1 ] =1
( з\ м„ - м/
(11)
Ясно, что сходимость а5 к а* следует при 5 . Учет неуправляемых, но контролируемых входных переменных ^ в приведенном выше примере, не вызывает затруднений.
Об одной особенности моделирования «трубчатых» процессов. Приведем следующий пример, имеющий отношение к идентификации безынерционной системы. Рассмотрим следующий простой частный случай. Пусть объект описывается уравнением:
х(м) = /(м1,м2,м3), (12)
где трехмерный вектор м = (м1, м2, м3) е Я3 является
входной переменной, а х е Я1 - выходная переменная. Традиционный путь построения модели процесса, описываемого (12), состоит в определении класса параметрических зависимостей х(и) = /(и1;и2,и3,а) и последующей оценки параметров а тем или иным способом по выборке наблюдений (мi, xi), г = 1,5, где 5 - объем выборки. Проанализируем этот пример с разных точек зрения. Пусть компоненты вектора входных переменных м = (м1, м2, м3) стохастически никак не связаны, т. е. независимы. В этом случае естественно использовать обычный традиционный прием, описанный выше. Теперь предположим, что объективно компоненты вектора входных переменных функционально связаны, например,
М2 = Ф1 (М\ ) , Мз =ф2(М2) = Ф2 (ф1 (М] )). (13)
Такая ситуация рассматривалась в [9]. Естественно, исследователь не знает о существовании зависимостей (13). В противном случае можно было бы сделать подстановку (13) в (12) и получить следующую зависимость х уже от одной переменной м1 вида
х(м) = / (мъ ф^мД Ф2 (ф1 (М\ ))). (14)
Таким образом, зависимость (12) в приведенных выше условиях может быть сведена к одномерной зависимости х от м1 . В случае, если зависимость м3 от м2 объективно отсутствует, то (12) легко приводится к виду
х(м) = / (мх, ф1(м1), м3), (15)
т. е. к двумерной зависимости х от м1, м3. Отсюда можно заключить, что при наличии функциональной зависимости между компонентами вектора м мы получаем зависимость х от м , в данном случае одно-, двух-, трехмерные. Подчеркнем еще раз, что о наличии функциональных зависимостей между компонентами вектора входных переменных исследователю не известно. Просто мы проанализировали случай «если бы...». А теперь проанализируем наиболее интересный случай, имеющий непосредственное отношение к Я-процессам [1]. Пусть м3 и м2, хотя и неиз-
вестным образом, но стохастически связаны. Подчеркнем - стохастически, а не функционально. Вернемся еще раз к анализу того, что произошло. Во-первых, если компоненты вектора м независимы, то исследуемый процесс описывается функцией трех переменных. Если две компоненты вектора входных переменных м связаны функциональной зависимостью, то процесс описывается функцией двух переменных. Наконец, если две переменные связаны стохастически, то процесс описывается функцией более чем двух переменных, но менее чем трех?! Можно считать, что мы приходим к зависимости от дробного числа переменных и, следовательно, к пространству дробной размерности. Например, Б. Мондельброт в [10] замечает: «Кровеносная система человека -пульсирующая, живая - имеет размерность 2.7». Дробная размерность пространств, по видимому, впервые была отмечена в работах Хаусдорфа и Бези-ковича.
Рассмотрим следующую ситуацию. Из простоты соображений, пусть интересующий нас процесс описывается (12).
В случае стохастической зависимости между переменными м2(м1), м3(м1) по имеющимся в наличии обучающим выборкам можно вычислить квадратичную ошибку прогноза м25 (м1), м35 (м1). Здесь м25 (м1), м35 (м1) есть непараметрические оценки (6) [3]:
821 = Х(м2 - м25 (м1))2
"2 '
831 =Х(м3 - м35 (м1))/ст23 .
(16)
г =1
«Силу» стохастической связи X между двумя произвольными переменными можно, например, вычислить по формуле:
Х = 1 -5 , (17)
где 5 может быть равно 521 либо 531.
Отсюда видно, что самая сильная стохастическая связь (функциональная) равна 1, отсутствие связи -при X = 0, а при стохастической зависимости между входными переменными 0 < X < 1. Последнее легко прослеживается из следующих рассуждений: если переменные м1, м2, м3 независимы, то (12) не более чем функция 3-х переменных. Если м1, м2 связаны функционально, то в результате подстановки в (12) мы приходим к функции 2-х переменных. Если функционально связаны все переменные м1, м2, м3, то в результате подстановки приходим к функции 1-й переменной. Наконец, если м1, м2, м3 связаны стохастически, то размерность пространства естественно уменьшается и может лежать в интервале от 2 до 3 либо от 2 до 4. Это зависит от тесноты связи, т. е., например, от значения X.
Если в более общем случае такого рода процессы интерпретировать как функции многих переменных, то изменчивость этой функции во времени может быть, например, показана на нижеследующей цепочке соотношений, действующих во времени:
г=1
X = / (/, их, и2, из, и4, и5, , u7, u8, и10) - Т1
X = /(/, иь и2, из, и4, и5, , u7, u8, u9, и10) - Т2
X =/(¿, и1, и2, из, и4, и5, , и7, и8, и9, ихо) - Т3
X =/(', щ, и2, из, и4, и 5, , u7, и8> u9, и10) - Т4
X =/с, и2, из, и4, и5, и !б , u7, u8, , и10) - Т5
X =м , и 2, и3, и4, и5, и и7, и8, , и10) - Т6
X =м , и 2, u3, u4, и5у u6, u7, , , и10) - Т7
X =м , и2, из, , и5, иб, и7, и8, , ию) - Т8
X =/(', щ, и2, из, и4, и5, иб, и7, и8, , ию) - Т9
Поясним наши обозначения. Наиболее темным цветом (их) обозначены переменные, которые оказывают самое сильное влияние на х (возможно, функциональная зависимость). Менее темное обозначение (И1) говорит о более слабом влиянии переменной на х (возможно, стохастическая зависимость), более слабое влияние на х оказывают и1 и щ. Т, где 1 = 1,9, -интервалы существования соответствующих зависимостей. Таким образом, в реально действующих процессах подобного рода роли значения переменных существенно изменчивы. Из приведенных выше зависимостей видно, что некоторые переменные могут утрачивать свое значение, а некоторые утрачивают, а потом восстанавливаются, а некоторые новые переменные появляются впервые, как например и6, и7.
Если сохранить математический «облик» интерпретации функции многих переменных как точку многомерного пространства, то мы приходим к наличию пространства дробной размерности ЬХ . Вычис-
г-Л
ление размерности Ь можно осуществить, например, так:
и-1
Шш = (п +1)-^Х
1,1+1
(19)
где п - размерность вектора и; 1+1 означает «силу»
стохастической связи между и, и и1+1.
В принципе, могут быть предложены и другие схемы вычисления размерности пространства. Например,
п-1
Шш Ь1Х = (п+1)-ХХ
1,1+1
(20)
1 Здесь речь идет о теории упругости.
меньшие по величине "физически важные" члены, дало много включаемых в теорию катастроф результатов».
При наличии нескольких выходных переменных процесса он может быть представлен следующим образом.
На рис. 2 показаны два Я-процесса по каналам « х1 - и », « х2 - и ». Как видно, эти Я-процессы не пересекаются в пространстве входных-выходных переменных. Из этого следует, что оба задающие воздей-
* *
ствия х1 и х2 не могут быть достигнуты одновременно при управлении этой системой по и. На рис. 3 показаны пересекающиеся Я-процессы по тем же самым каналам [12].
Рис. 2. Непересекающиеся Я-процессы
Из рис. 3 видно, что у двух Я-процессов есть совместная подобласть 0н (и, х) = 0.н (и, х) п Пн (и, х).
Именно в этой области Я следует искать задающие
**
воздействия х1 и х2. Для этого сформируем из исходной обучающей выборки |хг-, и1, 1 = 1,элементы, принадлежащие области 0н (и,х). Эту выборку обозначим через | X-, Uj, ] = 1, л0, л0 << . На ее основании строится оценка задающего воздействия:
где Х11+1 - зависимость всех компонент вектора и от
одной компоненты и1.
При достаточно внимательном анализе разложения функций в ряды уместно вспомнить фразу В. И. Арнольда из замечательной книги «Теория катастроф» [11]: «Вычисление в этих прикладных1 исследованиях обычно проводилось без общей теории за счет правильного отбрасывания одних членов ряда Тейлора и оставления других, наиболее важных. Из физиков, особенно систематически применявших теорию катастроф до ее возникновения, стоит особо выделить Л. Д. Ландау. В его руках искусство отбрасывать "несущественные" члены ряда Тейлора, сохраняя
х- =
л0 к X х- П Ф 1=1 г=1 г г и : - и: 3 1 1 п П Ф ) ?=1
СЩ V л
л0 к ( ХП Ф 1=1 г=1 V г г \ и : - Щ J 1 п П Ф ?=1 V
сЩ л /
^-и? 1
и? -и?
V,-Л^. (21)
Непараметрический алгоритм дуального управления многомерной безынерционной системой примем в виде [13; 14]:
( *<1>
и? =
X и^ Ф
]=1
П Ф
?=1
? -и? ^
X Ф
-=1
( *<1> * х* - х-
ПФ
?=1
Ч?-и? ^
(22)
1=1
1=1
о
о
Рис. 3. Пересекающиеся Я-процессы
Второе значение управляющего воздействия может быть оценено с учетом уже найденного первого:
u<2> =
( *<2>
Х«г2>ф
3=1
Пф
q=1
Л.<1> „,<1>Ч\
( *<2>
ХФ
3=1
Пф
q=1
fHq л
Ф
ф us - uJ
cu s
V
Г uf> - ufl
,(23)
G f *<k> *
xs - Xj
Х^Ф
j=1
uf> = n f.
ПФ
q=i
HS -H
Пф
fu<r>-u<r>Л
s J
(24)
Хф xf>- x* )ПФ j=i v ' q=i
fH -Hq Л k
cH
Пф
fu<r>- и<гЛ
J
Как видно из (24), все последующие значения управления вычисляются с учетом предыдущих. Таким образом, алгоритмы управления многомерными .-процессами представляют собой цепочку соотношений (21)-(24).
Рис. 4. Матричная система управления .-процессами
Матричная система управления многомерным безынерционным объектом показана на рис. 4, где представлены следующие обозначения: УУ - управляющие устройства, УУ" блок расчета управляющих воздействий, УУ" - блок расчета задающих воздействий. Такими образом, вышеприведенная схема иллюстрирует функционирование системы управления многомерным безынерционным .-процессом при фиксированном неуправляемом входном воздействии
Заключение. Итог статьи состоит в анализе особенностей, возникающих при моделировании процессов «трубчатой» структуры, которая имеет место всегда, если компоненты вектора входных переменных процесса стохастически зависимы. В этом случае традиционно [15] используемые модели статических систем с запаздыванием неприменимы или могут приводить к значительным ошибкам. Наиболее интересным является тот факт, что мы приходим к необходимости введения пространства дробной размерности. Безусловно, важным является факт исчезновения и появления роли некоторых входных переменных в различные периоды времени на значения выходных переменных процесса, что тесно связано не столько с пространством дробной размерности, сколько с пространством изменяющейся размерности.
При управлении многомерными .-процессами сначала необходимо определить соответствующие задающие воздействие для выходных переменных, а уже затем использовать непараметрические алгоритмы управления, которые названы матричным
Am -регулятором. Это является существенной особенностью построения управляющей системы в отличие от традиционных алгоритмов управления многомерными безынерционными системами. В настоящей статье все .-процессы, показанные на рис. 1-3, из соображений простоты визуализации рассматриваются в трехмерном пространстве и в «застывшем» виде. На самом деле, приведенные .-модели и алгоритмы управления функционируют, когда реальные процессы находятся в движении.
Библиографические ссылки
1. Medvedev A. V. Nonparametric approximation in adaptive systems theory // Works of Applied Methods of Statical Analysis. Simulation and Statistical Inference. Novosibirsk : SFU, 2011. С. 195-212.
2. Медведев А. В. Анализ данных в задаче идентификации // Компьютерный анализ данных и моделирование. Минск : БГУ, 1995. Т. 2. С. 201-206.
3. Медведев А. В. Теория непараметрических систем. Моделирование // Вестник СибГАУ. 2010. Вып. 4(30), С. 4-9.
4. Медведев А. В. .-модели для безынерционных систем с запаздыванием // Вестник СибГАУ. 2012. Вып. 5(45), С. 84-89.
5. Кошкин Г. М., Пивен И. Г. Непараметрическая идентификация стохастических объектов : науч. изд. Хабаровск : Российская академия наук, Дальневосточное отделение, 2009. 336 с.
G
G
G
6. Васильев В., Добровидов А., Кошкин Г. Непараметрическое оценивание функционалов от распределений стационарных последовательностей. М. : Наука, 2004. 512 с.
7. Катковник В. Я. Непараметрическая идентификация и сглаживание данных: метод локальной аппроксимации. М. : Главная редакция физико-математической литературы, 1985. 336 с.
8. Хардле В. Прикладная непараметрическая регрессия. М. : Мир, 1993. 349 с.
9. Айвазян С. А., Енюков И. С., Мешалкин Л. Д. Прикладная статистика, исследование зависимостей. М. : Финансы и статистика, 1985. 488 с.
10. Мондельброт Б. Фрактальная геометрия природы. М. ; Ижевск : Ижевский институт компьютерных исследований : НИЦ «Регулярная и хаотическая динамика», 2010. 656 с.
11. Арнольд В. И. Теория катастроф. М. : Наука, 1990. 128 с.
12. Медведев А. В., Михов Е. Д. О компьютерном исследовании Н-моделей // Вестник СибГАУ. 2014. № 3(55). C. 107-113.
13. Медведев А. В. Непараметрические системы адаптации. Новосибирск : Наука, 1983. 173 с.
14. Медведев А. В. Основы теории адаптивных систем / СибГАУ. Красноярск, 2015. 525 с.
15. Методы классической и современной теории автоматического управления : учебник. В 5 т. Т. 3. Синтез регуляторов систем автоматического управления. 2-е изд., перераб. и доп. / под ред. К. А. Пупкова и Н. Д. Егупова. М. : Издательство МГТУ им. Н. Э. Баумана, 2004. 616 с.
References
1. Medvedev A. V. Nonparametric approximation in adaptive systems theory. Works of Applied Methods of Statical Analysis. Simulation and Statistical Inference, Novosibirsk, SFU Publ., 2011, P. 195-212.
2. Medvedev A. V. [Data analysis in the identification problem]. Kompyuternyy analiz dannykh modelirovaniya. Minsk, BGU Publ., 1995, P. 201-206 (In Russ.).
3. Medvedev A. V. [The theory of non-parametric systems. Modeling]. Vestnik SibGAU. 2010, No. 4(30), P. 4-9 (In Russ.).
4. Medvedev A. V. [H-models for non-inertia systems with a delay]. VestnikSibGAU. 2012, No. 5(45), P. 84-89 (In Russ.).
5. Koshkin G. M., Piven I. G. Neprametricheskaya identifikatsiya stokhasticheskikh ob"ektov. [Non-parametric
identification of stochastic objects]. Habarovsk, Rossiy-skaya akademiya nauk, Dal'nevostochnoe otdelenie Publ., 2009, 336 a
6. Vasil'ev V., Dobrovidov A., Koshkin G. Neparametricheskoe otsenivanie funktsionalov ot raspre-deleniy statsionarnykh posledovatel 'nostey. [Nonparametric Estimation of Functional of distributions of stationary sequences]. Moscow, Nauka Publ., 2004, 512 p.
7. Katkovnik V. Ya. Neparametricheskaya identifikatsiya i sglazhivanie dannykh: metod lokal'noj approk-simatsii. [Non-parametric identification and data smoothing: local approximation method]. Moscow, Glavnaya redakciya fiziko-matematicheskoy literatury Publ., 1985, 336 p.
8. Hardle V. Prikladnaya neparametricheskaya regressiya [Applied nonparametric regression]. Moscow, Mir Publ., 1993, 349 p.
9. Ajvazyan S. A., Enyukov I. S., Meshalkin L. D. Prikladnaya statistika, issledovanie zavisimostey. [Applied Statistics, the study of dependencies]. Moscow, Finansy i statistika Publ., 1985, 488 p.
10. Mondelbrot B. Fraktalnaya geometriya prirody. [Fractal Geometry of Nature]. Moscow - Izhevsk, Izhevskiy institut kompyuternykh issledovaniy, NITS "Regulyarnaya i haoticheskaya dinamika" Publ., 2010, 656 p.
11. Arnold V. I. Teoriya katastrof [Catastrophe Theory]. Moscow, Nauka Publ., 1990, 128 p.
12. Medvedev A. V., Mihov E. D. [On a computer study of ^-models]. Vestnik SibGAU. 2014, No. 3 (55), P. 107-113 (In Russ.).
13. Medvedev A. V. Neparametricheskie sistemy adaptatsii. [Non-parametric adaptation system]. Novosibirsk, Nauka Publ., 1983, 173 p.
14. Medvedev A. V. Osnovy teorii adaptivnykh system. [Basic theory of adaptive systems]. Krasnoyarsk, Izd. SibGAU Publ., 2015, 525 p.
15. Metody klassicheskoy i sovremennoy teorii av-tomaticheskogo upravleniya. T. 3. Sintez regulyatorov sistem avtomaticheskogo upravleniya [Methods of classical and modern control theory. Vol. 3. Synthesis of regulators of automatic control systems]. Ed. K. A. Pupkova, N. D. Egupova. Moscow, Izdatel'stvo MGTU im. N. E. Baumana Publ., 2004, 616 p.
© Медведев А. В., 2016