求解刚性微分方程-梯形法则+逆向微分公式
ode23tb
是TR-BDF2的实现,TR-BDF2是隐式龙格-库塔公式,其第一阶段为梯形规则步,第二阶段为二阶倒向微分公式。通过构造,采用相同的迭代矩阵对两个阶段进行求解。就像ode23s
而且ode23t
,这个求解器可能比ode15s
用于粗公差问题[1],[2].
[1] Bank, R. E., W. C. Coughran, Jr., W. Fichtner, E. Grosse, D. Rose,和R. Smith,“硅器件和电路的瞬态模拟”,IEEE反式。计算机辅助设计, 4(1985),第436-451页。
[2] Shampine, L. F.和M. E. Hosea,“TR-BDF2的分析与实现”,应用数值数学20, 1996年。