用matlab做单摆,用Matlab软件求单摆的运动

用M atlab 软件求单摆的运动

!

罗

颖，罗兴垅

(赣南师范学院物理与电子信息科学系，江西赣州341000)

摘

要：直接从单摆的运动微分方程出发，应用M atlab 软件，

求出了线性单摆阻尼振动的解析解和非线性单摆无阻尼振动、阻尼振动的数值解；画出了单摆的振动曲线和相图.

关键词：单摆；阻尼；非线性；数值解；M atlab 软件中图分类号：0439

文献标识码：A

文章编号：1004-8332(2005)03-0090-03

1

线性单摆的无阻尼振动

一根不会伸缩、长度为l 的细线，

上端固定(或一根刚性轻杠，上端与无摩擦的铰链相连)，下端悬挂一质量为m 的小球就构成一单摆.当摆角 很小时，si n "0，

单摆作简谐振动，令s =l ，其运动方程为s =A cos ( 0t + 0)(1)式中 0=g ／#

l ，称为单摆的固有圆频率，其周期为T 0=2 ／ 0=2 l ／#

g (2)式(1)

对时间t 求一阶导数，得小球的速度O =-A 0si n ( 0t + 0)(3)

由式(1)与式(3)

得s ()A 2

+O A ()

2

=1(4)

上式表明单摆的速度-位移曲线图

(通常