In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied tho...In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.展开更多
Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of...Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.展开更多
We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, ...We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.展开更多
The symmetric positive definite solutions of matrix equations (AX,XB)=(C,D) and AXB=C are considered in this paper. Necessary and sufficient conditions for the matrix equations to have symmetric positive de...The symmetric positive definite solutions of matrix equations (AX,XB)=(C,D) and AXB=C are considered in this paper. Necessary and sufficient conditions for the matrix equations to have symmetric positive definite solutions are derived using the singular value and the generalized singular value decompositions. The expressions for the general symmetric positive definite solutions are given when certain conditions hold.展开更多
Generally unitary solution to the system of martix equations over the quaternion field [X mA ns =B ns ,X nn C nt =D nt ] is considered. A necessary and sufficient condition for the existence o...Generally unitary solution to the system of martix equations over the quaternion field [X mA ns =B ns ,X nn C nt =D nt ] is considered. A necessary and sufficient condition for the existence of and the expression for the generally unitary solution of the system are derived.展开更多
Minor self conjugate (msc) and skewpositive semidefinite (ssd) solutions to the system of matrix equations over skew fields [A mn X nn =A mn ,B sn X nn =O sn ] are considered. Necessary and su...Minor self conjugate (msc) and skewpositive semidefinite (ssd) solutions to the system of matrix equations over skew fields [A mn X nn =A mn ,B sn X nn =O sn ] are considered. Necessary and sufficient conditions for the existence of and the expressions for the msc solutions and the ssd solutions are obtained for the system.展开更多
In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solv...In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solvability and the general expression of the solution are obtained. Moreover, the related optimal approximation problem to a given matrix over the solution set is solved.展开更多
Dehghan and Hajarian, [4], investigated the matrix equations A^TXB+B^TX^TA = C and A^TXB + B^TXA = C providing inequalities for the determinant of the solutions of these equations. In the same paper, the authors pre...Dehghan and Hajarian, [4], investigated the matrix equations A^TXB+B^TX^TA = C and A^TXB + B^TXA = C providing inequalities for the determinant of the solutions of these equations. In the same paper, the authors presented a lower bound for the product of the eigenvalues of the solutions to these matrix equations. Inspired by their work, we give some generalizations of Dehghan and Hajarian results. Using the theory of the numerical ranges, we present an inequality involving the trace of C when A, B, X are normal matrices satisfying A^T B = BA^T.展开更多
We derive necessary and sufficient conditions for the existence and an expression of the (anti)reflexive solution with respect to the nontrivial generalized reflection matrix P to the system of complex matrix equati...We derive necessary and sufficient conditions for the existence and an expression of the (anti)reflexive solution with respect to the nontrivial generalized reflection matrix P to the system of complex matrix equations AX = B and XC = D. The explicit solutions of the approximation problem min x∈Ф ||X - E||F was given, where E is a given complex matrix and Ф is the set of all reflexive (or antireflexive) solutions of the system mentioned above, and ||·|| is the Frobenius norm. Furthermore, it was pointed that some results in a recent paper are special cases of this paper.展开更多
Some new oscillation criteria are established for the second-order matrix differential system (r(t)Z'(t))' + p(t)Z'(t) + Q(t)F(Z'(t))G(Z(t)) = 0, t ≥ to 〉 0, are different from most known ...Some new oscillation criteria are established for the second-order matrix differential system (r(t)Z'(t))' + p(t)Z'(t) + Q(t)F(Z'(t))G(Z(t)) = 0, t ≥ to 〉 0, are different from most known ones in the sense that they are based on the information only on a sequence of subintervals of [to,∞), rather than on the whole half-line. The results weaken the condition of Q(t) and generalize some well-known results of Wong (1999) to nonlinear matrix differential equation.展开更多
In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermit...In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.展开更多
Necessary and sufficient conditions are given for the existence of the general solution, the centrosymmetric solution, and the centroskewsymmetric solution to a system of linear matrix equations over an arbitrary skew...Necessary and sufficient conditions are given for the existence of the general solution, the centrosymmetric solution, and the centroskewsymmetric solution to a system of linear matrix equations over an arbitrary skew field. The representations of such the solutions of the system are also derived.展开更多
In this paper, solutions to the generalized Sylvester matrix equations AX -XF = BY and MXN -X = TY with A, M ∈ R^n×n, B, T ∈ Rn×r, F, N ∈ R^p×p and the matrices N, F being in companion form, are est...In this paper, solutions to the generalized Sylvester matrix equations AX -XF = BY and MXN -X = TY with A, M ∈ R^n×n, B, T ∈ Rn×r, F, N ∈ R^p×p and the matrices N, F being in companion form, are established by a singular value decomposition of a matrix with dimensions n × (n + pr). The algorithm proposed in this paper for the euqation AX - XF = BY does not require the controllability of matrix pair (A, B) and the restriction that A, F do not have common eigenvalues. Since singular value decomposition is adopted, the algorithm is numerically stable and may provide great convenience to the computation of the solution to these equations, and can perform important functions in many design problems in control systems theory.展开更多
In this paper, the maximal and minimal ranks of the solution to a system of matrix equations over H, the real quaternion algebra, were derived. A previous known result could be regarded as a special case of the new re...In this paper, the maximal and minimal ranks of the solution to a system of matrix equations over H, the real quaternion algebra, were derived. A previous known result could be regarded as a special case of the new result.展开更多
In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
We consider matrix integrable fifth-order mKdV equations via a kind of group reductions of the Ablowitz–Kaup–Newell–Segur matrix spectral problems. Based on properties of eigenvalue and adjoint eigenvalue problems,...We consider matrix integrable fifth-order mKdV equations via a kind of group reductions of the Ablowitz–Kaup–Newell–Segur matrix spectral problems. Based on properties of eigenvalue and adjoint eigenvalue problems, we solve the corresponding Riemann–Hilbert problems, where eigenvalues could equal adjoint eigenvalues, and construct their soliton solutions, when there are zero reflection coefficients. Illustrative examples of scalar and two-component integrable fifthorder mKdV equations are given.展开更多
In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared wi...In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.展开更多
An AOR(Accelerated Over-Relaxation)iterative method is suggested by introducing one more parameter than SOR(Successive Over-Relaxation)method for solving coupled Lyapunov matrix equations(CLMEs)that come from continuo...An AOR(Accelerated Over-Relaxation)iterative method is suggested by introducing one more parameter than SOR(Successive Over-Relaxation)method for solving coupled Lyapunov matrix equations(CLMEs)that come from continuous-time Markovian jump linear systems.The proposed algorithm improves the convergence rate,which can be seen from the given illustrative examples.The comprehensive theoretical analysis of convergence and optimal parameter needs further investigation.展开更多
In this paper, we introduce a Hermite operational matrix collocation method for solving higher-order linear complex differential equations in rectangular or elliptic domains. We show that based on a linear algebra the...In this paper, we introduce a Hermite operational matrix collocation method for solving higher-order linear complex differential equations in rectangular or elliptic domains. We show that based on a linear algebra theorem, the use of different polynomials such as Hermite, Bessel and Taylor in polynomial collocation methods for solving differential equations leads to an equal solution, and the difference in the numerical results arises from the difference in the coefficient matrix of final linear systems of equations. Some numerical examples will also be given.展开更多
Let P∈C^(n×n)be a Hermitian and{k+1}-potent matrix,i.e.,P^(k+1)=P=P^(*),where(·)^(*)stands for the conjugate transpose of a matrix.A matrix X∈C^(n×n)is called{P,k+1}-reflexive(anti-reflexive)if PXP=X(...Let P∈C^(n×n)be a Hermitian and{k+1}-potent matrix,i.e.,P^(k+1)=P=P^(*),where(·)^(*)stands for the conjugate transpose of a matrix.A matrix X∈C^(n×n)is called{P,k+1}-reflexive(anti-reflexive)if PXP=X(P XP=-X).The system of matrix equations AX=C,XB=D subject to{P,k+1}-reflexive and anti-reflexive constraints are studied by converting into two simpler cases:k=1 and k=2,the least squares solution and the associated optimal approximation problem are also considered.展开更多
文摘In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.
基金Supported by the National Natural Science Foundation of China(12001395)the special fund for Science and Technology Innovation Teams of Shanxi Province(202204051002018)+1 种基金Research Project Supported by Shanxi Scholarship Council of China(2022-169)Graduate Education Innovation Project of Taiyuan Normal University(SYYJSYC-2314)。
文摘Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.
基金Project supported by the National Natural Science Foundation of China(Grant No.11471262)the National Basic Research Program of China(Grant No.2012CB025904)the State Key Laboratory of Science and Engineering Computing and the Center for High Performance Computing of Northwestern Polytechnical University,China
文摘We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.
文摘The symmetric positive definite solutions of matrix equations (AX,XB)=(C,D) and AXB=C are considered in this paper. Necessary and sufficient conditions for the matrix equations to have symmetric positive definite solutions are derived using the singular value and the generalized singular value decompositions. The expressions for the general symmetric positive definite solutions are given when certain conditions hold.
文摘Generally unitary solution to the system of martix equations over the quaternion field [X mA ns =B ns ,X nn C nt =D nt ] is considered. A necessary and sufficient condition for the existence of and the expression for the generally unitary solution of the system are derived.
文摘Minor self conjugate (msc) and skewpositive semidefinite (ssd) solutions to the system of matrix equations over skew fields [A mn X nn =A mn ,B sn X nn =O sn ] are considered. Necessary and sufficient conditions for the existence of and the expressions for the msc solutions and the ssd solutions are obtained for the system.
基金supported by National Natural Science Foundation of China (10571047)and by Scientific Research Fund of Hunan Provincial Education Department of China Grant(06C235)+1 种基金by Central South University of Forestry and Technology (06Y017)by Specialized Research Fund for the Doctoral Program of Higher Education (20060532014)
文摘In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solvability and the general expression of the solution are obtained. Moreover, the related optimal approximation problem to a given matrix over the solution set is solved.
基金partially supported by FCT(Portugal)with national funds through Centro de Matemática da Universidade de Trás-os-Montes e Alto Douro(PEst-OE/MAT/UI4080/2014)
文摘Dehghan and Hajarian, [4], investigated the matrix equations A^TXB+B^TX^TA = C and A^TXB + B^TXA = C providing inequalities for the determinant of the solutions of these equations. In the same paper, the authors presented a lower bound for the product of the eigenvalues of the solutions to these matrix equations. Inspired by their work, we give some generalizations of Dehghan and Hajarian results. Using the theory of the numerical ranges, we present an inequality involving the trace of C when A, B, X are normal matrices satisfying A^T B = BA^T.
基金supported by the National Natural Science Foundation of China (Grant No.60672160)
文摘We derive necessary and sufficient conditions for the existence and an expression of the (anti)reflexive solution with respect to the nontrivial generalized reflection matrix P to the system of complex matrix equations AX = B and XC = D. The explicit solutions of the approximation problem min x∈Ф ||X - E||F was given, where E is a given complex matrix and Ф is the set of all reflexive (or antireflexive) solutions of the system mentioned above, and ||·|| is the Frobenius norm. Furthermore, it was pointed that some results in a recent paper are special cases of this paper.
基金Supported by NECC and NSF of Shandong Proyilice,China(Y2005A06).
文摘Some new oscillation criteria are established for the second-order matrix differential system (r(t)Z'(t))' + p(t)Z'(t) + Q(t)F(Z'(t))G(Z(t)) = 0, t ≥ to 〉 0, are different from most known ones in the sense that they are based on the information only on a sequence of subintervals of [to,∞), rather than on the whole half-line. The results weaken the condition of Q(t) and generalize some well-known results of Wong (1999) to nonlinear matrix differential equation.
文摘In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.
基金Supported by the National Natural Science Foundation of China(10471085)
文摘Necessary and sufficient conditions are given for the existence of the general solution, the centrosymmetric solution, and the centroskewsymmetric solution to a system of linear matrix equations over an arbitrary skew field. The representations of such the solutions of the system are also derived.
基金This work was supported by the Chinese Outstanding Youth Foundation(No.69925308)Program for Changjiang Scholars and Innovative ResearchTeam in University.
文摘In this paper, solutions to the generalized Sylvester matrix equations AX -XF = BY and MXN -X = TY with A, M ∈ R^n×n, B, T ∈ Rn×r, F, N ∈ R^p×p and the matrices N, F being in companion form, are established by a singular value decomposition of a matrix with dimensions n × (n + pr). The algorithm proposed in this paper for the euqation AX - XF = BY does not require the controllability of matrix pair (A, B) and the restriction that A, F do not have common eigenvalues. Since singular value decomposition is adopted, the algorithm is numerically stable and may provide great convenience to the computation of the solution to these equations, and can perform important functions in many design problems in control systems theory.
基金Project supported by the National Natural Science Foundation of China (Grant No.60672160)
文摘In this paper, the maximal and minimal ranks of the solution to a system of matrix equations over H, the real quaternion algebra, were derived. A previous known result could be regarded as a special case of the new result.
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.
基金supported in part by the National Natural Science Foundation of China (Grant Nos. 11975145, 11972291, and 51771083)the Ministry of Science and Technology of China (Grant No. G2021016032L)the Natural Science Foundation for Colleges and Universities in Jiangsu Province, China (Grant No. 17 KJB 110020)。
文摘We consider matrix integrable fifth-order mKdV equations via a kind of group reductions of the Ablowitz–Kaup–Newell–Segur matrix spectral problems. Based on properties of eigenvalue and adjoint eigenvalue problems, we solve the corresponding Riemann–Hilbert problems, where eigenvalues could equal adjoint eigenvalues, and construct their soliton solutions, when there are zero reflection coefficients. Illustrative examples of scalar and two-component integrable fifthorder mKdV equations are given.
基金Project supported by the National Natural Science Foundation of China (Grant No.10271074), and the Special Funds for Major Specialities of Shanghai Education Commission (Grant No.J50101)
文摘In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.
基金Supported by Key Scientific Research Project of Colleges and Universities in Henan Province of China(Grant No.20B110012)。
文摘An AOR(Accelerated Over-Relaxation)iterative method is suggested by introducing one more parameter than SOR(Successive Over-Relaxation)method for solving coupled Lyapunov matrix equations(CLMEs)that come from continuous-time Markovian jump linear systems.The proposed algorithm improves the convergence rate,which can be seen from the given illustrative examples.The comprehensive theoretical analysis of convergence and optimal parameter needs further investigation.
文摘In this paper, we introduce a Hermite operational matrix collocation method for solving higher-order linear complex differential equations in rectangular or elliptic domains. We show that based on a linear algebra theorem, the use of different polynomials such as Hermite, Bessel and Taylor in polynomial collocation methods for solving differential equations leads to an equal solution, and the difference in the numerical results arises from the difference in the coefficient matrix of final linear systems of equations. Some numerical examples will also be given.
基金Supported by the Education Department Foundation of Hebei Province(QN2015218)Supported by the Natural Science Foundation of Hebei Province(A2015403050)
文摘Let P∈C^(n×n)be a Hermitian and{k+1}-potent matrix,i.e.,P^(k+1)=P=P^(*),where(·)^(*)stands for the conjugate transpose of a matrix.A matrix X∈C^(n×n)is called{P,k+1}-reflexive(anti-reflexive)if PXP=X(P XP=-X).The system of matrix equations AX=C,XB=D subject to{P,k+1}-reflexive and anti-reflexive constraints are studied by converting into two simpler cases:k=1 and k=2,the least squares solution and the associated optimal approximation problem are also considered.