利用Matlab进行弹簧振子运动状态的仿真（打印版）.pdf

利用Matlab进行弹簧振子运动状态的仿真第1页共8页

利用Matlab进行弹簧振子运动状态的仿真

汪家培（安庆师范学院物理与电气工程学院安徽安庆246001）

指导教师：张杰

摘要:简谐振动是最基本的振动，弹簧振子是最典型的简谐振动。本文从牛顿第二定律出发推导

了弹簧振子的运动学方程，从而进一步给出由两相同轻质弹簧和一个质量为m的小球组成的非线性弹

簧振子系统的运动学方程，利用Matlab建立了该系统的仿真模型，研究了在不同条件下系统的的运动

图像、能量图和相图，详细分析了此非线性系统的特点。

关键词：弹簧振子，能量图，相图，运动图，计算机仿真

1.引言

[1-2]

普通物理力学中讲振动时，常以水平弹簧振子和单摆为例作为问题的切入点，它们只是非线性

振动的近似(胡克定律有一定的适用范围)。所谓非线性振动是指恢复力与位移不成线性比例或阻尼力与

速度不成线性比例的系统的振动。一般说，线性振动只适用于小运动范围，超过此范围，就变成非线性

振动。非线性系统的运动微分方程是非线性的，不能用叠加原理求解。方程中不显含时间的非线性系统

称为非线性自治系统；显含时间的称为非线性非自治系统。保守非线性自治系统的自由振动仍是周期性

的，但其周期依赖于振幅。对于渐硬弹簧，振幅越大，周期越短；对于渐软弹簧，振幅越大，周期越长。

非保守非线性自治系统具有非线性阻尼，阻尼系数随运动而变化，因而有可能在某个中间振幅下等效阻

尼为零，从而能把外界非振动性能量转变为振动激励而建立起稳定的自激振动（简称自振）。弦乐器和

钟表是常见的自振系统。周期地改变系统的某个参量而激起系统的大幅振动称参变激发非线性问题是当

前物理研究中的热门问题，非线性光学、非线性声学等已成为专门的学科分支。因此在普通物理力学教

学中，适当的介绍一些非线性振动及处理方法，对学生后继课程的学习及今后的研究工作都是有益的。

本文用一个小球、两个弹簧设计了一个单自由度非线性振动装置，并用计算机数值计算出小球的运动图、

能量图和相图，对非线性振动的特点进行了较为系统全面的分析。

2.弹簧振子系统概述

弹簧振子是一种简谐振子。它的典型结构是由一

个一端固定，质量可以忽略的轻弹簧和连在它另一端

（自由端）与一个跨在气垫导轨上的滑块相连，如图

（1）。这样的，以及与此类似的系统称为弹簧振子。

将滑块自弹簧自由伸展时的O点移动一小位移至A

点，然后释放。可观察到滑块作变加速运动，由A点

通过O点到达B点，继而又经过O点回到A点。并

且再度开始与上述过程完全相同的运动，如此往返。

不考虑作用于滑块的空气阻力及气垫的摩擦。将滑块

视作质点，弹簧自由伸展时质点的位置是平衡位置，

以此为坐标原点建立坐标系Ox。X表示质点的位置坐

标，又等于相对与原点的位移，也是弹簧的伸长（压

缩）量。x很小是，力f与x之间成线性关系，即f=-kx，

k是弹簧劲度系数，弹簧弹性力是线性回复力，弹簧

图（1）弹簧振子系统

振子做简谐振动。

根据牛顿第二定律有[7]

利用Matlab进行弹簧振子运动状态的仿真第2页共8页

2

dx

mkx（1）

dt2

k2

用m除上式两端，并令0，上式可写作

m