• Bulletin of the Korean Mathematical Society
• /
• v.60 no.5
• /
• pp.1181-1199
• /
• 2023

For solving complex symmetric positive definite linear systems, we propose a single step real-valued (SSR) iterative method, which does not involve the complex arithmetic. The upper bound on the spectral radius of the iteration matrix of the SSR method is given and its convergence properties are analyzed. In addition, the quasi-optimal parameter which minimizes the upper bound for the spectral radius of the proposed method is computed. Finally, numerical experiments are given to demonstrate the effectiveness and robustness of the propose methods.

• Journal of the Korean Mathematical Society
• /
• v.55 no.1
• /
• pp.73-82
• /
• 2018

In this article we consider certain types of nonlinear matrix equations including the stochastic rational Riccati equation and show the existence and uniqueness of the positive definite solution by using Bhaskar-Lakshmikantham's coupled fixed point theorem.

• Journal of applied mathematics & informatics
• /
• v.5 no.2
• /
• pp.383-392
• /
• 1998

This paper proposes Generalized Stationary Iterative called GSI method. It is shown that the existing stationary iterative methods are special cases of GSI method. Convergence properties of this method are provided and their numerical experiments for linear systems with symmetric positive definite matrix are also provided.

• Journal of applied mathematics & informatics
• /
• v.31 no.3_4
• /
• pp.545-558
• /
• 2013

In this paper, the formulas for calculating the extremal ranks and inertias of the Hermitian least squares solutions to matrix equation AX = B are established. In particular, the necessary and sufficient conditions for the existences of the positive and nonnegative definite solutions to this matrix equation are given. Meanwhile, the least squares problem of the above matrix equation with Hermitian R-symmetric and R-skew symmetric constraints are also investigated.

• Journal of the Korean Data and Information Science Society
• /
• v.14 no.4
• /
• pp.1007-1012
• /
• 2003

A recursive procedure for finding the Cholesky root of the inverse of sample covariance matrix, leading to a direct solution for the inverse of a positive definite matrix, is developed using the likelihood equation for the maximum likelihood estimation of the Cholesky root under normality assumptions. An example of the Hilbert matrix is considered for an illustration of the procedure.

• Communications for Statistical Applications and Methods
• /
• v.24 no.1
• /
• pp.81-96
• /
• 2017

Cumulative logit random effects models are typically used to analyze longitudinal ordinal data. The random effects covariance matrix is used in the models to demonstrate both subject-specific and time variations. The covariance matrix may also be homogeneous; however, the structure of the covariance matrix is assumed to be homoscedastic and restricted because the matrix is high-dimensional and should be positive definite. To satisfy these restrictions two Cholesky decomposition methods were proposed in linear (mixed) models for the random effects precision matrix and the random effects covariance matrix, respectively: modified Cholesky and moving average Cholesky decompositions. In this paper, we use these two methods to model the random effects precision matrix and the random effects covariance matrix in cumulative logit random effects models for longitudinal ordinal data. The methods are illustrated by a lung cancer data set.

• 제어로봇시스템학회:학술대회논문집
• /
• 2003.10a
• /
• pp.2004-2009
• /
• 2003

A robust optimal control design is proposed in this study for rigid robotic systems under the unknown load and the other uncertainties. The uncertainties are quadratically bounded for some positive definite matrix. Iterative method finding the Q weighting matrix is shown. Computer simulations have been done for a weight-lifting operation of a two-link manipulator and the result of the simulation shows that the proposed algorithm is very effective for a robust control of robotic systems.

• Communications for Statistical Applications and Methods
• /
• v.25 no.1
• /
• pp.61-70
• /
• 2018

Modeling of the random effects covariance matrix in generalized linear mixed models (GLMMs) is an issue in analysis of longitudinal categorical data because the covariance matrix can be high-dimensional and its estimate must satisfy positive-definiteness. To satisfy these constraints, we consider the autoregressive and moving average Cholesky decomposition (ARMACD) to model the covariance matrix. The ARMACD creates a more flexible decomposition of the covariance matrix that provides generalized autoregressive parameters, generalized moving average parameters, and innovation variances. In this paper, we analyze longitudinal count data with overdispersion using GLMMs. We propose negative binomial loglinear mixed models to analyze longitudinal count data and we also present modeling of the random effects covariance matrix using the ARMACD. Epilepsy data are analyzed using our proposed model.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.455-458
• /
• 2003

An optimal control approach to robust control design is proposed in this study for rigid robotic systems under the unknown load and the other uncertainties. The uncertainties are quadratically bounded for some positive definite matrix. Iterative method to find the matrix is shown. Simulations arc made for a weight-lifting operation of a two-link manipulator and the robust control performance of robotic systems by the proposed algorithm is remarkable.

• Communications of the Korean Mathematical Society
• /
• v.9 no.3
• /
• pp.753-765
• /
• 1994

The linearly constrained quadratic programming(QP) considered is : $$ min f(x) = c^T x + \frac{1}{2}x^T Hx $$ $$ (1) subject to A^T x \geq b,$$ where $c,x \in R^n, b \in R^m, H \in R^{n \times n)}$, symmetric, and $A \in R^{n \times n}$. If there are bounds on x, these are included in the matrix $A^T$. The Hessian matrix H may be positive definite or negative semi-difinite. For large problems H and the constraint matrix A are assumed to be sparse.