线性定常系统非齐次两点边值问题的扩展精细积分方法

提出了一种求解非齐次线性两点边值问题的高精度和高稳定的扩展精细积分方法（EPIM）.首先引入了区段量（即区段矩阵和区段向量）来离散非齐次线性微分方程,建立了非齐次两点边值问题基于区段量的求解框架.在该框架下,不同区段的区段量可以并行计算,整体代数方程组的集成不依赖于边界条件.然后引入区段响应矩阵来处理两点边值问题的非齐次项,导出了多项式函数、指数函数、正/余弦函数及其组合函数形式的非齐次项对应的区段响应矩阵的加法定理,结合增量存储技术提出了EPIM.对具有上述函数形式的非齐次项,该方法可以得到计算机上的精确解,一般形式的非齐次项则利用上述函数近似求解.最后通过两个具有刚性特征的数值算例验证了该方法的高精度和高稳定性.     

正规升列在参数代数方程组求解中的应用

本文提出一种利用多项式系统的正规零点分解的算法来求解代数方程组以及带有参数的代数方程组的方法．对于给定的的代数方城组,通过正规分解,可以得到一组具有三角形式的分解．根据这种三角形式,我们可以给出代数方程组的所有解．而对于带有参数方程组,将给出方程组有解时参数需满足的条件．进一步,对于给定的参数值,正规分解中得到三角形式仍然保持,通过求解三角形式的方程组从而得出原参数方程组的解．     

病态代数方程求解的一种改进精细积分法

通过在病态代数方程精细积分法的基础上增加一个迭代改善算法,建立了病态代数方程求解的改进精细积分法．该方法进一步提高了病态代数方程精细积分法的精度和效率,具有良好的应用前景．算例证明了该方法在病态代数方程求解中的有效性．     

一种求解变系数渗压固结微分方程的数值方法

对变系数渗压固结微分方程的求解过程进行了深入研究,提出一种精细积分半解析数值方法,首先对渗压固结微分方程在空间离散,建立起对于时间的常微分方程组,然后对时程积分,利用矩阵指数函数可以在计算机字长范围内精确计算的特点,得到精密解答.当指数函数用Taylor展开式的一阶近似替代时,精细积分转化为差分方程.用matlab语言编写程序进行求解,得到孔隙比在固结过程中的分布规律,并通过模型试验进行了验证.     

带常系数的Cauchy型奇异积分方程的快速方法

Petrov-Galerkin 方法是研究Cauchy型奇异积分方程的最基本的数值方法. 用此方法离散积分方程可得一系数矩阵是稠密的线性方程组. 如果方程组的阶比较大, 则求解此方程组所需要的计算复杂度则会变得很大. 因此, 发展此类方程的快速数值算法就变成了必然. 该文将就对带常系数的Cauchy型奇异积分方程给出一种快速数值方法. 首先用一稀疏矩阵来代替稠密系数矩阵, 其次用数值积分公式离散上述方程组得到其完全离散的形式,然后用多层扩充方法求解此完全离散的线性方程组. 证明经过上述过程得到方程组的逼进解仍然保持了最优阶, 并且整个过程所需要的计算复杂度是拟线性的. 最后通过数值实验证明结论.     

基于模块脉冲函数的非线性随机It?-Volterra积分方程数值解（英文）

为求解非线性随机It?o-Volterra积分方程,本文介绍了一种基于模块脉冲函数的有效数值方法.运用模块脉冲函数的积分算子矩阵将非线性随机积分方程转化为代数方程.通过误差分析,证明该方法收敛速度良好.最后,利用实例验证了此方法的有效性.     

热弹性动力学耦合问题的微分代数方法

对一般的热机械问题提出了一种有效的数值方法,并对二维的热弹性问题进行了测试．该方法的基本思路是将描述热机械耦合问题的偏微分方程进行降阶,使之成为一组微分代数方程,应力应变关系被写成代数方程．所得到的微分代数系统采用全隐式的向后差分公式进行求解．对该方法进行了详细的说明．为了验证该方法的有效性,将其应用于一个动态非耦合的热弹性问题的求解和一个耦合的二维热弹性问题的求解．     

基于模块脉冲函数的非线性随机It\^{o}-Volterra积分方程数值解[英文]

为求解非线性随机It\^{o}-Volterra积分方程, 本文介绍了一种基于模块脉冲函数的有效数值方法. 运用模块脉冲函数的积分算子矩阵将非线性随机积分方程转化为代数方程. 通过误差分析, 证明该方法收敛速度良好. 最后, 利用实例验证了此方法的有效性.     

预条件广义极小残余新算法

研究Krylov子空间广义极小残余算法(GMRES(m))的基本理论，给出GMRES(m)算法透代求解所满足的代数方程组．深入探讨算法的收敛性与方程组系数矩阵的密切关系，提出一种改进GMRES(m)算法收敛性的新的预条件方法，并作出相关论证．     

具对称半正定矩阵线代数方程组的直接法

有限元方法应用于椭圆型第二类达值问题所得到的线代数方程组的系数矩阵都是对称半定的。在土木结构、造船、航空、机械零件的应力分析中常常遇到。目前,在组合结构大型程序中,纷纷采用各种形式的直接法,可是遇到具奇异半定矩阵的情形,那些算法与程序会失效。本文将对文献[3]的方法进行修改,提出了求解这类线代数方程的直接法,这种方法无需行列变换。同样可与变带宽的分块直接法和波前直接法配合使用,克     

Bernstein polynomial basis for numerical solution of boundary value problems

The purpose of this paper is to propose a computational method for the approximate solution of linear and nonlinear two-point boundary value problems. In order to approximate the solution, the expansions in terms of the Bernstein polynomial basis have been used. The properties of the Bernstein polynomial basis and the corresponding operational matrices of integration and product are utilized to reduce the given boundary value problem to a system of algebraic equations for the unknown expansion coefficients of the solution. On this approach, the problem can be solved as a system of algebraic equations. By considering a special case of the problem, an error analysis is given for the approximate solution obtained by the present method. At last, five examples are examined in order to illustrate the efficiency of the proposed method.     

病态代数系统求解的精细迭代方法

提出了病态代数系统求解的精细迭代方法．首先利用一个小参数对病态矩阵加以改良,将原病态系统的求解问题转化为该改良系统的求解问题．然后利用精细积分法给出了改良矩阵求逆的高精度方法．该方法具有高精度、高效率的优点,且对改良参数的适应性较好,具有良好的应用前景．理论和数值分析证明了该方法的有效性．     

Optimal Control of One-Dimensional Partial Differential Algebraic Equations with Applications

We present an approach to compute optimal control functions in dynamic models based on one-dimensional partial differential algebraic equations (PDAE). By using the method of lines, the PDAE is transformed into a large system of usually stiff ordinary differential algebraic equations and integrated by standard methods. The resulting nonlinear programming problem is solved by the sequential quadratic programming code NLPQL. Optimal control functions are approximated by piecewise constant, piecewise linear or bang-bang functions. Three different types of cost functions can be formulated. The underlying model structure is quite flexible. We allow break points for model changes, disjoint integration areas with respect to spatial variable, arbitrary boundary and transition conditions, coupled ordinary and algebraic differential equations, algebraic equations in time and space variables, and dynamic constraints for control and state variables. The PDAE is discretized by difference formulae, polynomial approximations with arbitrary degrees, and by special update formulae in case of hyperbolic equations. Two application problems are outlined in detail. We present a model for optimal control of transdermal diffusion of drugs, where the diffusion speed is controlled by an electric field, and a model for the optimal control of the input feed of an acetylene reactor given in form of a distributed parameter system.     

Laguerre series direct method for variational problems

A direct method for solving variational problems via Laguerre series is presented. First, an operational matrix for the integration of Laguerre polynomials is introduced. The variational problems are reduced to the solution of algebraic equations. An illustrative example is given.     

Free vibration of non-uniform Euler–Bernoulli beam under various supporting conditions using Chebyshev wavelet collocation method

This paper proposes operational matrix of rth integration of Chebyshev wavelets. A general procedure of this matrix is given. Operational matrix of rth integration is taken as rth power of operational matrix of first integration in literature. But, this study removes this disadvantage of Chebyshev wavelets method. Free vibration problems of non-uniform Euler–Bernoulli beam under various supporting conditions are investigated by using Chebyshev Wavelet Collocation Method. The proposed method is based on the approximation by the truncated Chebyshev wavelet series. A homogeneous system of linear algebraic equations has been obtained by using the Chebyshev collocation points. The determinant of coefficients matrix is equated to the zero for nontrivial solution of homogeneous system of linear algebraic equations. Hence, we can obtain ith natural frequencies of the beam and the coefficients of the approximate solution of Chebyshev wavelet series that satisfied differential equation and boundary conditions. Mode shapes functions corresponding to the natural frequencies can be obtained by normalizing of approximate solutions. The computed results well fit with the analytical and numerical results as in the literature. These calculations demonstrate that the accuracy of the Chebyshev wavelet collocation method is quite good even for small number of grid points.     

Bernoulli wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations

In this paper, a new numerical method for solving fractional differential equations is presented. The fractional derivative is described in the Caputo sense. The method is based upon Bernoulli wavelet approximations. The Bernoulli wavelet is first presented. An operational matrix of fractional order integration is derived and is utilized to reduce the initial and boundary value problems to system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique.     

Newton-like methods for two-point boundary value problems

A class of Newton-like methods for discrete two-point boundary value problems is constructed from the sum equation formulation of the problem. Each step of the Newton-like method can be described as first solving a system of linear algebraic equations. The solution vector of this system gives boundary values to a number of discrete boundary value problems which can be solved explicitly.     

Method of successive approximations for three-dimensional laminar boundary layer problems (locally self-similar case) : PMM vol. 37, n6, 1973, pp. 974–983

We present an analytical method for the computation of problems of incompressible boundary layer theory based on an application of the method of successive approximations. The system of equations is reduced to a form suitable for integration. Parameters characterizing the external flow and the body geometry are contained only in the coefficients of the system and do not enter into the boundary conditions. The transformed momentum equations are integrated across the boundary layer from a current value to infinity with the boundary conditions taken into account. If the integration is made from zero to infinity, then the equations pass over into the Kármán relations. Integrating the system of equations a second time, using the boundary conditions at the wall, we obtain a system of nonlinear integro-differential equations. To solve this system of equations we apply the method of successive approximations. To satisfy the boundary Conditions at infinity we introduce, at each step of the iterations, unknown “governing” functions. From the conditions at the outer side of the boundary layer we obtain additional equations for their determination. With the iterational algorithm formulated in this way, the boundary conditions, both on the body and at the outer side of the boundary layer; are satisfied automatically.We consider a locally self-similar approximation. In this case, relative to the “governing” functions, we obtain an algebraic system of equations. We write out the solution in the first approximation. The results obtained in the first approximation are compared with the results of finite-difference computations for a wide range of problems. The results obtained in this paper are compared with those obtained in [1] for the flow in the neighborhood of a stagnation point. An indication is given of the nonuniqueness of the solutions of the three-dimensional boundary layer equations.