Now showing 1 - 5 of 5
  • Some of the metrics are blocked by your 
    Item type:Publication,
    Gradient iterative method with optimal convergent factor for solving a generalized sylvester matrix equation with applications to diffusion equations
    (2020-10-01)
    Boonruangkan, Nunthakarn
    ;
    We introduce a gradient iterative scheme with an optimal convergent factor for solving a generalized Sylvester matrix equation ∑<sup>p</sup>i=1<sup>A</sup> i XB<inf>i</inf> = F, where A<inf>i</inf>, B<inf>i</inf> and F are conformable rectangular matrices. The iterative scheme is derived from the gradients of the squared norm-errors of the associated subsystems for the equation. The convergence analysis reveals that the sequence of approximated solutions converge to the exact solution for any initial value if and only if the convergent factor is chosen properly in terms of the spectral radius of the associated iteration matrix. We also discuss the convergent rate and error estimations. Moreover, we determine the fastest convergent factor so that the associated iteration matrix has the smallest spectral radius. Furthermore, we provide numerical examples to illustrate the capability and efficiency of this method. Finally, we apply the proposed scheme to discretized equations for boundary value problems involving convection and diffusion.
  • Some of the metrics are blocked by your 
    Item type:Publication,
    Modified jacobi-gradient iterative method for generalized sylvester matrix equation
    (2020-11-01)
    Sasaki, Nopparut
    ;
    We propose a new iterative method for solving a generalized Sylvester matrix equation A<inf>1</inf> XA<inf>2</inf> + A<inf>3</inf> XA<inf>4</inf> = E with given square matrices A<inf>1</inf>, A<inf>2</inf>, A<inf>3</inf>, A<inf>4</inf> and an unknown rectangular matrix X. The method aims to construct a sequence of approximated solutions converging to the exact solution, no matter the initial value is. We decompose the coefficient matrices to be the sum of its diagonal part and others. The recursive formula for the iteration is derived from the gradients of quadratic norm-error functions, together with the hierarchical identification principle. We find equivalent conditions on a convergent factor, relied on eigenvalues of the associated iteration matrix, so that the method is applicable as desired. The convergence rate and error estimation of the method are governed by the spectral norm of the related iteration matrix. Furthermore, we illustrate numerical examples of the proposed method to show its capability and efficacy, compared to recent gradient-based iterative methods.
  • Some of the metrics are blocked by your 
    Item type:Publication,
    The steepest descent of gradient-based iterative method for solving rectangular linear systems with an application to Poisson’s equation
    (2020-12-01)
    Kittisopaporn, Adisorn
    ;
    We introduce an effective iterative method for solving rectangular linear systems, based on gradients along with the steepest descent optimization. We show that the proposed method is applicable with any initial vectors as long as the coefficient matrix is of full column rank. Convergence analysis produces error estimates and the asymptotic convergence rate of the algorithm, which is governed by the term 1−κ−2, where κ is the condition number of the coefficient matrix. Moreover, we apply the proposed method to a sparse linear system arising from a discretization of the one-dimensional Poisson equation. Numerical simulations illustrate the capability and effectiveness of the proposed method in comparison to the well-known and recent methods.
  • Some of the metrics are blocked by your 
    Item type:Publication,
    Convergence analysis of gradient-based iterative algorithms for a class of rectangular Sylvester matrix equations based on Banach contraction principle
    (2021-12-01)
    Kittisopaporn, Adisorn
    ;
    ;
    Lewkeeratiyutkul, Wicharn
    We derive an iterative procedure for solving a generalized Sylvester matrix equation AXB+ CXD= E, where A, B, C, D, E are conforming rectangular matrices. Our algorithm is based on gradients and hierarchical identification principle. We convert the matrix iteration process to a first-order linear difference vector equation with matrix coefficient. The Banach contraction principle reveals that the sequence of approximated solutions converges to the exact solution for any initial matrix if and only if the convergence factor belongs to an open interval. The contraction principle also gives the convergence rate and the error analysis, governed by the spectral radius of the associated iteration matrix. We obtain the fastest convergence factor so that the spectral radius of the iteration matrix is minimized. In particular, we obtain iterative algorithms for the matrix equation AXB= C, the Sylvester equation, and the Kalman–Yakubovich equation. We give numerical experiments of the proposed algorithm to illustrate its applicability, effectiveness, and efficiency.
  • Some of the metrics are blocked by your 
    Item type:Publication,
    Convergence analysis of a gradient iterative algorithm with optimal convergence factor for a generalized sylvester-transpose matrix equation
    (2021-01-01)
    Boonruangkan, Nunthakarn
    ;
    Consider a generalized Sylvester-transpose matrix equation with rectangular coefficient matrices. Based on gradients and hierarchical identification principle, we derive an iterative algorithm to produce a sequence of approximated solutions with a reasonable stopping rule concerning a relative norm-error. A convergence analysis via Banach fixed-point theorem reveals the sequence converges to a unique solution of the matrix equation for any given initial matrix if and only if the convergence factor is chosen appropriately in a certain range. The performance of algorithm is theoretically analysed through the convergence rate and error estimations. The optimal convergence factor is chosen to attain the fastest asymptotic behaviour. Finally, numerical experiments are provided to illustrate the capability and efficiency of the proposed algorithm, compared to recent gradient-based iterative algorithms.