An algorithm for numerical solution of discrete Hamilton-Jacobi-Bellman equations is proposed. The method begins with a suitable initial guess value of the solution,then finds a suitable matrix to linearize the system...An algorithm for numerical solution of discrete Hamilton-Jacobi-Bellman equations is proposed. The method begins with a suitable initial guess value of the solution,then finds a suitable matrix to linearize the system and constructs an iteration algorithm to generate the monotone sequence. The convergence of the algorithm for nonlinear discrete Hamilton-Jacobi-Bellman equations is proved. Some numerical examples are presented to confirm the effciency of this algorithm.展开更多
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.展开更多
Vector quantization (VQ) is an important data compression method. The key of the encoding of VQ is to find the closest vector among N vectors for a feature vector. Many classical linear search algorithms take O(N)...Vector quantization (VQ) is an important data compression method. The key of the encoding of VQ is to find the closest vector among N vectors for a feature vector. Many classical linear search algorithms take O(N) steps of distance computing between two vectors. The quantum VQ iteration and corresponding quantum VQ encoding algorithm that takes O(√N) steps are presented in this paper. The unitary operation of distance computing can be performed on a number of vectors simultaneously because the quantum state exists in a superposition of states. The quantum VQ iteration comprises three oracles, by contrast many quantum algorithms have only one oracle, such as Shor's factorization algorithm and Grover's algorithm. Entanglement state is generated and used, by contrast the state in Grover's algorithm is not an entanglement state. The quantum VQ iteration is a rotation over subspace, by contrast the Grover iteration is a rotation over global space. The quantum VQ iteration extends the Grover iteration to the more complex search that requires more oracles. The method of the quantum VQ iteration is universal.展开更多
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.展开更多
A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is g...A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is given, which shows that the same order of accuracy can be achieved as solving the system directly in the fine mesh when h = H2. Both theoretical analysis and numerical experiments illustrate the efficiency of the algorithm for solving the coupled problem.展开更多
The large scale linear systems with M-matrices often appear in a wide variety of areas of physical,fluid dynamics and economic sciences.It is reported in[1]that the convergence rate of the IMGS method,with the precond...The large scale linear systems with M-matrices often appear in a wide variety of areas of physical,fluid dynamics and economic sciences.It is reported in[1]that the convergence rate of the IMGS method,with the preconditioner I+S_α,is superior to that of the basic SOR iterative method for the M-matrix.This paper considers the preconditioned Jacobi(PJ)method with the preconditioner P=I+S_α+S_β,and proves theoretically that the convergence rate of the PJ method is better than that of the basic AOR method.Numerical examples are provided to illustrate the main results obtained.展开更多
This paper presents a new nonstationary iterative method for solving non linear algebraic equations that does not require the use of any derivative. The study uses only the Newton’s divided differences of first and s...This paper presents a new nonstationary iterative method for solving non linear algebraic equations that does not require the use of any derivative. The study uses only the Newton’s divided differences of first and second orders instead of the derivatives of (1).展开更多
Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are increasingly used in modeling practical super diffusive problems in fluid flow, finance and others areas of...Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are increasingly used in modeling practical super diffusive problems in fluid flow, finance and others areas of application. This paper presents the analytical solutions of the space fractional diffusion equations by variational iteration method (VIM). By using initial conditions, the explicit solutions of the equations have been presented in the closed form. Two examples, the first one is one-dimensional and the second one is two-dimensional fractional diffusion equation, are presented to show the application of the present techniques. The present method performs extremely well in terms of efficiency and simplicity.展开更多
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.展开更多
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.展开更多
In this paper, we consider two extended model equations for shallow water waves. We use He’s variational iteration method (VIM) to solve them. It is proved that this method is a very good tool for shallow water wave ...In this paper, we consider two extended model equations for shallow water waves. We use He’s variational iteration method (VIM) to solve them. It is proved that this method is a very good tool for shallow water wave equations and the obtained solutions are shown graphically.展开更多
In this paper, a user friendly algorithm based on the variational iteration method (VIM) is proposed to solve singular integral equations with generalized Abel’s kernel. It is observed that an approximate solutions y...In this paper, a user friendly algorithm based on the variational iteration method (VIM) is proposed to solve singular integral equations with generalized Abel’s kernel. It is observed that an approximate solutions yn(x) converges to the exact solution irrespective of the initial choice y0 (x). Illustrative numerical examples are given to demonstrate the efficiency and simplicity of the method in solving these types of singular integral equations.展开更多
This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones a...This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.展开更多
文摘An algorithm for numerical solution of discrete Hamilton-Jacobi-Bellman equations is proposed. The method begins with a suitable initial guess value of the solution,then finds a suitable matrix to linearize the system and constructs an iteration algorithm to generate the monotone sequence. The convergence of the algorithm for nonlinear discrete Hamilton-Jacobi-Bellman equations is proved. Some numerical examples are presented to confirm the effciency of this algorithm.
文摘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.
文摘Vector quantization (VQ) is an important data compression method. The key of the encoding of VQ is to find the closest vector among N vectors for a feature vector. Many classical linear search algorithms take O(N) steps of distance computing between two vectors. The quantum VQ iteration and corresponding quantum VQ encoding algorithm that takes O(√N) steps are presented in this paper. The unitary operation of distance computing can be performed on a number of vectors simultaneously because the quantum state exists in a superposition of states. The quantum VQ iteration comprises three oracles, by contrast many quantum algorithms have only one oracle, such as Shor's factorization algorithm and Grover's algorithm. Entanglement state is generated and used, by contrast the state in Grover's algorithm is not an entanglement state. The quantum VQ iteration is a rotation over subspace, by contrast the Grover iteration is a rotation over global space. The quantum VQ iteration extends the Grover iteration to the more complex search that requires more oracles. The method of the quantum VQ iteration is universal.
文摘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.
基金supported by National Foundation of Natural Science(11471092,11326231)Zhejiang Provincial Natural Science Foundation of China(LZ13A010003)
文摘A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is given, which shows that the same order of accuracy can be achieved as solving the system directly in the fine mesh when h = H2. Both theoretical analysis and numerical experiments illustrate the efficiency of the algorithm for solving the coupled problem.
基金The work was supported by Xi'an Jiaotong University Natural Science Foundation, China.
文摘The large scale linear systems with M-matrices often appear in a wide variety of areas of physical,fluid dynamics and economic sciences.It is reported in[1]that the convergence rate of the IMGS method,with the preconditioner I+S_α,is superior to that of the basic SOR iterative method for the M-matrix.This paper considers the preconditioned Jacobi(PJ)method with the preconditioner P=I+S_α+S_β,and proves theoretically that the convergence rate of the PJ method is better than that of the basic AOR method.Numerical examples are provided to illustrate the main results obtained.
文摘This paper presents a new nonstationary iterative method for solving non linear algebraic equations that does not require the use of any derivative. The study uses only the Newton’s divided differences of first and second orders instead of the derivatives of (1).
文摘Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are increasingly used in modeling practical super diffusive problems in fluid flow, finance and others areas of application. This paper presents the analytical solutions of the space fractional diffusion equations by variational iteration method (VIM). By using initial conditions, the explicit solutions of the equations have been presented in the closed form. Two examples, the first one is one-dimensional and the second one is two-dimensional fractional diffusion equation, are presented to show the application of the present techniques. The present method performs extremely well in terms of efficiency and simplicity.
文摘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.
文摘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.
文摘In this paper, we consider two extended model equations for shallow water waves. We use He’s variational iteration method (VIM) to solve them. It is proved that this method is a very good tool for shallow water wave equations and the obtained solutions are shown graphically.
文摘In this paper, a user friendly algorithm based on the variational iteration method (VIM) is proposed to solve singular integral equations with generalized Abel’s kernel. It is observed that an approximate solutions yn(x) converges to the exact solution irrespective of the initial choice y0 (x). Illustrative numerical examples are given to demonstrate the efficiency and simplicity of the method in solving these types of singular integral equations.
文摘This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.