This communique is opted to study the approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices.We choose the geodesic distance betweenAHXXA an...This communique is opted to study the approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices.We choose the geodesic distance betweenAHXXA and P as the cost function,and put forward the Extended Hamiltonian algorithm(EHA)and Natural gradient algorithm(NGA)for the solution.Finally,several numerical experiments give you an idea about the effectiveness of the proposed algorithms.We also show the comparison between these two algorithms EHA and NGA.Obtained results are provided and analyzed graphically.We also conclude that the extended Hamiltonian algorithm has better convergence speed than the natural gradient algorithm,whereas the trajectory of the solution matrix is optimal in case of Natural gradient algorithm(NGA)as compared to Extended Hamiltonian Algorithm(EHA).The aim of this paper is to show that the Extended Hamiltonian algorithm(EHA)has superior convergence properties as compared to Natural gradient algorithm(NGA).Upto the best of author’s knowledge,no approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices is found so far in the literature.展开更多
In this paper, we focus on the construction of structure preserving algorithms for Birkhoffian systems, based on existing symplectic schemes for the Hamiltonian equations. The key of the method is to seek an invertibl...In this paper, we focus on the construction of structure preserving algorithms for Birkhoffian systems, based on existing symplectic schemes for the Hamiltonian equations. The key of the method is to seek an invertible transformation which drives the Birkhoffian equations reduce to the Hamiltonian equations. When there exists such a transformation, applying the corresponding inverse map to symplectic discretization of the Hamiltonian equations, then resulting difference schemes are verified to be Birkhoftian symplectic for the original Birkhoffian equations. To illustrate the operation process of the method, we construct several desirable algorithms for the linear damped oscillator and the single pendulum with linear dissipation respectively. All of them exhibit excellent numerical behavior, especially in preserving conserved quantities.展开更多
An efficient and stable structure preserving algorithm, which is a variant of the QR like (SR) algorithm due to Bunse-Gerstner and Mehrmann, is presented for computing the eigenvalues and stable invariant subspaces of...An efficient and stable structure preserving algorithm, which is a variant of the QR like (SR) algorithm due to Bunse-Gerstner and Mehrmann, is presented for computing the eigenvalues and stable invariant subspaces of a Hamiltonian matrix. In the algorithm two strategies are employed, one of which is called dis-unstabilization technique and the other is preprocessing technique. Together with them, a so-called ratio-reduction equation and a backtrack technique are introduced to avoid the instability and breakdown in the original algorithm. It is shown that the new algorithm can overcome the instability and breakdown at low cost. Numerical results have demonstrated that the algorithm is stable and can compute the eigenvalues to very high accuracy.展开更多
The 3-dimensional couple equations of magneto-electro-elastic structures are derived under Hamiltonian system based on the Hamilton principle. The problem of single sort of variables is converted into the problem of d...The 3-dimensional couple equations of magneto-electro-elastic structures are derived under Hamiltonian system based on the Hamilton principle. The problem of single sort of variables is converted into the problem of double sorts of variables, and the Hamilton canonical equations are established. The 3-dimensional problem of magneto-electro-elastic structure which is investigated in Euclidean space commonly is converted into symplectic system. At the same time the Lagrange system is converted into Hamiltonian system. As an example, the dynamic characteristics of the simply supported functionally graded magneto-electro-elastic material （FGMM） plate and pipe are investigated. Finally, the problem is solved by symplectic algorithm. The results show that the physical quantities of displacement, electric potential and magnetic potential etc. change continuously at the interfaces between layers under the transverse pressure while some other physical quantities such as the stress, electric and magnetic displacement are not continuous. The dynamic stiffness is increased by the piezoelectric effect while decreased by the piezomagnetic effect.展开更多
We present a quantum adiabatic algorithm for a set of quantum 2-satisfiability(Q2SAT)problem,which is a generalization of 2-satisfiability(2SAT)problem.For a Q2SAT problem,we construct the Hamiltonian which is similar...We present a quantum adiabatic algorithm for a set of quantum 2-satisfiability(Q2SAT)problem,which is a generalization of 2-satisfiability(2SAT)problem.For a Q2SAT problem,we construct the Hamiltonian which is similar to that of a Heisenberg chain.All the solutions of the given Q2SAT problem span the subspace of the degenerate ground states.The Hamiltonian is adiabatically evolved so that the system stays in the degenerate subspace.Our numerical results suggest that the time complexity of our algorithm is O(n^(3.9))for yielding non-trivial solutions for problems with the number of clauses m=dn(n-1)/2(d■0.1).We discuss the advantages of our algorithm over the known quantum and classical algorithms.展开更多
In this paper we give a formula for enumerating the equivalent classes of orderly labeled Hamiltonian graphs under group D. and two algorithms for constructing these equivalent classes and all nonisomorphic Hamiltonia...In this paper we give a formula for enumerating the equivalent classes of orderly labeled Hamiltonian graphs under group D. and two algorithms for constructing these equivalent classes and all nonisomorphic Hamiltonian graphs. Some computational results obtained by microcomputers are listed.展开更多
This paper develops a new approach to construct variational integrators. A simplified unconventional Hamilton＇s variational principle corresponding to initial value problems is proposed, which is convenient for appli...This paper develops a new approach to construct variational integrators. A simplified unconventional Hamilton＇s variational principle corresponding to initial value problems is proposed, which is convenient for applications. The displacement and mo- mentum are approximated with the same Lagrange interpolation. After the numerical integration and variational operation, the original problems are expressed as algebraic equations with the displacement and momentum at the interpolation points as unknown variables. Some particular variational integrators are derived. An optimal scheme of choosing initial values for the Newton-Raphson method is presented for the nonlinear dynamic system. In addition, specific examples show that the proposed integrators are symplectic when the interpolation point coincides with the numerical integration point, and both are Gaussian quadrature points. Meanwhile, compared with the same order symplectic Runge-Kutta methods, although the accuracy of the two methods is almost the same, the proposed integrators are much simpler and less computationally expensive.展开更多
In this paper, the dissipative and the forced terms of the Duffing equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, b...In this paper, the dissipative and the forced terms of the Duffing equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, based on the gradient- Hamiltonian decomposition theory of vector fields, by using splitting methods, this paper constructs structure-preserving algorithms （SPAs） for the Duffing equation. Then, according to the Liouville formula, it proves that the Jacobian matrix determinants of the SPAs are equal to that of the exact flow of the Duffing equation. However, considering the explicit Runge Kutta methods, this paper finds that there is an error term of order p＋l for the Jacobian matrix determinants. The volume evolution law of a given region in phase space is discussed for different algorithms, respectively. As a result, the sum of Lyapunov exponents is exactly invariable for the SPAs proposed in this paper. Finally, through numerical experiments, relative norm errors and absolute energy errors of phase trajectories of the SPAs and the Heun method （a second-order Runge-Kutta method） are compared. Computational results illustrate that the SPAs are evidently better than the Heun method when e is small or equal to zero.展开更多
We study the statistical properties of an ensemble of disordered 1D spatial spin-chains (SSCs) of certain length in the external field. On nodes of spin-chain lattice the recurrent equations and corresponding inequal-...We study the statistical properties of an ensemble of disordered 1D spatial spin-chains (SSCs) of certain length in the external field. On nodes of spin-chain lattice the recurrent equations and corresponding inequal-ity conditions are obtained for calculation of local minimum of a classical Hamiltonian. Using these equa-tions for simulation of a model of 1D spin-glass an original high-performance parallel algorithm is developed. Distributions of different parameters of unperturbed spin-glass are calculated. It is analytically proved and shown by numerical calculations that the distribution of the spin-spin interaction constant in the Heisenberg nearest-neighboring Hamiltonian model as opposed to the widely used Gauss-Edwards-Anderson distribu-tion satisfies the Lévy alpha-stable distribution law which does not have variance. We have studied critical properties of spin-glass depending on the external field amplitude and have shown that even at weak external fields in the system strong frustrations arise. It is shown that frustrations have a fractal character, they are self-similar and do not disappear at decreasing of calculations area scale. After averaging over the fractal structures the mean values of polarizations of the spin-glass on the scales of external field's space-time peri-ods are obtained. Similarly, Edwards-Anderson’s ordering parameter depending on the external field ampli-tude is calculated. It is shown that the mean values of polarizations and the ordering parameter depending on the external field demonstrate phase transitions of first-order.展开更多
文摘This communique is opted to study the approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices.We choose the geodesic distance betweenAHXXA and P as the cost function,and put forward the Extended Hamiltonian algorithm(EHA)and Natural gradient algorithm(NGA)for the solution.Finally,several numerical experiments give you an idea about the effectiveness of the proposed algorithms.We also show the comparison between these two algorithms EHA and NGA.Obtained results are provided and analyzed graphically.We also conclude that the extended Hamiltonian algorithm has better convergence speed than the natural gradient algorithm,whereas the trajectory of the solution matrix is optimal in case of Natural gradient algorithm(NGA)as compared to Extended Hamiltonian Algorithm(EHA).The aim of this paper is to show that the Extended Hamiltonian algorithm(EHA)has superior convergence properties as compared to Natural gradient algorithm(NGA).Upto the best of author’s knowledge,no approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices is found so far in the literature.
基金supported by the National Natural Science Foundation of China(Grant No.11272050)the Excellent Young Teachers Program of North China University of Technology(Grant No.XN132)the Construction Plan for Innovative Research Team of North China University of Technology(Grant No.XN129)
文摘In this paper, we focus on the construction of structure preserving algorithms for Birkhoffian systems, based on existing symplectic schemes for the Hamiltonian equations. The key of the method is to seek an invertible transformation which drives the Birkhoffian equations reduce to the Hamiltonian equations. When there exists such a transformation, applying the corresponding inverse map to symplectic discretization of the Hamiltonian equations, then resulting difference schemes are verified to be Birkhoftian symplectic for the original Birkhoffian equations. To illustrate the operation process of the method, we construct several desirable algorithms for the linear damped oscillator and the single pendulum with linear dissipation respectively. All of them exhibit excellent numerical behavior, especially in preserving conserved quantities.
文摘An efficient and stable structure preserving algorithm, which is a variant of the QR like (SR) algorithm due to Bunse-Gerstner and Mehrmann, is presented for computing the eigenvalues and stable invariant subspaces of a Hamiltonian matrix. In the algorithm two strategies are employed, one of which is called dis-unstabilization technique and the other is preprocessing technique. Together with them, a so-called ratio-reduction equation and a backtrack technique are introduced to avoid the instability and breakdown in the original algorithm. It is shown that the new algorithm can overcome the instability and breakdown at low cost. Numerical results have demonstrated that the algorithm is stable and can compute the eigenvalues to very high accuracy.
文摘The 3-dimensional couple equations of magneto-electro-elastic structures are derived under Hamiltonian system based on the Hamilton principle. The problem of single sort of variables is converted into the problem of double sorts of variables, and the Hamilton canonical equations are established. The 3-dimensional problem of magneto-electro-elastic structure which is investigated in Euclidean space commonly is converted into symplectic system. At the same time the Lagrange system is converted into Hamiltonian system. As an example, the dynamic characteristics of the simply supported functionally graded magneto-electro-elastic material （FGMM） plate and pipe are investigated. Finally, the problem is solved by symplectic algorithm. The results show that the physical quantities of displacement, electric potential and magnetic potential etc. change continuously at the interfaces between layers under the transverse pressure while some other physical quantities such as the stress, electric and magnetic displacement are not continuous. The dynamic stiffness is increased by the piezoelectric effect while decreased by the piezomagnetic effect.
基金Project supported by the National Key R&D Program of China(Grant Nos.2017YFA0303302 and 2018YFA0305602)the National Natural Science Foundation of China(Grant No.11921005)Shanghai Municipal Science and Technology Major Project,China(Grant No.2019SHZDZX01)。
文摘We present a quantum adiabatic algorithm for a set of quantum 2-satisfiability(Q2SAT)problem,which is a generalization of 2-satisfiability(2SAT)problem.For a Q2SAT problem,we construct the Hamiltonian which is similar to that of a Heisenberg chain.All the solutions of the given Q2SAT problem span the subspace of the degenerate ground states.The Hamiltonian is adiabatically evolved so that the system stays in the degenerate subspace.Our numerical results suggest that the time complexity of our algorithm is O(n^(3.9))for yielding non-trivial solutions for problems with the number of clauses m=dn(n-1)/2(d■0.1).We discuss the advantages of our algorithm over the known quantum and classical algorithms.
基金Project supported by National Natural Foundation of China
文摘In this paper we give a formula for enumerating the equivalent classes of orderly labeled Hamiltonian graphs under group D. and two algorithms for constructing these equivalent classes and all nonisomorphic Hamiltonian graphs. Some computational results obtained by microcomputers are listed.
基金Project supported by the National Natural Science Foundation of China(Nos.11172334 and11202247)the Fundamental Research Funds for the Central Universities(No.2013390003161292)
文摘This paper develops a new approach to construct variational integrators. A simplified unconventional Hamilton＇s variational principle corresponding to initial value problems is proposed, which is convenient for applications. The displacement and mo- mentum are approximated with the same Lagrange interpolation. After the numerical integration and variational operation, the original problems are expressed as algebraic equations with the displacement and momentum at the interpolation points as unknown variables. Some particular variational integrators are derived. An optimal scheme of choosing initial values for the Newton-Raphson method is presented for the nonlinear dynamic system. In addition, specific examples show that the proposed integrators are symplectic when the interpolation point coincides with the numerical integration point, and both are Gaussian quadrature points. Meanwhile, compared with the same order symplectic Runge-Kutta methods, although the accuracy of the two methods is almost the same, the proposed integrators are much simpler and less computationally expensive.
基金Project supported by the National Natural Science Foundation of China (Grant No 10572021)the Doctoral Programme Foundation of Institute of Higher Education of China (Grant No 20040007022)
文摘In this paper, the dissipative and the forced terms of the Duffing equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, based on the gradient- Hamiltonian decomposition theory of vector fields, by using splitting methods, this paper constructs structure-preserving algorithms （SPAs） for the Duffing equation. Then, according to the Liouville formula, it proves that the Jacobian matrix determinants of the SPAs are equal to that of the exact flow of the Duffing equation. However, considering the explicit Runge Kutta methods, this paper finds that there is an error term of order p＋l for the Jacobian matrix determinants. The volume evolution law of a given region in phase space is discussed for different algorithms, respectively. As a result, the sum of Lyapunov exponents is exactly invariable for the SPAs proposed in this paper. Finally, through numerical experiments, relative norm errors and absolute energy errors of phase trajectories of the SPAs and the Heun method （a second-order Runge-Kutta method） are compared. Computational results illustrate that the SPAs are evidently better than the Heun method when e is small or equal to zero.
文摘We study the statistical properties of an ensemble of disordered 1D spatial spin-chains (SSCs) of certain length in the external field. On nodes of spin-chain lattice the recurrent equations and corresponding inequal-ity conditions are obtained for calculation of local minimum of a classical Hamiltonian. Using these equa-tions for simulation of a model of 1D spin-glass an original high-performance parallel algorithm is developed. Distributions of different parameters of unperturbed spin-glass are calculated. It is analytically proved and shown by numerical calculations that the distribution of the spin-spin interaction constant in the Heisenberg nearest-neighboring Hamiltonian model as opposed to the widely used Gauss-Edwards-Anderson distribu-tion satisfies the Lévy alpha-stable distribution law which does not have variance. We have studied critical properties of spin-glass depending on the external field amplitude and have shown that even at weak external fields in the system strong frustrations arise. It is shown that frustrations have a fractal character, they are self-similar and do not disappear at decreasing of calculations area scale. After averaging over the fractal structures the mean values of polarizations of the spin-glass on the scales of external field's space-time peri-ods are obtained. Similarly, Edwards-Anderson’s ordering parameter depending on the external field ampli-tude is calculated. It is shown that the mean values of polarizations and the ordering parameter depending on the external field demonstrate phase transitions of first-order.
