我正在尝试使用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个就足够了。
如何改善这一点或创建另一个以获得更准确的结果?