Continuously differentiable radial basis functions (C<sup>∞</sup>-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are ...Continuously differentiable radial basis functions (C<sup>∞</sup>-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are full and can become very ill- conditioned. Similarly, the Hilbert and Vandermonde have full matrices and become ill-conditioned. The difference between a coefficient matrix generated by C<sup>∞</sup>-RBFs for partial differential or integral equations and Hilbert and Vandermonde systems is that C<sup>∞</sup>-RBFs are very sensitive to small changes in the adjustable parameters. These parameters affect the condition number and solution accuracy. The error terrain has many local and global maxima and minima. To find stable and accurate numerical solutions for full linear equation systems, this study proposes a hybrid combination of block Gaussian elimination (BGE) combined with arbitrary precision arithmetic (APA) to minimize the accumulation of rounding errors. In the future, this algorithm can execute faster using preconditioners and implemented on massively parallel computers.展开更多
A novel method based on ant colony optimization （ACO）, algorithm for solving the ill-conditioned linear systems of equations is proposed. ACO is a parallelized bionic optimization algorithm which is inspired from th...A novel method based on ant colony optimization （ACO）, algorithm for solving the ill-conditioned linear systems of equations is proposed. ACO is a parallelized bionic optimization algorithm which is inspired from the behavior of real ants. ACO algorithm is first introduced, a kind of positive feedback mechanism is adopted in ACO. Then, the solu- tion problem of linear systems of equations was reformulated as an unconstrained optimization problem for solution by an ACID algorithm. Finally, the ACID with other traditional methods is applied to solve a kind of multi-dimensional Hilbert ill-conditioned linear equations. The numerical results demonstrate that ACO is effective, robust and recommendable in solving ill-conditioned linear systems of equations.展开更多
In this paper a new ODE numerical integration method was successfully applied to solving nonlinear equations. The method is of same simplicity as fixed point iteration, but the efficiency has been significantly improv...In this paper a new ODE numerical integration method was successfully applied to solving nonlinear equations. The method is of same simplicity as fixed point iteration, but the efficiency has been significantly improved, so it is especially suitable for large scale systems. For Brown’s equations, an existing article reported that when the dimension of the equation N = 40, the subroutines they used could not give a solution, as compared with our method, we can easily solve this equation even when N = 100. Other two large equations have the dimension of N = 1000, all the existing available methods have great difficulties to handle them, however, our method proposed in this paper can deal with those tough equations without any difficulties. The sigularity and choosing initial values problems were also mentioned in this paper.展开更多
Let Fqbe the finite field,q=p^(k),with p being a prime and k being a positive integer.Let F_(q)^(*)be the multiplicative group of Fq,that is F_(q)^(*)=F_(q){0}.In this paper,by using the Jacobi sums and an analog of H...Let Fqbe the finite field,q=p^(k),with p being a prime and k being a positive integer.Let F_(q)^(*)be the multiplicative group of Fq,that is F_(q)^(*)=F_(q){0}.In this paper,by using the Jacobi sums and an analog of Hasse-Davenport theorem,an explicit formula for the number of solutions of cubic diagonal equation x_(1)^(3)+x_(2)^(3)+…+x_(n)^(3)=c over Fqis given,where c∈F_(q)^(*)and p≡1(mod 3).This extends earlier results.展开更多
This paper deals with the singular perturbation of the boundary value problem of the systems for quasi-linear ordinary differential equationswhere x,f, y , h, A, B and C all belong to Rn , and g is an n×n matrix ...This paper deals with the singular perturbation of the boundary value problem of the systems for quasi-linear ordinary differential equationswhere x,f, y , h, A, B and C all belong to Rn , and g is an n×n matrix function. Under suitable conditions we prove the existence of the solutions by diagonalization and the fixed point theorem and also estimate the remainder.展开更多
We present the existence of solution for a coupled system of fractional integro-differential equations. The differential operator is taken in the Caputo fractional sense. We combine the diagonalization method with Arz...We present the existence of solution for a coupled system of fractional integro-differential equations. The differential operator is taken in the Caputo fractional sense. We combine the diagonalization method with Arzela-Ascoli theorem to show a fixed point theorem of Schauder.展开更多
In Kronecker products works, matrices are sometimes regarded as vectors and vectors are sometimes made in to matrices. To be precise about these reshaping we use the vector and diagonal extraction operators. In the pr...In Kronecker products works, matrices are sometimes regarded as vectors and vectors are sometimes made in to matrices. To be precise about these reshaping we use the vector and diagonal extraction operators. In the present paper, the results are organized in the following ways. First, we formulate the coupled matrix linear least-squares problem and present the efficient solutions of this problem that arises in multistatic antenna array processing problem. Second, we extend the use of connection between the Hadamard (Kronecker) product and diagonal extraction (vector) operator in order to construct a computationally-efficient solution of non-homogeneous coupled matrix differential equations that useful in various applications. Finally, the analysis indicates that the Kronecker (Khatri-Rao) structure method can achieve good efficient while the Hadamard structure method achieve more efficient when the unknown matrices are diagonal.展开更多
Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate...Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate regularized matrix, and a singular decomposition method was used to generate regularization parameters. Numerical test results suggest that the regularized ambiguity float solution is more stable and reliable than the least-squares float solution. The mean square error matrix of the new method possesses a lower correlation than the variancecovariance matrix of the least-squares estimation. The size of the ambiguity search space is reduced and the search efficiency is improved. The success rate of the integer ambiguity searching process is improved significantly when the ambiguity resolution by using constraint equation method is used to determine the correct ambiguity integervector. The ambiguity resolution by using constraint equation method requires an initial input of the ambiguity float solution candidates which are obtained from the LAMBDA method in the new method. In addition, the observation time required to fix reliable integer ambiguities can he significantly reduced.展开更多
In this paper we first give an a priori estimate of maximum modulus ofsolutions for a class of systems of diagonally degenerate elliptic equations in the case of p > 2.
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.展开更多
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.展开更多
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.展开更多
In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved dire...In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy 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.展开更多
The Krein-Rutman theorem is vital in partial differential equations that are non-linear and provides evidence of the presence of several significant eigenvalues useful in topological degree calculations, stability ana...The Krein-Rutman theorem is vital in partial differential equations that are non-linear and provides evidence of the presence of several significant eigenvalues useful in topological degree calculations, stability analysis, and bifurcation theory. Schr?der’s equation which has been used extensively in studies of turbulence is an equation with a single independent variable suitable for encoding self-similarity. The concept of Hilbert spaces has been an inner product space frequently used due to its convenience in countless dimensional vector analysis. This paper is aimed at proving a number of solutions through the Krein-Rutman theorem in unitary spaces especially in Hilbert spaces. It has been certainly observed that the whole Krein-Rutman theorem system has a fairly stable scope, and has strong regular features, and many non-linear elliptic operators need the most ethical principles to satisfy the comparison policy.展开更多
In previous works, the theoretical and experimental deterministic scalar kinematic structures, the theoretical and experimental deterministic vector kinematic structures, the theoretical and experimental deterministic...In previous works, the theoretical and experimental deterministic scalar kinematic structures, the theoretical and experimental deterministic vector kinematic structures, the theoretical and experimental deterministic scalar dynamic structures, and the theoretical and experimental deterministic vector dynamic structures have been developed to compute the exact solution for deterministic chaos of the exponential pulsons and oscillons that is governed by the nonstationary three-dimensional Navier-Stokes equations. To explore properties of the kinetic energy, rectangular, diagonal, and triangular summations of a matrix of the kinetic energy and general terms of various sums have been used in the current paper to develop quantization of the kinetic energy of deterministic chaos. Nested structures of a cumulative energy pulson, an energy pulson of propagation, an internal energy oscillon, a diagonal energy oscillon, and an external energy oscillon have been established. In turn, the energy pulsons and oscillons include group pulsons of propagation, internal group oscillons, diagonal group oscillons, and external group oscillons. Sequentially, the group pulsons and oscillons contain wave pulsons of propagation, internal wave oscillons, diagonal wave oscillons, and external wave oscillons. Consecutively, the wave pulsons and oscillons are composed of elementary pulsons of propagation, internal elementary oscillons, diagonal elementary oscillons, and external elementary oscillons. Topology, periodicity, and integral properties of the exponential pulsons and oscillons have been studied using the novel method of the inhomogeneous Fourier expansions via eigenfunctions in coordinates and time. Symbolic computations of the exact expansions have been performed using the experimental and theoretical programming in Maple. Results of the symbolic computations have been justified by probe visualizations.展开更多
文摘Continuously differentiable radial basis functions (C<sup>∞</sup>-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are full and can become very ill- conditioned. Similarly, the Hilbert and Vandermonde have full matrices and become ill-conditioned. The difference between a coefficient matrix generated by C<sup>∞</sup>-RBFs for partial differential or integral equations and Hilbert and Vandermonde systems is that C<sup>∞</sup>-RBFs are very sensitive to small changes in the adjustable parameters. These parameters affect the condition number and solution accuracy. The error terrain has many local and global maxima and minima. To find stable and accurate numerical solutions for full linear equation systems, this study proposes a hybrid combination of block Gaussian elimination (BGE) combined with arbitrary precision arithmetic (APA) to minimize the accumulation of rounding errors. In the future, this algorithm can execute faster using preconditioners and implemented on massively parallel computers.
文摘A novel method based on ant colony optimization （ACO）, algorithm for solving the ill-conditioned linear systems of equations is proposed. ACO is a parallelized bionic optimization algorithm which is inspired from the behavior of real ants. ACO algorithm is first introduced, a kind of positive feedback mechanism is adopted in ACO. Then, the solu- tion problem of linear systems of equations was reformulated as an unconstrained optimization problem for solution by an ACID algorithm. Finally, the ACID with other traditional methods is applied to solve a kind of multi-dimensional Hilbert ill-conditioned linear equations. The numerical results demonstrate that ACO is effective, robust and recommendable in solving ill-conditioned linear systems of equations.
文摘In this paper a new ODE numerical integration method was successfully applied to solving nonlinear equations. The method is of same simplicity as fixed point iteration, but the efficiency has been significantly improved, so it is especially suitable for large scale systems. For Brown’s equations, an existing article reported that when the dimension of the equation N = 40, the subroutines they used could not give a solution, as compared with our method, we can easily solve this equation even when N = 100. Other two large equations have the dimension of N = 1000, all the existing available methods have great difficulties to handle them, however, our method proposed in this paper can deal with those tough equations without any difficulties. The sigularity and choosing initial values problems were also mentioned in this paper.
基金Supported by the Natural Science Foundation of Henan Province(232300420123)the National Natural Science Foundation of China(12026224)the Research Center of Mathematics and Applied Mathematics,Nanyang Institute of Technology。
文摘Let Fqbe the finite field,q=p^(k),with p being a prime and k being a positive integer.Let F_(q)^(*)be the multiplicative group of Fq,that is F_(q)^(*)=F_(q){0}.In this paper,by using the Jacobi sums and an analog of Hasse-Davenport theorem,an explicit formula for the number of solutions of cubic diagonal equation x_(1)^(3)+x_(2)^(3)+…+x_(n)^(3)=c over Fqis given,where c∈F_(q)^(*)and p≡1(mod 3).This extends earlier results.
文摘This paper deals with the singular perturbation of the boundary value problem of the systems for quasi-linear ordinary differential equationswhere x,f, y , h, A, B and C all belong to Rn , and g is an n×n matrix function. Under suitable conditions we prove the existence of the solutions by diagonalization and the fixed point theorem and also estimate the remainder.
文摘We present the existence of solution for a coupled system of fractional integro-differential equations. The differential operator is taken in the Caputo fractional sense. We combine the diagonalization method with Arzela-Ascoli theorem to show a fixed point theorem of Schauder.
文摘In Kronecker products works, matrices are sometimes regarded as vectors and vectors are sometimes made in to matrices. To be precise about these reshaping we use the vector and diagonal extraction operators. In the present paper, the results are organized in the following ways. First, we formulate the coupled matrix linear least-squares problem and present the efficient solutions of this problem that arises in multistatic antenna array processing problem. Second, we extend the use of connection between the Hadamard (Kronecker) product and diagonal extraction (vector) operator in order to construct a computationally-efficient solution of non-homogeneous coupled matrix differential equations that useful in various applications. Finally, the analysis indicates that the Kronecker (Khatri-Rao) structure method can achieve good efficient while the Hadamard structure method achieve more efficient when the unknown matrices are diagonal.
文摘Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate regularized matrix, and a singular decomposition method was used to generate regularization parameters. Numerical test results suggest that the regularized ambiguity float solution is more stable and reliable than the least-squares float solution. The mean square error matrix of the new method possesses a lower correlation than the variancecovariance matrix of the least-squares estimation. The size of the ambiguity search space is reduced and the search efficiency is improved. The success rate of the integer ambiguity searching process is improved significantly when the ambiguity resolution by using constraint equation method is used to determine the correct ambiguity integervector. The ambiguity resolution by using constraint equation method requires an initial input of the ambiguity float solution candidates which are obtained from the LAMBDA method in the new method. In addition, the observation time required to fix reliable integer ambiguities can he significantly reduced.
文摘In this paper we first give an a priori estimate of maximum modulus ofsolutions for a class of systems of diagonally degenerate elliptic equations in the case of p > 2.
文摘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.
基金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.
文摘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.
文摘In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy 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.
文摘The Krein-Rutman theorem is vital in partial differential equations that are non-linear and provides evidence of the presence of several significant eigenvalues useful in topological degree calculations, stability analysis, and bifurcation theory. Schr?der’s equation which has been used extensively in studies of turbulence is an equation with a single independent variable suitable for encoding self-similarity. The concept of Hilbert spaces has been an inner product space frequently used due to its convenience in countless dimensional vector analysis. This paper is aimed at proving a number of solutions through the Krein-Rutman theorem in unitary spaces especially in Hilbert spaces. It has been certainly observed that the whole Krein-Rutman theorem system has a fairly stable scope, and has strong regular features, and many non-linear elliptic operators need the most ethical principles to satisfy the comparison policy.
文摘In previous works, the theoretical and experimental deterministic scalar kinematic structures, the theoretical and experimental deterministic vector kinematic structures, the theoretical and experimental deterministic scalar dynamic structures, and the theoretical and experimental deterministic vector dynamic structures have been developed to compute the exact solution for deterministic chaos of the exponential pulsons and oscillons that is governed by the nonstationary three-dimensional Navier-Stokes equations. To explore properties of the kinetic energy, rectangular, diagonal, and triangular summations of a matrix of the kinetic energy and general terms of various sums have been used in the current paper to develop quantization of the kinetic energy of deterministic chaos. Nested structures of a cumulative energy pulson, an energy pulson of propagation, an internal energy oscillon, a diagonal energy oscillon, and an external energy oscillon have been established. In turn, the energy pulsons and oscillons include group pulsons of propagation, internal group oscillons, diagonal group oscillons, and external group oscillons. Sequentially, the group pulsons and oscillons contain wave pulsons of propagation, internal wave oscillons, diagonal wave oscillons, and external wave oscillons. Consecutively, the wave pulsons and oscillons are composed of elementary pulsons of propagation, internal elementary oscillons, diagonal elementary oscillons, and external elementary oscillons. Topology, periodicity, and integral properties of the exponential pulsons and oscillons have been studied using the novel method of the inhomogeneous Fourier expansions via eigenfunctions in coordinates and time. Symbolic computations of the exact expansions have been performed using the experimental and theoretical programming in Maple. Results of the symbolic computations have been justified by probe visualizations.
