The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identica...The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.展开更多
An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byu...An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byusing SU(2), SU(1,1) Lie algebraic method, respectively. Meanwhile, the eigenstates of the above two models are revealedto be SU(2), SU(1,1) coherent states, respectively. The relation between the usual Bogoliubov Valatin transformationand the algebraic method in a special case is also discussed.展开更多
An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byu...An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byusing SU(2), SU(1,1) Lie algebraic method, respectively. Meanwhile, the eigenstates of the above two models are revealedto be SU(2), SU(1,1) coherent states, respectively. The relation between the usual Bogoliubov Valatin transformationand the algebraic method in a special case is also discussed.展开更多
In this paper the method and technique of the diagonalization are employed to transform a vector second-order nonlinear system into two first-order approximate diagonalized systems. The existence and the asymptotic be...In this paper the method and technique of the diagonalization are employed to transform a vector second-order nonlinear system into two first-order approximate diagonalized systems. The existence and the asymptotic behavior of the solutions are obtained for a vector second-order nonlinear Robin problem of singular perturbation type.展开更多
The mathematical model of hydraulic drive unit of quadruped robot was built in this paper. According to the coupling characteristics between position control system and force control system, the decoupling control str...The mathematical model of hydraulic drive unit of quadruped robot was built in this paper. According to the coupling characteristics between position control system and force control system, the decoupling control strategy was realized based on diagonal matrix method in AMESim?. The results of simulation show that using diagonal matrix method can achieve the decoupling control effectively and it can achieve the decoupling control more effectively with the method of not offset pole-zero in the S coordinate. This research can provide theoretical basis for the application of test system of hydraulic drive unit.展开更多
Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the f...Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.展开更多
It is well known that the matrix equations play a significant role in engineering and applicable sciences. In this research article, a new modification of the homotopy perturbation method (HPM) will be proposed to obt...It is well known that the matrix equations play a significant role in engineering and applicable sciences. In this research article, a new modification of the homotopy perturbation method (HPM) will be proposed to obtain the approximated solution of the matrix equation in the form AX = B. Moreover, the conditions are deduced to check the convergence of the homotopy series. Numerical implementations are adapted to illustrate the properties of the modified method.展开更多
A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the ...A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the iteration parameters of the method are tuned in such a way that fast convergence to the value of corrector is achieved.展开更多
In this paper, a modified direct product method of scattering matrix （DPSM） was presented and the cal- culation formulawas derived as follows： φ（z）=∑n1/n!（Mz）nφ（O）and φ（z＋εj）=∑n1/n!（Mεj）nφ（z）,...In this paper, a modified direct product method of scattering matrix （DPSM） was presented and the cal- culation formulawas derived as follows： φ（z）=∑n1/n!（Mz）nφ（O）and φ（z＋εj）=∑n1/n!（Mεj）nφ（z）,whereM is the scattering matrix of which the dimension can be reduced by ＇Bethe potential method＇ drastically and therefore the calculation speed can be increased tremendously without losing accuracy very much. The results calculated with the DPSM method are in almost exact agreement with those calculated with BW method. However, the calculation speed for the modified DPSM method is approximately three times faster than that for the BW method. Furthermore, the DPSM is suitable for computing all types of ma- trices without requiring symmetry or conjugate symmetry.展开更多
In this work, the three dimensional Poisson’s equation in Cartesian coordinates with the Dirichlet’s boundary conditions in a cube is solved directly, by extending the method of Hockney. The Poisson equation is appr...In this work, the three dimensional Poisson’s equation in Cartesian coordinates with the Dirichlet’s boundary conditions in a cube is solved directly, by extending the method of Hockney. The Poisson equation is approximated by 19-points and 27-points fourth order finite difference approximation schemes and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The efficiency of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results. It is shown that 19-point formula produces comparable results with 27-point formula, though computational efforts are more in 27-point formula.展开更多
The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore...The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.展开更多
The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulat...The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,展开更多
The finite temperature Lanczos method(FTLM),which is an exact diagonalization method intensively used in quantum many-body calculations,is formulated in the framework of orthogonal polynomials and Gauss quadrature.The...The finite temperature Lanczos method(FTLM),which is an exact diagonalization method intensively used in quantum many-body calculations,is formulated in the framework of orthogonal polynomials and Gauss quadrature.The main idea is to reduce finite temperature static and dynamic quantities into weighted summations related to one-and twodimensional Gauss quadratures.Then lower order Gauss quadrature,which is generated from Lanczos iteration,can be applied to approximate the initial weighted summation.This framework fills the conceptual gap between FTLM and kernel polynomial method,and makes it easy to apply orthogonal polynomial techniques in the FTLM calculation.展开更多
Many numerical methods,such as tensor network approaches including density matrix renormalization group calculations,have been developed to calculate the extreme/ground states of quantum many-body systems.However,litt...Many numerical methods,such as tensor network approaches including density matrix renormalization group calculations,have been developed to calculate the extreme/ground states of quantum many-body systems.However,little attention has been paid to the central states,which are exponentially close to each other in terms of system size.We propose a delta-Davidson(DELDAV)method to efficiently find such interior(including the central)states in many-spin systems.The DELDAV method utilizes a delta filter in Chebyshev polynomial expansion combined with subspace diagonalization to overcome the nearly degenerate problem.Numerical experiments on Ising spin chain and spin glass shards show the correctness,efficiency,and robustness of the proposed method in finding the interior states as well as the ground states.The sought interior states may be employed to identify many-body localization phase,quantum chaos,and extremely long-time dynamical structure.展开更多
文摘The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.
文摘An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byusing SU(2), SU(1,1) Lie algebraic method, respectively. Meanwhile, the eigenstates of the above two models are revealedto be SU(2), SU(1,1) coherent states, respectively. The relation between the usual Bogoliubov Valatin transformationand the algebraic method in a special case is also discussed.
文摘An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byusing SU(2), SU(1,1) Lie algebraic method, respectively. Meanwhile, the eigenstates of the above two models are revealedto be SU(2), SU(1,1) coherent states, respectively. The relation between the usual Bogoliubov Valatin transformationand the algebraic method in a special case is also discussed.
文摘In this paper the method and technique of the diagonalization are employed to transform a vector second-order nonlinear system into two first-order approximate diagonalized systems. The existence and the asymptotic behavior of the solutions are obtained for a vector second-order nonlinear Robin problem of singular perturbation type.
文摘The mathematical model of hydraulic drive unit of quadruped robot was built in this paper. According to the coupling characteristics between position control system and force control system, the decoupling control strategy was realized based on diagonal matrix method in AMESim?. The results of simulation show that using diagonal matrix method can achieve the decoupling control effectively and it can achieve the decoupling control more effectively with the method of not offset pole-zero in the S coordinate. This research can provide theoretical basis for the application of test system of hydraulic drive unit.
文摘Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.
文摘It is well known that the matrix equations play a significant role in engineering and applicable sciences. In this research article, a new modification of the homotopy perturbation method (HPM) will be proposed to obtain the approximated solution of the matrix equation in the form AX = B. Moreover, the conditions are deduced to check the convergence of the homotopy series. Numerical implementations are adapted to illustrate the properties of the modified method.
文摘A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the iteration parameters of the method are tuned in such a way that fast convergence to the value of corrector is achieved.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11274263 and 11274264)
文摘In this paper, a modified direct product method of scattering matrix （DPSM） was presented and the cal- culation formulawas derived as follows： φ（z）=∑n1/n!（Mz）nφ（O）and φ（z＋εj）=∑n1/n!（Mεj）nφ（z）,whereM is the scattering matrix of which the dimension can be reduced by ＇Bethe potential method＇ drastically and therefore the calculation speed can be increased tremendously without losing accuracy very much. The results calculated with the DPSM method are in almost exact agreement with those calculated with BW method. However, the calculation speed for the modified DPSM method is approximately three times faster than that for the BW method. Furthermore, the DPSM is suitable for computing all types of ma- trices without requiring symmetry or conjugate symmetry.
文摘In this work, the three dimensional Poisson’s equation in Cartesian coordinates with the Dirichlet’s boundary conditions in a cube is solved directly, by extending the method of Hockney. The Poisson equation is approximated by 19-points and 27-points fourth order finite difference approximation schemes and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The efficiency of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results. It is shown that 19-point formula produces comparable results with 27-point formula, though computational efforts are more in 27-point formula.
基金supported by the National Natural Science Foundation of China(Nos. 41130418 and 40925013)the National Basic Research Program(973 Program)(No.2009CB219403)
文摘The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.
文摘The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,
基金supported by the National Natural Science Foundation of China(Grant Nos.11734002 and U1930402)。
文摘The finite temperature Lanczos method(FTLM),which is an exact diagonalization method intensively used in quantum many-body calculations,is formulated in the framework of orthogonal polynomials and Gauss quadrature.The main idea is to reduce finite temperature static and dynamic quantities into weighted summations related to one-and twodimensional Gauss quadratures.Then lower order Gauss quadrature,which is generated from Lanczos iteration,can be applied to approximate the initial weighted summation.This framework fills the conceptual gap between FTLM and kernel polynomial method,and makes it easy to apply orthogonal polynomial techniques in the FTLM calculation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.91836101,U1930201,and 11574239).
文摘Many numerical methods,such as tensor network approaches including density matrix renormalization group calculations,have been developed to calculate the extreme/ground states of quantum many-body systems.However,little attention has been paid to the central states,which are exponentially close to each other in terms of system size.We propose a delta-Davidson(DELDAV)method to efficiently find such interior(including the central)states in many-spin systems.The DELDAV method utilizes a delta filter in Chebyshev polynomial expansion combined with subspace diagonalization to overcome the nearly degenerate problem.Numerical experiments on Ising spin chain and spin glass shards show the correctness,efficiency,and robustness of the proposed method in finding the interior states as well as the ground states.The sought interior states may be employed to identify many-body localization phase,quantum chaos,and extremely long-time dynamical structure.
