Я пытаюсь создать свой собственный интегратор дискретного времени в Simulink, используя правило Bogacki Shampine.Общая формула для правила (когда оно является только функцией времени):
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).Результаты были близки друг к другу, но не слишком много.Также я не уверен, что я правильно создаю этот интегратор.Так как Богацкий Шампин 3-го порядка.Я думал, что должен был использовать задержку в 3 единицы, но мне было достаточно 2.
Как я могу улучшить это или создать еще одну, чтобы получить более точные результаты?