RU2528133C1 - Device for identifying lagrange dynamic systems based on iterative regularisation - Google Patents

Device for identifying lagrange dynamic systems based on iterative regularisation Download PDF

Info

Publication number
RU2528133C1
RU2528133C1 RU2013106875/08A RU2013106875A RU2528133C1 RU 2528133 C1 RU2528133 C1 RU 2528133C1 RU 2013106875/08 A RU2013106875/08 A RU 2013106875/08A RU 2013106875 A RU2013106875 A RU 2013106875A RU 2528133 C1 RU2528133 C1 RU 2528133C1
Authority
RU
Russia
Prior art keywords
formation
block
unit
information
information output
Prior art date
Application number
RU2013106875/08A
Other languages
Russian (ru)
Other versions
RU2013106875A (en
Inventor
Дмитрий Сергеевич Андрашитов
Андрей Алексеевич Костоглотов
Original Assignee
Российская Федерация, от имени которой выступает Министерство обороны РФ
Федеральное государственное казенное военное образовательное учреждение высшего профессионального образования Военная академия Ракетных войск стратегического назначения имени Петра Великого МО РФ
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Российская Федерация, от имени которой выступает Министерство обороны РФ, Федеральное государственное казенное военное образовательное учреждение высшего профессионального образования Военная академия Ракетных войск стратегического назначения имени Петра Великого МО РФ filed Critical Российская Федерация, от имени которой выступает Министерство обороны РФ
Priority to RU2013106875/08A priority Critical patent/RU2528133C1/en
Publication of RU2013106875A publication Critical patent/RU2013106875A/en
Application granted granted Critical
Publication of RU2528133C1 publication Critical patent/RU2528133C1/en

Links

Images

Landscapes

  • Complex Calculations (AREA)

Abstract

FIELD: information technology.
SUBSTANCE: device includes constant storage units, function generation units, derivative generation units, product generation units, difference generation units, summation units, transposition units and integration units. Said technical result is achieved via calculation and successive refinement of regularisation parameters.
EFFECT: high accuracy of identifying parameters of dynamic systems.
2 dwg

Description

Изобретение относится к области цифровой вычислительной техники и может быть использовано в автоматических и автоматизированных системах различного назначения для идентификации параметров.The invention relates to the field of digital computing and can be used in automatic and automated systems for various purposes to identify parameters.

Известно устройство идентификации параметров динамических систем, в основе которого лежит фильтр Калмана. Его эффективное функционирование возможно в случае априорной определенности законов распределения внешних воздействий, что является существенным недостатком.A device for identifying parameters of dynamic systems is known, which is based on a Kalman filter. Its effective functioning is possible in the case of a priori certainty of the laws of distribution of external influences, which is a significant drawback.

Наиболее близким по технической сущности к заявленному изобретению является устройство идентификации параметров динамических систем на основе вариационных принципов [1]. Однако его использование в условиях неизвестных воздействий не позволяет получить требуемую точность.The closest in technical essence to the claimed invention is a device for identifying parameters of dynamic systems based on variational principles [1]. However, its use in conditions of unknown effects does not allow to obtain the required accuracy.

Цель изобретения - повысить точность идентификации параметров динамических систем. Указанный технический результат достигается за счет устройства идентификации Лагранжевых динамических систем на основе итерационной регуляризации, которое содержит следующие блоки: первый, второй блоки хранения констант; первый, второй, третий блоки формирования функций f ( x ^ i , x ˙ ^ i , z ^ i )

Figure 00000001
, H ( x ^ i , ( z ^ i ) , t )
Figure 00000002
, (1-γi-1), соответственно; первый, второй, третий, четвертый, пятый блоки формирования производной; первый, второй, третий, четвертый, пятый, шестой, седьмой, восьмой, девятый, десятый, одиннадцатый блоки формирования произведения; первый, второй, третий блоки формирования разности; первый, второй, третий блоки формирования суммы; первый, второй блоки транспонирования; первый, второй, третий, четвертый, пятый, шестой блоки интегрирования; при этом на первый информационный вход первого блока хранения констант, который является входом устройства, поступает значение наблюдаемой величины; первый информационный выход первого блока хранения констант соединен с первым информационным входом первого блока формирования произведения; второй, третий и четвертый информационные выходы первого блока хранения констант соединены соответственно с первым, вторым и третьим информационными входами первого блока формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
; второй и четвертый информационные выходы первого блока хранения констант соединены соответственно с первым и вторым информационными входами второго блока формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
; пятый информационный выход первого блока хранения констант соединен со вторыми информационными входами первого, второго, третьего и шестого блоков интегрирования, а также с первыми информационными входами четвертого и пятого блоков интегрирования; шестой информационный выход первого блока хранения констант соединен со вторым информационным входом второго блока формирования разности; седьмой информационный выход первого блока хранения констант соединен с первым информационным входом третьего блока формирования функции (1-γi-1); восьмой информационный выход первого блока хранения констант соединен со вторым информационным входом третьего блока формирования функции (1-γi-1) и третьим информационным входом первого блока формирования суммы; девятый информационный выход первого блока хранения констант соединен со вторым информационным входом третьего блока формирования произведения; десятый информационный выход первого блока хранения констант соединен с первым и вторым информационными входами соответственно шестого и девятого блоков формирования произведения; одиннадцатый информационный выход первого блока хранения констант соединен с первым информационным входом второго блока транспонирования, а также с первыми информационными входами восьмого и девятого блоков формирования произведения; двенадцатый информационный выход первого блока хранения констант соединен с первым информационным входом седьмого блока формирования произведения; тринадцатый информационный выход первого блока хранения констант соединен с первым информационным входом третьего блока формирования разности; информационный выход первого блока формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
соединен со вторым информационнымThe purpose of the invention is to increase the accuracy of identification of parameters of dynamic systems. The specified technical result is achieved due to the identification device of Lagrange dynamical systems based on iterative regularization, which contains the following blocks: first, second blocks of constant storage; first, second, third blocks of the formation of functions f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
, H ( x ^ i , ( z ^ i ) , t )
Figure 00000002
, (1-γ i-1 ), respectively; first, second, third, fourth, fifth blocks of derivative formation; first, second, third, fourth, fifth, sixth, seventh, eighth, ninth, tenth, eleventh blocks of formation of a work; first, second, third blocks of difference formation; first, second, third blocks of the formation of the amount; first, second transpose blocks; first, second, third, fourth, fifth, sixth integration blocks; at the same time, the value of the observed value is supplied to the first information input of the first constant storage unit, which is the input of the device; the first information output of the first constant storage unit is connected to the first information input of the first product formation unit; the second, third and fourth information outputs of the first constant storage unit are connected respectively to the first, second and third information inputs of the first function forming unit f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
; the second and fourth information outputs of the first constant storage unit are connected respectively to the first and second information inputs of the second function generating unit H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
; the fifth information output of the first constant storage unit is connected to the second information inputs of the first, second, third and sixth integration units, as well as to the first information inputs of the fourth and fifth integration units; the sixth information output of the first constant storage unit is connected to the second information input of the second difference generating unit; the seventh information output of the first constant storage unit is connected to the first information input of the third function generating unit (1-γ i-1 ); the eighth information output of the first constant storage unit is connected to the second information input of the third function generating unit (1-γ i-1 ) and the third information input of the first sum generating unit; the ninth information output of the first constant storage unit is connected to the second information input of the third product generating unit; the tenth information output of the first constant storage unit is connected to the first and second information inputs of the sixth and ninth product formation blocks, respectively; the eleventh information output of the first constant storage unit is connected to the first information input of the second transpose unit, as well as to the first information inputs of the eighth and ninth product formation units; the twelfth information output of the first constant storage unit is connected to the first information input of the seventh product formation unit; the thirteenth information output of the first constant storage unit is connected to the first information input of the third difference generating unit; information output of the first block forming the function f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
connected to the second information

входом первого блока формирования разности, а также с первыми информационными входами второго, третьего и четвертого блоков формирования производной; информационный выход второго блока формирования функции H ( x ^ i , ( z ^ i ) , t )

Figure 00000003
соединен с первыми информационными входами первого блока транспонирования и второго блока формирования разности, информационный выход которого соединен с первым информационным входом второго блока формирования произведения, 1 информационный выход которого соединен с первыми информационными входами второго блока хранения констант и первого блока формирования суммы; информационный выход второго блока хранения констант соединен со вторым информационным входом первого блока формирования суммы, информационный выход которого соединен с первым информационным входом третьего блока формирования произведения, информационный выход которого соединен с первым информационным входом четвертого блока формирования произведения; информационный выход первого блока транспонирования соединен с первым информационным входом первого блока формирования производной, информационный выход которого соединен со вторым информационным входом четвертого блока формирования произведения, информационный выход которого соединен соответственно со вторым информационным входом первого блока формирования произведения и первыми информационными входами пятого блока формирования произведения и пятого блока формирования производной, информационный выход которого соединен с первым информационным входом десятого блока формирования произведения; информационный выход второго блока транспонирования соединен со вторым информационным входом шестого блока формирования произведения, информационный выход которого соединен соответственно со вторыми информационными входами пятого и десятого блоков формирования произведения, информационный выход которого соединен со вторым информационным входом одиннадцатого блока формирования произведения; информационный выход девятого блока формирования произведения соединен с первым информационным входом одиннадцатого блока формирования произведения, информационный выход которого соединен со вторым информационным входом третьего блока формирования разности, информационный выход которого соединен с первым информационным входом шестого блока интегрирования, информационный выход которого соединен с первым и вторым информационными входами соответственно шестого и девятого блоков формирования произведения; информационный выход первого блока формирования произведения соединен с первым информационным входом первого блока формирования разности, информационный выход которого соединен с первым информационным входом первого блока интегрирования, информационный выход которого соединен со вторым информационным входом первого блока формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
и первым информационным входом второго блока интегрирования, информационный выход которого соединен с первыми информационными входами первого и второго блоков формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
; информационный выход пятого блока формирования произведения соединен с первым информационным входом третьего блока интегрирования, информационный выход которого является выходом устройства, с него же снимается значение искомого вектора идентификации параметров и подается на второй и третий информационные входы соответственно второго и первого блоков формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
, и f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
; информационный выход второго блока формирования производной соединен со вторым информационным входом седьмого блока формирования произведения, информационный выход которого соединен с первым информационным входом второго блока формирования суммы; информационный выход третьего блока формирования производной соединен со вторым информационным входом восьмого блока формирования произведения, информационный выход которого соединен со вторым информационным входом второго блока формирования суммы, информационный выход которого соединен с первым информационным входом третьего блока формирования суммы; информационный выход четвертого блока формирования производной соединен со вторым информационным входом третьего блока формирования суммы, информационный выход которого соединен со вторым информационным входом пятого блока интегрирования, информационный выход которого соединен с первым информационным входом седьмого блока формирования произведения и вторым информационным входом четвертого блока интегрирования, информационный выход которого соединен с первым информационным входом второго блока транспонирования, а также с первыми информационными входами восьмого и девятого блоков формирования умножения.the input of the first block of the formation of the difference, as well as with the first information inputs of the second, third and fourth blocks of the formation of the derivative; information output of the second block forming the function H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
connected to the first information inputs of the first transpose unit and the second difference forming unit, the information output of which is connected to the first information input of the second product forming unit, 1 information output of which is connected to the first information inputs of the second constant storage unit and the first sum forming unit; the information output of the second constant storage unit is connected to the second information input of the first sum forming unit, the information output of which is connected to the first information input of the third product generating unit, the information output of which is connected to the first information input of the fourth product generating unit; the information output of the first transpose unit is connected to the first information input of the first derivative formation unit, the information output of which is connected to the second information input of the fourth product formation unit, the information output of which is connected respectively to the second information input of the first product formation unit and the first information inputs of the fifth product formation unit and the fifth derivative formation unit, the information output of which is connected to ne first information input of the tenth block of the formation of the work; the information output of the second transpose unit is connected to the second information input of the sixth product formation unit, the information output of which is connected respectively to the second information inputs of the fifth and tenth product formation units, the information output of which is connected to the second information input of the eleventh product formation unit; the information output of the ninth product formation unit is connected to the first information input of the eleventh product formation unit, the information output of which is connected to the second information input of the third difference forming unit, the information output of which is connected to the first information input of the sixth integration unit, the information output of which is connected to the first and second information the inputs of the sixth and ninth blocks of the formation of the work, respectively; the information output of the first product formation unit is connected to the first information input of the first difference formation unit, the information output of which is connected to the first information input of the first integration unit, the information output of which is connected to the second information input of the first function formation unit f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
and the first information input of the second integration unit, the information output of which is connected to the first information inputs of the first and second function formation units H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
; the information output of the fifth product formation unit is connected to the first information input of the third integration unit, the information output of which is the output of the device, the value of the sought parameter identification vector is taken from it and fed to the second and third information inputs of the second and first function formation units, respectively H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
, and f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
; the information output of the second derivative formation unit is connected to the second information input of the seventh product formation unit, the information output of which is connected to the first information input of the second sum formation unit; the information output of the third derivative formation unit is connected to the second information input of the eighth product formation unit, the information output of which is connected to the second information input of the second sum formation unit, the information output of which is connected to the first information input of the third sum formation unit; the information output of the fourth derivative formation unit is connected to the second information input of the third sum formation unit, the information output of which is connected to the second information input of the fifth integration unit, the information output of which is connected to the first information input of the seventh product formation unit and the second information input of the fourth integration unit, information output which is connected to the first information input of the second transpose block, as well as to the first and the data inputs of the eighth and ninth multiplication block formation.

Сущность изобретения поясняется чертежами, где на фиг.1 представлено устройство идентификации параметров Лагранжевых динамических систем на основе итерационной регуляризации, которое содержит:The invention is illustrated by drawings, where figure 1 shows a device for identifying parameters of Lagrange dynamic systems based on iterative regularization, which contains:

1.1 - первый блок хранения констант;1.1 - the first block storage constants;

1.2 - второй блок хранения констант;1.2 - the second block storage constants;

2.1- первый блок формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )

Figure 00000001
;2.1- first block of function formation f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
;

2.2 - второй блок формирования функции H ( x ^ i , ( z ^ i ) , t )

Figure 00000003
;2.2 - the second block of the formation of the function H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
;

2.3 - третий блок формирования функции (1-γi1);2.3 - the third block of the formation of the function (1-γ i1 );

3.1 - первый блок формирования производной;3.1 - the first block of derivative formation;

3.2 - второй блок формирования производной;3.2 - the second derivative formation unit;

3.3 - третий блок формирования производной;3.3 - the third block of derivative formation;

3.4 - четвертый блок формирования производной;3.4 - the fourth block of derivative formation;

3.5 - пятый блок формирования производной;3.5 - fifth block of derivative formation;

4.1 - первый блок формирования произведения;4.1 - the first block of the formation of the work;

4.2 - второй блок формирования произведения;4.2 - the second block of the formation of the work;

4.3 - третий блок формирования произведения;4.3 - the third block of the formation of the work;

4.4 - четвертый блок формирования произведения;4.4 - the fourth block of the formation of the work;

4.5 - пятый блок формирования произведения;4.5 - the fifth block of the formation of the work;

4.6 - шестой блок формирования произведения;4.6 - the sixth block of the formation of the work;

4.7 - седьмой блок формирования произведения;4.7 - the seventh block of the formation of the work;

4.8 - восьмой блок формирования произведения;4.8 - the eighth block of the formation of the work;

4.9 - девятый блок формирования произведения;4.9 - the ninth block of the formation of the work;

4.10 - десятый блок формирования произведения;4.10 - the tenth block of the formation of the work;

4.11 - одиннадцатый блок формирования произведения;4.11 - the eleventh block of the formation of the work;

5.1 - первый блок формирования разности;5.1 - the first block forming the difference;

5.2 - второй блок формирования разности;5.2 - the second block forming the difference;

5.3 - третий блок формирования разности;5.3 - the third block of the formation of the difference;

6.1 - первый блок формирования суммы;6.1 - the first block of the formation of the amount;

6.2 - второй блок формирования суммы;6.2 - the second block of the formation of the amount;

6.3 - третий блок формирования суммы;6.3 - the third block of the formation of the amount;

7.1 - первый блок транспонирования;7.1 - the first transpose block;

7.2 - второй блок транспонирования;7.2 - the second transpose block;

8.1 - первый блок интегрирования;8.1 - the first integration unit;

8.2 - второй блок интегрирования;8.2 - the second integration unit;

8.3 - третий блок интегрирования;8.3 - the third integration unit;

8.4 - четвертый блок интегрирования;8.4 - the fourth integration unit;

8.5 - пятый блок интегрирования;8.5 - the fifth integration unit;

8.6 - шестой блок интегрирования.8.6 is the sixth integration unit.

На фиг.2 представлен первый блок интегрирования:Figure 2 presents the first integration unit:

8.1.1 - блок формирования произведения;8.1.1 - the unit of formation of the work;

8.1.2 - блок формирования суммы;8.1.2 - the unit for the formation of the amount;

8.1.3 - линия задержки.8.1.3 - delay line.

Устройство идентификации Лагранжевых динамических систем на основе итерационной регуляризации работает следующим образом: на первый информационный вход первого блока 1.1 хранения констант (фиг.1) поступает значение наблюдаемой величины; первый информационный выход первого блока 1.1 хранения констант соединен с первым информационным входом первого блока 4.1 формирования произведения; второй, третий и четвертый информационные выходы первого блока 1.1 хранения констант соединены соответственно с первым, вторым и третьим информационными входами первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )

Figure 00000004
; второй и четвертый информационные выходы первого блока 1.1 хранения констант соединены соответственно с первым и вторым информационными входами второго блока 2.2 формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
; пятый информационный выход первого блока 1.1 хранения констант соединен со вторыми информационными входами первого блока 8.1, второго блока 8.2, третьего блока 8.3 и шестого блока 8.6 интегрирования, а также с первыми информационными входами четвертого блока 8.4 и пятого блока 8.5 интегрирования; шестой информационный выход первого блока 1.1 хранения констант соединен со вторым информационным входом второго блока 5.2 формирования разности; седьмой информационный выход первого блока 1.1 хранения констант соединен с первым информационным входом третьего блока 2.3 формирования функции (1-γi-1); восьмой информационный выход первого блока 1.1 хранения констант соединен со вторым информационным входом третьего блока 2.3 формирования функции (1-γi-1) и третьим информационным входом первого блока 6.1 формирования суммы; девятый информационный выход первого блока 1.1 хранения констант соединен со вторым информационным входом третьего блока 4.3 формирования произведения; десятый информационный выход первого блока 1.1 хранения констант соединен с первым и вторым информационными входами соответственно шестого блока 4.6 и девятого блока 4.9 формирования произведения; одиннадцатый информационный выход первого блока 1.1 хранения констант соединен с первым информационным входом второго блока 7.2 транспонирования, а также с первыми информационными входами восьмого блока 4.8 и девятого блока 4.9 формирования произведения; двенадцатый информационный выход первого блока 1.1 хранения констант соединен с первым информационным входом седьмого блока 4.7 формирования произведения; тринадцатый информационный выход первого блока 1:1 хранения констант соединен с первым информационным входом третьего блока 5.3 формирования разности; информационный выход первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
соединен со вторым информационным входом первого блока 5.1 формирования разности, а также с первыми информационными входами второго блока 3.2, третьего блока 3.3 и четвертого блока 3.4 формирования производной; информационный выход второго блока 2.2 формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
соединен с первыми информационными входами первого блока 7.1 транспонирования и второго блока 5.2 формирования разности, информационный выход которого соединен с первым информационным входом второго блока 4.2 формирования произведения, информационный выход которого соединен с первыми информационными входами второго блока 1.2 хранения констант и первого блока 6.1 формирования суммы; информационный выход второго блока 1.2 хранения констант соединен со вторым информационным входом первого блока 6.1 формирования суммы, информационный выход которого соединен с первым информационным входом третьего блока 4.3 формирования произведения, информационный выход которого соединен с первым информационным входом четвертого блока 4.4 формирования произведения; информационный выход первого блока 7.1 транспонирования соединен с первым информационным входом первого блока 3.1 формирования производной, информационный выход которого соединен со вторым информационным входом четвертого блока 4.4 формирования произведения, информационный выход которого соединен соответственно со вторым информационным входом первого блока 4.1 формирования произведения и первыми информационными входами пятого блока 4.5 формирования произведения и пятого блока 3.5 формирования производной, информационный выход которого соединен с первым информационным входом десятого блока 4.10 формирования произведения; информационный выход второго блока 7.2 транспонирования соединен со вторым информационным входом шестого блока 4.6 формирования произведения, информационный выход которого соединен соответственно со вторыми информационными входами пятого блока 4.5 и десятого блока 4.10 формирования произведения, информационный выход которого соединен со вторым информационным входом одиннадцатого блока 4.11 формирования произведения; информационный выход девятого блока 4.9 формирования произведения соединен с первым информационным входом одиннадцатого блока 4.11 формирования произведения, информационный выход которого соединен со вторым информационным входом третьего блока 5.3 формирования разности, информационный выход которого соединен с первым информационным входом шестого блока 8.6 интегрирования, информационный выход которого соединен с первым и вторым информационными входами соответственно шестого блока 4.6 и девятого блока 4.9 формирования произведения; информационный выход первого блока 4.1 формирования умножения соединен с первым информационным входом первого блока 5.1 формирования разности, информационный выход которого соединен с первым информационным входом первого блока 8.1 интегрирования, информационный выход которого соединен со вторым информационным входом первого блока 2.1 формирования функции f(xitihz.i) и первым информационным входом второго блока 8.2 интегрирования, информационный выход которого соединен с первыми информационными входами первого блока 2.1 и второго блока 2.2 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
и H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
, соответственно; информационный выход пятого блока 4.5 формирования произведения соединен с первым информационным входом третьего блока 8.3 интегрирования, с информационного выхода которого снимается значение искомого вектора идентификации параметров z и подается на второй и третий информационные входы соответственно второго блока 2.2 и первого блока 2.1 формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
и f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
соответственно; информационный выход второго блока 3.2 формирования производной соединен со вторым информационным входом седьмого блока 4.7 формирования произведения, информационный выход которого соединен с первым информационным входом второго блока 6.2 формирования суммы; информационный выход третьего блока 3.3 формирования производной соединен со вторым информационным входом восьмого блока 4.8 формирования произведения, информационный выход которого соединен со вторым информационным входом второго блока 6.2 формирования суммы, информационный выход которого соединен с первым информационным входом третьего блока 6.3 формирования суммы; информационный выход четвертого блока 3.4 формирования производной соединен со вторым информационным входом третьего блока 6.3 формирования суммы, информационный выход которого соединен со вторым информационным входом пятого блока 8.5 интегрирования, информационный выход которого соединен с первым информационным входом седьмого блока 4.7 формирования произведения и вторым информационным входом четвертого блока 8.4 интегрирования, информационный выход которого соединен с первым информационным входом второго блока 7.2 транспонирования, а также с первыми информационными входами восьмого блока 4.8 и девятого блока 4.9 формирования произведения.The identification device of Lagrangian dynamical systems based on iterative regularization works as follows: the value of the observed value is supplied to the first information input of the first constant storage unit 1.1 (Fig. 1); the first information output of the first constant storage unit 1.1 is connected to the first information input of the first product generating unit 4.1; the second, third and fourth information outputs of the first constant storage unit 1.1 are connected respectively to the first, second and third information inputs of the first function generating unit 2.1 f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000004
; the second and fourth information outputs of the first constant storage unit 1.1 are connected respectively to the first and second information inputs of the second function generating unit 2.2 H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
; the fifth information output of the first constant storage unit 1.1 is connected to the second information inputs of the first block 8.1, the second block 8.2, the third block 8.3 and the sixth integration block 8.6, as well as the first information inputs of the fourth block 8.4 and the fifth integration block 8.5; the sixth information output of the first constant storage unit 1.1 is connected to the second information input of the second difference generating unit 5.2; the seventh information output of the first constant storage unit 1.1 is connected to the first information input of the third function generation unit 2.3 (1-γ i-1 ); the eighth information output of the first constant storage unit 1.1 is connected to the second information input of the third function generation unit 2.3 (1-γ i-1 ) and the third information input of the first sum formation unit 6.1; the ninth information output of the first constant storage unit 1.1 is connected to the second information input of the third product generating unit 4.3; the tenth information output of the first constant storage unit 1.1 is connected to the first and second information inputs of the sixth block 4.6 and the ninth block 4.9 of product formation, respectively; the eleventh information output of the first constant storage unit 1.1 is connected to the first information input of the second transposition unit 7.2, as well as to the first information inputs of the eighth block 4.8 and the ninth block 4.9 of creating the product; the twelfth information output of the first constant storage unit 1.1 is connected to the first information input of the seventh product formation block 4.7; the thirteenth information output of the first constant storage unit 1: 1 is connected to the first information input of the third difference generating unit 5.3; information output of the first block 2.1 function formation f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
connected to the second information input of the first difference generating unit 5.1, as well as to the first information inputs of the second block 3.2, the third block 3.3 and the fourth derivative formation block 3.4; information output of the second block 2.2 function formation H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
connected to the first information inputs of the first transposition unit 7.1 and the second difference generating unit 5.2, the information output of which is connected to the first information input of the second product generating unit 4.2, the information output of which is connected to the first information inputs of the second constant storage unit 1.2 and the first sum forming unit 6.1; the information output of the second constant storage unit 1.2 is connected to the second information input of the first sum forming unit 6.1, the information output of which is connected to the first information input of the third work generating unit 4.3, the information output of which is connected to the first information input of the fourth work forming unit 4.4; the information output of the first transposition unit 7.1 is connected to the first information input of the first derivative formation unit 3.1, the information output of which is connected to the second information input of the fourth product formation unit 4.4, the information output of which is connected respectively to the second information input of the first product formation unit 4.1 and the first information inputs of the fifth block 4.5 formation of the product and the fifth block 3.5 of the formation of the derivative, the information output of which connected to the first information input of the tenth block 4.10 formation of the product; the information output of the second transposition unit 7.2 is connected to the second information input of the sixth creation unit 4.6, the information output of which is connected respectively to the second information inputs of the fifth block 4.5 and the tenth creation unit 4.10, the information output of which is connected to the second information input of the eleventh composition formation block 4.11; the information output of the ninth product generation unit 4.9 is connected to the first information input of the eleventh product formation unit 4.11, the information output of which is connected to the second information input of the third difference generating unit 5.3, the information output of which is connected to the first information input of the sixth integration unit 8.6, the information output of which is connected to the first and second information inputs, respectively, of the sixth block 4.6 and the ninth block 4.9 of the formation of the work; the information output of the first block 4.1 of the formation of multiplication is connected to the first information input of the first block 5.1 of the formation of the difference, the information output of which is connected to the first information input of the first block 8.1 of integration, the information output of which is connected to the second information input of the first block 2.1 of the formation of function f (xitihz.i) and the first information input of the second integration block 8.2, the information output of which is connected to the first information inputs of the first block 2.1 and the second block 2.2 of the form IAOD function f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
and H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
, respectively; the information output of the fifth block 4.5 of product formation is connected to the first information input of the third integration block 8.3, from the information output of which the desired parameter identification vector z is taken and fed to the second and third information inputs, respectively, of the second block 2.2 and the first block 2.1 of the function formation H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
and f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
respectively; the information output of the second derivative generating unit 3.2 is connected to the second information input of the seventh product generating unit 4.7, the information output of which is connected to the first information input of the second sum forming unit 6.2; the information output of the third derivative generating unit 3.3 is connected to the second information input of the eighth product generating unit 4.8, the information output of which is connected to the second information input of the second sum forming unit 6.2, the information output of which is connected to the first information input of the third sum forming unit 6.3; the information output of the fourth derivative formation block 3.4 is connected to the second information input of the third sum formation block 6.3, the information output of which is connected to the second information input of the fifth integration unit 8.5, the information output of which is connected to the first information input of the seventh product formation block 4.7 and the second information input of the fourth block 8.4 integration, the information output of which is connected to the first information input of the second transposition block 7.2, and ose from the first information input unit of the eighth and ninth block 4.8 4.9 product formation.

Информационный выход первого блока 5.1 формирования разности соединен с первым информационным входом блока 8.1.1 формирования произведения (фиг.2); пятый информационный выход первого блока 1.1 хранения констант соединен со вторым информационным входом первого 8.1 блока интегрирования, информационный выход которого соединен с первым информационным входом блока 8.1.2 формирования суммы; информационный выход линии задержки 8.1.3 соединен со вторым информационным входом блока 8.1.2 формирования суммы, с информационного выхода которого снимается значение первообразной и подается на первый вход линии задержки 8.1.3 и на второй информационный вход первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )

Figure 00000001
и первыйThe information output of the first block 5.1 of the formation of the difference is connected to the first information input of the block 8.1.1 of the formation of the product (figure 2); the fifth information output of the first constant storage unit 1.1 is connected to the second information input of the first 8.1 integration unit, the information output of which is connected to the first information input of the summing unit 8.1.2; the information output of the delay line 8.1.3 is connected to the second information input of the summing block 8.1.2, the antiderivative value is taken from the information output of which is fed to the first input of the delay line 8.1.3 and to the second information input of the first block 2.1 of the function formation f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
and first

информационный вход второго блока 8.2 интегрирования.information input of the second integration unit 8.2.

Блоки 8.2, 8.3, 8.4, 8.5, 8.6 интегрирования (фиг.1) имеют структуру и связи, аналогичные блоку 8.1, рассмотренному выше.Blocks 8.2, 8.3, 8.4, 8.5, 8.6 integration (figure 1) have a structure and relationships similar to block 8.1, discussed above.

Заявленное изобретение направлено на повышение эффективности идентификации параметров динамических систем, что весьма важно на всех этапах создания, экспериментальной обработки и эксплуатации объектов ракетно-космической, авиационной, корабельной и других видов техники.The claimed invention is aimed at increasing the efficiency of identification of parameters of dynamic systems, which is very important at all stages of the creation, experimental processing and operation of objects of rocket and space, aviation, ship and other types of equipment.

Приведенные результаты позволяют утверждать, что идентификация параметров на основе вариационных принципов с использованием итерационной регуляризации обеспечивает увеличение точности оценок в сравнении с известным фильтром Калмана.The above results suggest that the identification of parameters based on variational principles using iterative regularization provides an increase in the accuracy of estimates in comparison with the well-known Kalman filter.

Структура предлагаемого устройства определяется решением задачи, чья постановка приведена ниже.The structure of the proposed device is determined by the solution of the problem, whose statement is given below.

Принято, что динамическая система описывается следующим векторным уравнением дифференциальным уравнениемIt is accepted that a dynamical system is described by the following vector equation by the differential equation

x ¨ = f ( x , x ˙ , z ) , x ( 0 ) = x , x ˙ ( 0 ) = x , ( 1 )

Figure 00000005
x ¨ = f ( x , x ˙ , z ) , x ( 0 ) = x , x ˙ ( 0 ) = x , ( one )
Figure 00000005

где x∈Rn - вектор обобщенных координат,where x∈R n is the vector of generalized coordinates,

x R n

Figure 00000006
- вектор обобщенных скоростей, x R n
Figure 00000006
is the vector of generalized velocities,

z∈Rm - вектор неизвестных постоянных параметров, подлежащих идентификации,z∈R m is the vector of unknown constant parameters to be identified,

f - вектор-функция, непрерывная и дифференцируемая по совокупности аргументов,f is a vector function that is continuous and differentiable by the totality of arguments,

n, m - натуральные числа,n, m are natural numbers,

T∈[0,T].T∈ [0, T].

Уравнение наблюдения имеет видThe observation equation has the form

y = H ( x , t ) + n ( t ) , ( 2 )

Figure 00000007
y = H ( x , t ) + n ( t ) , ( 2 )
Figure 00000007

где y∈Rk - вектор наблюдения,where y∈R k is the observation vector,

H(x, f) - непрерывная вместе с частными производными вектор-функция,H (x, f) is a continuous vector-function with partial derivatives,

k - натуральное число,k is a natural number,

n(t) - вектор белого гауссовского шума.n (t) is the vector of white Gaussian noise.

Пусть динамика идентифицируемых параметров z определяется уравнениемLet the dynamics of identifiable parameters z be determined by the equation

z ˙ = η , z ( 0 ) = z , ( 3 )

Figure 00000008
z ˙ = η , z ( 0 ) = z , ( 3 )
Figure 00000008

где η∈Rm - вектор неизвестных неслучайных возмущений, удовлетворяющий требованиям физической реализуемости η ( t ) L 2 m [ 0, T ]

Figure 00000009
.where η∈R m is the vector of unknown nonrandom disturbances satisfying the requirements of physical realizability η ( t ) L 2 m [ 0 T ]
Figure 00000009
.

Пусть динамика исследуемой системы определяется принципом Гамильтона-Остроградского. В соответствии с нимLet the dynamics of the system under study be determined by the Hamilton-Ostrogradsky principle. According to him

δ ' W = δ S + 0 T δ ' A d t = 0 T [ δ L + Q δ x ] d t = 0, ( 4 )

Figure 00000010
δ '' W = δ S + 0 T δ '' A d t = 0 T [ δ L + Q δ x ] d t = 0 ( four )
Figure 00000010

где S - действие по Гамильтону на интервале времени [0, T],where S is the Hamilton action on the time interval [0, T],

L - кинетический потенциал,L is the kinetic potential,

А - работа вектора обобщенных внешних сил Q∈Rn,And - the work of the vector of generalized external forces Q∈R n ,

δx - вектор вариаций обобщенных координат, знак δ' обозначает бесконечно малую величину, которая не является вариацией.δx is the vector of variations of the generalized coordinates, the sign of δ 'denotes an infinitely small quantity, which is not a variation.

Ставится задача определения оценки z ^

Figure 00000011
вектора z из условия минимума функционала невязкиThe task of determining the assessment z ^
Figure 00000011
vector z from the minimum condition for the residual functional

J 1 = 1 2 0 T [ y H ( x ^ ( z ^ ) , t ) ] T N 1 [ y H ( x ^ ( z ^ ) , t ) ] d t min , ( 5 )

Figure 00000012
J one = one 2 0 T [ y - H ( x ^ ( z ^ ) , t ) ] T N - one [ y - H ( x ^ ( z ^ ) , t ) ] d t min , ( 5 )
Figure 00000012

где N∈Rm×Rm - матрица односторонней спектральной плотности шума наблюдения, характеризующая интенсивность помех в канале наблюдений.where N∈R m × R m is the matrix of one-sided spectral density of the observation noise, which characterizes the intensity of interference in the observation channel.

Задача идентификации параметров динамических систем (1)-(5) является обратной, некорректно поставленной по Адамару. Для ее решения используется метод регуляризации А.Н. Тихонова.The problem of identifying the parameters of dynamical systems (1) - (5) is the inverse, incorrectly posed by Hadamard. To solve it, the regularization method A.N. Tikhonov.

Это требует рассмотрения условий минимума сглаживающего функционалаThis requires consideration of the minimum smoothing functional conditions.

J α [ z , x ] = J 1 + μ W + α Ω [ η ] , ( 6 )

Figure 00000013
J α [ z , x ] = J one + μ W + α Ω [ η ] , ( 6 )
Figure 00000013

где Ω [ η ] = 1 2 0 T η T ( t ) η ( t ) d t

Figure 00000014
- стабилизирующий функционал,Where Ω [ η ] = one 2 0 T η T ( t ) η ( t ) d t
Figure 00000014
- stabilizing functional,

α - положительное число.α is a positive number.

Функционал Jα[z, х] неотрицательный, поэтому существует его нижняя грань, и если {αn} - убывающая последовательность положительных чисел, сходящихся к нулю (α→0), и соответствующая последовательность {ηn(t)}, согласно, сходится к η*(t) при условии, что числовой параметр α удовлетворяет требованиямThe functional J α [z, x] is non-negative, therefore there exists its lower bound, and if {α n } is a decreasing sequence of positive numbers converging to zero (α → 0), and the corresponding sequence {η n (t)}, according to, converges to η * (t) provided that the numerical parameter α satisfies the requirements

lim n α n = 0, lim n k = 1 n α k = , lim n k = 1 n α k < .

Figure 00000015
lim n α n = 0 lim n k = one n α k = , lim n k = one n α k < .
Figure 00000015

Для дифференцируемого выпуклого функционала Jα[z, x] минимум находится путем определения стационарной точки, в которой For a differentiable convex functional J α [z, x], the minimum is found by determining the stationary point at which

g r a d z J α [ z ^ , x ^ ] = g r a n d x J α [ z ^ , x ^ ] = 0, ( 7 )

Figure 00000016
g r a d z J α [ z ^ , x ^ ] = g r a n d x J α [ z ^ , x ^ ] = 0 ( 7 )
Figure 00000016

при ограничениях (1) - (3).under the constraints (1) - (3).

Для определения градиента в (7) воспользуемся игольчатым варьированием возмущения η(t) и асинхронным варьированием траектории x(t).To determine the gradient in (7), we use needle-like variation of the perturbation η (t) and asynchronous variation of the trajectory x (t).

Игольчатой вариацией возмущения будем называть функцию видаThe needle variation of the perturbation will be called a function of the form

η ε ( t ) = { v , t [ τ , τ + ε l ] , v * ( t ) , t [ τ , τ + ε l ] .

Figure 00000017
η ε ( t ) = { v , t [ τ , τ + ε l ] , v * ( t ) , t [ τ , τ + ε l ] .
Figure 00000017

где τ - заданная произвольная точка непрерывности функции η*(t),where τ is a given arbitrary point of continuity of the function η * (t),

v - некоторый постоянный вектор,v is some constant vector,

l -заданное положительное число,l is a given positive number,

ε - произвольное положительное число, такое, что τ+εl<T.ε is an arbitrary positive number such that τ + εl <T.

Разность ηε(t)-η*(t)=δη(t) будем называть игольчатой вариацией возмущения, которая порождает вариацию параметров δz(t), вызывающую вариацию траектории δx(t). Обозначим xε(t) - траекторию, соответствующую ηε(t), x*(t) - соответствующую Вариацию траектории вида δxc(t)=xε(t)-x*(t) будем называть синхронной, а вариациюThe difference η ε (t) -η * (t) = δη (t) will be called the needle-like variation of the perturbation, which generates the variation of the parameters δz (t), causing the variation of the trajectory δx (t). Denote by x ε (t) the trajectory corresponding to η ε (t), x * (t) the corresponding variation of the trajectory of the form δx c (t) = x ε (t) -x * (t) will be called synchronous, and the variation

δ x a ( t ) = x ε ( t + d t ) x * ( t ) = x z δ z = G δ z

Figure 00000018
δ x a ( t ) = x ε ( t + d t ) - x * ( t ) = x z δ z = G δ z
Figure 00000018

асинхронной согласно терминологии.asynchronous according to terminology.

Здесь G - матрица чувствительности системы (1) по вектору параметров z, удовлетворяющая уравнениюHere G is the sensitivity matrix of system (1) with respect to the vector of parameters z, satisfying the equation

G ¨ = f x ˙ G ˙ + f x G + f x , G ( 0 ) = G ˙ ( 0 ) = 0

Figure 00000019
G ¨ = f x ˙ G ˙ + f x G + f x , G ( 0 ) = G ˙ ( 0 ) = 0
Figure 00000019

Вариации δxc(t) и δxa(t) обладают следующими свойствамиThe variations δx c (t) and δx a (t) have the following properties

δxc(t)=0, t∈[0,τ+εl),δx c (t) = 0, t∈ [0, τ + εl),

δxc(t)=0, t∉[0,τ+εl),δx c (t) = 0, t∉ [0, τ + εl),

δxc(t)=δxa(t), t+τ+εl.δx c (t) = δx a (t), t + τ + εl.

Пусть η*(t) - возмущение, доставляющее минимум функционалу (6). Найдем приращение функционала (6), обусловленное игольчатой вариацией ηε(t). Заметим, что на интервале [τ,τ+εl] за счет скачка ηε(t) возникает асинхронная вариация δxa(t), а на интервале [τ+εl,T] возникает синхронная δxc(t), поскольку приращения обобщенных координатLet η * (t) be a perturbation delivering a minimum to functional (6). We find the increment of functional (6) due to the needle-shaped variation of η ε (t). Note that on the interval [τ, τ + εl] due to the jump η ε (t), an asynchronous variation δx a (t) occurs, and on the interval [τ + εl, T], the synchronous δx c (t) arises, since the increments of the generalized coordinates

определяются решением дифференциальных уравнений в вариациях при начальных условиях в момент времени t=τ+εl. Таким образом, имеем are determined by solving differential equations in variations under initial conditions at time t = τ + εl. Thus, we have

Δ J = τ τ + ε l { [ ( y H ( x , t ) ) T N 1 H x ] δ x a + α η T δ η + μ δ L a = Q δ x a } d t + τ + ε l T { [ ( y H ( x , t ) ) T N 1 H x ] δ x c + μ δ L c + Q δ x c } d t ( 8 )

Figure 00000020
Δ J = τ τ + ε l { [ - ( y - H ( x , t ) ) T N - one H x ] δ x a + α η T δ η + μ δ L a = Q δ x a } d t + τ + ε l T { [ - ( y - H ( x , t ) ) T N - one H x ] δ x c + μ δ L c + Q δ x c } d t ( 8 )
Figure 00000020

где δLa и δLc - вариации кинетического потенциала, связанные соответственное δxa(t) и δxc(Y)-where δL a and δL c are the kinetic potential variations related respectively by δx a (t) and δx c (Y) -

Учитывая, что z ˙ = η

Figure 00000021
, выполняя преобразования с учетомGiven that z ˙ = η
Figure 00000021
performing transformations considering

τ τ + ε l η T δ z ˙ d t = η T δ z | τ τ + ε l τ τ + ε l η ˙ T δ z d t = τ τ + ε l η ˙ T δ z d t , ( 9 )

Figure 00000022
τ τ + ε l η T δ z ˙ d t = η T δ z | | | τ τ + ε l - τ τ + ε l η ˙ T δ z d t = - τ τ + ε l η ˙ T δ z d t , ( 9 )
Figure 00000022

и поскольку δz(τ)=0, η(τ+εl)=0 получимand since δz (τ) = 0, η (τ + εl) = 0, we obtain

Δ J = τ τ + ε l [ ( y H ( x , t ) ) T N 1 H x G + α η ˙ T ] δ z d t + τ + ε l T ( y H ( x , t ) ) T N 1 H x μ [ x ¨ f ( x , x ˙ , z ) ] δ x c d t ( 10 )

Figure 00000023
Δ J = τ τ + ε l [ - ( y - H ( x , t ) ) T N - one H x G + α η ˙ T ] δ z d t + τ + ε l T - ( y - H ( x , t ) ) T N - one H x - μ [ x ¨ - f ( x , x ˙ , z ) ] δ x c d t ( 10 )
Figure 00000023

откуда следует, что в точке минимума функционалаwhence it follows that at the minimum point of the functional

g r a d z J α [ z ^ , x ^ ] = ( y H ( x ^ , τ ) ) T N 1 H x G + α η ˙ * T = 0 g r a d z J α [ z ^ , x ^ ] = ( y H ( x ^ , t ) ) T N 1 H x μ [ x ¨ ^ f ( x , x ˙ ^ , z ) ] = 0 ( 11 )

Figure 00000024
g r a d z J α [ z ^ , x ^ ] = - ( y - H ( x ^ , τ ) ) T N - one H x G + α η ˙ * T = 0 g r a d z J α [ z ^ , x ^ ] = - ( y - H ( x ^ , t ) ) T N - one H x - μ [ x ¨ ^ - f ( x , x ˙ ^ , z ) ] = 0 ( eleven )
Figure 00000024

Используя эти условия стационарности получимUsing these stationarity conditions, we obtain

η * ( τ ) = α 1 τ T N 1 G T H x [ y H ( x ^ ( z ^ , τ ) , t ) ] d t , x = f ( x ^ , x ˙ ^ , z ) μ 1 H x N 1 ( y H ( x ^ ( z ^ , τ ) t ) ) T . ( 12 )

Figure 00000025
η * ( τ ) = α - one τ T N - one G T H x [ y - H ( x ^ ( z ^ , τ ) , t ) ] d t , x = f ( x ^ , x ˙ ^ , z ) - μ - one H x N - one ( y - H ( x ^ ( z ^ , τ ) t ) ) T . ( 12 )
Figure 00000025

Найти точку [ z ^ , η * ]

Figure 00000026
путем решения данных уравнений в условиях некорректности исходной задачи достаточно сложно, поэтому широкое распространение получили методы последовательного приближения к стационарной точке. Рассмотрим метод простой итерацииFind point [ z ^ , η * ]
Figure 00000026
by solving these equations under conditions of incorrectness of the original problem, it is quite difficult, therefore, methods of successive approximation to a stationary point are widespread. Consider a simple iteration method

η k + 1 = η k α k g r a n d z J α [ z ^ k + 1 , η k ] . ( 13 )

Figure 00000027
η k + one = η k - α k g r a n d z J α [ z ^ k + one , η k ] . ( 13 )
Figure 00000027

Для сокращения записи вводится обозначениеTo shorten the entry, the notation

B ( z ^ , τ ) = N 1 G T H T x [ y H ( x ˙ ^ ( z ^ , τ ) , t ) ] .

Figure 00000028
B ( z ^ , τ ) = N - one G T H T x [ y - H ( x ˙ ^ ( z ^ , τ ) , t ) ] .
Figure 00000028

Тогда с учетом (11) алгоритм (13) может быть представлен в видеThen, taking into account (11), algorithm (13) can be represented as

η k + 1 ( τ ) = η k ( τ ) [ 1 α k ] = τ T B ( z k + 1 ( τ ) , t ) d t ,

Figure 00000029
η k + one ( τ ) = η k ( τ ) [ one - α k ] = τ T B ( z k + one ( τ ) , t ) d t ,
Figure 00000029

Приняв за нулевое приближение η0=0, итерационную последовательность в развернутой форме согласно запишем следующим образомTaking η 0 = 0 as the zeroth approximation, we write the iterative sequence in the expanded form as follows

η 0 ( τ ) = 0, z ^ ˙ = η 0 ( τ ) , x ^ ¨ 0 = f ( x ^ i , x ˙ ^ i , z ^ i ) μ 1 H x N 1 ( y H ( x ^ ( z ^ , τ ) , t ) ) T , η 1 ( τ ) = τ T B ( x ^ 1 , x ˙ ^ 1 , z ^ 1 ) μ 1 H x 1 N 1 ( y H ( x ^ 1 ( z ^ 1 , τ ) , t ) ) T , x ^ ¨ 2 = f ( x ^ i , x ˙ ^ i , z ^ i ) μ 1 H x 2 N 1 ( y H ( x ^ ( z ^ , τ ) , t ) ) ( 14 ) η k ( τ ) = τ T i = 1 k γ i k B ( x ^ i ( z ^ i , τ ) , t ) d t , z ^ ˙ = η k ( τ ) x ^ ¨ 0 = f ( x ^ i , x ^ ˙ i , z ^ i ) μ 1 H x N 1 ( y H ( x ^ ( z ^ , τ ) , t ) ) T , x ^ i ( t 0 ) = x ( 0 ) , x ^ ˙ i ( t 0 ) = x ˙ ( 0 ) z ^ i ( t 0 ) = z ( 0 ) , τ [ 0, T ] ,

Figure 00000030
η 0 ( τ ) = 0 z ^ ˙ = η 0 ( τ ) , x ^ ¨ 0 = f ( x ^ i , x ˙ ^ i , z ^ i ) - μ - one H x N - one ( y - H ( x ^ ( z ^ , τ ) , t ) ) T , η one ( τ ) = τ T B ( x ^ one , x ˙ ^ one , z ^ one ) - μ - one H x one N - one ( y - H ( x ^ one ( z ^ one , τ ) , t ) ) T , x ^ ¨ 2 = f ( x ^ i , x ˙ ^ i , z ^ i ) - μ - one H x 2 N - one ( y - H ( x ^ ( z ^ , τ ) , t ) ) ... ( fourteen ) η k ( τ ) = - τ T i = one k γ i k B ( x ^ i ( z ^ i , τ ) , t ) d t , z ^ ˙ = η k ( τ ) x ^ ¨ 0 = f ( x ^ i , x ^ ˙ i , z ^ i ) - μ - one H x N - one ( y - H ( x ^ ( z ^ , τ ) , t ) ) T , x ^ i ( t 0 ) = x ( 0 ) , x ^ ˙ i ( t 0 ) = x ˙ ( 0 ) z ^ i ( t 0 ) = z ( 0 ) , τ [ 0 T ] ,
Figure 00000030

где γ i k

Figure 00000031
определяется по правилу γ i k = [ 1 α i 1 ] [ 1 α k 1 ] , i = 1, k ¯ .
Figure 00000032
Where γ i k
Figure 00000031
determined by the rule γ i k = [ one - α i - one ] ... [ one - α k - one ] , i = one, k ¯ .
Figure 00000032

Пусть задача оценки z ^ k ( τ )

Figure 00000033
процесса (1)-(3) рассматривается для момента времени τ, когда желательно получить оценку, соответствующую точке T интервала наблюдения [0,T], а T увеличивается. Требуется получить z ^ k
Figure 00000034
как функцию от Г, доставляющую минимум критерию (5). Для решения данной задачи преобразуем уравнение (14) к виду двухточечной краевой задачи. Для этого продифференцируем правую часть уравнения (14) по параметру τ в результате преобразований с учетом правила Лейбница получимLet the estimation problem z ^ k ( τ )
Figure 00000033
of process (1) - (3) is considered for a point in time τ when it is desirable to obtain an estimate corresponding to the point T of the observation interval [0, T], and T increases. Required to get z ^ k
Figure 00000034
as a function of Γ, delivering a minimum to criterion (5). To solve this problem, we transform equation (14) to the form of a two-point boundary value problem. To do this, we differentiate the right-hand side of equation (14) with respect to the parameter τ as a result of transformations, taking into account the Leibniz rule, we obtain

d η k ( τ ) d τ = i = 1 k γ i k B ( τ , x ^ i ( z ^ i , ( τ ) ) ) . ( 15 )

Figure 00000035
d η k ( τ ) d τ = i = one k γ i k B ( τ , x ^ i ( z ^ i , ( τ ) ) ) . ( fifteen )
Figure 00000035

Полученная каждая k - тая двухточечная краевая задача выглядит следующим образомThe resulting every kth two-point boundary value problem is as follows

z ˙ ^ k = η k ( τ ) , η ˙ k ( τ ) = i = 1 k γ i k B ( τ , x ^ i ) , x ^ ¨ 0 = f ( x ^ i , x ^ ˙ i , z ^ i ) μ 1 H x N 1 ( y H ( x ^ ( z ^ , τ ) , t ) ) T n k ( T ) = 0, z k ( 0 ) = z 0 + P 0 η k ( 0 ) , ( 16 )

Figure 00000036
z ˙ ^ k = η k ( τ ) , η ˙ k ( τ ) = i = one k γ i k B ( τ , x ^ i ) , x ^ ¨ 0 = f ( x ^ i , x ^ ˙ i , z ^ i ) - μ - one H x N - one ( y - H ( x ^ ( z ^ , τ ) , t ) ) T n k ( T ) = 0 z k ( 0 ) = z 0 + P 0 η k ( 0 ) , ( 16 )
Figure 00000036

где P0 - некоторая матрица размера m×m, i = 1, k ¯

Figure 00000037
.where P 0 is some m × m matrix, i = one, k ¯
Figure 00000037
.

Применение метода инвариантного погружения позволяет получить итерационный алгоритм идентификации, который выглядит следующим образомApplication of the invariant immersion method allows one to obtain an iterative identification algorithm, which is as follows

z ˙ ^ 0 = P 0 G 0 T H T x ^ N 1 γ 0 1 ( y H ( x ^ 0 ( z ^ 0 ) , t ) ) , P ˙ 0 = I P 0 G 0 T x ^ { H T x N 1 γ 0 1 ( y H ( x ^ 0 ( z ^ 0 ) , t ) ) } G 0 P 0 , G ¨ 0 = f x ˙ ^ G ˙ 0 + f x ^ G 0 + f x ^ , x ¨ ^ 0 = f ( x ^ i , x ˙ ^ i , z ^ i ) μ 1 H T x N 1 ( y H ( x ^ 0 ( z ^ 0 ) , t ) ) , z ˙ ^ k = P k G k T H T x ^ N 1 i = 1 k γ i k ( y H ( x ^ ( z ^ ) , t ) ) , P ˙ k = I P k G k T x ^ { H T x ^ N 1 i = 1 k γ i k ( y H ( x ^ i ( z ^ i ) , t ) ) } G k P k , G ¨ k = f x ˙ ^ G ˙ k + f x ^ G k + f z ^ , G k ( 0 ) = G ˙ k ( 0 ) = 0, x ¨ ^ 0 = f ( x ^ k , x ˙ ^ k , z ^ k ) μ 1 H T x ^ N 1 ( y H ( x ^ i ( z ^ i ) , t ) ) , ( 18 ) P i ( t 0 ) = P ( 0 ) , x ^ i ( 0 ) = x ( 0 ) , x ˙ ^ i ( 0 ) = x ˙ ( 0 ) . i = 1, k ¯

Figure 00000038
z ˙ ^ 0 = P 0 G 0 T H T x ^ N - one γ 0 one ( y - H ( x ^ 0 ( z ^ 0 ) , t ) ) , P ˙ 0 = I - P 0 G 0 T x ^ { H T x N - one γ 0 one ( y - H ( x ^ 0 ( z ^ 0 ) , t ) ) } G 0 P 0 , G ¨ 0 = f x ˙ ^ G ˙ 0 + f x ^ G 0 + f x ^ , x ¨ ^ 0 = f ( x ^ i , x ˙ ^ i , z ^ i ) - μ - one H T x N - one ( y - H ( x ^ 0 ( z ^ 0 ) , t ) ) , ... z ˙ ^ k = P k G k T H T x ^ N - one i = one k γ i k ( y - H ( x ^ ( z ^ ) , t ) ) , P ˙ k = I - P k G k T x ^ { H T x ^ N - one i = one k γ i k ( y - H ( x ^ i ( z ^ i ) , t ) ) } G k P k , G ¨ k = f x ˙ ^ G ˙ k + f x ^ G k + f z ^ , G k ( 0 ) = G ˙ k ( 0 ) = 0 x ¨ ^ 0 = f ( x ^ k , x ˙ ^ k , z ^ k ) - μ - one H T x ^ N - one ( y - H ( x ^ i ( z ^ i ) , t ) ) , ( eighteen ) P i ( t 0 ) = P ( 0 ) , x ^ i ( 0 ) = x ( 0 ) , x ˙ ^ i ( 0 ) = x ˙ ( 0 ) . i = one, k ¯
Figure 00000038

где I - единичная матрица.where I is the identity matrix.

Первая группа уравнений (18) представляет собой рекуррентные уравнения последовательной идентификации, а каждое последующее уравнение оценки для к+1 использует в качестве входных параметров y, z ^ 1 , , z ^ k

Figure 00000039
и x ^ 1 , , x ^ k
Figure 00000040
, что позволяет получить более точные значения идентифицируемых параметров.The first group of equations (18) are recurrent equations of sequential identification, and each subsequent estimation equation for k + 1 uses y as input parameters, z ^ one , ... , z ^ k
Figure 00000039
and x ^ one , ... , x ^ k
Figure 00000040
, which allows to obtain more accurate values of the identified parameters.

Устройство работает следующим образом. В исходном состоянии на вход устройства подается у, а в первом блоке 1.1 хранения констант записываются значения µ-1, x ^ , x ˙ ^ , z ^

Figure 00000041
, Δt, γ, k, N-1, P, G, Ġ, I. Значение µ-1 с первого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход первого блока 4.1 формирования произведения. Значение x ^ k
Figure 00000042
со второго информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
и на первый информационный вход второго блока 2.2 формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
. Значение x ˙ ^ k
Figure 00000043
с третьего информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
. Значение z ^ k
Figure 00000044
с четвертого информационного выхода первого блока 1.1 хранения констант поступает на третий информационный вход первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
и на второй информационный вход второго блока 2.2 формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
. Значение Δt с пятого информационного выхода блока 1.1 хранения констант поступает на второй информационный вход первого 8,1, второго 8.2, третьего 8.3 и шестого 8.6 блоков интегрирования, а также на первый информационный вход четвертого 8.4 и пятого 8.5 блоков интегрирования. Значение у с шестого информационного выхода блока 1.1 хранения констант поступает на второй информационный вход второго блока 5.2 формирования разности. Значение γ i k
Figure 00000045
с седьмого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход третьего блока 2.3 формирования функции (1-γi-1). Значение k с восьмого информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход третьего блока 2.3 формирования функции (1-γi-1) и на третий информационный вход первого блока 6.1 формирования суммы. Значение N-1 с девятого информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход третьего блока 4.3 формирования произведения. Значение P ˙ k
Figure 00000046
с десятого информационного выхода первого блока 1.1 хранения констант поступает на первый и второй информационные входы соответственно шестого блока 4.6 и девятого блока 4.9 формирования произведения. Значение Ġk c одиннадцатого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход второго блока 7.2 транспонирования, а также на первые информационные входы восьмого блока 4.8 и девятого блока 4.9 формирования произведения. Значение Ġk с двенадцатого информационногоThe device operates as follows. In the initial state, y is supplied to the input of the device, and in the first constant storage unit 1.1, the values μ -1 are written, x ^ , x ˙ ^ , z ^
Figure 00000041
, Δt, γ, k, N -1 , P, G, Ġ, I. The value of µ -1 from the first information output of the first constant storage unit 1.1 is supplied to the first information input of the first work formation block 4.1. Value x ^ k
Figure 00000042
from the second information output of the first constant storage unit 1.1 is fed to the first information input of the first function generating unit 2.1 f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
and to the first information input of the second block 2.2 function formation H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
. Value x ˙ ^ k
Figure 00000043
from the third information output of the first constant storage unit 1.1, it enters the second information input of the first function generating unit 2.1 f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
. Value z ^ k
Figure 00000044
from the fourth information output of the first block 1.1 of storage of constants goes to the third information input of the first block 2.1 of the formation of the function f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
and to the second information input of the second block 2.2 of the formation of the function H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
. The value Δt from the fifth information output of the constant storage unit 1.1 is supplied to the second information input of the first 8.1, second 8.2, third 8.3 and sixth 8.6 integration blocks, as well as the first information input of the fourth 8.4 and fifth 8.5 integration blocks. The value of y from the sixth information output of the constant storage unit 1.1 is supplied to the second information input of the second difference generating unit 5.2. Value γ i k
Figure 00000045
from the seventh information output of the first constant storage unit 1.1, it is supplied to the first information input of the third function generation unit 2.3 (1-γ i-1 ). The value of k from the eighth information output of the first constant storage unit 1.1 is supplied to the second information input of the third function generation unit 2.3 (1-γ i-1 ) and to the third information input of the first sum formation unit 6.1. The value of N -1 from the ninth information output of the first block 1.1 storage of constants is fed to the second information input of the third block 4.3 of the formation of the product. Value P ˙ k
Figure 00000046
from the tenth information output of the first constant storage unit 1.1, it is supplied to the first and second information inputs of the sixth block 4.6 and the ninth block 4.9 of product formation, respectively. The value Ġ k c of the eleventh information output of the first constant storage unit 1.1 is supplied to the first information input of the second transposition block 7.2, as well as to the first information inputs of the eighth block 4.8 and the ninth block 4.9 of the product formation. Ġ k value from the twelfth informational

выхода первого блока 1.1 хранения констант поступает на первый информационный вход седьмого блока 4.7 формирования произведения. Значение I с тринадцатого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход третьего блока 5.3 формирования разности. Значение f ( x ^ k , x ˙ ^ k , z ^ k )

Figure 00000047
с информационного выхода первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000048
поступает на второй информационный вход первого блока 5.1 формирования разности, а также на первые информационные входы второго блока 3.2, третьего блока 3.3 и четвертого блока 3.4 формирования производной. Значение H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
с информационного выхода второго блока 2.2 формирования функции H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
, поступает на первые информационные входы первого блока 7.1 транспонирования и второго блока 5.2 формирования разности, с информационного выхода которого значение y H ( x ^ i ( z ^ i ) , t )
Figure 00000049
поступает на второй информационный вход второго блока 4.2 формирования произведения, с информационного выхода которого значение γ i k ( y H ( x ^ i ( z i ^ ) , t ) )
Figure 00000050
поступает на первые информационные входы второго блока 1.2 хранения констант и первого блока 6.1 формирования суммы. Значение γ i k ( y H ( x i ^ ( z i ^ ) , t ) )
Figure 00000051
с информационного выхода второго блока 1.2 хранения констант поступает на второй информационный вход первого блока 6.1 формирования суммы, с информационного выхода которого значение i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) )
Figure 00000052
поступает на первый информационный вход третьего блока 4.3 формирования произведения, с информационного выхода которого значение N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) )
Figure 00000053
поступает на первый информационный вход четвертого блока 4.4 формирования произведения. Значение HT с информационного выхода первого блока 7.1 транспонирования поступает на первый информационный вход первого блока 3.1 формирования производной, с информационного выхода которого значение H T x ^
Figure 00000054
поступает на второй информационный вход четвертого блока 4.4 формирования произведения, с информационного выхода которого значение H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) )
Figure 00000055
поступает соответственно на второй информационный вход первого блока 4.1 формирования произведения и на первые информационные входы пятого блока 4.5 формирования произведения и пятого блока 3.5 формирования производной, с информационного выхода которого значение x ^ { H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) ) }
Figure 00000056
поступает на первый информационный вход десятого блока 4.10 формирования произведения. Значение G k T
Figure 00000057
с информационного выхода второго блока 7.2 транспонирования поступает на второй информационный вход шестого блока 4.6 формирования произведения, с информационного выхода которого значение P k G k T
Figure 00000058
поступает на вторые информационные входы пятого блока 4.5 и десятого блока 4.10 формирования произведения, с информационного выхода которого значение P k G k T x ^ { H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) ) }
Figure 00000059
поступает на второй информационный вход одиннадцатого блока 4.11 формирования произведения. Значение GkPk информационного выхода девятого блока 4.9 формирования произведения поступает на первый информационный вход одиннадцатого блока 4.11 формирования произведения, с информационного выхода которого значение P k G k T x ^ { H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) ) } G k P k
Figure 00000060
поступает на второй информационный вход третьего блока 5.3 формирования разности, с информационного выхода которого значение I P k G k T x ^ { H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) ) } G k P k
Figure 00000061
поступает на первый информационный вход шестого блока 8.6 интегрирования, с информационного выхода которого значение Pk поступает на первый и второй информационные входы соответственно шестого блока 4.6 и девятого блока 4.9 формирования произведения. Значение μ 1 H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) )
Figure 00000062
с информационного выхода первого блока 4.1 формирования произведения поступает на первый информационный вход первого блока 5.1 формирования разности, с информационного выхода которого значение f ( x ^ k , x ˙ ^ k , z ^ k ) μ 1 H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) )
Figure 00000063
поступает на первый информационный вход первого блока 8.1 интегрирования, с информационного выхода которого значение x ˙ ^ k
Figure 00000064
поступает на второй информационный вход первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
и первый информационный вход второго блока 8.2 интегрирования, с информационного выхода которого значение x ^ i
Figure 00000066
поступает на первые информационные входы первого блока 2.1 и второго блока 2.2 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
и H ( x i ^ ( z i ^ ) , t )
Figure 00000067
соответственно. Значение P k G k T H T x ^ N 1 i = 1 k γ i k ( y H ( x i ^ ( z i ^ ) , t ) )
Figure 00000068
с информационного выхода пятого блока 4.5 формирования умножения поступает на первый информационный вход третьего блока 8.3 интегрирования, с информационного выхода которого снимается значение искомого вектора идентификации параметров z и подается на второй и третий информационные входы соответственно второго блока 2.2 и первого блока 2.1 формирования функции H ( x i ^ ( z i ^ ) , t )
Figure 00000067
и f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
соответственно. Значение f x ˙ ^
Figure 00000069
с информационного выхода второго блока 3.2 формирования производной поступает на второй информационный вход седьмого блока 4.7 формирования произведения, с информационного выхода которого значение f x ˙ ^ G ˙ k
Figure 00000070
поступает на первый информационный вход второго блока 6.2 формирования суммы. Значение f x ^
Figure 00000071
с информационного выхода третьего блока 3.3 формирования производной поступает на второй информационный вход восьмого блока 4.8 формирования произведения, с информационного выхода которого значение f x ^ G k
Figure 00000072
поступает на второй информационный вход второго блока 6.2 формирования суммы, с информационного выхода которого значение f x ˙ ^ G ˙ k + f x ^ G k
Figure 00000073
поступает на первый информационный вход третьего блока 6.3 формирования суммы. Значение f x ^
Figure 00000071
с информационного выхода четвертого блока 3.4 формирования производной поступает на второй информационный вход третьего блока 6.3 формирования суммы, с информационного выхода которого значение G ¨ k
Figure 00000074
поступает на второй информационный вход пятого блока 8.5 интегрирования, с информационного выхода которого значение Ġk поступает на первый информационный вход седьмого блока 4.7 формирования произведения и второй информационный вход четвертого блока 8.4 интегрирования, с информационного выхода которого значение Gk поступает на первый информационный вход второго блока 7.2 транспонирования, а также на первые информационные входы восьмого блока 4.8 и девятого блока 4.9 формирования произведения.the output of the first block 1.1 storage of constants is fed to the first information input of the seventh block 4.7 of the formation of the product. The value I from the thirteenth information output of the first constant storage unit 1.1 is supplied to the first information input of the third difference generating unit 5.3. Value f ( x ^ k , x ˙ ^ k , z ^ k )
Figure 00000047
from the information output of the first block 2.1 function formation f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000048
arrives at the second information input of the first block 5.1 of the formation of the difference, as well as the first information inputs of the second block 3.2, the third block 3.3 and the fourth block 3.4 of the formation of the derivative. Value H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
from the information output of the second block 2.2 function formation H ( x ^ i , ( z ^ i ) , t )
Figure 00000003
arrives at the first information inputs of the first transposition block 7.1 and the second difference generating block 5.2, from the information output of which the value y - H ( x ^ i ( z ^ i ) , t )
Figure 00000049
arrives at the second information input of the second block 4.2 of the formation of the work, from the information output of which the value γ i k ( y - H ( x ^ i ( z i ^ ) , t ) )
Figure 00000050
arrives at the first information inputs of the second block 1.2 of storage of constants and the first block 6.1 of the formation of the sum. Value γ i k ( y - H ( x i ^ ( z i ^ ) , t ) )
Figure 00000051
from the information output of the second block 1.2 of storage of constants goes to the second information input of the first block 6.1 of the formation of the amount from the information output of which i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) )
Figure 00000052
arrives at the first information input of the third block 4.3 of the formation of the work, from the information output of which the value N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) )
Figure 00000053
arrives at the first information input of the fourth block 4.4 of the formation of the work. The value of H T from the information output of the first transposition block 7.1 is supplied to the first information input of the first derivative formation block 3.1, from the information output of which H T x ^
Figure 00000054
arrives at the second information input of the fourth block 4.4 of the formation of the work, from the information output of which the value H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) )
Figure 00000055
arrives, respectively, at the second information input of the first block 4.1 of product formation and the first information inputs of the fifth block 4.5 of product formation and the fifth block 3.5 of derivative formation, from the information output of which x ^ { H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) ) }
Figure 00000056
arrives at the first information input of the tenth block 4.10 formation of the work. Value G k T
Figure 00000057
from the information output of the second block 7.2 transposition is fed to the second information input of the sixth block 4.6 of the formation of the product, from the information output of which P k G k T
Figure 00000058
arrives at the second information inputs of the fifth block 4.5 and the tenth block 4.10 of the formation of the work, from the information output of which the value P k G k T x ^ { H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) ) }
Figure 00000059
arrives at the second information input of the eleventh block 4.11 of the formation of the work. The value G k P k of the information output of the ninth block 4.9 of the creation of the product is supplied to the first information input of the eleventh block 4.11 of the formation of the work, from the information output of which the value P k G k T x ^ { H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) ) } G k P k
Figure 00000060
arrives at the second information input of the third block 5.3 of the formation of the difference, from the information output of which the value I - P k G k T x ^ { H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) ) } G k P k
Figure 00000061
arrives at the first information input of the sixth integration unit 8.6, from the information output of which the value of P k goes to the first and second information inputs of the sixth block 4.6 and the ninth block 4.9 of product formation, respectively. Value μ - one H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) )
Figure 00000062
from the information output of the first block 4.1 of the formation of the product goes to the first information input of the first block 5.1 of the formation of the difference, from the information output of which the value f ( x ^ k , x ˙ ^ k , z ^ k ) μ - one H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) )
Figure 00000063
arrives at the first information input of the first integration unit 8.1, from the information output of which the value x ˙ ^ k
Figure 00000064
arrives at the second information input of the first block 2.1 function formation f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
and the first information input of the second integration unit 8.2, from the information output of which the value x ^ i
Figure 00000066
arrives at the first information inputs of the first block 2.1 and the second block 2.2 of the formation of the function f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
and H ( x i ^ ( z i ^ ) , t )
Figure 00000067
respectively. Value P k G k T H T x ^ N - one i = one k γ i k ( y - H ( x i ^ ( z i ^ ) , t ) )
Figure 00000068
from the information output of the fifth block 4.5 of multiplication formation, it goes to the first information input of the third integration block 8.3, from the information output of which the value of the sought parameter identification vector z is taken and fed to the second and third information inputs of the second block 2.2 and the first block 2.1 of the function formation H ( x i ^ ( z i ^ ) , t )
Figure 00000067
and f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
respectively. Value f x ˙ ^
Figure 00000069
from the information output of the second block 3.2 formation of the derivative goes to the second information input of the seventh block 4.7 of the formation of the product, from the information output of which f x ˙ ^ G ˙ k
Figure 00000070
arrives at the first information input of the second block 6.2 of the formation of the amount. Value f x ^
Figure 00000071
from the information output of the third block 3.3 of the formation of the derivative is fed to the second information input of the eighth block 4.8 of the formation of the product, from the information output of which the value f x ^ G k
Figure 00000072
arrives at the second information input of the second block 6.2 of the formation of the amount, from the information output of which the value f x ˙ ^ G ˙ k + f x ^ G k
Figure 00000073
arrives at the first information input of the third block 6.3 of the formation of the amount. Value f x ^
Figure 00000071
from the information output of the fourth block 3.4 of the formation of the derivative is fed to the second information input of the third block 6.3 of the formation of the sum, from the information output of which the value G ¨ k
Figure 00000074
arrives at the second information input of the fifth integration block 8.5, from the information output of which the value Ġ k goes to the first information input of the seventh product formation block 4.7 and the second information input of the fourth integration block 8.4, from the information output of which the value G k goes to the first information input of the second block 7.2 transposition, as well as the first information inputs of the eighth block 4.8 and the ninth block 4.9 of the formation of the work.

Первый блок интегрирования 8.1 (фиг.2) работает следующим образом. С информационного выхода первого блока 5.1 формирования разности значение x ¨ ^

Figure 00000075
поступает на первый информационный вход блока 8.1.1 формирования произведения. Значение Δt с пятого информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход первого 8.1 блока интегрирования, с информационного выхода которого значение x ¨ Δ t ^
Figure 00000076
поступает на первый информационный вход блока 8.1.2 формирования суммы. Значение x ˙ ^ i 1
Figure 00000077
информационного выхода линии задержки 8.1.3 поступает на второй информационный вход блока 8.1.2 формирования суммы, с информационного выхода которого снимается значение первообразной x ˙ ^ i
Figure 00000078
и подается на первый вход линии задержки 8.1.3 и на второй информационный вход первого блока 2.1 формирования функции f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
и первый информационный вход второго блока 8.2 интегрирования. Остальные блоки интегрирования работают аналогично.The first integration unit 8.1 (figure 2) works as follows. From the information output of the first block 5.1 the formation of the difference value x ¨ ^
Figure 00000075
arrives at the first information input of the block 8.1.1 formation of the work. The value Δt from the fifth information output of the first constant storage unit 1.1 is supplied to the second information input of the first 8.1 integration unit, from the information output of which x ¨ Δ t ^
Figure 00000076
arrives at the first information input of the block 8.1.2 forming the amount. Value x ˙ ^ i - one
Figure 00000077
the information output of the delay line 8.1.3 goes to the second information input of the block 8.1.2 of the formation of the sum, from the information output of which the value of the antiderivative x ˙ ^ i
Figure 00000078
and is fed to the first input of the delay line 8.1.3 and to the second information input of the first block 2.1 of the formation of the function f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000065
and the first information input of the second integration unit 8.2. The remaining integration units work in a similar way.

ИСТОЧНИКИ ИНФОРМАЦИИINFORMATION SOURCES

1. http://wwwl.fips.ru/fips_servl/fips_servlet1.http: //wwwl.fips.ru/fips_servl/fips_servlet

2. Колесников А.А. Синергетическая теория управления. Таганрог: ТРТУ. -М.: Энергоатомиздат, 1994, с.117-124.2. Kolesnikov A.A. Synergetic control theory. Taganrog: TRTU. -M .: Energoatomizdat, 1994, p. 117-124.

3. Лурье А.И. Аналитическая механика. -М.: Гос.изд. физ. - мат.лит., 1961, с.642-648.3. Lurie A.I. Analytical mechanics. -M .: State publ. physical - Mat. lit., 1961, p. 642-648.

4. Сейдж Э.П., Меле Д.Л. Идентификация систем управления. -М.: Наука, 1974. с.204-212.4. Sage E.P., Mele D.L. Identification of management systems. -M .: Nauka, 1974.p.204-212.

5. Справочник по теории автоматического управления. / Под ред. А.А. Красовского. М.: Наука, 1987, с.347-360.5. Reference on the theory of automatic control. / Ed. A.A. Krasovsky. M .: Nauka, 1987, p. 347-360.

6. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука, 1986, с.15-18.6. Tikhonov A.N., Arsenin V.Ya. Methods for solving incorrect tasks. M .: Nauka, 1986, p. 15-18.

7. Фарина А., Студер Ф. Цифровая обработка радиолокационной информации. Сопровождение целей. - М.: Радио и связь, 1993, с.38-47.7. Farina A., Studer F. Digital processing of radar information. Tracking goals. - M.: Radio and Communications, 1993, p. 38-47.

Claims (1)

Устройство идентификации Лагранжевых динамических систем на основе итерационной регуляризации, содержащее первый блок хранения констант; первый, второй блоки формирования функций f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
, H ( x ^ i , ( z ^ i ) , t )
Figure 00000002
соответственно; первый, второй, третий, четвертый, пятый блоки формирования производной; первый, второй, третий, четвертый, пятый, шестой, седьмой, восьмой, девятый, десятый, одиннадцатый блоки формирования произведения; первый, второй, третий блоки формирования разности; первый, второй блоки формирования суммы; первый, второй блоки транспонирования; первый, второй, третий, четвертый, пятый, шестой блоки интегрирования; при этом первый информационный вход первого блока хранения констант является входом устройства, а информационный выход третьего блока интегрирования является выходом устройства, отличающееся тем, что в устройство введены второй блок хранения констант, третий блок формирования функции (1-γi-1), третий блок формирования суммы, причем информационный выход второго блока формирования произведения соединен с информационным входом второго блока хранения констант, информационный выход которого соединен со вторым входом первого блока формирования суммы; седьмой информационный выход первого блока хранения констант соединен с первым информационным входом третьего блока формирования функции (1-γi-1); восьмой информационный выход первого блока хранения констант соединен со вторым информационным входом третьего блока формирования функции (1-γi-1), выход которого соединен со вторым информационным входом второго блока формирования произведения; информационный выход второго блока формирования суммы соединен с первым информационным входом третьего блока формирования суммы; информационный выход четвертого блока формирования производной соединен со вторым информационным входом третьего блока формирования суммы, информационный выход которого соединен со вторым информационным входом пятого блока интегрирования.
An identification device for Lagrangian dynamical systems based on iterative regularization, comprising a first constant storage unit; first, second blocks of the formation of functions f ( x ^ i , x ˙ ^ i , z ^ i )
Figure 00000001
, H ( x ^ i , ( z ^ i ) , t )
Figure 00000002
respectively; first, second, third, fourth, fifth blocks of derivative formation; first, second, third, fourth, fifth, sixth, seventh, eighth, ninth, tenth, eleventh blocks of formation of a work; first, second, third blocks of difference formation; first, second blocks of the formation of the amount; first, second transpose blocks; first, second, third, fourth, fifth, sixth integration blocks; the first information input of the first constant storage unit is the input of the device, and the information output of the third integration unit is the output of the device, characterized in that the second constant storage unit, the third function formation unit (1-γ i-1 ), the third unit are introduced into the device the formation of the sum, and the information output of the second unit of formation of the product is connected to the information input of the second unit of storage of constants, the information output of which is connected to the second input of the first unit of formation amount of money; the seventh information output of the first constant storage unit is connected to the first information input of the third function generating unit (1-γ i-1 ); the eighth information output of the first constant storage unit is connected to the second information input of the third function formation unit (1-γ i-1 ), the output of which is connected to the second information input of the second product formation unit; the information output of the second sum forming unit is connected to the first information input of the third sum forming unit; the information output of the fourth derivative formation unit is connected to the second information input of the third sum formation unit, the information output of which is connected to the second information input of the fifth integration unit.
RU2013106875/08A 2013-02-15 2013-02-15 Device for identifying lagrange dynamic systems based on iterative regularisation RU2528133C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2013106875/08A RU2528133C1 (en) 2013-02-15 2013-02-15 Device for identifying lagrange dynamic systems based on iterative regularisation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2013106875/08A RU2528133C1 (en) 2013-02-15 2013-02-15 Device for identifying lagrange dynamic systems based on iterative regularisation

Publications (2)

Publication Number Publication Date
RU2013106875A RU2013106875A (en) 2014-08-20
RU2528133C1 true RU2528133C1 (en) 2014-09-10

Family

ID=51384445

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2013106875/08A RU2528133C1 (en) 2013-02-15 2013-02-15 Device for identifying lagrange dynamic systems based on iterative regularisation

Country Status (1)

Country Link
RU (1) RU2528133C1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2653939C1 (en) * 2016-12-27 2018-05-15 ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ КАЗЕННОЕ ВОЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ "Военная академия Ракетных войск стратегического назначения имени Петра Великого" МИНИСТЕРСТВА ОБОРОНЫ РОССИЙСКОЙ ФЕДЕРАЦИИ Device for assessment of state and identification of parameters of the dynamic system models
RU2657477C1 (en) * 2016-11-18 2018-06-14 Дмитрий Сергеевич Андрашитов Device for estimating parameters using a priori information in form of the integral of action

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2357357C2 (en) * 2007-01-31 2009-05-27 Андрей Александрович Костоглотов Digital intellectual recursive filter
RU2452080C1 (en) * 2011-06-08 2012-05-27 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Южно-Российский государственный университет экономики и сервиса" (ФГБОУ ВПО "ЮРГУЭС") Digital multi-iterative filter
RU2464615C1 (en) * 2011-07-21 2012-10-20 Андрей Александрович Костоглотов Device for identifying parameters of dynamic systems based on variation principles

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2357357C2 (en) * 2007-01-31 2009-05-27 Андрей Александрович Костоглотов Digital intellectual recursive filter
RU2452080C1 (en) * 2011-06-08 2012-05-27 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Южно-Российский государственный университет экономики и сервиса" (ФГБОУ ВПО "ЮРГУЭС") Digital multi-iterative filter
RU2464615C1 (en) * 2011-07-21 2012-10-20 Андрей Александрович Костоглотов Device for identifying parameters of dynamic systems based on variation principles

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2657477C1 (en) * 2016-11-18 2018-06-14 Дмитрий Сергеевич Андрашитов Device for estimating parameters using a priori information in form of the integral of action
RU2653939C1 (en) * 2016-12-27 2018-05-15 ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ КАЗЕННОЕ ВОЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ "Военная академия Ракетных войск стратегического назначения имени Петра Великого" МИНИСТЕРСТВА ОБОРОНЫ РОССИЙСКОЙ ФЕДЕРАЦИИ Device for assessment of state and identification of parameters of the dynamic system models

Also Published As

Publication number Publication date
RU2013106875A (en) 2014-08-20

Similar Documents

Publication Publication Date Title
Sharma et al. Trend analysis and change point techniques: a survey
Kuznetsov et al. Generalization bounds for time series prediction with non-stationary processes
US20170017732A1 (en) Parameter-dependent model-blending with multi-expert based machine learning and proxy sites
Pascoal et al. Non-parametric wave spectral estimation using vessel motions
Hekmati et al. About the convergence of POD and EPOD modes computed from CFD simulation
US20200202176A1 (en) Method for correlation filter based visual tracking
Meng et al. Estimation of chirp signals with time-varying amplitudes
Harley et al. Application of Mellin transform features for robust ultrasonic guided wave structural health monitoring
RU2528133C1 (en) Device for identifying lagrange dynamic systems based on iterative regularisation
Alford et al. Estimation of extreme slamming pressures using the non-uniform Fourier phase distributions of a design loads generator
De Brouwer et al. Latent convergent cross mapping
Fischer et al. Beamforming of aeroacoustic sources in the time domain: An investigation of the intermittency of the noise radiated by a forward-facing step
Grebenyuk et al. Digital image correlation with fast Fourier transform for large displacement measurement
Xiong et al. Multifractal spectrum distribution based on detrending moving average
Barnard et al. Differences between the CME fronts tracked by an expert, an automated algorithm, and the Solar Stormwatch project
Breakey et al. Comparison of metrics for the evaluation of similarity in acoustic pressure signals
Cho Search and parameter estimate in gravitational wave data analysis and the fisher matrix
Roncoroni et al. Efficient extraction of seismic reflection with Deep Learning
RU2653939C1 (en) Device for assessment of state and identification of parameters of the dynamic system models
O'Leary et al. Polynomial approximation: An alternative to windowing in Fourier analysis
Brault et al. Fast and robust pattern detection: Application to spherical bead location in holographic microscopy
Toporkov Study of electromagnetic scattering from randomly rough ocean-like surfaces using integral-equation-based numerical technique
Muthyala et al. Quantum Interference for Counting Clusters
Bărbulescu et al. ARIMA models for the analysis of the precipitation evolution
Hasanov et al. Fourier Collocation Algorithm for identification of a spacewise dependent source in wave equation from Neumann-type measured data

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20160216