基于Gauss全局径向基函数的近岸浅水变形波高数值计算新方法

应用Gauss全局径向基函数来模拟波浪浅水变形波高变化方程中的未知函数,经实例分析探讨得到了一种可用于求解该方程数值解的新方法,并将其计算结果与常用数值分析方法得到的数值解相互对比印证,证明了基于Gauss全局径向基函数法计算结果的正确性.经验证,Gauss径向基函数法的平均计算误差相比其他方法均要小,表明该方法拥有更高的计算精度.同时,根据Gauss全局径向基函数的逼近结果,得出了浅水变形波高变化微分方程数值解的拟合函数,在实际工程中,可以利用该拟合函数来代替原方程的解析解,研究成果可为求解近岸浅水区域波浪运动提供一种新思路.     

求解Bratu型方程的径向基函数逼近法

基于径向基函数可以逼近几乎所有函数的强大逼近功能,借鉴弹塑性静力学的处理方法,提出位移、速度、加速度联合插值的径向基函数表达式,结合MATLAB数值软件进行计算机编程,成功求解了Bratu型强非线性方程,并给出相应的相对误差.通过分析几种典型的算例,并将计算结果与一些现有的数值分析法得到的数值解进行对比,表明了该方法的可行性和精确性,为求解强非线性Bratu型方程提供了一种新思路.     

求解非线性Black-Scholes模型的自适应小波精细积分法

针对非线性Black-Scholes方程,基于quasi-Shannon小波函数给出了一种求解非线性偏微分方程的自适应多尺度小波精细积分法.该方法首先利用插值小波理论构造了用于逼近连续函数的多尺度小波插值算子,利用该算子可以将非线性Black-Scholes方程自适应离散为非线性常微分方程组;然后将用于求解常微分方程组的精细积分法和小波变换的动态过程相结合,并利用非线性处理技术(如同伦分析技术)可有效求解非线性Black-Scholes方程.数值结果表明了该方法在数值精度和计算效率方面的优越性.     

RBF-PU方法求解二维非局部扩散问题和近场动力学问题

采用单位分解径向基函数(radial basis function partition of unity,RBF-PU)方法,数值求解了二维非局部扩散问题和近场动力学问题。主要思想是对求解区域进行局部划分,在局部子区域上分别进行函数逼近,然后加权得到未知函数的全局逼近。这种基于方程强形式的径向基函数方法在求解非局部问题时,不需要处理网格与球形邻域求交的问题,避免了额外的一层积分计算,实施简便,计算量小。数值实验显示计算结果与解析解吻合较好,RBF-PU方法可以准确有效地求解非局部扩散方程和近场动力学方程。     

基于径向基函数逼近的非线性动力系统数值求解

径向基函数具有形式简单、各向同性等优点.将径向基函数逼近的思想与加权余量配点法相结合,借鉴边值问题的求解,构造了一种求解非线性动力系统初值问题的数值方法.分析了几种较为成熟的非线性动力系统数值求解方法的优缺点.给出了实际算例,与已有方法对比,表明该方法计算过程简单、收敛性好、计算精度高.     

球面混合插值的逼近性质

球面径向基函数(SBF)和多项式样条函数均为处理球面散乱数据的有效工具. 本文考虑由球面径向基函数与球面多项式函数组成的混合插值模型, 并利用最小二乘法求解该模型. 对于该插值模型, 首先, 给出带Bessel势的Sobolev空间中的Bernstein不等式, 然后利用该不等式建立逼近正定理,并进一步给出该插值工具的误差估计. 最后, 研究该插值方式(即利用最小二乘法求解混合插值模型)的稳定性.     

基于区间三次Hermite样条小波的Poisson方程数值求解方法

提出一种新的求解Poisson方程的小波有限元方法,采用区间三次Hermite样条小波基作为多尺度有限元插值基函数,并详细讨论了小波有限元提升框架.由于小波基按照给定的内积正交,可实现相应的多尺度嵌套逼近小波有限元求解方程,在不同尺度上的插值基之间完全解耦和部分解耦.数值算例表明在求解Poisson方程时,该方法具有高的效率和精度.     

轴对称Poisson方程的Trefftz有限元解法

将径向基函数应用到一类轴对称Poisson方程的数值求解中,提出了一种Trefftz有限元计算格式.非0右端项将问题的特解引入Trefftz单元域内场,致使单元刚度方程涉及区域积分.利用径向基函数对特解近似处理,可消除区域积分,从而保持Trefftz有限元法只含边界积分的优势.为获得特解,选取求解域内所有单元的节点和形心作为基本插值点,而在求解域之外构造一个虚拟边界,在其上布置一定数目的虚拟点作为额外插值点.数值算例验证了该方法的有效性和可行性.     

基于三次样条插值函数的非线性动力系统数值求解

三次样条插值函数具有良好的收敛性、稳定性与二阶光滑性.研究了借助三次样条插值函数构造的非线性动力系统数值求解方法,分析了该方法与已有的非线性动力系统数值求解方法的优缺点,刻画了误差估计且给出了数值算例.结果表明基于三次样条插值函数构造的数值方法比已有的方法收敛速度快、逼近精度高且能够很好地逼近非线性动力系统的解析解.     

非齐次弹性力学问题双互易边界元方法研究北大核心CSCD

基于弹性力学边界元方法理论,将边界元法与双互易法结合,采用指数型基函数对非齐次项进行插值得到双互易边界积分方程.将边界积分方程离散为代数方程组,利用已知边界条件和方程特解求解方程组,得出域内位移和边界面力.指数型基函数的形状参数是由插值点最近距离的最小值决定,采用这种形状参数变化方案,分析径向基函数(RBF)插值精度以及插值稳定性.再次将指数型基函数应用到双互易边界元法中,分析双互易边界元方法下计算精度及稳定性,验证了指数型插值函数作为双互易边界元方法的径向基函数解决弹性力学域内体力项问题的有效性.     

A symmetric integrated radial basis function method for solving differential equations

In this article, integrated radial basis functions (IRBFs) are used for Hermite interpolation in the solution of differential equations, resulting in a new meshless symmetric RBF method. Both global and local approximation‐based schemes are derived. For the latter, the focus is on the construction of compact approximation stencils, where a sparse system matrix and a high‐order accuracy can be achieved together. Cartesian‐grid‐based stencils are possible for problems defined on nonrectangular domains. Furthermore, the effects of the RBF width on the solution accuracy for a given grid size are fully explored with a reasonable computational cost. The proposed schemes are numerically verified in some elliptic boundary‐value problems governed by the Poisson and convection‐diffusion equations. High levels of the solution accuracy are obtained using relatively coarse discretisations.     

Numerical solution of the nonlinear Fredholm integral equations by positive definite functions

A numerical method for solving the nonlinear Fredholom integral equations is presented. The method is based on interpolation by radial basis functions (RBF) to approximate the solution of the Fredholm nonlinear integral equations. Several examples are given and numerical examples are presented to demonstrate the validity and applicability of the method.     

Differential quadrature method based on the highest derivative and its applications

The numerical differentiation is often used when dealing with the differential equations. Using the numerical differentiation, the differential equations can be transformed into algebraic equations. Then we can get the numerical solution from the algebraic equations. But the numerical differentiation process is very sensitive to even a small level of errors. In contrast, it is expected that on average the numerical integration process is much less sensitive to errors. In this paper, we provide a new method using the DQ method based on the interpolation of the highest derivative (DQIHD) for the differential equations. The original function is then obtained by integration. In this paper, the DQIHD method was applied to the buckling analysis of thin isotropic plates and Winkler plates, the numerical results agree well with the analytic solutions, and the results show that our method is of high accuracy, of good convergence with little computational efforts. And it is easy to deal with the boundary conditions.     

基于卷积神经网络模型数值求解双曲型偏微分方程的研究

人工神经网络近年来得到了快速发展,将此方法应用于数值求解偏微分方程是学者们关注的热点问题.相比于传统方法其具有应用范围广泛（即同一种模型可用于求解多种类型方程）、网格剖分条件要求低等优势,并且能够利用训练好的模型直接计算区域中任意点的数值.该文基于卷积神经网络模型,对传统有限体积法格式中的权重系数进行优化,以得到在粗粒度网格下具有较高精度的新数值格式,从而更适用于复杂问题的求解.该网络模型可以准确、有效地求解Burgers方程和level set方程,数值结果稳定,且具有较高数值精度.     

Bivariate Polynomial Natural Spline Interpolation to Scattered Data

By means of the theory of spline interpolation in Hilbert spaces, the bivariate polynomial natural spline interpolation to scattered data is constructed. The method can easily be carried out on a computer, and parallelly generalized to high dimensional cases as well. The results can be used for numerical integration in higher dimensions and numerical solution of partial differential equations, and so on.     

Numerical solution of singular IVPs of Lane–Emden type using a modified Legendre-spectral method

In this paper, a numerical method for singular initial value problems of the Lane–Emden type in the second-order ordinary differential equations is proposed. The method changes solving the equation to solving a Volterra integral equation. We have applied the improved Legendre-spectral method to solve Lane–Emden type equations. The Legendre Gauss points are used as collocation nodes and Lagrange interpolation is employed in Volterra term. The results reveal that the method is effective, simple and accurate.     

On choosing “optimal” shape parameters for RBF approximation

Many radial basis function (RBF) methods contain a free shape parameter that plays an important role for the accuracy of the method. In most papers the authors end up choosing this shape parameter by trial and error or some other ad hoc means. The method of cross validation has long been used in the statistics literature, and the special case of leave-one-out cross validation forms the basis of the algorithm for choosing an optimal value of the shape parameter proposed by Rippa in the setting of scattered data interpolation with RBFs. We discuss extensions of this approach that can be applied in the setting of iterated approximate moving least squares approximation of function value data and for RBF pseudo-spectral methods for the solution of partial differential equations. The former method can be viewed as an efficient alternative to ridge regression or smoothing spline approximation, while the latter forms an extension of the classical polynomial pseudo-spectral approach. Numerical experiments illustrating the use of our algorithms are included.     

Convergence of Runge-Kutta Methods for Delay Differential Equations

This paper deals with the adaptation of Runge—Kutta methods to the numerical solution of nonstiff initial value problems for delay differential equations. We consider the interpolation procedure that was proposed in In 't Hout [8], and prove the new and positive result that for any given Runge—Kutta method its adaptation to delay differential equations by means of this interpolation procedure has an order of convergence equal to min {p,q}, where p denotes the order of consistency of the Runge—Kutta method and q is the number of support points of the interpolation procedure.This revised version was published online in October 2005 with corrections to the Cover Date.     

Numerical solution of differential equations using Sinc method based on the interpolation of the highest derivatives

In general, we will use the numerical differentiation when dealing with the differential equations. Thus the differential equations can be transformed into algebraic equations and then we can get the numerical solutions. But as we all have known, the numerical differentiation process is very sensitive to even a small level of errors. In contrast it is expected that on average the numerical integration process is much less sensitive to errors. In this paper, based on the Sinc method we provide a new method using Sinc method incorporated with the double exponential transformation based on the interpolation of the highest derivatives (SIHD) for the differential equations. The error in the approximation of the solution is shown to converge at an exponential rate. The numerical results show that compared with the exiting results, our method is of high accuracy, of good convergence with little computational efforts. It is easy to treat nonhomogeneous mixed boundary condition for our method, which is unlike the traditional Sinc method.