Я пытаюсь использовать Ямпу для какой-то базовой симуляции системы, как в Симулинке. В этом случае я хочу имитировать систему spring и демпфера, представленную этим учебным пособием по simulink. Я написал следующие сигнальные функции для представления системы:
system = time >>> force >>> displacement
force = constant (m * g)
displacement = feedback (-) (velocity >>> integral) (gain $ k / m) 0
velocity = feedback (-) integral (gain $ c / m) 0
Где функция feedback
создает базовый цикл обратной связи и реализуется следующим образом:
feedback op a b b0 = loopPre b0 inner
where inner = arr (uncurry op) >>> a >>> (identity &&& b)
Oh и:
gain x = arr (*x)
С разумными положительными константами я получаю дико неустойчивую систему:
Есть ли что-то явно неправильное в том, как я создаю контуры обратной связи или применяю интеграцию?