Устройство идентификации лагранжевых динамических систем на основе итерационной регуляризации

Изобретение относится к области цифровой вычислительной техники и может быть использовано в автоматических и автоматизированных системах различного назначения для идентификации параметров. Техническим результатом является повышение точности идентификации параметров динамических систем. Устройство содержит блоки хранения констант, блоки формирования функций, блоки формирования производной, блоки формирования произведения, блоки формирования разности, блоки формирования суммы, блоки транспонирования, блоки интегрирования. Указанный технический результат достигается за счет вычисления и последовательного уточнения значений параметров регуляризации. 2 ил.

 

Изобретение относится к области цифровой вычислительной техники и может быть использовано в автоматических и автоматизированных системах различного назначения для идентификации параметров.

Известно устройство идентификации параметров динамических систем, в основе которого лежит фильтр Калмана. Его эффективное функционирование возможно в случае априорной определенности законов распределения внешних воздействий, что является существенным недостатком.

Наиболее близким по технической сущности к заявленному изобретению является устройство идентификации параметров динамических систем на основе вариационных принципов [1]. Однако его использование в условиях неизвестных воздействий не позволяет получить требуемую точность.

Цель изобретения - повысить точность идентификации параметров динамических систем. Указанный технический результат достигается за счет устройства идентификации Лагранжевых динамических систем на основе итерационной регуляризации, которое содержит следующие блоки: первый, второй блоки хранения констант; первый, второй, третий блоки формирования функций $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ , $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ , (1-γ_i-1), соответственно; первый, второй, третий, четвертый, пятый блоки формирования производной; первый, второй, третий, четвертый, пятый, шестой, седьмой, восьмой, девятый, десятый, одиннадцатый блоки формирования произведения; первый, второй, третий блоки формирования разности; первый, второй, третий блоки формирования суммы; первый, второй блоки транспонирования; первый, второй, третий, четвертый, пятый, шестой блоки интегрирования; при этом на первый информационный вход первого блока хранения констант, который является входом устройства, поступает значение наблюдаемой величины; первый информационный выход первого блока хранения констант соединен с первым информационным входом первого блока формирования произведения; второй, третий и четвертый информационные выходы первого блока хранения констант соединены соответственно с первым, вторым и третьим информационными входами первого блока формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ ; второй и четвертый информационные выходы первого блока хранения констант соединены соответственно с первым и вторым информационными входами второго блока формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ ; пятый информационный выход первого блока хранения констант соединен со вторыми информационными входами первого, второго, третьего и шестого блоков интегрирования, а также с первыми информационными входами четвертого и пятого блоков интегрирования; шестой информационный выход первого блока хранения констант соединен со вторым информационным входом второго блока формирования разности; седьмой информационный выход первого блока хранения констант соединен с первым информационным входом третьего блока формирования функции (1-γ_i-1); восьмой информационный выход первого блока хранения констант соединен со вторым информационным входом третьего блока формирования функции (1-γ_i-1) и третьим информационным входом первого блока формирования суммы; девятый информационный выход первого блока хранения констант соединен со вторым информационным входом третьего блока формирования произведения; десятый информационный выход первого блока хранения констант соединен с первым и вторым информационными входами соответственно шестого и девятого блоков формирования произведения; одиннадцатый информационный выход первого блока хранения констант соединен с первым информационным входом второго блока транспонирования, а также с первыми информационными входами восьмого и девятого блоков формирования произведения; двенадцатый информационный выход первого блока хранения констант соединен с первым информационным входом седьмого блока формирования произведения; тринадцатый информационный выход первого блока хранения констант соединен с первым информационным входом третьего блока формирования разности; информационный выход первого блока формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ соединен со вторым информационным

входом первого блока формирования разности, а также с первыми информационными входами второго, третьего и четвертого блоков формирования производной; информационный выход второго блока формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ соединен с первыми информационными входами первого блока транспонирования и второго блока формирования разности, информационный выход которого соединен с первым информационным входом второго блока формирования произведения, 1 информационный выход которого соединен с первыми информационными входами второго блока хранения констант и первого блока формирования суммы; информационный выход второго блока хранения констант соединен со вторым информационным входом первого блока формирования суммы, информационный выход которого соединен с первым информационным входом третьего блока формирования произведения, информационный выход которого соединен с первым информационным входом четвертого блока формирования произведения; информационный выход первого блока транспонирования соединен с первым информационным входом первого блока формирования производной, информационный выход которого соединен со вторым информационным входом четвертого блока формирования произведения, информационный выход которого соединен соответственно со вторым информационным входом первого блока формирования произведения и первыми информационными входами пятого блока формирования произведения и пятого блока формирования производной, информационный выход которого соединен с первым информационным входом десятого блока формирования произведения; информационный выход второго блока транспонирования соединен со вторым информационным входом шестого блока формирования произведения, информационный выход которого соединен соответственно со вторыми информационными входами пятого и десятого блоков формирования произведения, информационный выход которого соединен со вторым информационным входом одиннадцатого блока формирования произведения; информационный выход девятого блока формирования произведения соединен с первым информационным входом одиннадцатого блока формирования произведения, информационный выход которого соединен со вторым информационным входом третьего блока формирования разности, информационный выход которого соединен с первым информационным входом шестого блока интегрирования, информационный выход которого соединен с первым и вторым информационными входами соответственно шестого и девятого блоков формирования произведения; информационный выход первого блока формирования произведения соединен с первым информационным входом первого блока формирования разности, информационный выход которого соединен с первым информационным входом первого блока интегрирования, информационный выход которого соединен со вторым информационным входом первого блока формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и первым информационным входом второго блока интегрирования, информационный выход которого соединен с первыми информационными входами первого и второго блоков формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ ; информационный выход пятого блока формирования произведения соединен с первым информационным входом третьего блока интегрирования, информационный выход которого является выходом устройства, с него же снимается значение искомого вектора идентификации параметров и подается на второй и третий информационные входы соответственно второго и первого блоков формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ , и $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ ; информационный выход второго блока формирования производной соединен со вторым информационным входом седьмого блока формирования произведения, информационный выход которого соединен с первым информационным входом второго блока формирования суммы; информационный выход третьего блока формирования производной соединен со вторым информационным входом восьмого блока формирования произведения, информационный выход которого соединен со вторым информационным входом второго блока формирования суммы, информационный выход которого соединен с первым информационным входом третьего блока формирования суммы; информационный выход четвертого блока формирования производной соединен со вторым информационным входом третьего блока формирования суммы, информационный выход которого соединен со вторым информационным входом пятого блока интегрирования, информационный выход которого соединен с первым информационным входом седьмого блока формирования произведения и вторым информационным входом четвертого блока интегрирования, информационный выход которого соединен с первым информационным входом второго блока транспонирования, а также с первыми информационными входами восьмого и девятого блоков формирования умножения.

Сущность изобретения поясняется чертежами, где на фиг.1 представлено устройство идентификации параметров Лагранжевых динамических систем на основе итерационной регуляризации, которое содержит:

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

1.2 - второй блок хранения констант;

2.1- первый блок формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ ;

2.2 - второй блок формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ ;

2.3 - третий блок формирования функции (1-γ_i1);

3.1 - первый блок формирования производной;

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

3.3 - третий блок формирования производной;

3.4 - четвертый блок формирования производной;

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

4.1 - первый блок формирования произведения;

4.2 - второй блок формирования произведения;

4.3 - третий блок формирования произведения;

4.4 - четвертый блок формирования произведения;

4.5 - пятый блок формирования произведения;

4.6 - шестой блок формирования произведения;

4.7 - седьмой блок формирования произведения;

4.8 - восьмой блок формирования произведения;

4.9 - девятый блок формирования произведения;

4.10 - десятый блок формирования произведения;

4.11 - одиннадцатый блок формирования произведения;

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

5.2 - второй блок формирования разности;

5.3 - третий блок формирования разности;

6.1 - первый блок формирования суммы;

6.2 - второй блок формирования суммы;

6.3 - третий блок формирования суммы;

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

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

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

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

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

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

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

8.6 - шестой блок интегрирования.

На фиг.2 представлен первый блок интегрирования:

8.1.1 - блок формирования произведения;

8.1.2 - блок формирования суммы;

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

Устройство идентификации Лагранжевых динамических систем на основе итерационной регуляризации работает следующим образом: на первый информационный вход первого блока 1.1 хранения констант (фиг.1) поступает значение наблюдаемой величины; первый информационный выход первого блока 1.1 хранения констант соединен с первым информационным входом первого блока 4.1 формирования произведения; второй, третий и четвертый информационные выходы первого блока 1.1 хранения констант соединены соответственно с первым, вторым и третьим информационными входами первого блока 2.1 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ ; второй и четвертый информационные выходы первого блока 1.1 хранения констант соединены соответственно с первым и вторым информационными входами второго блока 2.2 формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ ; пятый информационный выход первого блока 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\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ соединен со вторым информационным входом первого блока 5.1 формирования разности, а также с первыми информационными входами второго блока 3.2, третьего блока 3.3 и четвертого блока 3.4 формирования производной; информационный выход второго блока 2.2 формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ соединен с первыми информационными входами первого блока 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\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ , соответственно; информационный выход пятого блока 4.5 формирования произведения соединен с первым информационным входом третьего блока 8.3 интегрирования, с информационного выхода которого снимается значение искомого вектора идентификации параметров z и подается на второй и третий информационные входы соответственно второго блока 2.2 и первого блока 2.1 формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ и $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ соответственно; информационный выход второго блока 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 формирования произведения.

Информационный выход первого блока 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\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и первый

информационный вход второго блока 8.2 интегрирования.

Блоки 8.2, 8.3, 8.4, 8.5, 8.6 интегрирования (фиг.1) имеют структуру и связи, аналогичные блоку 8.1, рассмотренному выше.

Заявленное изобретение направлено на повышение эффективности идентификации параметров динамических систем, что весьма важно на всех этапах создания, экспериментальной обработки и эксплуатации объектов ракетно-космической, авиационной, корабельной и других видов техники.

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

Структура предлагаемого устройства определяется решением задачи, чья постановка приведена ниже.

Принято, что динамическая система описывается следующим векторным уравнением дифференциальным уравнением

$$\ddot{x}=f\left(x,\dot{x},z\right),x\left(0\right)=x^{\circ },\dot{x}\left(0\right)=x^{\circ },(1)$$

где x∈Rⁿ - вектор обобщенных координат,

$$x\in R^n$$ - вектор обобщенных скоростей,

z∈R^m - вектор неизвестных постоянных параметров, подлежащих идентификации,

f - вектор-функция, непрерывная и дифференцируемая по совокупности аргументов,

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

T∈[0,T].

Уравнение наблюдения имеет вид

$$y=H\left(x,t\right)+n\left(t\right),(2)$$

где y∈R^k - вектор наблюдения,

H(x, f) - непрерывная вместе с частными производными вектор-функция,

k - натуральное число,

n(t) - вектор белого гауссовского шума.

Пусть динамика идентифицируемых параметров z определяется уравнением

$$\dot{z}=\eta ,z\left(0\right)=z^{\circ },(3)$$

где η∈R^m - вектор неизвестных неслучайных возмущений, удовлетворяющий требованиям физической реализуемости $$\eta \left(t\right)\in L_2^m\left[\mathrm{0,}T\right]$$ .

Пусть динамика исследуемой системы определяется принципом Гамильтона-Остроградского. В соответствии с ним

$$\delta \text{'}W=\delta S+{\displaystyle \underset{0}{\overset{T}{\int }}\delta \text{'}Adt={\displaystyle \underset{0}{\overset{T}{\int }}\left[\delta L+Q\delta x\right]}}dt=\mathrm{0,}(4)$$

где S - действие по Гамильтону на интервале времени [0, T],

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

А - работа вектора обобщенных внешних сил Q∈Rⁿ,

δx - вектор вариаций обобщенных координат, знак δ' обозначает бесконечно малую величину, которая не является вариацией.

Ставится задача определения оценки $$\hat{z}$$ вектора z из условия минимума функционала невязки

$$J_1=\frac{1}{2}{\displaystyle \underset{0}{\overset{T}{\int }}{\left[y-H\left(\hat{x}\left(\hat{z}\right),t\right)\right]}^T{N}^{-1}\left[y-H\left(\hat{x}\left(\hat{z}\right),t\right)\right]dt\to \min ,(5)}$$

где N∈R^m×R^m - матрица односторонней спектральной плотности шума наблюдения, характеризующая интенсивность помех в канале наблюдений.

Задача идентификации параметров динамических систем (1)-(5) является обратной, некорректно поставленной по Адамару. Для ее решения используется метод регуляризации А.Н. Тихонова.

Это требует рассмотрения условий минимума сглаживающего функционала

$$J^{\alpha }\left[z,x\right]=J_1+\mu W+\alpha \Omega \left[\eta \right],(6)$$

где $$\Omega \left[\eta \right]=\frac{1}{2}{\displaystyle \underset{0}{\overset{T}{\int }}{\eta }^T\left(t\right)\eta \left(t\right)dt}$$ - стабилизирующий функционал,

α - положительное число.

Функционал J^α[z, х] неотрицательный, поэтому существует его нижняя грань, и если {α_n} - убывающая последовательность положительных чисел, сходящихся к нулю (α→0), и соответствующая последовательность {η_n(t)}, согласно, сходится к η*(t) при условии, что числовой параметр α удовлетворяет требованиям

$$\underset{n\to \infty }{\lim }{\alpha }_n=\mathrm{0,}\underset{n\to \infty }{\lim }{\displaystyle \sum _{k=1}^n{\alpha }_k}=\infty ,\underset{n\to \infty }{\lim }{\displaystyle \sum _{k=1}^n{\alpha }_k}<\infty .$$

Для дифференцируемого выпуклого функционала J^α[z, x] минимум находится путем определения стационарной точки, в которой

$$gra{d}_z{J}^{\alpha }\left[\hat{z},\hat{x}\right]=gran{d}_x{J}^{\alpha }\left[\hat{z},\hat{x}\right]=\mathrm{0,}(7)$$

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

Для определения градиента в (7) воспользуемся игольчатым варьированием возмущения η(t) и асинхронным варьированием траектории x(t).

Игольчатой вариацией возмущения будем называть функцию вида

$${\eta }_{\epsilon }\left(t\right)=\{\begin{array}{l}v,t\in \left[\tau ,\tau +\epsilon l\right],\\ v*\left(t\right),t\notin \left[\tau ,\tau +\epsilon l\right].\end{array}$$

где τ - заданная произвольная точка непрерывности функции η*(t),

v - некоторый постоянный вектор,

l -заданное положительное число,

ε - произвольное положительное число, такое, что τ+εl<T.

Разность η_ε(t)-η*(t)=δη(t) будем называть игольчатой вариацией возмущения, которая порождает вариацию параметров δz(t), вызывающую вариацию траектории δx(t). Обозначим x_ε(t) - траекторию, соответствующую η_ε(t), x*(t) - соответствующую Вариацию траектории вида δx_c(t)=x_ε(t)-x*(t) будем называть синхронной, а вариацию

$$\delta x_a\left(t\right)=x_{\epsilon }\left(t+dt\right)-x*\left(t\right)=\frac{\partial x}{\partial z}\delta z=G\delta z$$

асинхронной согласно терминологии.

Здесь G - матрица чувствительности системы (1) по вектору параметров z, удовлетворяющая уравнению

$$\ddot{G}=\frac{\partial f}{\partial \dot{x}}\dot{G}+\frac{\partial f}{\partial x}G+\frac{\partial f}{\partial x},G\left(0\right)=\dot{G}\left(0\right)=0$$

Вариации δx_c(t) и δx_a(t) обладают следующими свойствами

δx_c(t)=0, t∈[0,τ+εl),

δx_c(t)=0, t∉[0,τ+εl),

δx_c(t)=δx_a(t), t+τ+εl.

Пусть η*(t) - возмущение, доставляющее минимум функционалу (6). Найдем приращение функционала (6), обусловленное игольчатой вариацией η_ε(t). Заметим, что на интервале [τ,τ+εl] за счет скачка η_ε(t) возникает асинхронная вариация δx_a(t), а на интервале [τ+εl,T] возникает синхронная δx_c(t), поскольку приращения обобщенных координат

определяются решением дифференциальных уравнений в вариациях при начальных условиях в момент времени t=τ+εl. Таким образом, имеем

$$\begin{array}{cc}\begin{array}{l}\Delta J={\displaystyle \underset{\tau }{\overset{\tau +\epsilon l}{\int }}\left\{\left[-{\left(y-H\left(x,t\right)\right)}^T{N}^{-1}\frac{\partial H}{\partial x}\right]\delta x_a+\alpha {\eta }^T\delta \eta +\mu \delta L_a=Q\delta x_a\right\}dt+}\\ {\displaystyle \underset{\tau +\epsilon l}{\overset{T}{\int }}\left\{\left[-{\left(y-H\left(x,t\right)\right)}^T{N}^{-1}\frac{\partial H}{\partial x}\right]\delta x_c+\mu \delta L_c+Q\delta x_c\right\}dt}\end{array}& (8)\end{array}$$

где δL_a и δL_c - вариации кинетического потенциала, связанные соответственное δx_a(t) и δx_c(Y)-

Учитывая, что $$\dot{z}=\eta$$ , выполняя преобразования с учетом

$${\displaystyle \underset{\tau }{\overset{\tau +\epsilon l}{\int }}{\eta }^T\delta \dot{z}dt={\eta }^T\delta z|{}_{\tau }^{\tau +\epsilon l}-}{\displaystyle \underset{\tau }{\overset{\tau +\epsilon l}{\int }}{\dot{\eta }}^T\delta zdt=-{\displaystyle \underset{\tau }{\overset{\tau +\epsilon l}{\int }}{\dot{\eta }}^T\delta zdt}},(9)$$

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

$$\begin{array}{cc}\begin{array}{l}\Delta J={\displaystyle \underset{\tau }{\overset{\tau +\epsilon l}{\int }}\left[-{\left(y-H\left(x,t\right)\right)}^T{N}^{-1}\frac{\partial H}{\partial x}G+\alpha {\dot{\eta }}^T\right]\delta zdt+}\\ {\displaystyle \underset{\tau +\epsilon l}{\overset{T}{\int }}-{\left(y-H\left(x,t\right)\right)}^T{N}^{-1}\frac{\partial H}{\partial x}-\mu \left[\ddot{x}-f\left(x,\dot{x},z\right)\right]\delta x_{c}dt}\end{array}& (10)\end{array}$$

откуда следует, что в точке минимума функционала

$$\begin{array}{cc}\begin{array}{l}gra{d}_z{J}^{\alpha }\left[\hat{z},\hat{x}\right]=-{\left(y-H\left(\hat{x},\tau \right)\right)}^T{N}^{-1}\frac{\partial H}{\partial x}G+\alpha \dot{\eta }{*}^T=0\\ gra{d}_z{J}^{\alpha }\left[\hat{z},\hat{x}\right]=-{\left(y-H\left(\hat{x},t\right)\right)}^T{N}^{-1}\frac{\partial H}{\partial x}-\mu \left[\hat{\ddot{x}}-f\left(x,\hat{\dot{x}},z\right)\right]=0\end{array}& (11)\end{array}$$

Используя эти условия стационарности получим

$$\begin{array}{cc}\begin{array}{l}\eta *\left(\tau \right)={\alpha }^{-1}{\displaystyle \underset{\tau }{\overset{T}{\int }}{N}^{-1}{G}^T\frac{\partial H}{\partial x}\left[y-H\left(\hat{x}\left(\hat{z},\tau \right),t\right)\right]dt,}\\ x=f\left(\hat{x},\hat{\dot{x}},z\right)-{\mu }^{-1}\frac{\partial H}{\partial x}{N}^{-1}{\left(y-H\left(\hat{x}\left(\hat{z},\tau \right)t\right)\right)}^T.\end{array}& (12)\end{array}$$

Найти точку $$\left[\hat{z},\eta *\right]$$ путем решения данных уравнений в условиях некорректности исходной задачи достаточно сложно, поэтому широкое распространение получили методы последовательного приближения к стационарной точке. Рассмотрим метод простой итерации

$${\eta }_{k+1}={\eta }_k-{\alpha }_{k}gran{d}_z{J}^{\alpha }\left[{\hat{z}}_{{}_{k+1}},{\eta }_k\right].(13)$$

Для сокращения записи вводится обозначение

$$B\left(\hat{z},\tau \right)=N^{-1}{G}^T\frac{\partial H^T}{\partial x}\left[y-H\left(\hat{\dot{x}}\left(\hat{z},\tau \right),t\right)\right].$$

Тогда с учетом (11) алгоритм (13) может быть представлен в виде

$${\eta }_{k+1}\left(\tau \right)={\eta }_k\left(\tau \right)\left[1-{\alpha }_k\right]={\displaystyle \underset{\tau }{\overset{T}{\int }}B}\left(z_{k+1}\left(\tau \right),t\right)dt,$$

Приняв за нулевое приближение η₀=0, итерационную последовательность в развернутой форме согласно запишем следующим образом

$$\begin{array}{l}{\eta }_0\left(\tau \right)=\mathrm{0,}\stackrel{\dot{^}}{z}={\eta }_0\left(\tau \right),\\ {\stackrel{\ddot{^}}{x}}_0=f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)-{\mu }^{-1}\frac{\partial H}{\partial x}{N}^{-1}{\left(y-H\left(\hat{x}\left(\hat{z},\tau \right),t\right)\right)}^T,\\ {\eta }_1\left(\tau \right)={\displaystyle \underset{\tau }{\overset{T}{\int }}B\left({\hat{x}}_1,{\hat{\dot{x}}}_1,{\hat{z}}_1\right)-{\mu }^{-1}\frac{\partial H}{\partial x_1}{N}^{-1}{\left(y-H\left({\hat{x}}_1\left({\hat{z}}_1,\tau \right),t\right)\right)}^T,}\\ {\stackrel{\ddot{^}}{x}}_2=f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)-{\mu }^{-1}\frac{\partial H}{\partial x_2}{N}^{-1}\left(y-H\left(\hat{x}\left(\hat{z},\tau \right),t\right)\right)\\ \dots (14)\\ {\eta }_k\left(\tau \right)=-{\displaystyle \underset{\tau }{\overset{T}{\int }}{\displaystyle \sum _{i=1}^k{\gamma }_i^{k}B\left({\hat{x}}_i\left({\hat{z}}_i,\tau \right),t\right)dt,\stackrel{\dot{^}}{z}={\eta }_k\left(\tau \right)}}\\ {\stackrel{\ddot{^}}{x}}_0=f\left({\hat{x}}_i,{\stackrel{\dot{^}}{x}}_i,{\hat{z}}_i\right)-{\mu }^{-1}\frac{\partial H}{\partial x}{N}^{-1}{\left(y-H\left(\hat{x}\left(\hat{z},\tau \right),t\right)\right)}^T,\\ {\hat{x}}_i\left(t_0\right)=x\left(0\right),{\stackrel{\dot{^}}{x}}_i\left(t_0\right)=\dot{x}\left(0\right){\hat{z}}_i\left(t_0\right)=z\left(0\right),\tau \in \left[\mathrm{0,}T\right],\end{array}$$

где $${\gamma }_i^k$$ определяется по правилу $${\gamma }_i^k=\left[1-{\alpha }_{i-1}\right]\cdot \dots \cdot \left[1-{\alpha }_{k-1}\right],i=\overline{\mathrm{1,}k}.$$

Пусть задача оценки $${\hat{z}}_k\left(\tau \right)$$ процесса (1)-(3) рассматривается для момента времени τ, когда желательно получить оценку, соответствующую точке T интервала наблюдения [0,T], а T увеличивается. Требуется получить $${\hat{z}}_k$$ как функцию от Г, доставляющую минимум критерию (5). Для решения данной задачи преобразуем уравнение (14) к виду двухточечной краевой задачи. Для этого продифференцируем правую часть уравнения (14) по параметру τ в результате преобразований с учетом правила Лейбница получим

$$\frac{d{\eta }_k\left(\tau \right)}{d\tau }={\displaystyle \sum _{i=1}^k{\gamma }_i^k}B\left(\tau ,{\hat{x}}_i\left({\hat{z}}_i,\left(\tau \right)\right)\right).(15)$$

Полученная каждая k - тая двухточечная краевая задача выглядит следующим образом

$$\begin{array}{cc}\begin{array}{l}\hat{\dot{z}}k={\eta }_k\left(\tau \right),\\ {\dot{\eta }}_k\left(\tau \right)={\displaystyle \sum _{i=1}^k{\gamma }_i^{k}B\left(\tau ,{\hat{x}}_i\right),}\\ {\stackrel{\ddot{^}}{x}}_0=f\left({\hat{x}}_i,{\stackrel{\dot{^}}{x}}_i,{\hat{z}}_i\right)-{\mu }^{-1}\frac{\partial H}{\partial x}{N}^{-1}{\left(y-H\left(\hat{x}\left(\hat{z},\tau \right),t\right)\right)}^T\\ n_k\left(T\right)=\mathrm{0,}{z}_k\left(0\right)=z^0+P^0{\eta }_k\left(0\right),\end{array}& (16)\end{array}$$

где P⁰ - некоторая матрица размера m×m, $$i=\overline{\mathrm{1,}k}$$ .

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

$$\begin{array}{l}{\hat{\dot{z}}}_0=P_0{G}_0^T\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\gamma }_0^1\left(y-H\left({\hat{x}}_0\left({\hat{z}}_0\right),t\right)\right),\\ {\dot{P}}_0=I-P_0{G}_0^T\frac{\partial }{\partial \hat{x}}\left\{\frac{\partial H^T}{\partial x}{N}^{-1}{\gamma }_0^1\left(y-H\left({\hat{x}}_0\left({\hat{z}}_0\right),t\right)\right)\right\}{G}_0{P}_0,\\ {\ddot{G}}_0=\frac{\partial f}{\partial \hat{\dot{x}}}{\dot{G}}_0+\frac{\partial f}{\partial \hat{x}}{G}_0+\frac{\partial f}{\partial \hat{x}},\\ {\hat{\ddot{x}}}_0=f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)-{\mu }^{-1}\frac{\partial H^T}{\partial x}{N}^{-1}\left(y-H\left({\hat{x}}_0\left({\hat{z}}_0\right),t\right)\right),\\ \dots \\ \begin{array}{cc}\begin{array}{l}{\hat{\dot{z}}}_k=P_k{G}_k^T\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\hat{x}\left(\hat{z}\right),t\right)\right)},\\ {\dot{P}}_k=I-P_k{G}_k^T\frac{\partial }{\partial \hat{x}}\left\{\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left({\hat{x}}_i\left({\hat{z}}_i\right),t\right)\right)}\right\}{G}_k{P}_k,\\ {\ddot{G}}_k=\frac{\partial f}{\partial \hat{\dot{x}}}{\dot{G}}_k+\frac{\partial f}{\partial \hat{x}}{G}_k+\frac{\partial f}{\partial \hat{z}},G_k\left(0\right)={\dot{G}}_k\left(0\right)=\mathrm{0,}\\ {\hat{\ddot{x}}}_0=f\left({\hat{x}}_k,{\hat{\dot{x}}}_k,{\hat{z}}_k\right)-{\mu }^{-1}\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}\left(y-H\left({\hat{x}}_i\left({\hat{z}}_i\right),t\right)\right),\end{array}& (18)\end{array}\\ P_i\left(t_0\right)=P\left(0\right),{\hat{x}}_i\left(0\right)=x\left(0\right),{\hat{\dot{x}}}_i\left(0\right)=\dot{x}\left(0\right).i=\overline{\mathrm{1,}k}\end{array}$$

где I - единичная матрица.

Первая группа уравнений (18) представляет собой рекуррентные уравнения последовательной идентификации, а каждое последующее уравнение оценки для к+1 использует в качестве входных параметров y, $${\hat{z}}_1,\dots ,{\hat{z}}_k$$ и $${\hat{x}}_1,\dots ,{\hat{x}}_k$$ , что позволяет получить более точные значения идентифицируемых параметров.

Устройство работает следующим образом. В исходном состоянии на вход устройства подается у, а в первом блоке 1.1 хранения констант записываются значения µ^-1, $$\hat{x},\hat{\dot{x}},\hat{z}$$ , Δt, γ, k, N^-1, P, G, Ġ, I. Значение µ^-1 с первого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход первого блока 4.1 формирования произведения. Значение $${\hat{x}}_k$$ со второго информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход первого блока 2.1 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и на первый информационный вход второго блока 2.2 формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ . Значение $${\hat{\dot{x}}}_k$$ с третьего информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход первого блока 2.1 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ . Значение $${\hat{z}}_k$$ с четвертого информационного выхода первого блока 1.1 хранения констант поступает на третий информационный вход первого блока 2.1 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и на второй информационный вход второго блока 2.2 формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ . Значение Δt с пятого информационного выхода блока 1.1 хранения констант поступает на второй информационный вход первого 8,1, второго 8.2, третьего 8.3 и шестого 8.6 блоков интегрирования, а также на первый информационный вход четвертого 8.4 и пятого 8.5 блоков интегрирования. Значение у с шестого информационного выхода блока 1.1 хранения констант поступает на второй информационный вход второго блока 5.2 формирования разности. Значение $${\gamma }_i^k$$ с седьмого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход третьего блока 2.3 формирования функции (1-γ_i-1). Значение k с восьмого информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход третьего блока 2.3 формирования функции (1-γ_i-1) и на третий информационный вход первого блока 6.1 формирования суммы. Значение N^-1 с девятого информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход третьего блока 4.3 формирования произведения. Значение $${\dot{P}}_k$$ с десятого информационного выхода первого блока 1.1 хранения констант поступает на первый и второй информационные входы соответственно шестого блока 4.6 и девятого блока 4.9 формирования произведения. Значение Ġ_k c одиннадцатого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход второго блока 7.2 транспонирования, а также на первые информационные входы восьмого блока 4.8 и девятого блока 4.9 формирования произведения. Значение Ġ_k с двенадцатого информационного

выхода первого блока 1.1 хранения констант поступает на первый информационный вход седьмого блока 4.7 формирования произведения. Значение I с тринадцатого информационного выхода первого блока 1.1 хранения констант поступает на первый информационный вход третьего блока 5.3 формирования разности. Значение $$f\left({\hat{x}}_k,{\hat{\dot{x}}}_k,{\hat{z}}_k\right)$$ с информационного выхода первого блока 2.1 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ поступает на второй информационный вход первого блока 5.1 формирования разности, а также на первые информационные входы второго блока 3.2, третьего блока 3.3 и четвертого блока 3.4 формирования производной. Значение $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ с информационного выхода второго блока 2.2 формирования функции $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ , поступает на первые информационные входы первого блока 7.1 транспонирования и второго блока 5.2 формирования разности, с информационного выхода которого значение $$y-H\left({\hat{x}}_i\left({\hat{z}}_i\right),t\right)$$ поступает на второй информационный вход второго блока 4.2 формирования произведения, с информационного выхода которого значение $${\gamma }_i^k\left(y-H\left({\hat{x}}_i\left(\widehat{z_i}\right),t\right)\right)$$ поступает на первые информационные входы второго блока 1.2 хранения констант и первого блока 6.1 формирования суммы. Значение $${\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)$$ с информационного выхода второго блока 1.2 хранения констант поступает на второй информационный вход первого блока 6.1 формирования суммы, с информационного выхода которого значение $${\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}$$ поступает на первый информационный вход третьего блока 4.3 формирования произведения, с информационного выхода которого значение $$N^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}$$ поступает на первый информационный вход четвертого блока 4.4 формирования произведения. Значение H^T с информационного выхода первого блока 7.1 транспонирования поступает на первый информационный вход первого блока 3.1 формирования производной, с информационного выхода которого значение $$\frac{\partial H^T}{\partial \hat{x}}$$ поступает на второй информационный вход четвертого блока 4.4 формирования произведения, с информационного выхода которого значение $$\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}$$ поступает соответственно на второй информационный вход первого блока 4.1 формирования произведения и на первые информационные входы пятого блока 4.5 формирования произведения и пятого блока 3.5 формирования производной, с информационного выхода которого значение $$\frac{\partial }{\partial \hat{x}}\left\{\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}\right\}$$ поступает на первый информационный вход десятого блока 4.10 формирования произведения. Значение $$G_k^T$$ с информационного выхода второго блока 7.2 транспонирования поступает на второй информационный вход шестого блока 4.6 формирования произведения, с информационного выхода которого значение $$P_k{G}_k^T$$ поступает на вторые информационные входы пятого блока 4.5 и десятого блока 4.10 формирования произведения, с информационного выхода которого значение $$P_k{G}_k^T\frac{\partial }{\partial \hat{x}}\left\{\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}\right\}$$ поступает на второй информационный вход одиннадцатого блока 4.11 формирования произведения. Значение G_kP_k информационного выхода девятого блока 4.9 формирования произведения поступает на первый информационный вход одиннадцатого блока 4.11 формирования произведения, с информационного выхода которого значение $$P_k{G}_k^T\frac{\partial }{\partial \hat{x}}\left\{\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}\right\}{G}_k{P}_k$$ поступает на второй информационный вход третьего блока 5.3 формирования разности, с информационного выхода которого значение $$I-P_k{G}_k^T\frac{\partial }{\partial \hat{x}}\left\{\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}\right\}{G}_k{P}_k$$ поступает на первый информационный вход шестого блока 8.6 интегрирования, с информационного выхода которого значение P_k поступает на первый и второй информационные входы соответственно шестого блока 4.6 и девятого блока 4.9 формирования произведения. Значение $${\mu }^{-1}\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}$$ с информационного выхода первого блока 4.1 формирования произведения поступает на первый информационный вход первого блока 5.1 формирования разности, с информационного выхода которого значение $$f\left({\hat{x}}_k,{\hat{\dot{x}}}_k,{\hat{z}}_k\right){\mu }^{-1}\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}$$ поступает на первый информационный вход первого блока 8.1 интегрирования, с информационного выхода которого значение $${\hat{\dot{x}}}_k$$ поступает на второй информационный вход первого блока 2.1 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и первый информационный вход второго блока 8.2 интегрирования, с информационного выхода которого значение $${\hat{x}}_i$$ поступает на первые информационные входы первого блока 2.1 и второго блока 2.2 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и $$H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)$$ соответственно. Значение $$P_k{G}_k^T\frac{\partial H^T}{\partial \hat{x}}{N}^{-1}{\displaystyle \sum _{i=1}^k{\gamma }_i^k\left(y-H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)\right)}$$ с информационного выхода пятого блока 4.5 формирования умножения поступает на первый информационный вход третьего блока 8.3 интегрирования, с информационного выхода которого снимается значение искомого вектора идентификации параметров z и подается на второй и третий информационные входы соответственно второго блока 2.2 и первого блока 2.1 формирования функции $$H\left(\widehat{x_i}\left(\widehat{z_i}\right),t\right)$$ и $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ соответственно. Значение $$\frac{\partial f}{\partial \hat{\dot{x}}}$$ с информационного выхода второго блока 3.2 формирования производной поступает на второй информационный вход седьмого блока 4.7 формирования произведения, с информационного выхода которого значение $$\frac{\partial f}{\partial \hat{\dot{x}}}{\dot{G}}_k$$ поступает на первый информационный вход второго блока 6.2 формирования суммы. Значение $$\frac{\partial f}{\partial \hat{x}}$$ с информационного выхода третьего блока 3.3 формирования производной поступает на второй информационный вход восьмого блока 4.8 формирования произведения, с информационного выхода которого значение $$\frac{\partial f}{\partial \hat{x}}{G}_k$$ поступает на второй информационный вход второго блока 6.2 формирования суммы, с информационного выхода которого значение $$\frac{\partial f}{\partial \hat{\dot{x}}}{\dot{G}}_k+\frac{\partial f}{\partial \hat{x}}{G}_k$$ поступает на первый информационный вход третьего блока 6.3 формирования суммы. Значение $$\frac{\partial f}{\partial \hat{x}}$$ с информационного выхода четвертого блока 3.4 формирования производной поступает на второй информационный вход третьего блока 6.3 формирования суммы, с информационного выхода которого значение $${\ddot{G}}_k$$ поступает на второй информационный вход пятого блока 8.5 интегрирования, с информационного выхода которого значение Ġ_k поступает на первый информационный вход седьмого блока 4.7 формирования произведения и второй информационный вход четвертого блока 8.4 интегрирования, с информационного выхода которого значение G_k поступает на первый информационный вход второго блока 7.2 транспонирования, а также на первые информационные входы восьмого блока 4.8 и девятого блока 4.9 формирования произведения.

Первый блок интегрирования 8.1 (фиг.2) работает следующим образом. С информационного выхода первого блока 5.1 формирования разности значение $$\hat{\ddot{x}}$$ поступает на первый информационный вход блока 8.1.1 формирования произведения. Значение Δt с пятого информационного выхода первого блока 1.1 хранения констант поступает на второй информационный вход первого 8.1 блока интегрирования, с информационного выхода которого значение $$\widehat{\ddot{x}\cdot \Delta t}$$ поступает на первый информационный вход блока 8.1.2 формирования суммы. Значение $${\hat{\dot{x}}}_{i-1}$$ информационного выхода линии задержки 8.1.3 поступает на второй информационный вход блока 8.1.2 формирования суммы, с информационного выхода которого снимается значение первообразной $${\hat{\dot{x}}}_i$$ и подается на первый вход линии задержки 8.1.3 и на второй информационный вход первого блока 2.1 формирования функции $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ и первый информационный вход второго блока 8.2 интегрирования. Остальные блоки интегрирования работают аналогично.

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

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

2. Колесников А.А. Синергетическая теория управления. Таганрог: ТРТУ. -М.: Энергоатомиздат, 1994, с.117-124.

3. Лурье А.И. Аналитическая механика. -М.: Гос.изд. физ. - мат.лит., 1961, с.642-648.

4. Сейдж Э.П., Меле Д.Л. Идентификация систем управления. -М.: Наука, 1974. с.204-212.

5. Справочник по теории автоматического управления. / Под ред. А.А. Красовского. М.: Наука, 1987, с.347-360.

6. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука, 1986, с.15-18.

7. Фарина А., Студер Ф. Цифровая обработка радиолокационной информации. Сопровождение целей. - М.: Радио и связь, 1993, с.38-47.

Устройство идентификации Лагранжевых динамических систем на основе итерационной регуляризации, содержащее первый блок хранения констант; первый, второй блоки формирования функций $$f\left({\hat{x}}_i,{\hat{\dot{x}}}_i,{\hat{z}}_i\right)$$ , $$H\left({\hat{x}}_i,\left({\hat{z}}_i\right),t\right)$$ соответственно; первый, второй, третий, четвертый, пятый блоки формирования производной; первый, второй, третий, четвертый, пятый, шестой, седьмой, восьмой, девятый, десятый, одиннадцатый блоки формирования произведения; первый, второй, третий блоки формирования разности; первый, второй блоки формирования суммы; первый, второй блоки транспонирования; первый, второй, третий, четвертый, пятый, шестой блоки интегрирования; при этом первый информационный вход первого блока хранения констант является входом устройства, а информационный выход третьего блока интегрирования является выходом устройства, отличающееся тем, что в устройство введены второй блок хранения констант, третий блок формирования функции (1-γ_i-1), третий блок формирования суммы, причем информационный выход второго блока формирования произведения соединен с информационным входом второго блока хранения констант, информационный выход которого соединен со вторым входом первого блока формирования суммы; седьмой информационный выход первого блока хранения констант соединен с первым информационным входом третьего блока формирования функции (1-γ_i-1); восьмой информационный выход первого блока хранения констант соединен со вторым информационным входом третьего блока формирования функции (1-γ_i-1), выход которого соединен со вторым информационным входом второго блока формирования произведения; информационный выход второго блока формирования суммы соединен с первым информационным входом третьего блока формирования суммы; информационный выход четвертого блока формирования производной соединен со вторым информационным входом третьего блока формирования суммы, информационный выход которого соединен со вторым информационным входом пятого блока интегрирования.