我正在尝试使用 Bogacki Shampine 规则在 Simulink 中创建自己的离散时间积分器。该规则的一般公式（当它只是时间的函数时）是：
y(n+1) = y(n) + (t/9)*(2*s1+3*s2+4s3)
在哪里：
s1 = x(n)
s2 = x(n+h/2)
s3 = x(n+3h/4)
这也等于：
y(n) = y(n-1) + (t/9)*(2*s1+3*s2+4s3) ;
在哪里：
s1 = x(n-1)
s2 = x(n-h/2)
s3 = x(n-h/4)
然后，我将结果与使用 ode3 (Bogacki Shampine) 的简单积分器模块进行了比较。结果彼此接近，但并不过分。我也不确定我是否以正确的方式创建了这个积分器。由于 Bogacki Shampine 是 3 阶。我认为我应该使用 3 个单位延迟，但 2 对我来说就足够了。
我该如何改进或创建另一个以获得更准确的结果？