Starting from general Jacobi polynomials we derive for the Ul-traspherical polynomials as their special case a set of related polynomials which can be extended to an orthogonal set of functions with interesting proper...Starting from general Jacobi polynomials we derive for the Ul-traspherical polynomials as their special case a set of related polynomials which can be extended to an orthogonal set of functions with interesting properties. It leads to an alternative definition of the Ultraspherical polynomials by a fixed integral operator in application to powers of the variable u in an analogous way as it is possible for Hermite polynomials. From this follows a generating function which is apparently known only for the Legendre and Chebyshev polynomials as their special case. Furthermore, we show that the Ultraspherical polynomials form a realization of the SU(1,1) Lie algebra with lowering and raising operators which we explicitly determine. By reordering of multiplication and differentiation operators we derive new operator identities for the whole set of Jacobi polynomials which may be applied to arbitrary functions and provide then function identities. In this way we derive a new “convolution identity” for Jacobi polynomials and compare it with a known convolution identity of different structure for Gegenbauer polynomials. In short form we establish the connection of Jacobi polynomials and their related orthonormalized functions to the eigensolution of the Schrödinger equation to Pöschl-Teller potentials.展开更多
In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear...In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear equations by expanding the approximate solution in terms of shifted Legendre polynomials with unknown coefficients. The performance of presented method has been compared with other methods, namely Sinc-Galerkin, quadratic spline collocation and LiuLin method. Numerical examples show better accuracy of the proposed method. Moreover, the computation cost decreases at least by a factor of 6 in this method.展开更多
This paper is devoted to studying the approximate solution of singular integral equations by means of Chebyshev polynomials. Some examples are presented to illustrate the method.
In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermit...In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.展开更多
This work is aim at providing a numerical technique for the Volterra integral equations using Galerkin method. For this purpose, an effective matrix formulation is proposed to solve linear Volterra integral equations ...This work is aim at providing a numerical technique for the Volterra integral equations using Galerkin method. For this purpose, an effective matrix formulation is proposed to solve linear Volterra integral equations of the first and second kind respectively using orthogonal polynomials as trial functions which are constructed in the interval [-1,1] with respect to the weight function w(x)=1+x<sup>2</sup>. The efficiency of the proposed method is tested on several numerical examples and compared with the analytic solutions available in the literature.展开更多
Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some eq...Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some equivalent system,,; are given. Using them the paper shows how to state properly initial and boundary conditions for these DAEs. The existence and uniqueness theory of the solution of the initial and boundary value problems for higher index DAEs are proposed.展开更多
In this paper, we propose a new variation of the Adomian polynomials, which we call the degenerate Adomian polynomials, for the power series solutions of nonlinear ordinary differential equations with nonseparable non...In this paper, we propose a new variation of the Adomian polynomials, which we call the degenerate Adomian polynomials, for the power series solutions of nonlinear ordinary differential equations with nonseparable nonlinearities. We establish efficient algorithms for the degenerate Adomian polynomials. Next we compare the results by the Adomian decomposition method using the classic Adomian polynomials with the results by the Rach-Adomian-Meyers modified decomposition method incorporating the degenerate Adomian polynomials, which itself has been shown to be a confluence of the Adomian decomposition method and the power series method. Convergence acceleration techniques including the diagonal Pade approximants are considered, and new numeric algorithms for the multistage decomposition are deduced using the degenerate Adomian polynomials. Our new technique provides a significant advantage for automated calculations when computing the power series form of the solution for nonlinear ordinary differential equations. Several expository examples are investigated to demonstrate its reliability and efficiency.展开更多
One of the solution techniques used for ordinary differential equations, partial and integral equations is the Elzaki Transform. This paper is an extension of Mamadu and Njoseh [1] numerical procedure (Elzaki transfor...One of the solution techniques used for ordinary differential equations, partial and integral equations is the Elzaki Transform. This paper is an extension of Mamadu and Njoseh [1] numerical procedure (Elzaki transform method (ETM)) for computing delay differential equations (DDEs). Here, a reconstructed Elzaki transform method (RETM) is proposed for the solution of DDEs where Mamadu-Njoseh polynomials are applied as basis functions in the approximation of the analytic solution. Using this strategy, a numerical illustration as in Ref.[1] is provided to the RETM as a basis for comparison to guarantee accuracy and consistency of the method. All numerical computations were performed with MAPLE 18 software.展开更多
The second order Euler-Lagrange equations are transformed to a set of first order differential/algebraic equations, which are then transformed to state equations by using local parameterization. The corresponding disc...The second order Euler-Lagrange equations are transformed to a set of first order differential/algebraic equations, which are then transformed to state equations by using local parameterization. The corresponding discretization method is presented, and the results can be used to implementation of various numerical integration methods. A numerical example is presented finally.展开更多
This research aims to solve Differential Algebraic Equation (DAE) problems in their original form, wherein both the differential and algebraic equations remain. The Newton or Newton-Broyden technique along with some i...This research aims to solve Differential Algebraic Equation (DAE) problems in their original form, wherein both the differential and algebraic equations remain. The Newton or Newton-Broyden technique along with some integrators such as the Runge-Kutta method is coupled together to solve the problems. Experiments show that the method developed in this paper is efficient, as it demonstrates that implementation of the method is not difficult, and such method is able to provide approximate solutions with ease within some desired accuracy standards.展开更多
In this paper,the Bernstein polynomials method is proposed for the numerical solution of Fredholm integro-differential-difference equation with variable coefficients and mixed conditions.This method is using a simple ...In this paper,the Bernstein polynomials method is proposed for the numerical solution of Fredholm integro-differential-difference equation with variable coefficients and mixed conditions.This method is using a simple computational manner to obtain a quite acceptable approximate solution.The main characteristic behind this method lies in the fact that,on the one hand,the problem will be reduced to a system of algebraic equations.On the other hand,the efficiency and accuracy of the Bernstein polynomials method for solving these equations are high.The existence and uniqueness of the solution have been proved.Moreover,an estimation of the error bound for this method will be shown by preparing some theorems.Finally,some numerical experiments are presented to show the excellent behavior and high accuracy of this algorithm in comparison with some other well-known methods.展开更多
How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linea...How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).展开更多
In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular inde...In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.展开更多
Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As...Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As far as we known, the nonnegativity of solutions of the nonlinear fractional differential-algebraic equations is still not studied. In this article, we investigate the nonnegativity of solutions of the equations. Firstly, we discuss the existence of nonnegative solutions of the equations, and then we show that the nonnegative solution can be approached by a monotone waveform relaxation sequence provided the initial iteration is chosen properly. The choice of initial iteration is critical and we give a method of finding it. Finally, we present an example to illustrate the efficiency of our method.展开更多
A certain variety of non-switched polynomials provides a uni-figure representation for a wide range of linear functional equations. This is properly adapted for the calculations. We reinterpret from this point of view...A certain variety of non-switched polynomials provides a uni-figure representation for a wide range of linear functional equations. This is properly adapted for the calculations. We reinterpret from this point of view a number of algorithms.展开更多
To solve Fredholm integral equations of the second kind, a generalized linear functional is introduced and a new function-valued Padé-type approximation is defined. By means of the power series expansion of the s...To solve Fredholm integral equations of the second kind, a generalized linear functional is introduced and a new function-valued Padé-type approximation is defined. By means of the power series expansion of the solution, this method can construct an approximate solution to solve the given integral equation. On the basis of the orthogonal polynomials, two useful determinant expressions of the numerator polynomial and the denominator polynomial for Padé-type approximation are explicitly given.展开更多
This is a study of the Durand-Kerner and Nourein methods for finding the roots of a given algebraic equation simultaneously. We consider the conditions under which the iterative methods fail. The numerical example is ...This is a study of the Durand-Kerner and Nourein methods for finding the roots of a given algebraic equation simultaneously. We consider the conditions under which the iterative methods fail. The numerical example is presented.展开更多
By looking at the situation when the coefficients Pj(z)(j=1,2,…,n-1)(or most of them) are exponential polynomials,we investigate the fact that all nontrivial solutions to higher order differential equations f((n))+Pn...By looking at the situation when the coefficients Pj(z)(j=1,2,…,n-1)(or most of them) are exponential polynomials,we investigate the fact that all nontrivial solutions to higher order differential equations f((n))+Pn-1(z)f((n-1))+…+P0(z)f=0 are of infinite order.An exponential polynomial coefficient plays a key role in these results.展开更多
Based on the dynamic equation, the performance functional and the system constraint equation of time-invariant discrete LQ control problem, the generalized Riccati equations of linear equality constraint system are ob...Based on the dynamic equation, the performance functional and the system constraint equation of time-invariant discrete LQ control problem, the generalized Riccati equations of linear equality constraint system are obtained according to the minimum principle, then a deep discussion about the above equations is given, and finally numerical example is shown in this paper.展开更多
文摘Starting from general Jacobi polynomials we derive for the Ul-traspherical polynomials as their special case a set of related polynomials which can be extended to an orthogonal set of functions with interesting properties. It leads to an alternative definition of the Ultraspherical polynomials by a fixed integral operator in application to powers of the variable u in an analogous way as it is possible for Hermite polynomials. From this follows a generating function which is apparently known only for the Legendre and Chebyshev polynomials as their special case. Furthermore, we show that the Ultraspherical polynomials form a realization of the SU(1,1) Lie algebra with lowering and raising operators which we explicitly determine. By reordering of multiplication and differentiation operators we derive new operator identities for the whole set of Jacobi polynomials which may be applied to arbitrary functions and provide then function identities. In this way we derive a new “convolution identity” for Jacobi polynomials and compare it with a known convolution identity of different structure for Gegenbauer polynomials. In short form we establish the connection of Jacobi polynomials and their related orthonormalized functions to the eigensolution of the Schrödinger equation to Pöschl-Teller potentials.
文摘In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear equations by expanding the approximate solution in terms of shifted Legendre polynomials with unknown coefficients. The performance of presented method has been compared with other methods, namely Sinc-Galerkin, quadratic spline collocation and LiuLin method. Numerical examples show better accuracy of the proposed method. Moreover, the computation cost decreases at least by a factor of 6 in this method.
文摘This paper is devoted to studying the approximate solution of singular integral equations by means of Chebyshev polynomials. Some examples are presented to illustrate the method.
文摘In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.
文摘This work is aim at providing a numerical technique for the Volterra integral equations using Galerkin method. For this purpose, an effective matrix formulation is proposed to solve linear Volterra integral equations of the first and second kind respectively using orthogonal polynomials as trial functions which are constructed in the interval [-1,1] with respect to the weight function w(x)=1+x<sup>2</sup>. The efficiency of the proposed method is tested on several numerical examples and compared with the analytic solutions available in the literature.
基金Project supported by the National Natural Science Foundation of China by Jiangsu Provincial Natural Science Foundation
文摘Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some equivalent system,,; are given. Using them the paper shows how to state properly initial and boundary conditions for these DAEs. The existence and uniqueness theory of the solution of the initial and boundary value problems for higher index DAEs are proposed.
文摘In this paper, we propose a new variation of the Adomian polynomials, which we call the degenerate Adomian polynomials, for the power series solutions of nonlinear ordinary differential equations with nonseparable nonlinearities. We establish efficient algorithms for the degenerate Adomian polynomials. Next we compare the results by the Adomian decomposition method using the classic Adomian polynomials with the results by the Rach-Adomian-Meyers modified decomposition method incorporating the degenerate Adomian polynomials, which itself has been shown to be a confluence of the Adomian decomposition method and the power series method. Convergence acceleration techniques including the diagonal Pade approximants are considered, and new numeric algorithms for the multistage decomposition are deduced using the degenerate Adomian polynomials. Our new technique provides a significant advantage for automated calculations when computing the power series form of the solution for nonlinear ordinary differential equations. Several expository examples are investigated to demonstrate its reliability and efficiency.
文摘One of the solution techniques used for ordinary differential equations, partial and integral equations is the Elzaki Transform. This paper is an extension of Mamadu and Njoseh [1] numerical procedure (Elzaki transform method (ETM)) for computing delay differential equations (DDEs). Here, a reconstructed Elzaki transform method (RETM) is proposed for the solution of DDEs where Mamadu-Njoseh polynomials are applied as basis functions in the approximation of the analytic solution. Using this strategy, a numerical illustration as in Ref.[1] is provided to the RETM as a basis for comparison to guarantee accuracy and consistency of the method. All numerical computations were performed with MAPLE 18 software.
文摘The second order Euler-Lagrange equations are transformed to a set of first order differential/algebraic equations, which are then transformed to state equations by using local parameterization. The corresponding discretization method is presented, and the results can be used to implementation of various numerical integration methods. A numerical example is presented finally.
文摘This research aims to solve Differential Algebraic Equation (DAE) problems in their original form, wherein both the differential and algebraic equations remain. The Newton or Newton-Broyden technique along with some integrators such as the Runge-Kutta method is coupled together to solve the problems. Experiments show that the method developed in this paper is efficient, as it demonstrates that implementation of the method is not difficult, and such method is able to provide approximate solutions with ease within some desired accuracy standards.
基金Supported by the Shiraz University of Technology,Shiraz,Iran.
文摘In this paper,the Bernstein polynomials method is proposed for the numerical solution of Fredholm integro-differential-difference equation with variable coefficients and mixed conditions.This method is using a simple computational manner to obtain a quite acceptable approximate solution.The main characteristic behind this method lies in the fact that,on the one hand,the problem will be reduced to a system of algebraic equations.On the other hand,the efficiency and accuracy of the Bernstein polynomials method for solving these equations are high.The existence and uniqueness of the solution have been proved.Moreover,an estimation of the error bound for this method will be shown by preparing some theorems.Finally,some numerical experiments are presented to show the excellent behavior and high accuracy of this algorithm in comparison with some other well-known methods.
基金support provided by the Ministry of Science and Technology,Taiwan,ROC under Contract No.MOST 110-2221-E-019-044.
文摘How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.
基金Project supported by the Foundation for the Authors of the National Excellent Doctoral Thesis Award of China (200720)
文摘In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.
基金supported by the Natural Science Foundation of China(NSFC)under grant 11501436Young Talent fund of University Association for Science and Technology in Shaanxi,China(20170701)
文摘Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As far as we known, the nonnegativity of solutions of the nonlinear fractional differential-algebraic equations is still not studied. In this article, we investigate the nonnegativity of solutions of the equations. Firstly, we discuss the existence of nonnegative solutions of the equations, and then we show that the nonnegative solution can be approached by a monotone waveform relaxation sequence provided the initial iteration is chosen properly. The choice of initial iteration is critical and we give a method of finding it. Finally, we present an example to illustrate the efficiency of our method.
文摘A certain variety of non-switched polynomials provides a uni-figure representation for a wide range of linear functional equations. This is properly adapted for the calculations. We reinterpret from this point of view a number of algorithms.
基金Project supported by the National Natural Science Foundation of China (No. 10271074)
文摘To solve Fredholm integral equations of the second kind, a generalized linear functional is introduced and a new function-valued Padé-type approximation is defined. By means of the power series expansion of the solution, this method can construct an approximate solution to solve the given integral equation. On the basis of the orthogonal polynomials, two useful determinant expressions of the numerator polynomial and the denominator polynomial for Padé-type approximation are explicitly given.
文摘This is a study of the Durand-Kerner and Nourein methods for finding the roots of a given algebraic equation simultaneously. We consider the conditions under which the iterative methods fail. The numerical example is presented.
基金supported partly by the National Natural Science Foundation of China(12171050,11871260)National Science Foundation of Guangdong Province(2018A030313508)。
文摘By looking at the situation when the coefficients Pj(z)(j=1,2,…,n-1)(or most of them) are exponential polynomials,we investigate the fact that all nontrivial solutions to higher order differential equations f((n))+Pn-1(z)f((n-1))+…+P0(z)f=0 are of infinite order.An exponential polynomial coefficient plays a key role in these results.
文摘Based on the dynamic equation, the performance functional and the system constraint equation of time-invariant discrete LQ control problem, the generalized Riccati equations of linear equality constraint system are obtained according to the minimum principle, then a deep discussion about the above equations is given, and finally numerical example is shown in this paper.
