悬崖高度估计模型及其检验与评价.doc

悬崖高度估计模型及其检验与评价 在生活中常常会遇到这样的问题：当我们站在悬崖顶上或者是一个高层建筑上观赏美丽风景时，突然想到要是能测量出悬崖的高度就好了。在这种情况下，我们不可能随身须待测量仪器，通常只会须带诸如手表之类的计时器，能否利用现有的工具来测量悬崖或高层建筑的高度呢？ 有一定物理学知识的人很容易想到用扔石头的办法来测量高度，即先分别记录下石头开始降落时刻和听到石头与地面撞击声音时刻，然后用这两个时刻之差来估计高度。具体办法如下： O假设：1、石头在降落过程中只受到空气阻力和重力的作用，空气阻力大小与石头的降落速度成正比。之所以这样假设是因为崖的高度一般都不是很高，石头的降落速度不会很快，根据物理学知识知，物体在空气中低速运动时所受到的空气阻力大小与运动速度成正比。 O 2、石头撞底的时刻与悬崖顶上的人听到撞底声音的时刻一致。由于声音在空气中的传播速度为330米/ 如图所示，以悬崖顶点为原点建立坐标系。记悬崖高度为H（m），石头开始降落的时刻为0 （s），悬崖顶上的人听到撞底声音的时刻为T （s），t时刻石头的坐标为，石头的质量为（kg），重力加速度为（），空气阻力的比列系数为（）。由牛顿运动定律的如下常微分方程初值问题： （1） 易得方程的解为： （2） 其中（1/s），将代入上式即得悬崖高度的估计式 （3） 只要确定了的值，由（3）式即可根据T值估计出悬崖的高度。但是，在模型投入使用前，我们应当先检验一下模型是否合理，即应当先讨论以下三个问题： 模型（3）与不考虑空气阻力所得模型是否一致？ 是否可以忽略空气的阻力？ 模型（3）对参数是否敏感？ 假设条件“忽略声音在空气中的传播时间”对模型结果的影响如何？ 问题1）的讨论 根据常识我们知道，模型（3）应当与不考虑空气阻力使得模型一致，显然T是k的函数，但是，由于T与k的解析表达式难以确定，因此，我们不能用在（3）中令k=0的方式来检验一致性，只能用求极限的方式来检验。注意到得 这说明模型（3）与不考虑空气阻力所得模型一致。 问题2）的讨论 由 （4） 知，随着T的增加，忽略空气阻力带来的误差不断增加，相对误差趋于无穷大，但是，当时，由得 （5） 即此时的相对误差约为。例如，取，当时，，米，米，相对误差为7%，与用（5）式估计的相对误差6.67%非常接近。 问题3）的讨论 由（3）得 将（3）式代入上式得 故有 （6） 当时，将（5）式代入（6）式得 （7） 另一方面，由（6）式知 故有 综上所述知，当时，模型（3）对参数k是稳定的，即当参数k增加1%时，模型结果减少，但是，当T充分大时模型（3）对参数k是不稳定的，即当时，参数k增加1%时，模型结果减少3%。取，当时，，米，米，k的相对误差为10%，H的相对误差为0.65%，与用（7）式估计的0.67%非常接近。 问题4）的讨论 记石头撞底声音在空气中的传播时间为，则石头坠落的实际时间为，注意到声音在空气中的传播速度为330米/秒，我们得如下关系式： （7） 其中悬崖的真实高度为。为了看清楚忽略声音的传播时间对模型结果的影响，我们取定参数值，对不同的T值用Matlab求解方程（7），得到悬崖的真实高度、模型（3）估计的高度和相应的误差： T(s) H(m) (s) (m) 绝对误差(m) 相对误差(％) 1.00 4.82 0.01 4.69 0.14 2.95 2.00 18.98 0.05 17.99 1.00 5.54 3.00 42