TY - JOUR

T1 - A scalar homotopy method for solving an Over/Under-determined system of non-linear algebraic equations

AU - Liu, Chein Shan

AU - Yeih, Weichung

AU - Kuo, Chung Lun

AU - Atluri, Satya N.

N1 - Copyright:
Copyright 2010 Elsevier B.V., All rights reserved.

PY - 2009

Y1 - 2009

N2 - Iterative algorithms for solving a system of nonlinear algebraic equations (NAEs): Fi(xj) = 0, i, j= 1,. . . ,n date back to the seminal work of Issac Newton. Nowadays a Newton-like algorithm is still the most popular one to solve the NAEs, due to the ease of its numerical implementation. However, this type of algorithm is sensitive to the initial guess of solution, and is expensive in terms of the computations of the Jacobian matrix ∂Fi= ∂×j and its inverse at each iterative step. In addition, the Newton-like methods restrict one to construct an iteration procedure for n-variables by using n-equations, which is not a necessary condition for the existence of a solution for underdetermined or overdetermined system of equations. In this paper, a natural system of first-order nonlinear Ordinary Differential Equations (ODEs) is derived from the given system of Nonlinear Algebraic Equations (NAEs), by introducing a scalar homotopy function gauging the total residual error of the system of equations. The iterative equations are obtained by numerically integrating the resultant ODEs, which does not need the inverse of ∂Fi=∂×j. The new method keeps the merit of homotopy method, such as the global convergence, but it does not involve the complicated computation of the inverse of the Jacobian matrix. Numerical examples given confirm that this Scalar Homotopy Method (SHM) is highly efficient to find the true solutions with residual errors being much smaller.

AB - Iterative algorithms for solving a system of nonlinear algebraic equations (NAEs): Fi(xj) = 0, i, j= 1,. . . ,n date back to the seminal work of Issac Newton. Nowadays a Newton-like algorithm is still the most popular one to solve the NAEs, due to the ease of its numerical implementation. However, this type of algorithm is sensitive to the initial guess of solution, and is expensive in terms of the computations of the Jacobian matrix ∂Fi= ∂×j and its inverse at each iterative step. In addition, the Newton-like methods restrict one to construct an iteration procedure for n-variables by using n-equations, which is not a necessary condition for the existence of a solution for underdetermined or overdetermined system of equations. In this paper, a natural system of first-order nonlinear Ordinary Differential Equations (ODEs) is derived from the given system of Nonlinear Algebraic Equations (NAEs), by introducing a scalar homotopy function gauging the total residual error of the system of equations. The iterative equations are obtained by numerically integrating the resultant ODEs, which does not need the inverse of ∂Fi=∂×j. The new method keeps the merit of homotopy method, such as the global convergence, but it does not involve the complicated computation of the inverse of the Jacobian matrix. Numerical examples given confirm that this Scalar Homotopy Method (SHM) is highly efficient to find the true solutions with residual errors being much smaller.

KW - Iterative method

KW - Nonlinear algebraic equations

KW - Ordinary differential equations

KW - Scalar homotopy method (SHM)

UR - http://www.scopus.com/inward/record.url?scp=77949806014&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:77949806014

VL - 53

SP - 47

EP - 71

JO - CMES - Computer Modeling in Engineering and Sciences

JF - CMES - Computer Modeling in Engineering and Sciences

SN - 1526-1492

IS - 1

ER -