具有增大稳定域的显式线性三阶四步方法

§1.引言 大步长、低精度实时仿真积分算法是根据某些应用领域尤其武器研制实际需要而提出来的。目的一般都是为了检验所设计的装置的操作性质。在微分方程(组)求近似解中采用大步长积分是为了加快积分速度,达到超时要求。而所求解的低精度是基于实时仿真系统的特点。     

一类线性隐式A稳定的单步法

1.引言 对于Stiff方程组初值问题的数值解法,Dahlquist在[1]中引进了 A稳定的概念,并且证明了显式的线性多步法(包括显式的Runge-Kutta方法)不可能是A稳定的.现在已经有许许多多隐式A稳定或Stiff稳定的方法,但绝大多数在数值解的过程中必须解由于隐式方法所产生的非线性方程组,而非线性方程组的求解过程往往又要采用Newton-Raphson迭代方法,因此需要计算方程y’=f(x,y)的右函数f(x,y)的Jacobi矩阵以及与此有关的逆矩阵.本文的主要思想是:既然在数值解过程中要计算f(x,y)的Jacobi矩阵,那么不妨在数值公式中明显的出现f(x,y)的一阶偏导数.我们将A稳定公式     

显式线性三步二阶方法的稳定性讨论

1 引 言 关于常微分方程初值问题 y′=f(t,y), t∈[t_0,T]R, y(t_0)=y_0, y,f∈R~m,m≥1, (1.1)存在许多数值方法,显式线性多步法,特别是Adams—Bashforth方法,由于其具有计算格式     

一类A(α)-稳定降阶的Gear方法

<正>1引言线性多步法是求解常微分方程(组)的初值问题:{y'=f(x,y),a≤x≤b y(a)=y0的重要经典算法之一~([1-29]),其中具有A(α)-稳定性,可用于解刚性常微分方程的线性多步法的文献有:Gear于1968年在[9]中给出了一类Gear(BDF)方法;杨大地等于2008年在[17]中给出了一类YL方法,它的绝对稳定域比同阶Gear方法略大一些;Cash于1980年在[15]中给出了一种广义的向后差分方法;孙耿等于1981年在[20]中给出了一类变更的     

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

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

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

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

具有A(α)-稳定的Adams-Moulton公式类型

给出了含参数的3阶3步法的A(α)-稳定的Adams-Moulton类型公式族.同时求出了公式的精确分数形式的系数,阶数和局部截断误差主项系数,计算出了它们的幅角α,最后用对比数值实验验证了公式是稳定的,并且适合于求解刚性常微分方程.     

具有大稳定域的线性多步方法

§1.引言 解常微分方程初值问题:的线性k步方法为 sum from j=0 to k (α_jy_(n j)=h sum from j=0 to k (β_jf_(n j),(2)其中α_0~2 β_0~2≠0,α_k≠0.当β_k≠0时,(2)为隐式k步法;当β_k=0时,(2)为显式k步法. 若将(2)应用于单个方程 y′=λy,Reλ<0,则得差分方程 ρ(E)y_n=μσ(E)y_(?),μ=λh,     

解线性不适定方程隐式迭代法的,γ-步拟残差准则

本文给出了解不适定算子方程隐式迭代法后验选取步数的一类准则,称为,ｒ－步拟残差准则,证明了它们总导致最佳收敛阶,这类准则包含著名的 Ｍｏｒｏｚｏｖ残差准则和 Ｇｆｒｅｒｅｒ 准则作为特例．     

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

考虑带常延迟的延迟积分微分方程线性系统零解的渐近稳定性，本文采用拉格朗日插值的线性多步方法，探讨了系统数值方法的线性稳定性。证明了所有A-稳定且强零-稳定的Pouzet型线性多步方法能够保持原线性系统的延迟不依赖稳定性。     

关于具有增大绝对稳定域的线性多步方法

考虑初值问题 其中, 其Jacobi矩阵为~~     

显式及对角隐式Runge-Kutta方法的非线性稳定性

1.引言 过去,stiff常微分方程初值问题数值方法的稳定性研究,主要集中于讨论方法的稳定区域,目标是针对线性自治系统的.最近十年,直接针对非线性非自治系统的理论研究,所谓非线性稳定性分析,才逐渐发展起来.     

解线性不适定方程隐式迭代法的r-步拟残差准则

本文给出了解不适定算子方程隐式迭代法后验选取步数的一类准则,称为r-步拟残差准则,证明了它们总导致最佳收敛阶.这类准则包含著名的Morozov残差准则和Gfrerer准则作为特例.     

一类A－稳定或L－稳定的经济隐式单块法

一类Ａ－稳定或Ｌ－稳定的经济隐式单块法赵双锁，张国凤（兰州大学数学系）ＡＣＬＡＳＳＯＦＡ－ＳＴＡＢＬＥＯＲＬ－ＳＴＡＢＬＥＥＣＯＮＯＭＩＣＡＬＩＭＰＬＩＣＩＴＳＩＮＧＬＥ－ＢＬＯＣＫＭＥＴＨＯＤＳ￥ＺｈａｏＳｈｕａｎｇ－ｓｕｏ；ＺｈａｎｇＧｕｏ－ｆｅ...     

一类A－稳定或L－稳定的经济隐式单块法

一类Ａ－稳定或Ｌ－稳定的经济隐式单块法赵双锁，张国凤（兰州大学数学系）ＡＣＬＡＳＳＯＦＡ－ＳＴＡＢＬＥＯＲＬ－ＳＴＡＢＬＥＥＣＯＮＯＭＩＣＡＬＩＭＰＬＩＣＩＴＳＩＮＧＬＥ－ＢＬＯＣＫＭＥＴＨＯＤＳ￥ＺｈａｏＳｈｕａｎｇ－ｓｕｏ；ＺｈａｎｇＧｕｏ－ｆｅ...     

线性随机时滞微分方程裂步Milstein方法的收敛性(英文)

给出一个新的求解线性随机时滞微分方程的显式分裂步长Milstein格式.运用ItoTaylor展开式证明该格式相对于已有的求解随机时滞微分方程的分裂步长方法而言具有更好的收敛性.数值实验验证了理论分析的正确性.     

解Stiff常微分方程组初值问题的线性隐式方法

对于Stiff常微分方程组初值问题的数值解,人们为了保证数值解过程误差传播的有界性,经常使用的方法之一是隐式的线性多步法.而在解由隐式线性多步法所产生的非线性方程组时,总是采用Newton-Raphson迭代方法.为此就要给出适当的预估式和计算     

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

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