output怎么用_用MatCont分析非线性动力学方程组（二）

在上一篇文章《用MatCont分析非线性动力学方程组（一）》中，笔者介绍了如何下载、安装与启动MatCont，这款基于Matlab开发的非线性动力学分析工具包。并且以最基本的Pitchfork Bifurcation为例，具体介绍了MatCont的使用方法。还没熟悉MatCont最基本环境的小伙伴们，可以随以下传送门先熟悉以下各种操作界面。

RiskSeeker Xu：用MatCont分析非线性动力学方程组（一）​zhuanlan.zhihu.com

本文中，我们将用MatCont来研究一个物理上非常经典的模型，The Duffing Equation，或称杜芬振子。它是在基本线型震动模型的基础上，增加了一个关于回复力的三次项。其应用极其广泛，从力学到电磁学，很多线性模型都是杜芬振子模型忽略掉三次项后的线性化结果。最简单的钟摆震动，是一个很标准的杜芬振子（含负三次项）；小振幅情况下，忽略掉三次项，就可以得到标准的线型简谐振动。

1. 受迫杜芬振子的动力学模型

初始条件

。

1.1. Exact Solutions

虽然有很非常多的论文给出了Duffing Equation的Exact Solution（多数基于雅可比椭圆函数），但是据笔者所知，目前还没有被当做绝对经典的Exact Solution。

1.2. Asymptotic Solutions

一般情况下，都是把由扰动法得到的渐近解作为经典解析解，具体方法是：由典型的Lindstedt Poincare扰动法，或Two Time-scales扰动法，结合Harmonic Balance法，加之考虑有正damping的存在，可比较方便地得到系统达到稳态时的渐进解。

Frequancy-Response：