自激振荡系统matlab仿真课程设计.doc

课程 计算物理和MATLAB课程设计 题目 自激振动系统的MATLAB仿真 专业 姓名 学号 主要内容、基本要求、主要参考资料等 主要内容： 研究范?德?波耳(Van der pol)方程所描述的非线性有阻尼的自激振动系统，其中是一个小的正的参量，是常数。下面简称范?德?波耳方程为VDP方程。在VDP方程中，增加外驱动力项所得到的方程称强迫VDP方程，其中外驱动力的振幅、角频率分别是和。试研究强迫VDP方程的行为。 基本要求： 1.演示VDP方程所描述的系统在非线性能源供给下，从任意初始条件出发都能产生稳定的周期性运动。 2.采用庞加莱映像，演示强迫VDP方程在不同参数下所存在四种吸引子，即周期1吸引子、周期2吸引子、不变环面吸引子和奇怪吸引子。 3.对于强迫VDP方程，在v和w为定值条件下，逐渐增大μ值，将出现周期倍分岔和混沌现象。 主要参考资料： [1] Steven E. Koonin, 秦克诚译. 计算物理学. 北京：高等教育出版社1993. [2] 马文淦等. 计算物理学. 合肥：中国科学技术大学出版社1992. [3] 张志涌. 精通MATLAB6.5. 北京：北京航空航天大学出版社，2003. 我们知道,对线性阻尼振动系统,严格的周期运动只能由受周期性驱动力作用的受迫振动产生;而对非线性系统,有一种自激振动系统,在非振动的能源供给下,它能产生严格的周期运动,这是人们十分感兴趣的现象。自激系统是一个非线性的有阻尼的振动系统,在振动过程中伴随有能量损耗,但系统存在一种机制,使能量能够由非振动能源通过系统本身的反馈调节,及时适量地得到补充,从而产生一个稳定的不衰减的周期运动,这样的振动称为自激振动[1]。 自激振动现象是一种普遍现象。如钟摆、弦乐器以及人的心脏的周期性跳动。活塞发动机的周期性运动等都是利用这种现象来建立不衰减的周期运动；但有些自激振动是十分有害的，这些现象应该设法避免。 1.2 范·德·波耳（Van der pol）方程 在此课程设计中，我们主要研究范?德?波耳(Van der pol)方程： （1-1） 所描述的非线性有阻尼的自激振动系统，其中是一个小的正的参量，是常数。下面简称范?德?波耳方程为VDP方程。 在VDP方程中，增加外驱动力项所得到的方程： （1-2） 称强迫VDP方程，其中外驱动力的振幅、角频率分别是和。 VDP方程所描述的系统在非线性能源供给下，从任意初始条件出发都能产生稳定的周期性运动。 而对于强迫VDP方程，在V和为 定值条件下，逐渐增大，将出现周期倍分岔和混浊现象。 第2章 自激振荡系统的周期性运动 2.1 课题解析 自激系统是一个非线性有阻尼的振动系统，在运动过程中伴随有能量损耗，但系统存在一种机制，使能量能够由非振动的能源通过系统本身的反馈 调节，及时适量地得到补充，从而产生一个稳定的不衰减的周期运动，这样的振动称为自激振动。 对VDP方程，可从机械振动角度理解，是阻尼系数，它是变化的，如果， 则阻尼系数为正，系统将受阻尼，能量将逐渐减少，但如果， 则发生负阻尼，意味着不仅不消耗系统的能量，反而给系统提供能量。此系统能通过自动的反馈调节，使得在一个振动过程中，补充的能量正好等于消耗的能量，从而系统作稳定的周期振动。 2.2 MATLAB程序设计及演示 取方程(1-1)中的，， ，0.66,，0.85，1.08这四个初始条件进行编程，程序详见附录-程序（1）。 2.3 MATLAB演示结果和分析 根据图3.1显示的结果,我们做出以下结论:给出任一初始条件，通过计算机数值求解可以证明它的相轨道都将趋向于一条闭合曲线，这一条闭合曲线，成为极限环，极限环以外的相轨道向里盘旋，而极限环以内的相轨道则向外盘旋，都趋向极限环，说明不论初始情况如何，系统最终都到达以极限环描述的周期性运动[2]。 图3.1 四种初始条件下的轨道 第3章 强迫VDP方程在不同参数下的四种吸引子 下面研究强迫VDP方程的行为，我们同时采用时间历程图、相图、庞加莱映像图来研究系统在不同参数条件下的动力学行为，可以看到存在不同的吸引子，即周期1吸引子，周期2吸引子，不变环面吸引子和奇怪吸引子。 3.1庞加莱映射 为了更清楚地了解运动的形态，庞加莱对连续运动的轨迹用一个截面（叫庞加莱截面）将其横截，那么根