布鲁塞尔振子(英语:Brusselator)也称布鲁塞尔方程是一组模拟自催化反应的非线性微分方程:[1] 布鲁塞尔振子模拟 d ( x ( t ) ) d t = a − b ∗ x ( t ) + x ( t ) 2 ∗ y ( t ) − x ( t ) {\displaystyle {\frac {d(x(t))}{dt}}=a-b*x(t)+x(t)^{2}*y(t)-x(t)} d ( y ( t ) ) d t = b ∗ x ( t ) − x ( t ) 2 ∗ y ( t ) {\displaystyle {\frac {d(y(t))}{dt}}=b*x(t)-x(t)^{2}*y(t)} 数值解 利用龙格-库塔法可求得布鲁塞尔振子的数值解,并利用Maple绘图[1] 布鲁塞尔振子图 布鲁塞尔振子极限环场矢图 布鲁塞尔振子 3D 图 参考文献 [1]Richard H. Enns George C. McCGuire, Nonlinear Physics with Maple, p199, Birkhauser,1997 外部链接