北京科技大学数学大物理组.docVIP

《Matlab数学实验》大作业 研究课题?? 轻杆单摆振动的周期和运动规律 学 院 数理学院 班 级?? 物理1101 组 长?? 周慧斌??? 组 员? 童 鑫 周 涛 2012年12月 ?? 应用题目及背景 单摆是质点振动系统的一种，是最简单的摆。绕一个悬点来回摆动的物体，都称为摆，但其周期一般和物体的形状、大小及密度的分布有关。但若把尺寸很小的质块悬于一端固定的长度且不能伸长的轻杆（质量可忽略不计）上，把质块拉离平衡位置，放手后质块往复振动，此即轻杆单摆的振动问题。当轻杆和过悬点铅垂线所成角度小于10°时，可视为质点的振动周期 T只和当地的重力加速度g有关，而和质块的质量、形状和振幅的大小都无关系，其运动状态可用简谐振动公式表示，称为单摆或数学摆。如果振动的角度大于 10°，则振动的周期将随振幅的增加而变大，就不成为单摆了。如摆球的尺寸相当大，轻杆的质量不能忽略，就成为复摆（物理摆），周期就和摆球的尺寸有关了。 本文将从严格的物理学角度，应用Matlab软件对各种角度下的轻杆类单摆的周期和运动规律进行模拟分析。此项研究对物理学其它模型的研究，各种科研活动中精密仪器的研制具有重要指导意义。 此项研究解决的主要问题如下： 设想一长为l的轻杆，连接一个质量为m的摆球，形成一个单摆，不计摩擦。建立单摆的周期与角振幅关系的数学模型，并用matlab进行曲线模拟. 编程演示单摆振动的动画，比较单摆振动和简谐振动的规律。 当单摆角振幅的度数为1°到7 °时(间隔为1 °)，将单摆运动的角位置和角速度与简谐振动进行比较。当单摆角振幅的度数为30 °到150 °时(间隔为30 °)，另加179 ° ，同样进行比较。 数学模型的建立 如图所示，设角位置为θ，摆锤的运动方程为 即： 在小角度的情况下，sinθ ≈ θ，可得 ω0为圆频率 振动周期为 可见：在小角振动的情况下，单摆的周期与角振幅无关，这称为单摆的等时性。 摆锤的角速度为ω = dθ/dt，因此 可得 积分得 当t = 0时，ω = 0，θ = θm，可得C = -gcosθm/l。因此角速度大小为 单摆的周期为： 对于任何角振幅θm，通过数值积分和符号积分都能计算周期。 利用半角公式可得 设 并设ksinx = sin(θ/2)，因此 第一类完全椭圆积分定义为 其中 将周期的椭圆积分公式按二项式展开得 可得用无穷级数表示的单摆周期 如果只取常数项，可得单摆小角度的周期T0。如果取前两个正弦项，则得 周期可表示为 其中. 对于一定的角振幅，根据精度在单摆周期的无穷级数中去有限项，计算周期的近似值。 【程序代码】 clear theta=0:179; thm=theta*pi/180; T0=ellipke(sin(thm/2).^2)*2/pi; figure plot(theta,T0,r-) fs=16; title(单摆的周期与角振幅的关系,FontSize,fs) xlabel(\it\theta\rm_m/\circ,FontSize,fs) ylabel(\itT/T\rm_0,FontSize,fs) grid on legend(第一类完全椭圆积分) text(0,2,\itT\rm_0=2\pi(\itl/g\rm)^{1/2},FontSize,fs) 轻杆单摆运动的动画模拟 【程序代码】 clear thetam=input(输入单摆角振幅：); thm=thetam*pi/180; T=ellipke(sin(thm/2)^2)*2/pi; t=linspace(0,T*2*pi); options.RelTol=1e-6; [t,TH]=o