[论文解读] Formulas for the amplitude of the van der Pol limit cycle

José L. López, Saied Abbasbandy+1 more|ArXiv.org|Jun 10, 2008

Fractional Differential Equations Solutions参考文献 17被引用 5

一句话总结

本文提出了一种新颖的解析近似方法，用于计算范德波尔极限环在整个非线性参数 ε 范围内的振幅，采用同伦分析法（HAM）。通过在相平面（x, ẋ）中推导递推公式，实现了振幅误差低于 0.1% 的精度，这是首个从弱非线性到强非线性区域均有效的解析解。

ABSTRACT

The limit cycle of the van der Pol oscillator, $\ddot{x}+ ε(x^2-1) \dot{x} + x =0$, is studied in the plane $(x,\dot{x})$ by applying the homotopy analysis method. A recursive set of formulas that approximate the amplitude and form of this limit cycle for the whole range of the parameter $ε$ is obtained. These formulas generate the amplitude with an error less than 0.1%. To our knowledge, this is the first time where an analytical approximation of the amplitude of the van der Pol limit cycle, with validity from the weakly up to the strongly nonlinear regime, is given.

研究动机与目标

• 开发一种适用于非线性参数 ε 全范围（从弱到强非线性）的范德波尔极限环振幅的解析近似方法。
• 克服在中间区域（0 ≪ ε ≪ ∞）缺乏解析结果的问题，该区域现有摄动法与渐近方法均失效。
• 基于同伦分析法（HAM）构建一个递推的解析框架，以同时获得极限环的振幅与形状。
• 实现振幅的高精度近似，误差低于 0.1%，优于现有数值与渐近方法。
• 推测所推导的二阶公式为真实振幅的上界，从而增强其在理论与实际应用中的价值。

提出的方法

• 直接在相平面（x, ẋ）中应用同伦分析法（HAM），而非时间域，从而简化对极限环形状与振幅的分析。
• 在同伦参数 h 中采用递推摄动方案，生成极限环形状的逐级近似 y₀(x), y₁(x), y₂(x)。
• 通过满足周期性与非线性恢复力 f(x) = x² − 1 的积分条件，推导出振幅修正项 a₀, a₁, a₂。
• 引入参数化 h(ε) = −bε/(c + ε²)，其中 b 与 c 经调优以最小化误差，实现对所有 ε 的全局有效性。
• 利用符号计算（通过 Mathematica）在每一阶次计算 yₙ(x) 与 aₙ(h) 的复杂初等函数表达式。
• 通过将所得振幅公式与高精度 Runge-Kutta 数值解在 ε ∈ (0, 50) 范围内对比，验证了该方法的有效性。

实验结果

研究问题

• RQ1 是否存在一个单一的解析公式，能够对从弱到强非线性范围内的范德波尔极限环振幅进行近似？
• RQ2 当标准摄动法与渐近方法在中间区域失效时，解析近似的振幅可达到多高的精度？
• RQ3 同伦分析法是否能够不仅生成振幅，还能在相平面中生成极限环的完整形状？
• RQ4 能否构建一个递推的解析框架，系统性地提高振幅精度，同时保持对所有 ε 的有效性？
• RQ5 所得公式是否可推测为真实振幅的上界，从而增强其理论与实际价值？

主要发现

• 本文推导出一阶解析振幅近似 a_R1(ε)，在所有 ε > 0 范围内误差低于 0.1%。
• 二阶公式 a_R2(ε) 将误差进一步降低至 0.05% 以下，显著优于现有解析与数值近似方法。
• 推测二阶公式 a_R2(ε) 为真实振幅的上界，与已知实验上界 2.02342 一致。
• 基于 HAM 的递推方法成功捕捉了当 ε 增大时，极限环从近似圆形向弛豫振荡形态的过渡。
• 该方法不仅生成振幅，还完整重构了极限环的形状，其中 y₀(x) + y₁(x) + y₂(x) 逐步提高近似精度。
• 所推导的公式具有全局有效性，而标准摄动级数（仅在 ε → 0 时有效）或渐近展开（仅在 ε → ∞ 时有效）则不具备此特性。