In this paper, a numerical solution of nonlinear partial differential equation, Benjamin-Bona-Mahony (BBM) and Cahn-Hilliard equation is presented by using Adomain Decomposition Method (ADM) and Variational Iteration ...In this paper, a numerical solution of nonlinear partial differential equation, Benjamin-Bona-Mahony (BBM) and Cahn-Hilliard equation is presented by using Adomain Decomposition Method (ADM) and Variational Iteration Method (VIM). The results reveal that the two methods are very effective, simple and very close to the exact solution.展开更多
This paper focuses on propagating perturbed two-body motion using orbital elements combined with a novel integration technique.While previous studies show that Modified Chebyshev Picard Iteration(MCPI)is a powerful to...This paper focuses on propagating perturbed two-body motion using orbital elements combined with a novel integration technique.While previous studies show that Modified Chebyshev Picard Iteration(MCPI)is a powerful tool used to propagate position and velocity,the present results show that using orbital elements to propagate the state vector reduces the number of MCPI iterations and nodes required,which is especially useful for reducing the computation time when including computationally-intensive calculations such as Spherical Harmonic gravity,and it also converges for>5.5x as many revolutions using a single segment when compared with cartesian propagation.Results for the Classical Orbital Elements and the Modified Equinoctial Orbital Elements(the latter provides singularity-free solutions)show that state propagation using these variables is inherently well-suited to the propagation method chosen.Additional benefits are achieved using a segmentation scheme,while future expansion to the two-point boundary value problem is expected to increase the domain of convergence compared with the cartesian case.MCPI is an iterative numerical method used to solve linear and nonlinear,ordinary differential equations(ODEs).It is a fusion of orthogonal Chebyshev function approximation with Picard iteration that approximates a long-arc trajectory at every iteration.Previous studies have shown that it outperforms the state of the practice numerical integrators of ODEs in a serial computing environment;since MCPI is inherently massively parallelizable,this capability is expected to increase the computational efficiency of the method presented.展开更多
Contrary to the opinion about approximation nature of a simple-iteration method, the exact solution of a system of linear algebraic equations (SLAE) in a finite number of iterations with a stationary matrix is demonst...Contrary to the opinion about approximation nature of a simple-iteration method, the exact solution of a system of linear algebraic equations (SLAE) in a finite number of iterations with a stationary matrix is demonstrated. We present a theorem and its proof that confirms the possibility to obtain the finite process and imposes the requirement for the matrix of SLAE. This matrix must be unipotent, i.e. all its eigenvalues to be equal to 1. An example of transformation of SLAE given analytically to the form with a unipotent matrix is presented. It is shown that splitting the unipotent matrix into identity and nilpotent ones results in Cramer’s analytical formulas in a finite number of iterations.展开更多
This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zer...This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.展开更多
In this work we will consider asynchronous iteration algorithms. As is well known in multiprocessor computers the parallel application of iterative methods often shows poor scaling and less optimal parallel efficiency...In this work we will consider asynchronous iteration algorithms. As is well known in multiprocessor computers the parallel application of iterative methods often shows poor scaling and less optimal parallel efficiency. The ordinary iterative asynchronous method often has much better parallel efficiency as they almost never need to wait to communicate between possessors. We will study probabilistic approach in asynchronous iteration algorithms and present a mathematical description of this computational process to the multiprocessor environment. The result of our simple numerical experiments shows a convergence and efficiency of asynchronous iterative processes for considered nonlinear problems.展开更多
We introduce a general iterative method for a finite family of generalized asymptotically quasi- nonexpansive mappings in a hyperbolic space and study its strong convergence. The new iterative method includes multi-st...We introduce a general iterative method for a finite family of generalized asymptotically quasi- nonexpansive mappings in a hyperbolic space and study its strong convergence. The new iterative method includes multi-step iterative method of Khan et al. [1] as a special case. Our results are new in hyperbolic spaces and generalize many known results in Banach spaces and CAT(0) spaces, simultaneously.展开更多
In this paper, variational iteration method and He-Laplace method are used to solve the nonlinear ordinary and partial differential equations. Laplace transformation with the homotopy perturbation method is called He-...In this paper, variational iteration method and He-Laplace method are used to solve the nonlinear ordinary and partial differential equations. Laplace transformation with the homotopy perturbation method is called He-Laplace method. A comparison is made among variational iteration method and He-Laplace. It is shown that, in He-Laplace method, the nonlinear terms of differential equation can be easily handled by the use of He’s polynomials and provides better results.展开更多
By benchmarking with the iteration of drilling technology,fracturing technology and well placement mode for shale oil and gas development in the United States and considering the geological characteristics and develop...By benchmarking with the iteration of drilling technology,fracturing technology and well placement mode for shale oil and gas development in the United States and considering the geological characteristics and development difficulties of shale oil in the Jiyang continental rift lake basin,East China,the development technology system suitable for the geological characteristics of shale oil in continental rift lake basins has been primarily formed through innovation and iteration of the development,drilling and fracturing technologies.The technology system supports the rapid growth of shale oil production and reduces the development investment cost.By comparing it with the shale oil development technology in the United States,the prospect of the shale oil development technology iteration in continental rift lake basins is proposed.It is suggested to continuously strengthen the overall three-dimensional development,improve the precision level of engineering technology,upgrade the engineering technical indicator system,accelerate the intelligent optimization of engineering equipment,explore the application of complex structure wells,form a whole-process integrated quality management system from design to implementation,and constantly innovate the concept and technology of shale oil development,so as to promote the realization of extensive,beneficial and high-quality development of shale oil in continental rift lake basins.展开更多
To reduce inter-symbol-interference (ISI) in underwater acoustic (UWA) communication systems, a method based on LDPC-QPSK joint iteration and Walsh-m composite sequence is proposed in this paper. The method is intende...To reduce inter-symbol-interference (ISI) in underwater acoustic (UWA) communication systems, a method based on LDPC-QPSK joint iteration and Walsh-m composite sequence is proposed in this paper. The method is intended for use in long-range and low signal-to-noise ratio (SNR) UWA communications. At the transmitter, Walsh-m composite sequence is introduced to resist multipath effect. At the receiver, a soft-input soft-output (SISO) module is implemented in a joint iterative process between QPSK demodulator and LDPC decoder. This method is demonstrated in three types of UWA channel models: positive, negative and invariable sound velocity gradients channels. It is shown that through contrastive simulation experiments, this method is more efficient than conventional methods based on independent decoding and demodulation. After two rounds of joint iteration, the proposed method can obtain 2.5 dB over conventional method at BER of 10-5. Numerical results verify that the proposed method is a good candidate for long-range underwater acoustic communication systems.展开更多
Minimum mean square error(MMSE) detection algorithm can achieve nearly optimal performance when the number of antennas at the base station(BS) is large enough compared to the number of users. But the traditional MMSE ...Minimum mean square error(MMSE) detection algorithm can achieve nearly optimal performance when the number of antennas at the base station(BS) is large enough compared to the number of users. But the traditional MMSE involves complicated matrix inversion. In this paper, we propose a modified MMSE algorithm which exploits the channel characteristics occurring in massive multiple-input multipleoutput(MIMO) channels and the relaxation iteration(RI) method to avoid the matrix inversion. A proper initial solution is given to accelerate the convergence speed. In addition, we point out that the channel estimation scheme used in channel hardening-exploiting message passing(CHEMP) receiver is very appropriate for our proposed detection algorithm. Simulation results verify that the proposed algorithm can achieve very close performance of the traditional MMSE algorithm with a small number of iterations.展开更多
In this paper, we investigate an accelerated version of the discrete-time Jacobi waveform relaxation iteration method. Based on the well known Chebyshev polynomial theory, we show that significant speed up can be achi...In this paper, we investigate an accelerated version of the discrete-time Jacobi waveform relaxation iteration method. Based on the well known Chebyshev polynomial theory, we show that significant speed up can be achieved by taking linear combinations of earlier iterates. The convergence and convergence speed of the new iterative method are presented and it is shown that the convergence speed of the new iterative method is sharper than that of the Jacobi method but blunter than that of the optimal SOR method. Moreover, at every iteration the new iterative method needs almost equal computation work and memory storage with the Jacobi method, and more important it can completely exploit the particular advantages of the Jacobi method in the sense of parallelism. We validate our theoretical conclusions with numerical experiments.展开更多
In this paper, Aitken’s extrapolation normally applied to convergent fixed point iteration is extended to extrapolate the solution of a divergent iteration. In addition, higher order Aitken extrapolation is introduce...In this paper, Aitken’s extrapolation normally applied to convergent fixed point iteration is extended to extrapolate the solution of a divergent iteration. In addition, higher order Aitken extrapolation is introduced that enables successive decomposition of high Eigen values of the iteration matrix to enable convergence. While extrapolation of a convergent fixed point iteration using a geometric series sum is a known form of Aitken acceleration, it is shown that in this paper, the same formula can be used to estimate the solution of sets of linear equations from diverging Gauss-Seidel iterations. In both convergent and divergent iterations, the ratios of differences among the consecutive values of iteration eventually form a convergent (divergent) series with a factor equal to the largest Eigen value of the iteration matrix. Higher order Aitken extrapolation is shown to eliminate the influence of dominant Eigen values of the iteration matrix in successive order until the iteration is determined by the lowest possible Eigen values. For the convergent part of the Gauss-Seidel iteration, further acceleration is made possible by coupling of the extrapolation technique with the successive over relaxation (SOR) method. Application examples from both convergent and divergent iterations have been provided. Coupling of the extrapolation with the SOR technique is also illustrated for a steady state two dimensional heat flow problem which was solved using MATLAB programming.展开更多
By using the asymptotic iteration method, we have calculated numerically the eigenvalues En of the hyperbolic single wave potential which is introduced by H. Bahlouli, and A. D. Alhaidari. They found a new approach (t...By using the asymptotic iteration method, we have calculated numerically the eigenvalues En of the hyperbolic single wave potential which is introduced by H. Bahlouli, and A. D. Alhaidari. They found a new approach (the “potential parameter” approach) which has been adopted for this eigenvalues problem. For a fixed energy, the problem is solvable for a set of values of the potential parameters (the “parameter spectrum”). This paper will introduce a related work to complete the goal of finding the eigenvalues, the Schr?dinger equation with hyperbolic single wave potential is solved by using asymptotic iteration method. It is found that asymptotically this method gives accurate results for arbitrary parameters, V0, γ, and λ.展开更多
This article presents the Parametric Iteration Method (PIM) for finding optimal control and its corresponding trajectory of linear systems. Without any discretization or transformation, PIM provides a sequence of func...This article presents the Parametric Iteration Method (PIM) for finding optimal control and its corresponding trajectory of linear systems. Without any discretization or transformation, PIM provides a sequence of functions which converges to the exact solution of problem. Our emphasis will be on an auxiliary parameter which directly affects on the rate of convergence. Comparison of PIM and the Variational Iteration Method (VIM) is given to show the preference of PIM over VIM. Numerical results are given for several test examples to demonstrate the applicability and efficiency of the method.展开更多
The purpose of initial orbit determination,especially in the case of angles-only data for observation,is to obtain an initial estimate that is close enough to the true orbit to enable subsequent precision orbit determ...The purpose of initial orbit determination,especially in the case of angles-only data for observation,is to obtain an initial estimate that is close enough to the true orbit to enable subsequent precision orbit determination processing to be successful.However,the classical angles-only initial orbit determination methods cannot deal with the observation data whose Earth-central angle is larger than 360°.In this paper,an improved double r-iteration initial orbit determination method to deal with the above case is presented to monitor geosynchronous Earth orbit objects for a spacebased surveillance system.Simulation results indicate that the improved double r-iteration method is feasible,and the accuracy of the obtained initial orbit meets the requirements of re-acquiring the object.展开更多
This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODE...This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.展开更多
For an upper bound of the spectral radius of the QHSS (quasi Hermitian and skew-Hermitian splitting) iteration matrix which can also bound the contraction factor of the QHSS iteration method,we give its minimum point ...For an upper bound of the spectral radius of the QHSS (quasi Hermitian and skew-Hermitian splitting) iteration matrix which can also bound the contraction factor of the QHSS iteration method,we give its minimum point under the conditions which guarantee that the upper bound is strictly less than one. This provides a good choice of the involved iteration parameters,so that the convergence rate of the QHSS iteration method can be significantly improved.展开更多
文摘In this paper, a numerical solution of nonlinear partial differential equation, Benjamin-Bona-Mahony (BBM) and Cahn-Hilliard equation is presented by using Adomain Decomposition Method (ADM) and Variational Iteration Method (VIM). The results reveal that the two methods are very effective, simple and very close to the exact solution.
文摘This paper focuses on propagating perturbed two-body motion using orbital elements combined with a novel integration technique.While previous studies show that Modified Chebyshev Picard Iteration(MCPI)is a powerful tool used to propagate position and velocity,the present results show that using orbital elements to propagate the state vector reduces the number of MCPI iterations and nodes required,which is especially useful for reducing the computation time when including computationally-intensive calculations such as Spherical Harmonic gravity,and it also converges for>5.5x as many revolutions using a single segment when compared with cartesian propagation.Results for the Classical Orbital Elements and the Modified Equinoctial Orbital Elements(the latter provides singularity-free solutions)show that state propagation using these variables is inherently well-suited to the propagation method chosen.Additional benefits are achieved using a segmentation scheme,while future expansion to the two-point boundary value problem is expected to increase the domain of convergence compared with the cartesian case.MCPI is an iterative numerical method used to solve linear and nonlinear,ordinary differential equations(ODEs).It is a fusion of orthogonal Chebyshev function approximation with Picard iteration that approximates a long-arc trajectory at every iteration.Previous studies have shown that it outperforms the state of the practice numerical integrators of ODEs in a serial computing environment;since MCPI is inherently massively parallelizable,this capability is expected to increase the computational efficiency of the method presented.
文摘Contrary to the opinion about approximation nature of a simple-iteration method, the exact solution of a system of linear algebraic equations (SLAE) in a finite number of iterations with a stationary matrix is demonstrated. We present a theorem and its proof that confirms the possibility to obtain the finite process and imposes the requirement for the matrix of SLAE. This matrix must be unipotent, i.e. all its eigenvalues to be equal to 1. An example of transformation of SLAE given analytically to the form with a unipotent matrix is presented. It is shown that splitting the unipotent matrix into identity and nilpotent ones results in Cramer’s analytical formulas in a finite number of iterations.
基金supported by the Development of airborne gravity gradiometer(No.2017YFC0601601)open subject of Key Laboratory of Petroleum Resources Research,Institute of Geology and Geophysics,Chinese Academy of Sciences(No.KLOR2018-8)
文摘This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.
文摘In this work we will consider asynchronous iteration algorithms. As is well known in multiprocessor computers the parallel application of iterative methods often shows poor scaling and less optimal parallel efficiency. The ordinary iterative asynchronous method often has much better parallel efficiency as they almost never need to wait to communicate between possessors. We will study probabilistic approach in asynchronous iteration algorithms and present a mathematical description of this computational process to the multiprocessor environment. The result of our simple numerical experiments shows a convergence and efficiency of asynchronous iterative processes for considered nonlinear problems.
文摘We introduce a general iterative method for a finite family of generalized asymptotically quasi- nonexpansive mappings in a hyperbolic space and study its strong convergence. The new iterative method includes multi-step iterative method of Khan et al. [1] as a special case. Our results are new in hyperbolic spaces and generalize many known results in Banach spaces and CAT(0) spaces, simultaneously.
文摘In this paper, variational iteration method and He-Laplace method are used to solve the nonlinear ordinary and partial differential equations. Laplace transformation with the homotopy perturbation method is called He-Laplace method. A comparison is made among variational iteration method and He-Laplace. It is shown that, in He-Laplace method, the nonlinear terms of differential equation can be easily handled by the use of He’s polynomials and provides better results.
基金Supported by the Strategic Research and Technical Consultation Project of Sinopec Science and Technology CommissionSinopec Major Science and Technology Project(P22037)。
文摘By benchmarking with the iteration of drilling technology,fracturing technology and well placement mode for shale oil and gas development in the United States and considering the geological characteristics and development difficulties of shale oil in the Jiyang continental rift lake basin,East China,the development technology system suitable for the geological characteristics of shale oil in continental rift lake basins has been primarily formed through innovation and iteration of the development,drilling and fracturing technologies.The technology system supports the rapid growth of shale oil production and reduces the development investment cost.By comparing it with the shale oil development technology in the United States,the prospect of the shale oil development technology iteration in continental rift lake basins is proposed.It is suggested to continuously strengthen the overall three-dimensional development,improve the precision level of engineering technology,upgrade the engineering technical indicator system,accelerate the intelligent optimization of engineering equipment,explore the application of complex structure wells,form a whole-process integrated quality management system from design to implementation,and constantly innovate the concept and technology of shale oil development,so as to promote the realization of extensive,beneficial and high-quality development of shale oil in continental rift lake basins.
基金Sponsored by the Fundamental Research Funds for the Central Universities(Grant No.HEUCF120814)
文摘To reduce inter-symbol-interference (ISI) in underwater acoustic (UWA) communication systems, a method based on LDPC-QPSK joint iteration and Walsh-m composite sequence is proposed in this paper. The method is intended for use in long-range and low signal-to-noise ratio (SNR) UWA communications. At the transmitter, Walsh-m composite sequence is introduced to resist multipath effect. At the receiver, a soft-input soft-output (SISO) module is implemented in a joint iterative process between QPSK demodulator and LDPC decoder. This method is demonstrated in three types of UWA channel models: positive, negative and invariable sound velocity gradients channels. It is shown that through contrastive simulation experiments, this method is more efficient than conventional methods based on independent decoding and demodulation. After two rounds of joint iteration, the proposed method can obtain 2.5 dB over conventional method at BER of 10-5. Numerical results verify that the proposed method is a good candidate for long-range underwater acoustic communication systems.
基金supported by the National Hightech R&D Program of China(2014AA01A704)the Natural Science Foundation of China(61201135)111 Project(B08038)
文摘Minimum mean square error(MMSE) detection algorithm can achieve nearly optimal performance when the number of antennas at the base station(BS) is large enough compared to the number of users. But the traditional MMSE involves complicated matrix inversion. In this paper, we propose a modified MMSE algorithm which exploits the channel characteristics occurring in massive multiple-input multipleoutput(MIMO) channels and the relaxation iteration(RI) method to avoid the matrix inversion. A proper initial solution is given to accelerate the convergence speed. In addition, we point out that the channel estimation scheme used in channel hardening-exploiting message passing(CHEMP) receiver is very appropriate for our proposed detection algorithm. Simulation results verify that the proposed algorithm can achieve very close performance of the traditional MMSE algorithm with a small number of iterations.
文摘In this paper, we investigate an accelerated version of the discrete-time Jacobi waveform relaxation iteration method. Based on the well known Chebyshev polynomial theory, we show that significant speed up can be achieved by taking linear combinations of earlier iterates. The convergence and convergence speed of the new iterative method are presented and it is shown that the convergence speed of the new iterative method is sharper than that of the Jacobi method but blunter than that of the optimal SOR method. Moreover, at every iteration the new iterative method needs almost equal computation work and memory storage with the Jacobi method, and more important it can completely exploit the particular advantages of the Jacobi method in the sense of parallelism. We validate our theoretical conclusions with numerical experiments.
文摘In this paper, Aitken’s extrapolation normally applied to convergent fixed point iteration is extended to extrapolate the solution of a divergent iteration. In addition, higher order Aitken extrapolation is introduced that enables successive decomposition of high Eigen values of the iteration matrix to enable convergence. While extrapolation of a convergent fixed point iteration using a geometric series sum is a known form of Aitken acceleration, it is shown that in this paper, the same formula can be used to estimate the solution of sets of linear equations from diverging Gauss-Seidel iterations. In both convergent and divergent iterations, the ratios of differences among the consecutive values of iteration eventually form a convergent (divergent) series with a factor equal to the largest Eigen value of the iteration matrix. Higher order Aitken extrapolation is shown to eliminate the influence of dominant Eigen values of the iteration matrix in successive order until the iteration is determined by the lowest possible Eigen values. For the convergent part of the Gauss-Seidel iteration, further acceleration is made possible by coupling of the extrapolation technique with the successive over relaxation (SOR) method. Application examples from both convergent and divergent iterations have been provided. Coupling of the extrapolation with the SOR technique is also illustrated for a steady state two dimensional heat flow problem which was solved using MATLAB programming.
文摘By using the asymptotic iteration method, we have calculated numerically the eigenvalues En of the hyperbolic single wave potential which is introduced by H. Bahlouli, and A. D. Alhaidari. They found a new approach (the “potential parameter” approach) which has been adopted for this eigenvalues problem. For a fixed energy, the problem is solvable for a set of values of the potential parameters (the “parameter spectrum”). This paper will introduce a related work to complete the goal of finding the eigenvalues, the Schr?dinger equation with hyperbolic single wave potential is solved by using asymptotic iteration method. It is found that asymptotically this method gives accurate results for arbitrary parameters, V0, γ, and λ.
文摘This article presents the Parametric Iteration Method (PIM) for finding optimal control and its corresponding trajectory of linear systems. Without any discretization or transformation, PIM provides a sequence of functions which converges to the exact solution of problem. Our emphasis will be on an auxiliary parameter which directly affects on the rate of convergence. Comparison of PIM and the Variational Iteration Method (VIM) is given to show the preference of PIM over VIM. Numerical results are given for several test examples to demonstrate the applicability and efficiency of the method.
文摘The purpose of initial orbit determination,especially in the case of angles-only data for observation,is to obtain an initial estimate that is close enough to the true orbit to enable subsequent precision orbit determination processing to be successful.However,the classical angles-only initial orbit determination methods cannot deal with the observation data whose Earth-central angle is larger than 360°.In this paper,an improved double r-iteration initial orbit determination method to deal with the above case is presented to monitor geosynchronous Earth orbit objects for a spacebased surveillance system.Simulation results indicate that the improved double r-iteration method is feasible,and the accuracy of the obtained initial orbit meets the requirements of re-acquiring the object.
文摘This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.
基金the National Natural Science Foundation (No.11671393),China.
文摘For an upper bound of the spectral radius of the QHSS (quasi Hermitian and skew-Hermitian splitting) iteration matrix which can also bound the contraction factor of the QHSS iteration method,we give its minimum point under the conditions which guarantee that the upper bound is strictly less than one. This provides a good choice of the involved iteration parameters,so that the convergence rate of the QHSS iteration method can be significantly improved.