算法工程师老王
-
告别数值发散 - 在Max/MSP gen~中运用RK4方法精确模拟洛伦兹吸引子
玩Max/MSP，特别是gen~的朋友，可能都尝试过模拟一些有趣的动态系统，比如经典的洛伦兹吸引子（Lorenz Attractor）。用简单的欧拉法（Euler method）快速搞个原型出来爽一下是挺方便，但当你开始追求更高的精度，或者在较低采样率（比如你想节省CPU资源时）、系统参数比较极端（临界混沌边缘）的情况下，欧拉法那点儿可怜的精度和稳定性问题就暴露无遗了，搞不好数值直接就飞了。 这时候，就该轮到更高级的数值积分方法出场了。今天咱们就来聊聊怎么在gen~环境里，用大名鼎鼎的四阶龙格-库塔法（RK4）来更精确、更稳定地模拟像洛伦兹吸引子这样的由微分方程定义的动态系...