In this paper,we propose a symmetric difference data enhancement physics-informed neural network(SDE-PINN)to study soliton solutions for discrete nonlinear lattice equations(NLEs).By considering known and unknown symm...In this paper,we propose a symmetric difference data enhancement physics-informed neural network(SDE-PINN)to study soliton solutions for discrete nonlinear lattice equations(NLEs).By considering known and unknown symmetric points,numerical simulations are conducted to one-soliton and two-soliton solutions of a discrete KdV equation,as well as a one-soliton solution of a discrete Toda lattice equation.Compared with the existing discrete deep learning approach,the numerical results reveal that within the specified spatiotemporal domain,the prediction accuracy by SDE-PINN is excellent regardless of the interior or extrapolation prediction,with a significant reduction in training time.The proposed data enhancement technique and symmetric structure development provides a new perspective for the deep learning approach to solve discrete NLEs.The newly proposed SDE-PINN can also be applied to solve continuous nonlinear equations and other discrete NLEs numerically.展开更多
In this article, the authors establish some new nonlinear difference inequalities in two independent variables, which generalize some existing results and can be used as handy tools in the study of qualitative as well...In this article, the authors establish some new nonlinear difference inequalities in two independent variables, which generalize some existing results and can be used as handy tools in the study of qualitative as well as quantitative properties of solutions of certain classes of difference equations.展开更多
The difference discrete system of Euler-beam with arbitrary supports was constructed by using the two order central difference formulas. This system is equivalent to the spring-mass-rigidrod model. By using the theory...The difference discrete system of Euler-beam with arbitrary supports was constructed by using the two order central difference formulas. This system is equivalent to the spring-mass-rigidrod model. By using the theory of oscillatory matrix, the signoscillatory property of stiffness matrices of this system was proved, and the necessary and sufficient condition for the system to be positive was obtained completely.展开更多
The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ...The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ability of the discrete forms expressing to the element functions was talked about. In discrete operator difference method, the displacements of the elements can be reproduced exactly in the discrete forms whether the displacements are conforming or not. According to this point, discrete operator difference method is a method with good performance.展开更多
A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative ...A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative position,velocity and attitude of two unmanned aerial vehicles (UAVs).The second-order divided difference filter which makes use of multidimensional interpolation formulations to approximate the nonlinear transformations could achieve more accurate estimation and faster convergence from inaccurate initial conditions than standard extended Kalman filter.The filter formulation is based on relative motion equations.The global attitude parameterization is given by quarternion,while a generalized three-dimensional attitude representation is used to define the local attitude error.Simulation results are shown to compare the performance of the second-order divided difference filter with a standard extended Kalman filter approach.展开更多
In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
In this paper, combining the idea of difference method and finite element method, we construct a difference scheme for a self-adjoint problem in conservation form. Its solution uniformly converges to that of the origi...In this paper, combining the idea of difference method and finite element method, we construct a difference scheme for a self-adjoint problem in conservation form. Its solution uniformly converges to that of the original differential equation problem with order h3.展开更多
When calculating the sampled-date representation of nonlinear systems second-order hold(SOH) assumption can be applied to improving the precision of the discretization results. This paper proposes a discretization met...When calculating the sampled-date representation of nonlinear systems second-order hold(SOH) assumption can be applied to improving the precision of the discretization results. This paper proposes a discretization method based on Taylor series and the SOH assumption for the nonlinear systems with the time delayed non-affine input. The mathematical structure of the proposed discretization method is explored. This proposed discretization method can provide a precise and finite dimensional discretization model for the nonlinear time-delayed non-affine system by keeping the truncation order of the Taylor series. The performance of the proposed discretization method is evaluated by doing the simulation using a nonlinear system with the time-delayed non-affine input.Different input signals, time-delay values and sampling periods are considered in the simulation to investigate the proposed method.The simulation results demonstrate that the proposed method is practical and easy for time-delayed nonlinear non-affine systems.The comparison between SOH assumption with first-order hold(FOH) and zero-order hold(ZOH) assumptions is given to show the advantages of the proposed method.展开更多
In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the ...In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the discretization of p.d.e. via the preconditioned conjugate gradient method. For the second-order partial differential equations with Dirichlel boundary conditions, we prove that the condition number of the preconditioned system is O(1) while the condition number of the original system is O(m 2) Here m is the number of interior gridpoints in each direction. Such condition number produces a linear convergence rale.展开更多
In this paper, we use our method to solve the extended Lotka-Volterra equation and discrete KdV equation. With the help of Maple, we obtain a number of exact solutions to the two equations including soliton solutions ...In this paper, we use our method to solve the extended Lotka-Volterra equation and discrete KdV equation. With the help of Maple, we obtain a number of exact solutions to the two equations including soliton solutions presented by hyperbolic functions of sinh and cosh, periodic solutions presented by trigonometric functions of sin and cos, and rational solutions. This method can be used to solve some other nonlinear difference-differential equations.展开更多
In order to study discrete nonconservative system,Hamilton's principle within fractional difference operators of Riemann-Liouville type is given. Discrete Lagrange equations of the nonconservative system as well a...In order to study discrete nonconservative system,Hamilton's principle within fractional difference operators of Riemann-Liouville type is given. Discrete Lagrange equations of the nonconservative system as well as the nonconservative system with dynamic constraint are established within fractional difference operators of Riemann-Liouville type from the view of time scales. Firstly,time scale calculus and fractional calculus are reviewed.Secondly,with the help of the properties of time scale calculus,discrete Lagrange equation of the nonconservative system within fractional difference operators of Riemann-Liouville type is presented. Thirdly,using the Lagrange multipliers,discrete Lagrange equation of the nonconservative system with dynamic constraint is also established.Then two special cases are discussed. Finally,two examples are devoted to illustrate the results.展开更多
Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the ...Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to展开更多
The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux ...The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.展开更多
We introduce the Euler-Lagrange cohomology to study the symplectic and multisymplectic structures and their preserving properties in finite and infinite dimensional Lagrangian systems respectively. We also explore the...We introduce the Euler-Lagrange cohomology to study the symplectic and multisymplectic structures and their preserving properties in finite and infinite dimensional Lagrangian systems respectively. We also explore their certain difference discrete counterparts in the relevant regularly discretized finite and infinite dimensional Lagrangian systems by means of the difference discrete variational principle with the difference being regarded as an entire geometric object and the noncommutative differential calculus on regular lattice. In order to show that in all these cases the symplectic and multisymplectic preserving properties do not necessarily depend on the relevant Euler-Lagrange equations, the Euler-Lagrange cohomological concepts and content in the configuration space are employed.展开更多
In this paper, we present a method to solve difference differential equation（s）. As an example, we apply this method to discrete KdV equation and Ablowitz-Ladik lattice equation. As a result, many exact solutions ar...In this paper, we present a method to solve difference differential equation（s）. As an example, we apply this method to discrete KdV equation and Ablowitz-Ladik lattice equation. As a result, many exact solutions are obtained with the help of Maple including soliton solutions presented by hyperbolic functions sinh and cosh, periodic solutions presented by sin and cos and rational solutions. This method can also be used to other nonlinear difference-differential equation（s）.展开更多
Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditi...Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.展开更多
A mathematical model was developed for layout optimization of truss structures with discrete variables subjected to dynamic stress, dynamic displacement and dynamic stability constraints. By using the quasi-static met...A mathematical model was developed for layout optimization of truss structures with discrete variables subjected to dynamic stress, dynamic displacement and dynamic stability constraints. By using the quasi-static method, the mathematical model of structure optimization under dynamic stress, dynamic displacement and dynamic stability constraints were transformed into one subjected to static stress, displacement and stability constraints. The optimization procedures include two levels, i.e., the topology optimization and the shape optimization. In each level, the comprehensive algorithm was used and the relative difference quotients of two kinds of variables were used to search the optimum solution. A comparison between the optimum results of model with stability constraints and the optimum results of model without stability constraint was given. And that shows the stability constraints have a great effect on the optimum solutions.展开更多
A definition is introduced about traveling waves of 2-1 dimension lattice difference equations. Discrete heat equation is introduced and a discussion is given for the existence of traveling waves. The theory of travel...A definition is introduced about traveling waves of 2-1 dimension lattice difference equations. Discrete heat equation is introduced and a discussion is given for the existence of traveling waves. The theory of traveling waves is extended on 2-1 dimension lattice difference equations. As an application, an example is presented to illustrate the main results.展开更多
This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and th...This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart （CR） element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.展开更多
The exploration of the memristor model in the discrete domain is a fascinating hotspot.The electromagnetic induction on neurons has also begun to be simulated by some discrete memristors.However,most of the current in...The exploration of the memristor model in the discrete domain is a fascinating hotspot.The electromagnetic induction on neurons has also begun to be simulated by some discrete memristors.However,most of the current investigations are based on the integer-order discrete memristor,and there are relatively few studies on the form of fractional order.In this paper,a new fractional-order discrete memristor model with prominent nonlinearity is constructed based on the Caputo fractional-order difference operator.Furthermore,the dynamical behaviors of the Rulkov neuron under electromagnetic radiation are simulated by introducing the proposed discrete memristor.The integer-order and fractional-order peculiarities of the system are analyzed through the bifurcation graph,the Lyapunov exponential spectrum,and the iterative graph.The results demonstrate that the fractional-order system has more abundant dynamics than the integer one,such as hyper-chaos,multi-stable and transient chaos.In addition,the complexity of the system in the fractional form is evaluated by the means of the spectral entropy complexity algorithm and consequences show that it is affected by the order of the fractional system.The feature of fractional difference lays the foundation for further research and application of the discrete memristor and the neuron map in the future.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.12071042)the Beijing Natural Science Foundation(Grant No.1202004)Promoting the Development of University Classification-Student Innovation and Entrepreneurship Training Programme(Grant No.5112410857)。
文摘In this paper,we propose a symmetric difference data enhancement physics-informed neural network(SDE-PINN)to study soliton solutions for discrete nonlinear lattice equations(NLEs).By considering known and unknown symmetric points,numerical simulations are conducted to one-soliton and two-soliton solutions of a discrete KdV equation,as well as a one-soliton solution of a discrete Toda lattice equation.Compared with the existing discrete deep learning approach,the numerical results reveal that within the specified spatiotemporal domain,the prediction accuracy by SDE-PINN is excellent regardless of the interior or extrapolation prediction,with a significant reduction in training time.The proposed data enhancement technique and symmetric structure development provides a new perspective for the deep learning approach to solve discrete NLEs.The newly proposed SDE-PINN can also be applied to solve continuous nonlinear equations and other discrete NLEs numerically.
基金a HKU Seed grant the Research Grants Council of the Hong Kong SAR(HKU7016/07P)
文摘In this article, the authors establish some new nonlinear difference inequalities in two independent variables, which generalize some existing results and can be used as handy tools in the study of qualitative as well as quantitative properties of solutions of certain classes of difference equations.
基金Project supported by the National Natural Science Foundation of China (No.60034010)
文摘The difference discrete system of Euler-beam with arbitrary supports was constructed by using the two order central difference formulas. This system is equivalent to the spring-mass-rigidrod model. By using the theory of oscillatory matrix, the signoscillatory property of stiffness matrices of this system was proved, and the necessary and sufficient condition for the system to be positive was obtained completely.
文摘The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ability of the discrete forms expressing to the element functions was talked about. In discrete operator difference method, the displacements of the elements can be reproduced exactly in the discrete forms whether the displacements are conforming or not. According to this point, discrete operator difference method is a method with good performance.
基金Sponsored by the Aerospace Technology Innovation Funding(Grant No. CASC0209)
文摘A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative position,velocity and attitude of two unmanned aerial vehicles (UAVs).The second-order divided difference filter which makes use of multidimensional interpolation formulations to approximate the nonlinear transformations could achieve more accurate estimation and faster convergence from inaccurate initial conditions than standard extended Kalman filter.The filter formulation is based on relative motion equations.The global attitude parameterization is given by quarternion,while a generalized three-dimensional attitude representation is used to define the local attitude error.Simulation results are shown to compare the performance of the second-order divided difference filter with a standard extended Kalman filter approach.
文摘In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
文摘In this paper, combining the idea of difference method and finite element method, we construct a difference scheme for a self-adjoint problem in conservation form. Its solution uniformly converges to that of the original differential equation problem with order h3.
基金supported by Jiangsu Province University Natural Science Research Project(No.13KJB510003)Jiangsu Province Research and Development Institute of Marine Resources Science and Technology Open Fund Project(No.JSIMR11B05)
文摘When calculating the sampled-date representation of nonlinear systems second-order hold(SOH) assumption can be applied to improving the precision of the discretization results. This paper proposes a discretization method based on Taylor series and the SOH assumption for the nonlinear systems with the time delayed non-affine input. The mathematical structure of the proposed discretization method is explored. This proposed discretization method can provide a precise and finite dimensional discretization model for the nonlinear time-delayed non-affine system by keeping the truncation order of the Taylor series. The performance of the proposed discretization method is evaluated by doing the simulation using a nonlinear system with the time-delayed non-affine input.Different input signals, time-delay values and sampling periods are considered in the simulation to investigate the proposed method.The simulation results demonstrate that the proposed method is practical and easy for time-delayed nonlinear non-affine systems.The comparison between SOH assumption with first-order hold(FOH) and zero-order hold(ZOH) assumptions is given to show the advantages of the proposed method.
文摘In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the discretization of p.d.e. via the preconditioned conjugate gradient method. For the second-order partial differential equations with Dirichlel boundary conditions, we prove that the condition number of the preconditioned system is O(1) while the condition number of the original system is O(m 2) Here m is the number of interior gridpoints in each direction. Such condition number produces a linear convergence rale.
文摘In this paper, we use our method to solve the extended Lotka-Volterra equation and discrete KdV equation. With the help of Maple, we obtain a number of exact solutions to the two equations including soliton solutions presented by hyperbolic functions of sinh and cosh, periodic solutions presented by trigonometric functions of sin and cos, and rational solutions. This method can be used to solve some other nonlinear difference-differential equations.
基金supported by the National Natural Science Foundation of China(Nos.11802193, 11572212,11272227)the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (18KJB130005)+1 种基金the Science Research Foundation of Suzhou University of Science and Technology(331812137)Natural Science Foundation of Suzhou University of Science and Technology
文摘In order to study discrete nonconservative system,Hamilton's principle within fractional difference operators of Riemann-Liouville type is given. Discrete Lagrange equations of the nonconservative system as well as the nonconservative system with dynamic constraint are established within fractional difference operators of Riemann-Liouville type from the view of time scales. Firstly,time scale calculus and fractional calculus are reviewed.Secondly,with the help of the properties of time scale calculus,discrete Lagrange equation of the nonconservative system within fractional difference operators of Riemann-Liouville type is presented. Thirdly,using the Lagrange multipliers,discrete Lagrange equation of the nonconservative system with dynamic constraint is also established.Then two special cases are discussed. Finally,two examples are devoted to illustrate the results.
文摘Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to
基金the National Key R&D Program of China(No.2018YFE0180900).
文摘The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.
文摘We introduce the Euler-Lagrange cohomology to study the symplectic and multisymplectic structures and their preserving properties in finite and infinite dimensional Lagrangian systems respectively. We also explore their certain difference discrete counterparts in the relevant regularly discretized finite and infinite dimensional Lagrangian systems by means of the difference discrete variational principle with the difference being regarded as an entire geometric object and the noncommutative differential calculus on regular lattice. In order to show that in all these cases the symplectic and multisymplectic preserving properties do not necessarily depend on the relevant Euler-Lagrange equations, the Euler-Lagrange cohomological concepts and content in the configuration space are employed.
基金The project supported by the State Key Basic Research Program of China under Grant No 2004CB318000
文摘In this paper, we present a method to solve difference differential equation（s）. As an example, we apply this method to discrete KdV equation and Ablowitz-Ladik lattice equation. As a result, many exact solutions are obtained with the help of Maple including soliton solutions presented by hyperbolic functions sinh and cosh, periodic solutions presented by sin and cos and rational solutions. This method can also be used to other nonlinear difference-differential equation（s）.
基金supported by National Natural Science Foundation of China(11101244,11271231)National Tackling Key Problems Program(20050200069)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.
基金Project supported by the National Natural Science Foundation of China (Nos. 10002005 and 10421002)the Natural Science Foundation of Tianjin (No.02360081)the Education Committee Foundation of Tianjin (No.20022104)the Program for Changjiang Scholars and Innovative Research Team in University of China and the 211 Foundation of Dalian University of Technology
文摘A mathematical model was developed for layout optimization of truss structures with discrete variables subjected to dynamic stress, dynamic displacement and dynamic stability constraints. By using the quasi-static method, the mathematical model of structure optimization under dynamic stress, dynamic displacement and dynamic stability constraints were transformed into one subjected to static stress, displacement and stability constraints. The optimization procedures include two levels, i.e., the topology optimization and the shape optimization. In each level, the comprehensive algorithm was used and the relative difference quotients of two kinds of variables were used to search the optimum solution. A comparison between the optimum results of model with stability constraints and the optimum results of model without stability constraint was given. And that shows the stability constraints have a great effect on the optimum solutions.
基金Supported by the National Natural Science Foundation of China（Ill61049）
文摘A definition is introduced about traveling waves of 2-1 dimension lattice difference equations. Discrete heat equation is introduced and a discussion is given for the existence of traveling waves. The theory of traveling waves is extended on 2-1 dimension lattice difference equations. As an application, an example is presented to illustrate the main results.
基金supported by the National Natural Science Foundation of China(Nos.11271273 and 11271298)
文摘This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart （CR） element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.
基金supported by the Major Research Plan of the National Natural Science Foundation of China(Grant No.91964108)the National Natural Science Foundation of China(Grant No.61971185)the Natural Science Foundation of Hunan Province,China(Grant No.2020JJ4218).
文摘The exploration of the memristor model in the discrete domain is a fascinating hotspot.The electromagnetic induction on neurons has also begun to be simulated by some discrete memristors.However,most of the current investigations are based on the integer-order discrete memristor,and there are relatively few studies on the form of fractional order.In this paper,a new fractional-order discrete memristor model with prominent nonlinearity is constructed based on the Caputo fractional-order difference operator.Furthermore,the dynamical behaviors of the Rulkov neuron under electromagnetic radiation are simulated by introducing the proposed discrete memristor.The integer-order and fractional-order peculiarities of the system are analyzed through the bifurcation graph,the Lyapunov exponential spectrum,and the iterative graph.The results demonstrate that the fractional-order system has more abundant dynamics than the integer one,such as hyper-chaos,multi-stable and transient chaos.In addition,the complexity of the system in the fractional form is evaluated by the means of the spectral entropy complexity algorithm and consequences show that it is affected by the order of the fractional system.The feature of fractional difference lays the foundation for further research and application of the discrete memristor and the neuron map in the future.
