Abstract:
The paper presents a new method for calculating the values of derivatives in the LD factorization of parametrized matrices, based on the direct procedure for the modified weighted Gram–Schmidt orthogonalization.
The need for calculating the values of derivatives in matrix orthogonal transformations arises in the theory of perturbations and control, in differential geometry, in solving problems such as the Lyapunov exponential calculation, the problems of automatic differentiation, the calculation of the numerical solution of the matrix differential Riccati equation, the calculation of high-order derivatives in the optimal input design. In the theory of parameter identification of mathematical models of discrete linear stochastic systems, such problems are solved by developing numerically effective algorithms for finding the solution of the matrix difference Riccati sensitivity equation.
In this paper, we have posed and solved a new problem of calculating the values of derivatives. Lemma 1 represents the main theoretical result. The practical result is the computational algorithm 2. The software implementation of the algorithm allows us to calculate the values of derivatives of the parametrized matrices that are the result of a direct procedure of the LD factorization quickly and with high accuracy. It is not necessary to calculate the values of derivatives of the matrix of weighted orthogonal transformation. The algorithm has a simple structure and does not contain complex operations of symbolic or numerical differentiation. Only one inversion of the triangular matrix and simple matrix operations of addition and multiplication are required.
Two numerical examples are considered that show the operability and numerical efficiency of the proposed algorithm 2.
The results obtained in this paper will be used to construct new classes of adaptive LD filters in the area of parameter identification of mathematical models of discrete linear stochastic systems.
Citation:
J. V. Tsyganova, A. V. Tsyganov, “On the computation of derivatives within LD factorization of parametrized matrices”, Bulletin of Irkutsk State University. Series Mathematics, 23 (2018), 64–79
\Bibitem{TsyTsy18}
\by J.~V.~Tsyganova, A.~V.~Tsyganov
\paper On the computation of derivatives within LD factorization of parametrized matrices
\jour Bulletin of Irkutsk State University. Series Mathematics
\yr 2018
\vol 23
\pages 64--79
\mathnet{http://mi.mathnet.ru/iigum331}
\crossref{https://doi.org/10.26516/1997-7670.2018.23.64}
Linking options:
https://www.mathnet.ru/eng/iigum331
https://www.mathnet.ru/eng/iigum/v23/p64
This publication is cited in the following 3 articles:
Andrey Tsyganov, Julia Tsyganova, “Algorithmic Differentiation of the MWGS-Based Arrays for Computing the Information Matrix Sensitivity Equations within the Problem of Parameter Identification”, Mathematics, 10:1 (2022), 126
Julia V. Tsyganova, Maria V. Kulikova, Andrey V. Tsyganov, “A general approach for designing the MWGS-based information-form Kalman filtering methods”, European Journal of Control, 56 (2020), 86
Innokentiy V. Semushin, Julia V. Tsyganova, Andrey V. Tsyganov, “Numerically Efficient LD-computations for the Auxiliary Performance Index Based Control Optimization under Uncertainties”, IFAC-PapersOnLine, 51:32 (2018), 568