多步Runge—Kutta方法的代数稳定性

本文获得了多步Runge—是Kutta方法代数稳定的一系列必要充分条件,其中多数结果可视为关于Radau I A、Radau ⅡA及Gauss型Runge—Kutta方法已有结果的推广。     

多步隐式Runge—Kutta方法的稳定性分析

本文对多步隐式Runge—Kutta方法进行数值稳定性分析。给出了广义压缩性及弱广义压缩性的概念,并导出了多步隐式Runge—Kutta方法为广义压缩的代数条件,最后还给出了数值例子。     

刚性Volterra泛函微分方程算法理论及高效算法

首先介绍刚性Volterra泛函微分方程的稳定性理论及其数值方法的B理论。这项工作为刚性延迟微分方程、刚性积分微分方程以及其它各种类型的刚性泛函微分方程的研究提供了统一的理论基础。其次以该理论为指针推荐高效算法，其中包括向后Euler方法、二阶BDF方法、并行多值混合方法及实特征值多步Runge—Kutta法。     

非线性Volterra延迟积分微分方程Runge-Kutta方法的散逸性

将(k,l)-代数稳定的Runge.Kutta方法应用于非线性沃尔泰拉延迟积分微分方程,讨论了该方法的数值散选性,证明了该方法具有有限维和无限维散逸性.     

实时仿真的数学模型及实时Runge-Kutta算法的收敛性分析

本文讨论了实时仿真的数学模型,和实时Runge—Kutta算法的收敛性分析,给出了补尝阶与实时算法收敛阶之间关系的理论证明。     

一类求解刚性方程的并行Runge—kutta方法

随着计算机趋于微型、巨型和网络化,开展适用于求解刚性方程的并行Runge—kutta方法的研究,已越来越引起从事常微数值解研究的科研人员的注意,我们在这方面进行了一些探讨。本文利用2级3阶的R—K方法构造了一系列适合于2台并行处理机并行计算且具有稳定性很好的4阶R—K方法。(所举例子分别是A—稳定和L—稳定的4阶方法)。     

两相流重粒子运动的非线性系统方程组的计算

本文阐述两相流重粒子运动的非线性常微分方程组实际上代表一个系统的方程组,讨论了常用的解法如Runge—Kutta法等和它的结果。本文采取变步长用有限分析法把非线性项线性化,构造一个迭代序列,逐步迭代修正非线性项,求解一般形式的非线性系统方程;并验证迭代逼近解的收缩运算条件;给出计算重粒子运动轨迹和速度的算例,结果分析和说明应用的实例。     

基于威胁特征点的启发式航迹规划方法

针对现有启发式航迹规划方法在决策变量、启发函数的选取和航迹表示方法等方面存在的不足,提出了一种基于威胁特征点的航迹规划新方法。该方法以导弹运动姿态角作为决策变量,通过引入特征点构造新的启发函数,以导弹经过空间的威胁度累积值、需用过载等为代价,结合导弹性能约束,采用龙格-库塔(Runge Kutta, R-K)积分得到平滑航迹。最后,以倾斜转弯面对称远程超音速导弹为例,通过攻防仿真实验,证明该方法规划生成的航迹满足各种约束条件且具有较高突防概率。     

一类并行块隐式方法

本文针对多处理机系统构造了一类并行块隐式方法。具体地给出了此类方法的一个具有两个进程和二阶精度的并行计算公式,并证明了该计算公式具有A稳定性,数值结果表明该计算公式对求解刚性常微分方程是有效的。     

基于对角隐式Runge-Kutta公式的无约束优化方法

通过把一个无约束优化问题转化为一个等价的常微分方程,利用二阶半对角隐式Runge Kutta公式构造了求解无约束优化问题的LRKOPT算法。LRKOPT算法具有与IMPBOT方法相似的数值特性,但LRKOPT算法可以看成是最速下降方向与牛顿法方向的非线性组合,而IMPBOT方法为它们两者之间的线性组合。在目标函数为一致凸函数的假设条件下,证明了LRKOPT方法的具有全局收敛和局部超线性收敛性。数值结果表明LRKOPT方法具有很好的数值稳定性并且LRKOPT方法的计算效率优于IMPBOT方法。     

数值求解常微分方程初值问题的一类并行算法

本文给出了一类数值求解常数微分方程初值问题的并行算法，该类并行算法适用于ＭＩＭＤ型多处理机系统，具有良好的收敛性和数值稳定性，此类并行算法对Ｍｉｒａｎｋｅｒ和Ｌｉｎｉｇｅｒ１９６７年提出的一种构造思想做了圆满的解闷。     

求解常微分方程初值问题的并行块隐式Runge—Kutta方法

本文针对多处理机系统构造了一类并行块隐式Runge-Kutta方法。在S=2的情况下,给出了几个具有三阶精度的并行计算公式,并证明了这类公式具有A稳定性,数值结果表明该计算公式对求解刚性常微分方程是有效的。     

一类并行隐式Runge-Kutta公式

本文针对多处理机系统构造了一类并行隐式Runge-Kutta公式,对2级Runge-Kutta公式给出具有4阶精度的公式族,并证明了它们的收敛性,进行稳定性分析。数值例子表明,该公式可以有效地数值求解较广泛类型的常微分方程初值问题。     

A Class of Parallel Implicit Runge-Kutta Formulas

A class of parallel implicit Runge-Kutta formulas is constructed for multiprocessor system. A family of parallel implicit two-stage fourth order Runge-Kutta formulas is given. For these formulas, the convergence is proved and the stability analysis is given. The numerical examples demonstrate that these formulas can solve an extensive class of initial value problems for the ordinary differential equations.     

Combination Method for Parallel Computation in ODEs

Ｃｏｍｂｉｎａｔｉｏｎ　Ｍｅｔｈｏｄ　ｆｏｒ　Ｐａｒａｌｌｅｌ　Ｃｏｍｐｕｔａｔｉｏｎ　ｉｎ　ＯＤＥｓＳｏｎｇＸｉａｏｑｉｕ（ＢｅｉｊｉｎｇＩｎｓｔｉｔｕｔｅｏｆＣｏｍｐｕｔｅｒＡｐｐｌｉｃａｔｉｏｎａｎｄＳｉｍｕｌａｔｉｏｎＴｅｃｈｎｏｌｏｇｙ）Ａ...     

A Class of Explicit Parallel Multistep Runge-Kutta Methods

In this paper, a rather general class of explicit parallel multistep Runge-Kutta methods is constructed for solving initial value problem of ordinary differential equations. Also, the corresponding convergence and stability are analysed. Several parallel computational formulae are given. The numerical experiments, including accuracy, speedup, and efficiency tests show that the methods are efficient.     

一类秩2微分代数系统的并行Runge-Kutta方法

对多处理机系统构造了一类秩 2微分代数系统的并行Runge -Kutta方法。对该类方法给出了阶条件 ,并且研究了收敛性理论。还研究了Runge -Kutta解的存在性和唯一性。已经构造了一系列使定理 4中的假定成立的具体公式 ,并在飞行器的轨道仿真中应用 ,也提出了一些要进一步研究的问题。     

A Class of Parallel Runge-Kutta Methods for Differential-Algebraic Systems of Index 2

1-INTHODUCTIONWeconsiderthedifferential-algebraicsystemwherefandgaresufficientIydifferentiable.Inorderthat(1a)isanindex2problemwesupposethatinaneighbourhoodoftheexactsolution.wefurtherassumethattheinitialva1uesyo,zoareconsistentwith(l),i.e.,(yo,zo)satisfies(lb)andThesolutionof(1a)isdenotedby(y(t),z(t)).Problemsoftheform(l)arefrequentlyencounteredinpractice.Forexamp1e,problem(l)hasthetypicalstructureofacontrolproblemwherethez-variableactsascontrolparameterwhichforcesthesolutionofthediffer…     

具有最大稳定域的实时RK4公式

实时仿真通常要求计算机在一个动态循环中与外部硬件相互作用。外部数据必须从输入信号中采集，并融入到数值积分算法中以计算出状态变量。同样，从积分过程中求得的数据也必须及时被输出。因此，一般希望能选择较大的积分步长，以保证能有足够的计算时间。首先将实时四阶龙格-库塔公式的稳定域问题化为一个约束求极大值问题，利用优化方法，找到了具有最大稳定域的实时四阶龙格-库塔公式应满足的条件。然后，根据局部截断误差与相关系数的关系，将其化为一个约束求极小值问题，并最终导出了一个新的实时RK4公式。该公式不仅具有最大的稳定域，而且局部截断误差函数值也小于其它已知的实时RK4公式。仿真结果表明，该公式是有效可行的。