非线性振动 matlab,用MATLAB分析非线性弹簧振子的振动

在偏离线性关系的恢复力作用下的振动系统的振动称为非线性振动. 非线性恢复力可用多项式表示为 f=- kx- k0x2 - k1x3+…，其中系数 k、k0、k1 均为常数．在许多实际情况中f=- kx- k1x3，式中的立方项保证恢复力 f 对正负位移仍有同一数值，因而振动对 x=0 是对称的.若 k1＞0，则对给定的位移来说恢复力比线性情况要大. 若这恢复力由弹簧产生，则该弹簧称为非线性硬弹簧. 若 k1＜0，恢复力就比线性情况小，该弹簧就称为非线性软弹簧. 非线性振动方程的求解通常有解析法和数值解法两种．解的解析表达虽然有利于进行理论分析，但求解过程烦琐，而且也只能在一些特殊情况下才能进行.随着计算机的飞速发展，发挥计算机强大的数值计算和图形模拟功能，已成为研究非线性问题的一种重要手段.本文借助 MATLAB 软件分析对称非线性硬弹簧与软弹簧的振动特性，并将其与线性弹簧振子进行比较. 1 势能曲线 对称非线性弹簧振子的势能函数为 Ep=- x 0 乙fdx=- x 0 乙(- kx- k1x3)dx= 1 2 kx2+ 1 4 k1x4 (1) 分别取 k1=0 和±0.008，根据(1)式可分别画出线性弹簧、非线性硬弹簧和软弹簧的势能曲线如图 1 所示.其中线性弹簧的势能曲线为抛物线. 图 1 显示，硬弹簧和软弹簧的势能曲线仍然是对称的. 硬弹簧的势能曲线比线性弹簧的势能曲线陡，曲线的斜率 较大，因而振子振动的速度较大，频率相应较高.软弹簧则相反，曲线的斜率较小，振子振动的速度较小，频率相应较低. 同时，软弹簧的势能存在极大值，说明对于软弹簧只有当总能量小于某一临界值时，才能作周期性的振动． 2 振动曲线与相图 由牛顿第二定律得，非线性弹簧振子的振动方程为mx 咬=- kx- k1x3 (2)这是一个二阶非线性常微分方程，求解析解困难，可以 借助 MATLAB 软件进行数值求解.设 y1=x， y2=x 觶，将(2)式化 为一阶微分方程 dy1 dt =y2 dy2 dt =- k m y1- k1 m y1 乙乙乙乙乙乙乙乙乙乙乙 3 (3) 设 m=1、k=1、x0=5，分别取 k1=0 和±0.008，利用 MATLAB 微分方程解算指令解微分方程组(3)，并绘出弹簧振子的x- t、v- x图，如图 2、图 3. 用 MATLAB 分析非线性弹簧振子的振动 谢善娟 (宁德师范高等专科学校物理系，福建宁德 352100) 摘 要：利用 MATLAB 软件进行数值计算，画出对称非线性弹簧振子的特性曲线，并对其运动进行分析. 关键词：非线性；弹簧振子；数值计算；MATLAB 中图分类号：O332 文献标识码： A 文章编号：1673- 260X(2009)04- 0007- 02 Vol. 25 No. 4 Apr. 2009 第 25 卷 第 4 期 2009 年 4 月 赤 峰 学 院 学 报( 自 然 科 学 版) Journal ofChifengUniversity(Natural Science Edition) 图 1 势能曲线 图 2 振动曲线 图 3 相图 7· · 由图 3 可知，相图 v- x 为一对称封闭曲线，表明对称非线性弹簧振子做对称周期运动.对于硬弹簧 k1＞0，振动速度大于线性弹簧. 同理对于软弹簧 k1＜0，振动速度小于软弹簧.由图 2 可知，非线性弹簧振子的振动曲线与线性振子的振动曲线略有不同，硬弹簧的振动周期小于线性弹簧，相反软弹簧的振动周期大于线性弹簧.这与上面由势能曲线定性讨论的结果一致. 利用 M