УДК 004.942
DOI: 10.17277/vestnik.2016.04.pp.559-564
IDENTIFICATION'S SOFTWARE MODULE OF THE DYNAMICS MATHEMATICAL MODELS OF MULTIDIMENSIONAL OBJECTS
T. Yu. Dorokhova1, E. Yu. Vyguzova2
Departments: "Designing of Radio Electronic and Microprocessor Systems " (1), "Foreign Languages" (2), TSTU, Tambov, Russia; ele-vigozova@mail.ru
Keywords: controlled object; identification of mathematical models; mathematical modeling; software module.
Abstract: The article considers the identification problem of dynamics' mathematical models of multidimensional objects, the algorithm of parametrical identification is shown in the example model one zone of a multidimensional object as a double aperiodic link.
The considered algorithm of parametrical identification of dynamics' mathematical model of one multi-chamber dryer drum-band type section in the developed software module is implemented which allows to check and to assess the adequacy of the obtained model.
The mathematical modelling problem of the controlled object is one of the key components in the synthesis of automatic systems [1]. The most practical objects have multiple inputs and multiple outputs, it means they belong to the class of multidimensional (MIMO - Multiple Input Multiple Output), so the identification problem of the dynamics mathematical models of MIMO-objects is relevant.
The identification problem is formulated as follows: according to the results of the measured values of input and output variables of the object should be obtained are optimal in some sense or the largest permissible error model suitable for the solution of management tasks [3, 6]. The General identification problem of dynamic model of the controlled object can be divided into a number of specific tasks, such as:
- identification of the model structure;
- estimation of model parameters;
- checking the adequacy of the obtained model.
For the identification of the model structure we will consider the simplified scheme of multidimensional object (Fig. 1).
i
U2
n—1
I Output^
Zone 1
Zone 2
Zone n - 1
Zone n
u
1
n
l
n
Fig. 1. The simplified scheme of multidimensional object
On this scheme: z,, i = 1, n - phase coordinate i-zone's; u, - the operating influence li, ri - the revolting influences from the next zones at the left and from the right.
Based on constructive and technological reasons each zone is considered as the local multistage object with the concentrated parameters described by system of the differential equations with explosive right part [2]. The i-zone's model is written as:
M, : Zi =
A1z1 (t) + B]u] (t) + C1/1 (t) + DV (t), z1 e [z0, zK);
A2z2 (t) + B2u2 (t) + C2/2 (t) + D2r2 (t), z2 e [z02, zK2);
(1)
,mi mi , Dmi mi , 1mi /.\ , nmi mi / mi r mi
a, 1 zi 1 (t) + Bi 1 ui 1 (t) + Ci 1 li 1 (t) + Dt 1 ri 1 (t), z, 1 e [zo 1, z,
where A/, Bj, Cj, Dj, j = 1, mi, i = 1, n - are the matrixes of parameters for j-stage
of i-zone; z/- a vector of phase coordinates; zj0, zfK - thej-stage border.
The dynamics of all objects can be presented as a set of zone's models
M = {M1, M2, ..., Mn}, (2)
where M■, i = 1, n - a model of i-zone.
Each stage in (1) can be presented to one of dynamics basic models [4]:
- aperiodic link (A);
- double integrator (DI);
- real double integrator (RDI);
- double aperiodic link (DA), etc.
The experiment on real object is made to definite the types of zone's models, number of stages and an assessment of parameters. As a result of the solution of structural identification problem optimum structures of dynamics models of MIMO object zones are allocated, one of the main choice criteria is the minimum of the maximum value of an absolute error of model.
The model parameters identification is the following stage after carrying out structural identification. Let's consider algorithm of parametrical identification on the example of model of one zone of multidimensional object in the form of a double aperiodic link:
i z1 = z2 (t);
¿2 = (t) + «2z2 (t) + bu (t) + cl (t) + dr (t).
The parametrical identification algorithm of dynamics model of one zone of multidimensional object consists of the following stages:
1. A choice of basic data for identification (for the set type of model six values of the phase coordinates operating and the revolting influences in different time points get out):
z1 (t), z2 (t), u1 (t), l1 (t), r1 (t), i = 16. (3)
2. The solution of the equations system
[ Ko z[ (t) + K1z2 (t) + K 2^ (t) + K3I1 (t) + K4 r1 (t) = z1+1 (t);
j . . . , (4)
[K5 z{ (t) + K6z2 (t) + Knu (t) + Kg/1 (t) + K9r1 (t) = z2+1 (t).
3. The calculation of dynamics model parameters:
K5
a =—; a? = 1 Kj 2
K6 - K0 . b = K7.
Ki
Ki
c = -
K8
Ki
d = K9 Kj
(5)
4. Checking the adequacy of the obtained model. To check the adequacy of the obtained model imitating modeling is carried out. The values of phase coordinates are calculated by the formulas:
I zm (i +1) = F1Z1m(i) + G1Zm(i) + H1u(i) + V1l(i) + W1r(i); [z2m (i +1) = F2 Zim (i) + G2 Z2m (i) + H2u(i) + V2I (i) + W? r(i),
(6)
where
a2 I a?
X = T +V T + a';
X? = TT "JT? + a
Ft =
a1 (e
aA e 1 - e 2
)
Gj =
Hi =
X1 -X? b X2e^°' -X1e
X,St X? St e 1 - e 2
X1 -X2
~X2St -X? + X1
X1 -X? X1X2
c X2eX1St -X1eX ?St - X? + X1
X1 -X2 X1X2
d X2eX1St - X1eX?t - - X? + X1
Fj =
G? =
H? =
= X1eX?St -X2eXjSt,
X1 -X2
XjX
V? =-
W? =
1A2
X1 -X2
X1eXiSt - X2eX?St X1 -X2 b (eXiSt - eX?St) X1 -X2 ;(eXiSt - eX?St) X1 -X2 d (eXiSt - eX?St) X1 -X2
According to the result of calculation are plotted phase coordinates and calculated the absolute error, based on the value which a conclusion about the adequacy of the
obtained dynamics model.
The program module in MathCAD system for realization above considered algorithm of identification is developed. As an example we will consider the mathematical dynamics model identification of one multi-chamber dryer drum-band type section [5, 7].
The basic data for dynamics model identification of multi-chamber dryer drum-band type section are given in the Table j.
The results of mathematical dynamics model identification of multi-chamber dryer drum-band type section and the check of its adequacy are given in the figure 2.
Table 1
Basic data for identification
№ T zx(t) z? (t) l(t) r(t) u (t)
1 0 20,0 0 0 0
2 2 27,4 4,626866 - 2,1 - 0,1
3 4 36,6 4,328358 - 3,2 0
100
4 6 44,0 3,283582 - 2,9 0
5 8 49,9 2,686567 - 3,0 - 0,1
6 10 54,1 1,940299 - 2,6 0,3
<jl ON Ы
Значения коэффициентов имитационной модели:
ся ся 2
Расчет значений фазовых координат
U := 100 Д, , 0 г0 := 0 ^ := -2.1 ь, := 0.1 Ц := -3.2 ^ := 0 ^ := -29 г, := 0 ^ , -3 rg := -0.1 ^ := -2.Í ^ := 0.3
Zlirig := 20 Z2m^ := 0
td te
H "I н
н
о
н
о C¡
и
Н ся н с;
Построение графиков фазовых координат
s := 0,2.. N
Zlittj
Расчет абсолютной ошибки модели
i:= 0,2.. 10 Максимально допустимое значение ошибки
Eclop := 3.2
i
Emax := тая(Е) Emax = 1.083
SR := return "Абсолютная ошибка превышает допустимое значение" if Emax > Edop
|retum "Абсолютная ошибка не превышает допустимого значения" if Emaj: < Edop
Результат расчета: SR = "Абсолютная ошибка не преввплает допустимого значения"
Значения параметров модели динамики:
Fig. 2. The results of mathematical dynamics model identification of one multi-chamber dryer drum-band type section
By the results of the software module operation it is possible to draw a conclusion that the obtained dynamics model is adequate as the absolute error of model doesn't exceed accepted value.
References
1. Starikov A.V [Parametric identification of linear static controlled objects], Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta [Bulletin of the Samara State Technical University], 2004, no. 27, pp. 74-77. (In Russ.)
2. Bessonov A.N., Zagashvili Yu.V., Markelov A.S. Metody i sredstva identifikatsii dinamicheskikh ob"ektov [Methods and means of dynamic objects identification], Leningrad: Energoatomizdat, 1989, 300 p. (In Russ.)
3. Muromtsev D.Yu. Metody i algoritmy sinteza energosberegayushchego upravleniya tekhnologicheskimi ob"ektami [Methods and synthesis algorithms of energy saving management of technological objects], Tambov, Moscow, St. Petersburg, Baku; Vienna: Nobelistika, 2005, 202 p. (In Russ.)
4. Muromtsev Yu.L., Orlova L.P. Mikroprotsessornye sistemy energosberegayushchego upravleniya [Microprocessor systems of energy saving management], Tambov: Izdatel'stvo Tambovskogo gosudarstvennogo universiteta, 2001, 80 p. (In Russ.)
5. Gribkov, A.N. [Management information system of the dynamic modes in multisection drying installations], Informatika i sistemy upravleniya [Informatics and control systems], 2009, no. 2(20), pp. 123-129. (In Russ.)
6. Gribkov A.N., Kurkin A.I. [Identification of mathematical models of a drum dryer' s dynamics on the set of States of the operation], Voprosy sovremennoi nauki i praktiki. Universitet imeni V.I. Vernadskogo [Problems of Contemporary Science and Practice. Vernadsky University], 2014, special issue (52), pp. 32-36. (In Russ., Abstract in Eng.)
7. Muromtsev D.Yu., Gribkov A.N., Churikov A.A., Pasechnikov I.I. [Methodological aspects of building software-analytical complex design of energy-efficient control systems], Transactions of the Tambov State Technical University, 2015, vol. 21, no. 4, pp. 542-547. doi: 10.17277/vestnik.2015.04.pp.542-547 (In Russ., Abstract in Eng.)
Программный модуль идентификации математических моделей динамики многомерных объектов
Т. Ю. Дорохова1, Е. Ю. Выгузова2
Кафедры: «Конструирование радиоэлектронных и микропроцессорных систем» (1), «Иностранные языки» (2), ФГБОУВО «ТГТУ», г. Тамбов, Россия; ele-vigozova@mail.ru
Ключевые слова: идентификация математических моделей; математическое моделирование; объект управления; программный модуль.
Аннотация: Рассмотрена проблема идентификации динамики математических моделей многомерных объектов, где алгоритм параметрической идентификации показан на примере модели одной зоны многомерного объекта в виде двойного апериодического звена.
Реализуется рассматриваемый алгоритм параметрической идентификации математической модели динамики одной секции многокамерной сушильной установки вальце-ленточного типа в разработанном программном модуле, который позволяет производить проверку и давать оценку адекватности полученной модели.
Список литературы
1. Стариков, А. В. Параметрическая идентификация линейных статических объектов управления / А. В. Стариков // Вестн. Самар. гос. техн. университета. -2004. - № 27. - С. 74 - 77.
2. Бессонов, А. Н. Методы и средства идентификации динамических объектов / А. Н. Бессонов, Ю. В. Загашвили, А. С. Маркелов. - Л. : Энергоатом-издат, 1989. - 300 с.
3. Муромцев, Д. Ю. Методы и алгоритмы синтеза энергосберегающего управления технологическими объектами : монография / Д. Ю. Муромцев. -Тамбов ; М. ; СПб. ; Баку ; Вена : Нобелистика, 2005. - 202 с.
4. Муромцев, Ю. Л. Микропроцессорные системы энергосберегающего управления : учеб. пособие / Ю. Л. Муромцев, Л. П. Орлова. - Тамбов : Изд-во Тамб. гос. техн. ун-та, 2001. - 80 с.
5. Грибков, А. Н. Информационно-управляющая система динамическими режимами в многосекционных сушильных установках / А. Н. Грибков // Информатика и системы управления. - 2009. - № 2(20). - 2009. - С. 123 - 129.
6. Грибков, А. Н. Идентификация математических моделей динамики барабанной сушильной установки на множестве состояний функционирования / А. Н. Грибков, И. А. Куркин // Вопр. соврем. науки и практики. Университет им. В. И. Вернадского. - 2014. - Спец. вып. (52). - С. 32 - 36.
7. Методологические аспекты построения программно-аналитического комплекса проектирования систем энергосберегающего управления / Д. Ю. Муромцев [и др.] // Вестн. Тамб. гос. техн. ун-та. - 2015. - Т. 21, № 4. - С. 542 - 547. ао1: 10.17277/уе$1тк.2015.04.рр.542-547
Programmmodul der Identifizierung der mathematischen Modelle der Dynamik der multidimensionalen Objekte
Zusammenfassung: Es ist das Problem der Identifizierung der Dynamik der mathematischen Modelle der multidimensionalen Objekte betrachtet. Es ist der Algorithmus der parametrischen Identifizierung am Beispiel des Modells einer Zone des multidimensionalen Objektes (der mehrzelligen Trockenanlage des Walzebandtyps) in Form des doppelten aperiodischen Gliedes dargelegt.
Module logiciel de l'identification des modèles mathématiques de la dynamique des objets multidimensionnels
Résumé: Est considéré le problème de l'identification de la dynamique des modèles mathématiques des objets multidimensionnels. Est présenté l'algorithme de l'identification paramétrique à l'exemple du modèle d'une zone de l'objet multidimensionnel (installation multicaméra du séchage du type rectifieuse à bande) sous la forme d'un double maillon apériodique.
Авторы: Дорохова Татьяна Юрьевна - кандидат педагогических наук, доцент кафедры «Конструирование радиоэлектронных и микропроцессорных систем»; Выгузова Елена Юрьевна - кандидат филологических наук, доцент кафедры «Иностранные языки», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.
Рецензент: Селиванова Зоя Михайловна - доктор технических наук, профессор кафедры «Конструирование радиоэлектронных и микропроцессорных систем», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.