质点振动
振动是日常生活和工程技术中常见的一种运动形式。利用常系数线性微分方程的理论来讨论有关自由振动和强迫振动的相关问题。利用MATLAB数学软件大致可分四类情况：（1）无阻尼自由振动情况；（2）有阻尼自由振动；（3）无阻尼强迫振动；（4）有阻尼强迫振动求其数值解和解析解；1、MATLAB软件求解微分方程解析解的命令“dsolve”1求通解的命令格式：（’微分方程’’自变量’）注：微分方程在输入时，一阶导数y’应输入Dy，y’’应输入D2y等，D应大写。（I）无阻尼自由振动情况常见的数学摆的无阻尼微小振动方程代码如下：t0pi502piy2si
3t2plottyb
2151050051152
0
1
2
3
4
5
6
7
（III）无阻尼强迫振动用ＭＡＴＬＡＢ操作的代码如下：t00110y3exp01tsi
5t4plottyk
f2，有阻尼自由振动由无阻尼振动的通解可以看出，无阻尼振动是按照正弦规律运动的，摆动似乎可以无限期的进行下去，但事实上，空气从在阻力，在运动时，我们必须把空气阻力考虑在内，所以我们得到有阻尼摆动方程为：记um2
glw2这里
w是正常数，所以：ydsolveD2y2
Dyw2y0t（443）解得yC3expt
w
w12C2expt
w
w12（1）小阻尼情形
w时，方程（443）的通解为：yexp
tc1cosw1tc2si
w1t和前面无阻尼的情形一样，可以把上式的通解改写为一下形式：yAexp
tsi
w1tQ445这里的A，Q为任意常数。用matlab操作得到：
t00110y3exp01tsi
5t4plottyk
如图：
f3
2
1
0
1
2
3
0
1
2
3
4
5
6
7
8
9
10
由（445）可见，摆动的运动不是周期的，振动的幅度随着时间的增加而不断减小。（2）大阻尼情形：你w时；r2r10方程（443）的通解：Yc1expr1tc2expr2t这里的c1c2为任意常数；（3）临界的情形：即，
w的情形，方程（443）的通解解为yexp
tc1c2t这里的，c1c2为任意常数；由MAYLAB绘制图像得：t001100yyexp02t0502tplottyyr
05040302010010203
0
10
20
30
40
50
60
70
80
90
100
3有阻尼自由振动无阻尼自由振动和有阻尼自由振动都属于自由振动，它对应于一个二阶常数齐次线性微分方程。当一个振动系统还经常受到一个外力作用时，这种振动称为强迫振动。
＝Ａｓｉｎ（ｔ＋）＋Ｈ／（２－ｐ２）ｓｉｎ（ｐｔ）
f取Ａ＝２；＝５；ｐ＝３；t0pi502piy2si
5t219si
3tplottyk
25215105005115225
0
1
2
3
4
5
6
7
4有r