On some inequalities for the gamma and psi functions

We present new inequalities for the gamma and psi functions, and we provide new classes of completely monotonic, star-shaped, and super-additive functions which are related to and .

  

任意流场中稀疏颗粒运动方程及其性质

通过研究和比较稀疏刚性颗凿相对流体运动时所受到的各种人艇力，并考虑在不同流动条件下各作用修正问题从而得到了任意流场中稀疏颗粒运动方程的一般形式，然后分析了该方程的数学性质，采用积分变换方法对该方程的基本形式进行了理论解析，最后探讨了闰物理性质对其运动规律的影响，以及几种典型流场中稀疏颗粒的运动特性，得到了一些有意义的结论  

The Fourier-series method for inverting transforms of probability distributions

This paper reviews the Fourier-series method for calculating cumulative distribution functions (cdf's) and probability mass functions (pmf's) by numerically inverting characteristic functions, Laplace transforms and generating functions. Some variants of the Fourier-series method are remarkably easy to use, requiring programs of less than fifty lines. The Fourier-series method can be interpreted as numerically integrating a standard inversion integral by means of the trapezoidal rule. The same formula is obtained by using the Fourier series of an associated periodic function constructed by aliasing; this explains the name of the method. This Fourier analysis applies to the inversion problem because the Fourier coefficients are just values of the transform. The mathematical centerpiece of the Fourier-series method is the Poisson summation formula, which identifies the discretization error associated with the trapezoidal rule and thus helps bound it. The greatest difficulty is approximately calculating the infinite series obtained from the inversion integral. Within this framework, lattice cdf's can be calculated from generating functions by finite sums without truncation. For other cdf's, an appropriate truncation of the infinite series can be determined from the transform based on estimates or bounds. For Laplace transforms, the numerical integration can be made to produce a nearly alternating series, so that the convergence can be accelerated by techniques such as Euler summation. Alternatively, the cdf can be perturbed slightly by convolution smoothing or windowing to produce a truncation error bound independent of the original cdf. Although error bounds can be determined, an effective approach is to use two different methods without elaborate error analysis. For this purpose, we also describe two methods for inverting Laplace transforms based on the Post-Widder inversion formula. The overall procedure is illustrated by several queueing examples.  

有限区间上的分数阶扩散-波方程定解问题与Laplace变换

求解了如下的分数阶扩散-波方程定解问题0Dαtu=2ux2,00,0<α≤2,u(0,t;α)=0,u(1,t;α)=θ(t),u(x,0+;α)=0,当1<α≤2时,还有ut(x,0+;α)=0.其中θ(t)是Heaviside单位阶跃函数,0Dαt为关于时间t的α阶Caputo分数阶导数算子,u=u(x,t;α)为时间t的因果函数(即t<0时恒为零的函数).利用Laplace变换的复围道积分反演和离散化反演及FoxH函数理论,给出在计算上对大的t和小的t分别适用的解的表达式.  

时间分数阶扩散-反应方程

1、引言 当前，对含有非整数阶导数和积分的方程的研究正引起越来越多学者的关注，这类分数阶导数和积分将广泛应用于科学和工程的各个领域．  

A parallel method for time discretization of parabolic equations based on Laplace transformation and quadrature

We consider the discretization in time of an inhomogeneous parabolicequation in a Banach space setting, using a representation ofthe solution as an integral along a smooth curve in the complexleft half-plane which, after transformation to a finite interval,is then evaluated to high accuracy by a quadrature rule. Thisreduces the problem to a finite set of elliptic equations withcomplex coefficients, which may be solved in parallel. The paperis a further development of earlier work by the authors, wherewe treated the homogeneous equation in a Hilbert space framework.Special attention is given here to the treatment of the forcingterm. The method is combined with finite-element discretizationin spatial variables.  

Time discretization of an evolution equation via Laplace transforms

Following earlier work by Sheen, Sloan, and Thomée concerningparabolic equations we study the discretization in time of aVolterra type integro-differential equation in which the integraloperator is a convolution of a weakly singular function andan elliptic differential operator in space. The time discretizationis accomplished by using a modified Laplace transform in timeto represent the solution as an integral along a smooth curveextending into the left half of the complex plane, which isthen evaluated by quadrature. This reduces the problem to afinite set of elliptic equations with complex coefficients,which may be solved in parallel. Stability and error boundsof high order are derived for two different choices of the quadraturerule. The method is combined with finite-element discretizationin the spatial variables.  

On the expected discounted penalty function at ruin of a surplus process with interest

In this paper, we study the expected value of a discounted penalty function at ruin of the classical surplus process modified by the inclusion of interest on the surplus. The ‘penalty’ is simply a function of the surplus immediately prior to ruin and the deficit at ruin. An integral equation for the expected value is derived, while the exact solution is given when the initial surplus is zero. Dickson’s [Insurance: Mathematics and Economics 11 (1992) 191] formulae for the distribution of the surplus immediately prior to ruin in the classical surplus process are generalised to our modified surplus process.  

基于直接数值积分的Laplace逆变换方法的比较研究

为了探讨各种数值积分方法,如梯形公式、Simpson法、Gauss积分方法和振荡函数积分方法等,在数值Laplace逆变换中的应用效果,本文进行了基于各种离散数值积分公式的Laplace逆变换方法的比较研究,涉及到24种方法,针对Davies和Martin的16个考题,给出了数值比较结果,得出了一些新的结论。  

超可微函数空间ε_*和D_*中的乘法和卷积运算

利用Fourier-Laplace变换对ω-超可微函数空间ε_*(R~N)和ω-试验函数空间D*(R~N)中的乘法和卷积运算进行了讨论,并且证明了D(R~N)是D*(R~N)的乘子空间,在卷积意义下D(R~N)是ε_*(R~N)的乘子空间,且在D*(R~N)中Parseval等式成立.