数值模拟混沌现象及频谱分析-复旦大学物理教学试验中心Fudan.PDFVIP

近代物理实验 非线性物理——混沌 数值模拟混沌现象及频谱分析 07300190082 朱海宁 07300190035 董伟男 指导教师：姚红英 （复旦大学物理系教学实验中心） 引言：随着科学技术的发展，物理学家对各种与物理有关的神奇自然现象研究越来越深，尤 其对非线性物理，这一继量子力学和相对论之后又一证明牛顿力学局限性的理论的研究，也 在物理界兴起了一阵浪潮。因此，增加对线性物理的学习和认识，也应该成为我们学习内容 中必不可少的一部分。 摘要：非线性物理与自然界的奇异现象有着密不可分的关系，也是现代物理关心的问题。蔡 氏电路是一个简单的非线性电路，其非线性元件的各线性分段函数又使其在具有混沌特征的 同时又有可分析性，本文主要介绍如何利用自己编写的软件进行对蔡氏电路的混沌现象的图 像展示和进行频谱分析。并且与实际测量的结果进行对比证明数值模拟的可靠性。 关键字：vb matlab 蔡氏电路 数值模拟 频谱分析 非线性 基本原理： Matlab 数值模拟 由于蔡氏电路的方程为非线性方程，我 们可以去分析他它的解析特性，但很难算出 具体的解析表达形式，因此我们采取了用 matlab 编写程序对上式的微分方程进行了 求解，使用ode45 这个函数对逼近它的解 程序设计的主要流程图如下： 图一 蔡氏电路（如图一所示）的主要元件有 可调电阻R （电路方程中以电导G=1/R 做参 数，以下方程求解过程都用G 来表示，而涉 及实验的内容采用R 表示）、电容C 和C 、 1 2 电感L 以及非线性负阻Nr 。它的运行状态 可以用以下方程组来描述：  dU1 C G(U U ) g (U )  1 dt 2 1 1   dU2 C G(U U ) I 2 1 2 L  dt  dIL L U  2  dt 非线性元件Nr 满足的 IV 特性曲线满足一 下关系式： 我们在子程序模块中改变我们所用的参数  4 就可以实现R 调制和C1 调制的情况。 I 7.558310 U U 1.06  4 4 I 3.947710 U 5.7377610 1.06 U 12.6  I 0.0037U 0.04593 U 1.06 自编Vb 程序界面展示  由于外界无源，因此方程初始量为 为了更加形象的表达混沌图像随着 R X(0)= （0,0,0 ） 值变换而变化的整个过程，我们设计了一个 基于之前的matlab 数值模拟的vb 程序。在 1 / 4