一类显式的k阶线性k步法基本公式

本文给出了一类比Adams-Bashforth方法的局部截断误差主项系数小和绝对稳定区间大的显式k阶线性k步法基本公式.作者求出了公式的分数形式的精确系数,阶数和局部截断误差主项系数,给出了3-9步公式的绝对稳定区间,构造了由新公式的4阶显式公式和一个同阶隐式基本公式组合而成的特殊预估-校正方法,它的绝对稳定区间大于预估公式而且等于校正公式, 比著名的Adams-Bashforth-Moulton预估校正方法的绝对稳定区间大, 最后用数值试验对结果进行了验证,适合于求解常微分方程初值问题.     

一类A(α)稳定的k阶线性k步法公式

本文给出了一类与Gear方法类似的k阶线性k步法隐式公式.作者还求出了公式的分数形式的系数,阶数和局部截断误差主项系数,并验证了2-6步公式都具有A(α)稳定的,计算出了它们的幅角α.最后用对比数值实验验证了公式确实是稳定的,并且适合于求解刚性常微分方程.     

一类A(α)稳定的k阶线性k步法公式

本文给出了一类与Gear方法类似的κ阶线性κ步法隐式公式.作者还求出了公式的分数形式的系数,阶数和局部截断误差主项系数,并验证了2-6步公式都具有A(α)稳定的,计算出了它们的幅角α.最后用对比数值实验验证了公式确实是稳定的,并且适合于求解刚性常微分方程.     

用加权平均方法构造新的隐式线性多步法公式

在已知的线性多步法公式中,用两个较适合的线性多步法进行加权平均就能构造出一系列新的隐式线性多步法公式,而且其中有些公式可能具有较好的性质,如稳定域增大.从而使得解刚性方程时,可以根据对稳定域与截断误差不同的需求来选择公式,以达到在适合的稳定域下,截断误差最小.经过数值试验验证,本文举出的实例中用加权平均方法构造出的有些新公式的稳定域大于原来两个公式任一个的稳定域,可应用于求解常微分方程初值问题的刚性问题.     

一类stiff稳定的线性多步法

§1.引言 在常微分方程初值问题的数值方法中,线性多步法是最简单、使用最广泛的方法之一.但由于现存的线性多步方法的绝对稳定区域较小,以致在解刚性(Stiff)微分方程中受到很大限制.本文在BDF方法及[2]的基础上增加二个修正项,构造一类修正BDF的线性多步法,具有较大的绝对稳定区域.其结果如下:此类修正方法的阶与同步数的BDF方法的阶一致,其绝对稳定区域与低二阶的BDF方法大致相同,甚至更好,并给出了参数的取值范围.     

一类A(α)稳定的κ阶线性κ步法公式

本文给出了一类与Gear方法类似的κ阶线性κ步法隐式公式.作者还求出了公式的分数形式的系数,阶数和局部截断误差主项系数,并验证了2-6步公式都具有A(α)稳定的,计算出了它们的幅角α.最后用对比数值实验验证了公式确实是稳定的,并且适合于求解刚性常微分方程.     

微分求积法的一类线性多步法

利用时域微分求积法和非等距网格,构造了一类A(α)-稳定或有限区间稳定的线性多步法.根据Dahlquist等价性定理,新的线性多步公式是收敛的.理论证明了新的隐式线性多步法公式是A(α)-稳定的或有限区间稳定的.通过数值实验的对比,表明了新的线性多步法比现有的线性多步法具有更好的数值性能.     

线性Volterra多延迟积分微分方程线性多步法的GPm稳定性

讨论了多步法求解线性Volterra多延迟积分微分方程数值方法的GPm稳定.证明了对任给的步长h>0,A-稳定的线性多步法保持原线性系统的渐近稳定性,从而是GPm稳定.     

一类线性多步法关于变延迟非线性中立型微分方程的渐近稳定性

1引言中立型微分方程广泛出现于生物学、物理学及工程技术等诸多领域.数值求解中立型微分方程时,数值方法的稳定性研究具有无容置疑的重要性,其中渐近稳定性的研究是其重要组成部分.对于线性中立型延迟微分方程,渐近稳定性研究已有许多重要结果,如文献[1,2,3,4,5,6]等.对于非线性中立型变延迟微分方程,数值方法的稳定性研究近几年才有进展.2000年,Bellen等在文献[7]中讨论了Runge-Kutta法求解一类特殊的中立型延迟微分     

非线性延迟积分微分方程线性多步法的渐近稳定性

本文研究求解非线性延迟积分微分方程的线性多步法的渐近稳定性,其中积分部分采用复化梯形公式计算,结果表明:在问题真解渐近稳定的条件下,A-稳定的线性多步法也是渐近稳定的.     

预估-校正方法的绝对稳定性讨论

预估-校正方法,即PECE方法,常被用于求解常微分方程的初值问题.而一般文献中常只讨论了单个线性多步法公式的稳定性问题,很少涉及由一个显式公式和一个隐式公式组合而成的PECE方法的稳定性.本文应用根轨迹法和对分法讨论了常用的PECE方法的稳定性,求出了一些常用PECE方法的组合公式的绝对稳定区间和绝对稳定区域,并用数值...     

THE STABILITY OF LINEAR MULTISTEP METHODS FOR LINEAR SYSTEMS OF NEUTRAL DIFFERENTIAL EQUATIONS

1. IntroductionFor a large class of electrical ~ks containing IOssless transmission lines the describingequations can be reduced to a system of neutral differential equations.This paper deals with the numerical solution of initial value problems for syste…     

CONJUGATE-SYMPLECTICITY OF LINEAR MULTISTEP METHODS

For the numerical treatment of Hamiltonian differential equations, symplectic integrators are the most suitable choice, and methods that are conjugate to a symplectic integrator share the same good long-time behavior. This note characterizes linear multistep methods whose underlying one-step method is conjugate to a symplectic integrator. The bounded- hess of parasitic solution components is not addressed.     

D-CONVERGENCE AND STABILITY OF A CLASS OF LINEAR MULTISTEP METHODS FOR NONLINEAR DDES

1. IntroductionConsider the following nonlinear delay problem{:;\f>>:v{t(tf,?,,<,>3,<'~">>,:: 5:,3Ti:,,,,, [l:::;where y: R - C",T > 0 is a delay term, f: [t.,T] x CN x CN - CN and W(t):[to -- T, tol - CN denotes a given initial function. Thoroughout this paper 9 the problem(1.1) is supposed to have a unique solution y(t), which satisfies11 y(')(t) 115 Mi, t e [to ~ T,T]here norm 11. 11 is defined by 11 x II'=< xgx > (Vx E C"), and Mi > 0 are someconstants.Definition 1.1.[1] The clas…     

DELAY-DEPENDENT TREATMENT OF LINEAR MULTISTEP METHODS FOR NEUTRAL DELAY DIFFERENTIAL EQUATIONS

This paper deals with a delay-dependent treatment of linear multistep methods for neutral delay differential equations y'(t) = ay(t) + by(t - τ) + cy'(t - τ), t > 0, y(t) = g(t), -τ≤ t ≤ 0, a,b andc ∈ R. The necessary condition for linear multistep methods to be Nτ(0)-stable is given. It is shown that the trapezoidal rule is Nτ(0)-compatible. Figures of stability region for some linear multistep methods are depicted.     

一类线性抛物型方程组的交替方向多步法及其理论分析

Efficient multistep procedure for time-stepping Galerkin method in which we use an alternating direction preconditioned iterative methods for approximately solving the linear equations arising at each timestep in a discrete Galerkin method for a class of linear parabolic systems is derived and analyzed. The optimal order error estimate is obtained. Numerical experiments show that the method has the characteristics of high efficiency and high accuracy.     

求解变延迟微分方程的一类线性多步方法的收缩性

This paper presents a sufficient condition on the contractivity of theoretical solution for a class of nonlinear systems of delay differential equations with many variable delays(MDDEs), which is weak,compared with the sufficient condition of previous articles.In addition,it discusses the numerical stability properties of a class of special linear nmltistep methods for this class nonlinear problems.And it is pointed out that not only the backwm‘d Euler method but also this class of linear multistep methods are GRNm-stable if linear interpolation is used.     

求解非线性中立型延迟微分方程一类线性多步方法的收敛性

本文致力于带有Lagrang插值的一类线性多步法求解非线性中立型延迟微分方程的误差分析.证明了一个p′阶的线性多步方法配上一个q阶的Lagrang插值导致一个minf[p′,q 1]阶的E-(或EB-)收敛的非线性中立型延迟微分方程数值方法.