ANCIENT CHINESE ALGORITHM: THE YING BUZU SHU (METHOD OF SURPLUS AND DEFICIENCY)VS NEWTON ITERATION METHOD

Air exploratory discussion of an ancient Chinese algorithm, the Ying Buzu Shu, in about 2nd century BC, known as the rule of double false position in the West is given. In addition to pointing out that the rule of double false position is actually a translation version of the ancient Chinese algorithm, a comparison with well-known Newton iteration method is also made. If derivative is introduced, the ancient Chinese algorithm reduces to the Newton method. A modification of the ancient Chinese algorithm is also proposed, and some of applications to nonlinear oscillators are illustrated.

3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.

Efficient Fast Independent Component Analysis Algorithm with Fifth-Order Convergence

Independent component analysis （ICA） is the primary statistical method for solving the problems of blind source separation. The fast ICA is a famous and excellent algorithm and its contrast function is optimized by the quadratic convergence of Newton iteration method. In order to improve the convergence speed and the separation precision of the fast ICA, an improved fast ICA algorithm is presented. The algorithm introduces an efficient Newton＇s iterative method with fifth-order convergence for optimizing the contrast function and gives the detail derivation process and the corresponding condition. The experimental results demonstrate that the convergence speed and the separation precision of the improved algorithm are better than that of the fast ICA.

ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS

This paper discusses the accelerating of nonlinear parabolic equations. Two iterative methods for solving the implicit scheme new nonlinear iterative methods named by the implicit-explicit quasi-Newton （IEQN） method and the derivative free implicit-explicit quasi-Newton （DFIEQN） method are introduced, in which the resulting linear equations from the linearization can preserve the parabolic characteristics of the original partial differential equations. It is proved that the iterative sequence of the iteration method can converge to the solution of the implicit scheme quadratically. Moreover, compared with the Jacobian Free Newton-Krylov （JFNK） method, the DFIEQN method has some advantages, e.g., its implementation is easy, and it gives a linear algebraic system with an explicit coefficient matrix, so that the linear （inner） iteration is not restricted to the Krylov method. Computational results by the IEQN, DFIEQN, JFNK and Picard iteration methods are presented in confirmation of the theory and comparison of the performance of these methods.

Numerical Approximation of a Nonlinear 3D Heat Radiation Problem

In this paper,we are concerned with the numerical approximation of a steady-state heat radiation problem with a nonlinear Stefan-Boltzmann boundary condition in IR^(3).We first derive an equivalent minimization problem and then present a finite element analysis to the solution of such a minimization problem.Moreover,we apply the Newton iterative method for solving the nonlinear equation resulting from the minimization problem.A numerical example is given to illustrate theoretical results.

New Algorithm for Fault Superimposed Quantities Based on Superimposed Network

A new algorithm for fault superimposed quantity(FSIQ)is presented and analyzed.The network equations are built up by combining fault superimposed networks(FSIN)with the boundary conditions of FSIQ at the fault point and are solved with the Newton iterative method.The algorithm has clear physical meaning and does not require an intermediate procedure to derive FSIQ.The algorithm is implemented by computer programming,and the results of calculations show that the algorithm is fast and accurate.The method can be used not only to calculate FSIQ in the complex power systems with simple or multiple faults,but also to analyze and evaluate the performance of the protective relays and automatic devices based on FSIQ.