• Title/Summary/Keyword: symmetric positive definite matrix

Search Result 34, Processing Time 0.036 seconds

Stabilizing Solutions of Algebraic Matrix riccati Equations in TEX>$H_\infty$ Control Problems

  • Kano, Hiroyuki;Nishimura, Toshimitsu
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 1994.10a
    • /
    • pp.364-368
    • /
    • 1994
  • Algebraic matrix Riccati equations of the form, FP+PF$^{T}$ -PRP+Q=0. are analyzed with reference to the stability of closed-loop system F-PR. Here F, R and Q are n * n real matrices with R=R$^{T}$ and Q=Q$^{T}$ .geq.0 (nonnegative-definite). Such equations have been playing key roles in optimal control and filtering problems with R .geq. 0. and also in the solutions of in H$_{\infty}$ control problems with R taking the form R=H$_{1}$$^{T}$ H$_{1}$-H$_{2}$$^{T}$ H$_{2}$. In both cases an existence of stabilizing solution, i.e. the solution yielding asymptotically stable closed-loop system, is an important problem. First, we briefly review the typical results when R is of definite form, namely either R .geq. 0 as in LQG problems or R .leq. 0. They constitute two extrence cases of Riccati to the cases H$_{2}$=0 and H$_{1}$=0. Necessary and sufficient conditions are shown for the existence of nonnegative-definite or positive-definite stabilizing solution. Secondly, we focus our attention on more general case where R is only assumed to be symmetric, which obviously includes the case for H$_{\infty}$ control problems. Here, necessary conditions are established for the existence of nonnegative-definite or positive-definite stabilizing solutions. The results are established by employing consistently the so-called algebraic method based on an eigenvalue problem of a Hamiltonian matrix.x.ix.x.

  • PDF

GENERALIZED STATIONARY ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS

  • Yun, Jae-Heon;Kim, Sang-Wook
    • 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.

THETA SUMS OF HIGHER INDEX

  • Yang, Jae-Hyun
    • Bulletin of the Korean Mathematical Society
    • /
    • v.53 no.6
    • /
    • pp.1893-1908
    • /
    • 2016
  • In this paper, we obtain some behaviours of theta sums of higher index for the $Schr{\ddot{o}}dinger$-Weil representation of the Jacobi group associated with a positive definite symmetric real matrix of degree m.

FOCK REPRESENTATIONS OF THE NEISENBERG GROUP $H_R^(G,H)$

  • Yang, Jae-Hyun
    • Journal of the Korean Mathematical Society
    • /
    • v.34 no.2
    • /
    • pp.345-370
    • /
    • 1997
  • In this paper, we introduce the Fock representation $U^{F, M}$ of the Heisenberg group $H_R^(g, h)$ associated with a positive definite symmetric half-integral matrix $M$ of degree h and prove that $U^{F, M}$ is unitarily equivalent to the Schrodinger representation of index $M$.

  • PDF

ITERATIVE METHODS FOR LARGE-SCALE CONVEX QUADRATIC AND CONCAVE PROGRAMS

  • Oh, Se-Young
    • 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.

  • PDF

SPECTRAL ANALYSIS OF THE MGSS PRECONDITIONER FOR SINGULAR SADDLE POINT PROBLEMS

  • RAHIMIAN, MARYAM;SALKUYEH, DAVOD KHOJASTEH
    • Journal of applied mathematics & informatics
    • /
    • v.38 no.1_2
    • /
    • pp.175-187
    • /
    • 2020
  • Recently Salkuyeh and Rahimian in (Comput. Math. Appl. 74 (2017) 2940-2949) proposed a modification of the generalized shift-splitting (MGSS) method for solving singular saddle point problems. In this paper, we present the spectral analysis of the MGSS preconditioner when it is applied to precondition the singular saddle point problems with the (1, 1) block being symmetric. Some eigenvalue bounds for the spectrum of the preconditioned matrix are given. We show that all the real eigenvalues of the preconditioned matrix are in a positive interval and all nonzero eigenvalues having nonzero imaginary part are contained in an intersection of two circles.

An Efficient Implementation of the Supernodal Multifrontal Method (초마디 멀티프런탈 방법의 효율적인 구현)

  • 박찬규;박순달
    • Korean Management Science Review
    • /
    • v.19 no.2
    • /
    • pp.155-168
    • /
    • 2002
  • In this paper, some efficient implementation techniques for the multifrontal method, which can be used to compute the Cholesky factor of a symmetric positive definite matrix, are presented. In order to use the cache effect in the cache-based computer architecture, a hybrid method for factorizing a frontal matrix is considered. This hybrid method uses the column Cholesky method and the submatrix Cholesky method alternatively. Experiments show that the hybrid method speeds up the performance of the supernodal multifrontal method by 5%~10%, and it is superior to the Cholesky method in some problems with dense columns or large frontal matrices.

A GENERALIZATION OF LOCAL SYMMETRIC AND SKEW-SYMMETRIC SPLITTING ITERATION METHODS FOR GENERALIZED SADDLE POINT PROBLEMS

  • Li, Jian-Lei;Luo, Dang;Zhang, Zhi-Jiang
    • Journal of applied mathematics & informatics
    • /
    • v.29 no.5_6
    • /
    • pp.1167-1178
    • /
    • 2011
  • In this paper, we further investigate the local Hermitian and skew-Hermitian splitting (LHSS) iteration method and the modified LHSS (MLHSS) iteration method for solving generalized nonsymmetric saddle point problems with nonzero (2,2) blocks. When A is non-symmetric positive definite, the convergence conditions are obtained, which generalize some results of Jiang and Cao [M.-Q. Jiang and Y. Cao, On local Hermitian and Skew-Hermitian splitting iteration methods for generalized saddle point problems, J. Comput. Appl. Math., 2009(231): 973-982] for the generalized saddle point problems to generalized nonsymmetric saddle point problems with nonzero (2,2) blocks. Numerical experiments show the effectiveness of the iterative methods.

A Study on the Application of Conjugate Gradient Method in Nonlinear Magnetic Field Analysis by FEM. (유한요소법에 의한 비선형 자계 해석에 공액 구배법 적응 연구)

  • 임달호;신흥교
    • The Transactions of the Korean Institute of Electrical Engineers
    • /
    • v.39 no.1
    • /
    • pp.22-28
    • /
    • 1990
  • This paper is a study on the reduction of computation time in case of nonlinear magnetic field analysis by finite element method and Newton-Raphson method. For the purpose, the nonlinear convergence equation is computed by the conjugate gradient method which is known to be applicable to symmetric positive definite matrix equations only. As the results, we can not prove mathematically that the system Jacobian is positive definite, but when we applied this method, the diverging case did not occur. And the computation time is reduced by 25-55% and 15-45% in comparison with the case of direct and successive over-relaxation method, respectively. Therefore, we proved the utility of conjugate gradient method.

  • PDF

COMPARISONS OF PARALLEL PRECONDITIONERS FOR THE COMPUTATION OF SMALLEST GENERALIZED EIGENVALUE

  • Ma, Sang-Back;Jang, Ho-Jong;Cho, Jae-Young
    • Journal of applied mathematics & informatics
    • /
    • v.11 no.1_2
    • /
    • pp.305-316
    • /
    • 2003
  • Recently, an iterative algorithm for finding the interior eigenvalues of a definite matrix by CG-type method has been proposed. This method compares to the inverse power method. The given matrices A, and B are assumed to be large and sparse, and SPD( Symmetric Positive Definite) The CG scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for large sparse eigenproblems for smallest eigenvalue. Also, it is very amenable to parallel computations, like the CG method for the linear systems. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. But for parallel computations we need to find an efficient parallel preconditioner. Our candidates we ILU(0) in the wave-front order, ILU(0) in the multi-coloring order, Point-SSOR(Symmetric Successive Overrelaxation), and Multi-Color Block SSOR preconditioner. Wavefront order is a simple way to increase parallelism in the natural order, and Multi-coloring realizes a parallelism of order(N), where N is the order of the matrix. Another choice is the Multi-Color Block SSOR(Symmetric Successive OverRelaxation) preconditioning. Block SSOR is a symmetric preconditioner which is expected to minimize the interprocessor communication due to the blocking. We implemented the results on the CRAY-T3E with 128 nodes. The MPI (Message Passing Interface) library was adopted for the interprocessor communications. The test problem was drawn from the discretizations of partial differential equations by finite difference methods. The results show that for small number of processors Multi-Color ILU(0) has the best performance, while for large number of processors Multi-Color Block SSOR performs the best.