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 -