线性Hamilton系统边值问题的保辛数值方法

以Hamilton系统的正则变换和生成函数为基础研究线性时变Hamilton系统边值问题的保辛数值求解算法.根据第二类生成函数系数矩阵与状态传递矩阵的关系,构造了生成函数系数矩阵的区段合并递推算法,并进一步将递推算法推广到线性非齐次边值问题中;然后利用生成函数的性质将边值问题转化为初值问题,最后采用初值问题的保辛算法求解以达到整个Hamilton系统保辛的目的.数值算例表明该方法能够有效地求解线性齐次与非齐次问题,并能很好地保持Hamilton系统的固有特性.     

一类二阶常系数非齐次线性微分方程及边值问题的解法

利用非齐次方程通解方法和Green函数法给出了非齐次项为点源函数的二阶常系数线性常微分方程及边值问题的求解方法和公式.然后以渗流力学一类具体问题为例进行了论证.结果表明这两种方法在本质上是一致的,所得到的结果是相互吻合的.该点源解可用于分析相关边值问题,并可用来求解具有一般非齐次项的微分方程及相关定解问题.     

一类变系数非齐次调和方程边值问题的级数解

利用由三角级数和幂级数复合构成的函数项级数的有关性质,得到了一类变系数非齐次调和方程边值问题的级数解.使变系数非齐次调和方程边值问题的求解有了新的进展.     

矩阵损失函数下回归系数矩阵的非齐次线性估计的可容许性

在二次矩阵损失函数下研究了协方差矩阵未知的多元线性模型中回归系数矩阵的可估线性函数的矩阵非齐次线性估计的可容许性,给出了矩阵非齐次线性估计在线性估计类中可容许的一个充要条件.     

应用第二类Chebyshev小波求解高阶多点边值问题的数值解（英文）

本文研究了一类高阶多点边值问题的数值解法问题.利用第二类Chebyhsev小波及其积分算子矩阵,将线性与非线性高阶常微分方程多点边值问题转化为代数方程组进行求解.通过与现有文献算法结果的比较,说明了该算法求解高阶多点边值问题的准确性与有效性.扩展了高阶多点边值问题的数值求解方法.     

分片Bernstein多项式的样条配点法求解四阶微分方程

本文对四阶微分方程边值问题,给出一种基于分片Bernstein多项式的样条配点法求解,该格式构造过程容易理解,形成的线代数方程组系数矩阵稀疏,可用迭代法求解.数值实验表明,该方法可有效求解一般四阶线性微分方程边值问题,结合非等距配置点亦可用于求解含小参数的扰动问题.     

多元线性模型中回归系数的线性可容许估计

基于Zellner的平衡损失的思想,本文提出了矩阵形式的平衡损失函数,并在该损失函数下讨论了多元回归系数线性估计的可容许性.给出了六种不同形式的可容许定义,证明了这六种容许性在齐次和非齐次线性估计类中是一致的,且得到了其共同的可容许估计的充要条件.     

多元线性模型中回归系数矩阵非齐次线性估计的可容许性

该文研究了协方差矩阵未知的多元线性模型中,二次矩阵损失函数下回归系数矩阵可估线性函数的非齐次线性估计的可容许性.不需正态分布的假设,作者给出矩阵非齐次线性估计在线性估计类中可容许的充要条件;在正态分布的假设下,作者给出矩阵非齐次线性估计在一切估计组成的估计类中可容许的充分条件.     

矩阵最小二乘问题的迭代求解及其在机器翻译中的应用

研究一类线性矩阵方程最小二乘问题的迭代法求解,利用目标函数与矩阵迹之间的关系构造了矩阵形式的"梯度"下降法迭代格式,推广了向量形式的经典"梯度"下降法,并引入了两个矩阵之间的弱正交性来刻画迭代修正量的特点.作为本文算法的应用,给出了机器翻译优化问题的一种迭代求解格式.     

二阶线性齐次微分方程边值问题相似构造解的应用及Matlab图版分析

在二阶线性齐次微分方程边值问题相似构造解式的基础上,首先利用相似构造法求解Bessel方程和变型的Bessel方程边值问题的解,然后建立了均质油藏的渗流规律的数学模型,再将均质油藏的渗流数学模型转换成变型的Bessel方程的边值问题,利用二阶线性齐次微分方程边值问题的相似构造法求解均质储层渗流的数学模型.最后通过Matlab编程进行图版分析,展示实例的函数解.这将极大地方便试进分析软件的编制,也提高了石油工作者的效率.     

Solving boundary value problems,integral, and integro-differential equations using Gegenbauer integration matrices

We introduce a hybrid Gegenbauer (ultraspherical) integration method (HGIM) for solving boundary value problems (BVPs), integral and integro-differential equations. The proposed approach recasts the original problems into their integral formulations, which are then discretized into linear systems of algebraic equations using Gegenbauer integration matrices (GIMs). The resulting linear systems are well-conditioned and can be easily solved using standard linear system solvers. A study on the error bounds of the proposed method is presented, and the spectral convergence is proven for two-point BVPs (TPBVPs). Comparisons with other competitive methods in the recent literature are included. The proposed method results in an efficient algorithm, and spectral accuracy is verified using eight test examples addressing the aforementioned classes of problems. The proposed method can be applied on a broad range of mathematical problems while producing highly accurate results. The developed numerical scheme provides a viable alternative to other solution methods when high-order approximations are required using only a relatively small number of solution nodes.     

Numerical treatment of singularly perturbed two point boundary value problems

We propose a method for numerically solving linear singularly perturbed two point boundary value problems in ordinary differential equations with a boundary layer on the left end of the underlying interval. This is a practical method and can be easily implemented on a computer. The original problem is divided into inner and outer region differential equation systems. The reduced problem is solved to obtain the terminal boundary condition. Then, a new inner region problem is created and solved as a two point boundary value problem (TPBVP). In turn, the outer region problem is also solved as a TPBVP. Both these TPBVPs are efficiently treated by employing a slightly modified classical finite difference scheme coupled with discrete invariant imbedding algorithm to obtain the numerical solutions. The stability of some recurrence relations involved in the algorithm is investigated. The proposed method is iterative on the terminal point. Some numerical examples are included, and the computational results are compared with exact solutions. It is observed that the accuracy predicted can always be achieved with very little computational effort.     

Finite integration method for partial differential equations

A finite integration method is proposed in this paper to deal with partial differential equations in which the finite integration matrices of the first order are constructed by using both standard integral algorithm and radial basis functions interpolation respectively. These matrices of first order can directly be used to obtain finite integration matrices of higher order. Combining with the Laplace transform technique, the finite integration method is extended to solve time dependent partial differential equations. The accuracy of both the finite integration method and finite difference method are demonstrated with several examples. It has been observed that the finite integration method using either radial basis function or simple linear approximation gives a much higher degree of accuracy than the traditional finite difference method.     

A numerical method for solving boundary value problems for fractional differential equations

A numerical scheme, based on the Haar wavelet operational matrices of integration for solving linear two-point and multi-point boundary value problems for fractional differential equations is presented. The operational matrices are utilized to reduce the fractional differential equation to system of algebraic equations. Numerical examples are provided to demonstrate the accuracy and efficiency and simplicity of the method.     

求解奇异摄动边值问题的精细积分法

提出了一种求解一端有边界层的奇异摄动边值问题的精细方法．首先将求解区域均匀离散,由状态参量在相邻节点间的精细积分关系式确定一组代数方程,并将其写成矩阵形式．代入边界条件后,该代数方程组的系数矩阵可化为块三对角形式,针对这一特性,给出了一种高效递推消元方法．由于在离散过程中,精细积分关系式不会引入离散误差,故所提出的方法具有极高的精度．数值算例充分证明了所提出方法的有效性．     

A meshless method based on boundary integral equations and radial basis functions for biharmonic-type problems

This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed.     

An efficient parallel processing optimal control scheme for a class of nonlinear composite systems

This article presents an efficient parallel processing approach for solving the optimal control problem of nonlinear composite systems. In this approach, the original high-order coupled nonlinear two-point boundary value problem (TPBVP) derived from the Pontryagin's maximum principle is first transformed into a sequence of lower-order decoupled linear time-invariant TPBVPs. Then, an optimal control law which consists of both feedback and forward terms is achieved by using the modal series method for the derived sequence. The feedback term specified by local states of each subsystem is determined by solving a matrix Riccati differential equation. The forward term for each subsystem derived from its local information is an infinite sum of adjoint vectors. The convergence analysis and parallel processing capability of the proposed approach are also provided. To achieve an accurate feedforward-feedback suboptimal control, we apply a fast iterative algorithm with low computational effort. Finally, some comparative results are included to illustrate the effectiveness of the proposed approach.     

结构动力方程求解的改进-精细Runge-Kutta方法

在已有精细Runge-Kutta(龙格-库塔)方法的基础上,考虑了状态空间方程非齐次项的特点和外荷载的特殊性,提出了求解结构动力方程的改进精细Runge-Kutta方法.通过对矩阵进行分块计算,在利用原有精细Runge-Kutta方法高精度的同时进一步提高了计算效率,有利于大型结构的长时间仿真.将改进精细Runge-Kutta方法应用于结构动力方程求解,为其求解提供一种新方法.数值算例表明了改进方法的正确性和有效性.     

A High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

The dual reciprocity boundary elements method for the linear and nonlinear two‐dimensional time‐fractional partial differential equations

In this paper, we apply the dual reciprocity boundary elements method for the numerical solution of two‐dimensional linear and nonlinear time‐fractional modified anomalous subdiffusion equations and time‐fractional convection–diffusion equation. The fractional derivative of problems is described in the Riemann–Liouville and Caputo senses. We employ the linear radial basis function for interpolation of the nonlinear, inhomogeneous and time derivative terms. This method is improved by using a predictor–corrector scheme to overcome the nonlinearity which appears in the nonlinear problems under consideration. The accuracy and efficiency of the proposed schemes are checked by five test problems. The proposed method is employed for solving some examples in two dimensions on unit square and also in complex regions to demonstrate the efficiency of the new technique. Copyright © 2016 John Wiley & Sons, Ltd.