Matrix Riccati Equations in Optimal Control

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.

Numerical Methods for a Class of Quadratic Matrix Equations

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.

Solution of the Matrix Second Semi-Tensor Product Equation A ∘ l X ∘ l B=C

In this paper, the solution of the matrix second semi-tensor product equation A∘lX∘lB=Cis studied. Firstly, the solvability of the matrix-vector second semi-tensor product equation is investigated. At the same time, the compatibility conditions, the sufficient and necessary conditions and the specific solution methods for the matrix solution are given. Secondly, we further consider the solvability of the second semi-tensor product equation of the matrix. For each part, several examples are given to illustrate the validity of the results.

An Extended Numerical Method by Stancu Polynomials for Solution of Integro-Differential Equations Arising in Oscillating Magnetic Fields

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.

Matrix integrable fifth-order mKdV equations and their soliton solutions

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.

Hermite Positive Definite Solution of the Quaternion Matrix Equation Xm + B*XB = C

This paper discusses the necessary and sufficient conditions for the existence of Hermite positive definite solutions of the quaternion matrix equation Xm+ B*XB = C (m > 0) and its iterative solution method. According to the characteristics of the coefficient matrix, a corresponding algebraic equation system is ingeniously constructed, and by discussing the equation system’s solvability, the matrix equation’s existence interval is obtained. Based on the characteristics of the coefficient matrix, some necessary and sufficient conditions for the existence of Hermitian positive definite solutions of the matrix equation are derived. Then, the upper and lower bounds of the positive actual solutions are estimated by using matrix inequalities. Four iteration formats are constructed according to the given conditions and existence intervals, and their convergence is proven. The selection method for the initial matrix is also provided. Finally, using the complexification operator of quaternion matrices, an equivalent iteration on the complex field is established to solve the equation in the Matlab environment. Two numerical examples are used to test the effectiveness and feasibility of the given method. .

The Least Squares {P,Q,k+1}-Reflexive Solution to a Matrix Equation

Let P∈C^(m×m)and Q∈C^(n×n)be Hermitian and{k+1}-potent matrices,i.e.,P k+1=P=P∗,Qk+1=Q=Q∗,where(·)∗stands for the conjugate transpose of a matrix.A matrix X∈C m×n is called{P,Q,k+1}-reflexive(anti-reflexive)if P XQ=X(P XQ=−X).In this paper,the least squares solution of the matrix equation AXB=C subject to{P,Q,k+1}-reflexive and anti-reflexive constraints are studied by converting into two simpler cases:k=1 and k=2.

Second-order two-scale analysis and numerical algorithms for the hyperbolic–parabolic equations with rapidly oscillating coefficients

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.

SYMMETRIC POSITIVE DEFINITE SOLUTIONS OF MATRIX EQUATIONS (AX,XB)=(C,D) AND AXB=C

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 a System of Matrix Equations over the Quaternion Field

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 and Skewpositive Semidefinite Solutions to a System of Matrix Equations over Skew Fields

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.

Improved Microwave Ocean Emissivity and Reflectivity Models Derived from Two-Scale Roughness Theory

The Geometrical Optics(GO)approach and the FAST Emissivity Model(FASTEM)are widely used to estimate the surface radiative components in atmospheric radiative transfer simulations,but their applications are limited in specific conditions.In this study,a two-scale reflectivity model(TSRM)and a two-scale emissivity model(TSEM)are developed from the two-scale roughness theory.Unlike GO which only computes six non-zero elements in the reflectivity matrix,The TSRM includes 16 elements of Stokes reflectivity matrix which are important for improving radiative transfer simulation accuracy in a scattering atmosphere.It covers the frequency range from L-to W-bands.The dependences of all TSRM elements on zenith angle,wind speed,and frequency are derived and analyzed in details.For a set of downwelling radiances in microwave frequencies,the reflected upwelling brightness temperature(BTs)are calculated from both TSRM and GO and compared for analyzing their discrepancies.The TSRM not only includes the effects of GO but also accounts for the small-scale Bragg scattering effect in an order of several degrees in Kelvins in brightness temperature.Also,the third and fourth components of the Stokes vector can only be produced from the TSRM.For the emitted radiation,BT differences in vertical polarization between a TSEM and FASTEM are generally less than 5 K when the satellite zenith angle is less than 40°,whereas those for the horizontal component can be quite significant,greater than 20 K.

THE GENERALIZED REFLEXIVE SOLUTION FOR A CLASS OF MATRIX EQUATIONS (AX-B,XC=D)

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.

SOME FURTHER NOTES ON THE MATRIX EQUATIONS A^TXB+B^TX^TA=C AND A^TXB+B^TXA=C

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.

Reflexive solution to a system of matrix equations

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.

NEW OSCILLATION CRITERIA FOR SECOND-ORDER NONLINEAR MATRIX DIFFERENTIAL EQUATIONS

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.

HERMITE MATRIX POLYNOMIALS AND SECOND ORDER MATRIX DIFFERENTIAL EQUATIONS

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.

The General and Centro(Kewsy) Symmetric Solutions to a System of Matrix Equations over an Arbitrary Skew Field

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.

Least square method based on Haar wavelet to solve multi-dimensional stochastic Ito-Volterra integral equations

This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.

Solutions to the generalized Sylvester matrixequations by a singular value decomposition

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.