李润岐宋馥璇朱家楠赵兴政
目录实验背景... 2 实验目的... 2 文献综述... 2 理想单摆周期公式推导... 2 空气阻力对小球周期的影响... 4 实验设计方法... 5 变量的初步选择... 5 响应变量... 5 控制变量... 5 固定变量... 6 干扰参数... 6 实验器械... 7 实验设计... 8 变量输入... 8 实验实施... 8 数据分析... 10 主效应图和交互作用图... 10 主效应图... 10 交互作用图... 11 一般线性模型... 11 以 m θ l 为控制变量,20T 为反应变量的一般线性模型... 11 以 m θ lnl 为控制变量,ln(20T) 为反应变量的一般线性模型... 12 对模型的验证... 15 参考文献... 15 1
实验背景 单摆 : 质点振动系统的一种, 是最简单的摆 绕一个悬点来回摆动的物体, 都称为摆, 但其周期一般和物体的形状 大小及密度的分布有关 但若把尺寸很小的质块悬于一端固定的长度为 l 且不能伸长的细绳上, 把质块拉离平衡位臵, 使细绳和过悬点铅垂线所成角度小于 5, 放手后质块往复振动, 可视为质点的振动, 其周期 T 只和摆长 l 及当地的重力加速度 g 有关, 即而和质点的质量 形状和振幅的大小都无明显关系, 其运动状态可用简谐振动公式表示, 称为单摆或数学摆 如果振动的角度大于 10, 则振动的周期将随振幅的增加而变大, 就不服从简单的简谐振动规律了 此时单摆的摆动的周期需要由高等数学来做推导 实验目的 本实验的目的在于探究单摆周期经验公式所反映的变量关系的正确性 ; 并发掘在不是完全理想条件下 ( 如存在空气阻力等 ) 影响单摆周期的其他因素, 如摆角大小 摆球质量 摆球大小等 ; 并分析这些新增变量对单摆周期的具体影响 从而给出在非理想条件下, 影响单摆周期长短的一般性结论 文献综述 理想单摆周期公式推导 当单摆在摆角小于 5 的条件下, 单摆运动近似为简谐运动, 其周期为 2
关于单摆运动近似为简谐运动的过程如下 : M = - m * g * l * Sin x. 其中 m 为质量,g 是重力加速度,l 是摆长,x 是摆角 我们希望得到摆角 x 的关于时间的函数, 来描述单摆运动 由力矩与角加速度的关系不难得到, M = J * β. 其中 J = m * l^2 是单摆的转动惯量,β = x''( 摆角关于时间的 2 阶导数 ) 是角加速度 于是化简得到 x'' * l = g * Sin x 我们对上式适当地选择比例系数, 就可以把常数 l 与 g 约去, 再移项就得到化简了的运动方程 x'' + Sin x = 0 因为单摆的运动方程 ( 微分方程 ) 是 x'' + Sin x = 0 (1) 而标准的简谐振动 ( 如弹簧振子 ) 则是 x'' + x = 0 (2) 我们知道 (1) 式是一个非线性微分方程, 而 (2) 式是一个线性微分方程 所以严格地说上面的 (1) 式描述的单摆的运动并不是简谐运动 3
不过, 在 x 比较小时, 近似地有 Sin x x ( 这里取的是弧度制 即当 x 0 时有 Sin x / x = o(1) ) 因而此时 (1) 式就变为 (2) 式, 单摆的非线性的运动被线性地近似为简谐运动 然后说一下为什么是 10 由于 Sin x x 这个近似公式只在角度比较小的时候成立 ( 这一个可以从正弦函数的在原点附近的图象近似看出 ), 所以只有在小角度下 (1) 式化作 (2) 式才是合理的 事实上 5 0.087266 弧度,Sin 5 0.087155, 二者相差只有千分之一点几, 是十分接近的 在低精度的实验中, 这种系统误差可以忽略不计 ( 因为实验操作中的偶然误差就比它大 ) 但如果换成 25, 误差高达百分之三, 就不宜再看成是简谐振动了 由于正弦函数的性质, 这个近似是角度越小, 越精确, 角度越大越不精确 如果角度很大 ( 比如 60 度处, 误差高达 17%), 就完全不能说它是简谐振动了 当单摆大于 5 时, 有关 T 的推导过程较为复杂, 这里只记录有关结论, 相应推导过程请参看附件 单摆周期推导 此时, 单摆的周期服从下式 : T = T0 (1+ sin2 + sin4 + ), 其中 T0 = 2π 从这一公式也可以看出, 当摆角增大时, 单摆的周期 T 是变长的 空气阻力对小球周期的影响 根据查阅的文献, 当小球速度不大时, 其受到的空气阻力正比于小球速度 v, 正比于小球半径 R 而小球质量正比于 R3 因此与回复力方向相反的空气阻力 4
所产生的加速度应反比于 R2 因此在同种材质的条件下, 小球的半径越大, 振 幅在空气阻力作用下衰减得也越快 实验设计方法 变量的初步选择 响应变量单摆摆动 20 个周期的时间和 因为希望反应变量能够准确而且尽可能明显地反应微小的差异, 我们通过多个周期的时间和来放大单摆周期的变化 反应变量单位测量方法与真实周期的关系 单摆 20 次周期 t 秒 (s) 测量单摆 20 次连续摆动的时间 从单摆第 1 次经过最低点时开始, 第 21 次经过最低点时停止 操作员 1 名, 每次由该操作员将摆球移动至相应位臵放开 测量员 3 名, 每次测量取 3 人平均值 由于空气阻力 摩擦力等噪音的存在,20 次摆动的周期一定是逐渐减小的 但是由于单摆单次摆动的周期很难精确测量, 且 20 次摆动的周期差异不大, 因此平均值对真实周期仍然具有很强的描述性 控制变量 控制变量主要是连续变量, 主要针对我们需要研究的影响因素 具体参见下 表 : 控制变量 单位 测量方法 变量取值 对反应变量的预期影响 摆球质量 m 克 (g) 精确测量 7.4; 8.4; 27.5 m 对 T 影响甚微 摆长 l 米 (m) 精确测量 0.20; 0.30; 0.40; 0.50; 0.60 l 与 T 正相关 摆角 θ 度 ( ) 精确测量 10; 20; 30; 40; 50 θ 对 T 有较小影响 5
在后面的实验中, 我们发现 : 由于空气阻力的存在, 摆球质量并不能较好地 反应摆球这一变量的影响 因此, 将这一变量细化为摆球密度和摆球体积两个离 散变量 控制变量 测量方法 变量取值 对反应变量的预期影响 摆球密度 精确测量 塑料 ; 钢 密度越大 t 越大 摆球体积 精确测量 直径 1.25cm; 1.85cm; 2.65cm 体积越大 t 越大 固定变量 我们通过固定某些对反应变量影响及其微小的变量来降低实验难度 如下表 : 固定变量 单位 一般水平 控制方法 对反应变量的影响 室温 t C 25 固定实验时间与地点 影响不大 重力加速度 g m/s 2 9.80 固定实验地点 影响不大 气体波动 \ \ 关门关窗 影响较大 干扰参数 实验中有些变量难以控制或无法侦测, 只能尽量降低其对实验结果的影响 如下表 : 干扰参数 应对策略 对实验结果的影响 操作员 Repetition 较大 单摆与挂钩的摩擦力 使用打磨得较光滑的挂钩 较小 6
实验器械为了确保单摆固定点的稳定, 并尽量避免锥形摆动的形成, 我们利用宿舍墙壁作为背景, 粘贴挂钩作固定点, 并在墙壁上粘贴相应的带有长度和角度标尺的自制背景来方便摆长 摆角的调节 当摆的摆动不在同一平面内即形成锥摆时, 摆球会碰撞墙壁 摆绳的一段穿过摆球以固定, 另一端绕过挂钩被粘贴在墙上, 以方便摆长的变换 真实的实验场景如下图 : 7
实验设计 变量输入 变量名称 控制水平 摆长 0.2m; 0.3m; 0.4m; 0.5m; 0.6m 摆角 10 ;20 ;30 ;40 ;50 摆球质量 7.4g; 8.4g; 27.5g 摆球体积 摆球密度 半径 1.25cm; 1.85cm; 2.65cm 塑料 ; 钢 实验实施摆球共有 4 种 : 小钢球 大钢球 小塑料球 大塑料球, 摆长分为 5 档, 摆角分为 5 档 将每种摆球的摆长摆角搭配所得的 20 个周期的时间总和 t 分别按序测量, 汇总得到 100 个实验数据, 并详细记录每个实验数据所对应的变量 实验数据记录如下 : 摆长 (m) 摆角 小钢球 (8.4g)t 大钢球 (27.5g)t 小塑料球 (7.4g)t 大塑料球 (8.4g)t 0.2 10 17.97 18 18.02 17.94 0.2 20 18.09 18.13 18.06 18.12 0.2 30 18.15 18.22 18.09 18.34 0.2 40 18.34 18.5 18.19 18.44 0.2 50 18.53 18.59 18.31 18.56 0.3 10 22 22 22.03 22 0.3 20 22.1 22.1 21.94 22.12 0.3 30 22.31 22.39 22.03 22.37 0.3 40 22.46 22.47 22.22 22.6 8
0.3 50 22.69 22.66 22.19 22.66 0.4 10 25.31 25.41 25.37 25.42 0.4 20 25.5 25.59 25.41 25.47 0.4 30 25.72 25.75 25.56 25.71 0.4 40 25.87 26.09 25.68 25.97 0.4 50 26.1 26.22 25.66 25.94 0.5 10 28.41 28.44 28.4 28.5 0.5 20 28.53 28.5 28.44 28.57 0.5 30 28.68 28.88 28.53 28.69 0.5 40 28.94 29.19 28.63 28.9 0.5 50 29.06 29.34 28.63 28.97 0.6 10 31 31.13 31.06 31.07 0.6 20 31.12 31.28 31.09 31.15 0.6 30 31.47 31.5 31.13 31.44 0.6 40 31.59 31.88 31.34 31.59 0.6 50 32 32.16 31.28 31.85 9
数据分析 主效应图和交互作用图 主效应图 从主效应图中可以看出, 摆长对于周期的影响最显著, 而摆角和质量对周期 略有影响, 但效果不明显 10
交互作用图 从交互作用图中可以看出, 三个变量之间基本不存在交互作用, 因此在建模 过程中不考虑交互作用 一般线性模型 以 m θ l 为控制变量, 20T 为反应变量的一般线性模型 由 Minitab 对数据拟合生成的残差图如下 : 11
由残差图可以看出, 残差的正态性并不理想, 具有一定的周期性, 因此, 以 m l 为控制变量,20T 为反应变量生成的一般线性模型并不能令人满意 为了使残差更符合正态性, 我们对 20T 和 l 取对数, 又进行了一次一般线性模型的拟合 以 m θ lnl 为控制变量, ln(20t) 为反应变量的一般线性模型 第一次拟合 由 Minitab 对数据拟合生成的残差图如下 : 12
在残差图中, 我们可以看到, 左下角有几个异常值, 经过检查, 发现这几个异常值都是在摆角较大时产生的, 而摆角较大时, 振幅衰减现象严重, 造成观测值与拟合值相比偏小, 这与残差图所展现的结果是一致的 因此, 去掉异常值后, 对数据进行下一步拟合 而 Minitab 的输出中, 三个控制变量的 P 值均为 0, 说明三个控制变量与反应变量之间是存在相关性的 去掉异常值后进行拟合 去掉 3 个异常值后, 由 Minitab 对数据拟合生成的残差图如下 : 13
由残差图可以看出, 残差的正态性是比较令人满意的 此时 Minitab 输出的 结果为 : ln( 周期 ) 的方差分析, 在检验中使用调整的 SS 来源 自由度 Seq SS Adj SS Adj MS F P 质量 1 0.0020 0.0005 0.0005 29.61 0.000 摆角 ( ) 1 0.0026 0.0086 0.0086 464.31 0.000 ln( 摆长 ) 1 3.6137 3.6137 3.6137 195331.28 0.000 误差 93 0.0017 0.0017 0.0000 合计 96 3.6200 S = 0.00430120 R-Sq = 99.95% R-Sq( 调整 ) = 99.95% 项系数系数标准误 T P 常量 3.67735 0.00163 2250.13 0.000 质量 0.000280 0.000051 5.44 0.000 摆角 ( ) 0.000677 0.000031 21.55 0.000 ln( 摆长 ) 0.496143 0.001123 441.96 0.000 R 2 的值为 99.95%, 也显示了模型是比较良好的 于是, 我们可以得到周期 与摆球质量 摆角 摆长的关系为 : 14
ln(20t) = 3.67735 + 0.00028m + 0.000677θ + 0.496143lnl 对模型的验证 l 考虑最简单的单摆周期公式 T 2π, 可以推导 : g ln(20t ) ln(40 ) 1 1 ln l 2 2 0.5ln l c 1 ln 2 l 1 ln 2 g ln g ln(40 ) 令重力加速度 g = 9.8m/s 2 1, 可以得到 c ln g ln(40 ) 3. 69 即: 2 ln(20t) = 3.69 + 0.5lnl 与模型所得比较, 非常接近, 因此可以判定, 所建模型符合理论 参考文献 单摆的周期与摆角的关系 www.docin.com 对单摆振动周期的讨论 卢伟辉 www.docin.com 在空气阻力下的抛物运动 www.docin.com 15