• Journal of applied mathematics & informatics
• /
• v.26 no.3_4
• /
• pp.471-479
• /
• 2008

Let $R\;{\in}\;C^{m{\times}m}$ and $S\;{\in}\;C^{n{\times}n}$ be nontrivial unitary involutions, i.e., $R^*\;=\;R\;=\;R^{-1}\;{\neq}\;I_m$ and $S^*\;=\;S\;=\;S^{-1}\;{\neq}\;I_m$. We say that $G\;{\in}\;C^{m{\times}n}$ is a generalized reflexive matrix if RGS = G. The set of all m ${\times}$ n generalized reflexive matrices is denoted by $GRC^{m{\times}n}$. In this paper, an efficient method for the least squares solution $X\;{\in}\;GRC^{m{\times}n}$ of the matrix equation AXB = D with arbitrary coefficient matrices $A\;{\in}\;C^{p{\times}m}$, $B\;{\in}\;C^{n{\times}q}$and the right-hand side $D\;{\in}\;C^{p{\times}q}$ is developed based on the canonical correlation decomposition(CCD) and, an explicit formula for the general solution is presented.

• Journal of applied mathematics & informatics
• /
• v.27 no.5_6
• /
• pp.1061-1071
• /
• 2009

For real generalized reflexive matrices R, S, i.e., $R^T$ = R, $R^2$ = I, $S^T$ = S, $S^2$ = I, we say that real matrix X is (R,S)-symmetric, if RXS = X. In this paper, an iterative algorithm is proposed to solve the least-squares problem of matrix equation AXB = C with (R,S)-symmetric X. Furthermore, the optimal approximation solution to given matrix $X_0$ is also derived by this iterative algorithm. Finally, given numerical example and its convergent curve show that this method is feasible and efficient.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.3
• /
• pp.57-65
• /
• 1990

The finite element menthod for the investigation of the static and dynamic stability of the plane framed structures subjected to non-conservative forces is presented. By using the Hermitian polynomial as the shape function, the geometric stiffness matrix, the load correction stiffness matrix for non-conservative forces, and the matrix equation of internal and external damping are derived. Then, a matrix equation of the motion for the non-conservative system is formulated and the critical divergence and flutter loads are determined from this equation.

• Journal of Institute of Control, Robotics and Systems
• /
• v.18 no.9
• /
• pp.806-811
• /
• 2012

In this paper Hartley functions are used to approximate the solutions of continuous time linear dynamical system. The Hartley function and its integral operational matrix are first presented, an efficient algorithm to solve the Stein equation is proposed. The algorithm is based on the compound matrix and the inverse of sum of matrices. Using the structure of the Hartley's integral operational matrix, the full order Stein equation should be solved in terms of the solutions of pure algebraic matrix equations, which reduces the computation time remarkably. Finally a numerical example is illustrated to demonstrate the validity of the proposed algorithm.

• Bulletin of the Korean Mathematical Society
• /
• v.55 no.2
• /
• pp.431-448
• /
• 2018

This paper is concerned with the positive definite solutions of the nonlinear matrix equation $X-A^*{\bar{X}}^{-1}A=Q$, where A, Q are given complex matrices with Q positive definite. We show that such a matrix equation always has a unique positive definite solution and if A is nonsingular, it also has a unique negative definite solution. Moreover, based on Sherman-Morrison-Woodbury formula, we derive elegant relationships between solutions of $X-A^*{\bar{X}}^{-1}A=I$ and the well-studied standard nonlinear matrix equation $Y+B^*Y^{-1}B=Q$, where B, Q are uniquely determined by A. Then several effective numerical algorithms for the unique positive definite solution of $X-A^*{\bar{X}}^{-1}A=Q$ with linear or quadratic convergence rate such as inverse-free fixed-point iteration, structure-preserving doubling algorithm, Newton algorithm are proposed. Numerical examples are presented to illustrate the effectiveness of all the theoretical results and the behavior of the considered algorithms.

• Journal of applied mathematics & informatics
• /
• v.12 no.1_2
• /
• pp.165-182
• /
• 2003

Toeplitz matrix method and the product Nystrom method are described for mixed Fredholm-Volterra singular integral equation of the second kind with Carleman Kernel and logarithmic kernel. The results are compared with the exact solution of the integral equation. The error of each method is calculated.

• Journal of IKEEE
• /
• v.23 no.2
• /
• pp.614-621
• /
• 2019

In this paper, we presents the integral equation-fast Fourier transform(IE-FFT) and block matrix preconditioner (BMP) to solve electromagnetic scattering problems of penetrable structures composed of dielectric or magnetic materials. IE-FFT can significantly improve the amount of calculation to solve the matrix equation constructed from the moment method(MoM). Moreover, the iterative method in conjunction with BMP can be significantly reduce the number of iterations required to solve the matrix equations which are constructed from electrically large structures. Numerical results show that IE-FFT and block matrix preconditioner can solve electromagnetic scattering problems for penetrable objects quickly and accurately.

• Journal of the Korean Mathematical Society
• /
• v.50 no.4
• /
• pp.755-770
• /
• 2013

One of the interesting nonlinear matrix equations is the quadratic matrix equation defined by $$Q(X)=AX^2+BX+C=0$$, where X is a $n{\times}n$ unknown real matrix, and A, B and C are $n{\times}n$ given matrices with real elements. Another one is the matrix polynomial $$P(X)=A_0X^m+A_1X^{m-1}+{\cdots}+A_m=0,\;X,\;A_i{\in}\mathbb{R}^{n{\times}n}$$. Newton's method is used to find the symmetric and bisymmetric solvents of the nonlinear matrix equations Q(X) and P(X). The method does not depend on the singularity of the Fr$\acute{e}$chet derivative. Finally, we give some numerical examples.

• International Journal of Control, Automation, and Systems
• /
• v.6 no.5
• /
• pp.776-784
• /
• 2008

In this paper, new upper matrix bounds for the solution of the continuous algebraic Riccati equation (CARE) are derived. Following the derivation of each bound, iterative algorithms are developed for obtaining sharper solution estimates. These bounds improve the restriction of the results proposed in a previous paper, and are more general. The proposed bounds are always calculated if the stabilizing solution of the CARE exists. Finally, numerical examples are given to demonstrate the effectiveness of the present schemes.

• Bulletin of the Korean Mathematical Society
• /
• v.54 no.1
• /
• pp.199-214
• /
• 2017

We consider fixed-point iterations constructed by simple transforming from a quadratic matrix equation to equivalent fixed-point equations and assume that the iterations are well-defined at some solutions. In that case, we suggest real valued functions. These functions provide radii at the solution, which guarantee the local convergence and the uniqueness of the solutions. Moreover, these radii obtained by simple calculations of some constants. We get the constants by arbitrary matrix norm for coefficient matrices and solution. In numerical experiments, the examples show that the functions give suitable boundaries which guarantee the local convergence and the uniqueness of the solutions for the given equations.