拉比诺维奇-法布里康特方程

拉比诺维奇-法布里康特方程(Rabinovich-Fabrikant equations)是 1979年苏联物理学家拉比诺维奇和法布里康特提出模拟非平衡介质自激波动的非线性常微分方程组 Rabinovich, Mikhail I.; Fabrikant, A. L. (1979). "Stochastic Self-Modulation of Waves in Nonequilibrium Media". Sov. Phys. JETP 50: 311. Bibcode:1979JETP...50..311R.：




其中 α, γ 是控制系统的参数.

Danca and ChenDanca指出由于拉比诺维奇-法布里康特方程包含平方项，因此比较难以分析，即便选择的参数相同，但由于求解微分方程组的步骤的不同也会导致不同的吸引子。

数值解
利用Maple中以龙格－库塔法rkf45为核心的软件包odeplot和plot、seq可以得出拉比诺维奇-法布里康特方程的数值解的3D动画图，以便观察拉比诺维奇-法布里康特系统随参数γ和时间t的变化：

参数值：α=1.1，γ=0.803..0.917，t=0...130

初始条件：x(0)=-1,y(0)=0,z(0)=0.5

在t<20时，系统表现为自激振动，当t>20,系统进入混沌态。

平衡点

拉比诺维奇-法布里康特系统具有5个双曲线平衡点，一个在原点，4个依赖于系统参数α   和γ ' ：  name="DancaChen"。





其中：



这些平衡点只存在于参数 α and γ > 0 的一些区域。

当α=1.1,γ=0.87
代人上式可得：
 [0,0,0]
  [.46748585798513339859, -2.3530123557983251267, .95430463972208895291]
  [-.46748585798513339859, -2.3530123557983251267, .95430463972208895291]
  [1.3347123182858183570, -.82414763460993508052, .62751354209609286530]
  [-1.3347123182858183570, -.82414763460993508052, .62751354209609286530]

γ = 0.87, α = 1.1
当 γ = 0.87 and α = 1.1，初始条件为(−1, 0, 0.5).Sprott The 关联维数 为 2.19 ± 0.01.Grassberger  李雅普诺夫指数, λ 约为 0.1981, 0, −0.6581  卡普兰 - 约克量纲, DKY ≈ 2.3010Sprott。

γ = 0.1


Danca and RomeraDanca指出当参数 γ = 0.1、 α = 0.98时，系统进入混沌状态，当 α'' = 0.14时，系统进入极限环。

利用Maple中以龙格－库塔法rkf45为核心的软件包odeplot和plot、seq可以得出拉比诺维奇-法布里康特方程的数值解的3D动画图，以便观察拉比诺维奇-法布里康特系统随参数α的变化：

参数值：γ=0.14，t=0...130,

初始条件：x(0)=-1,y(0)=0,z(0)=0.5

参考文献

外部链接
 Weisstein, Eric W. "Rabinovich–Fabrikant Equation."  From MathWorld—A Wolfram Web Resource.
Chaotics Models a more appropriate approach to the chaotic graph of the system "Rabinovich–Fabrikant Equation" 

Category:非线性常微分方程
Category:混沌理论