估计水塔水流量

《估计水塔水流量》由会员分享，可在线阅读，更多相关《估计水塔水流量（36页珍藏版）》请在金锄头文库上搜索。

1、 一、问题美国某州的各用水管理机构要求各社区提供以每小时多少 加仑计的用水率以及每天所用的总水量，但许多社区并没有测 量流入或流出当地水塔的水量的设备，他们只能代之以每小时 测量水塔中的水位，其精度在0.5%以内。更为重要的是，无论 什么时候，只要水塔中的水位下降到某一最底水位L时，水泵 就启动向水塔重新充水直至到某一最高水位H，但也无法得到 水泵的供水量的测量数据。因此，在水泵正在工作时，人们不 容易建立水塔中的水位与水泵工作时的用水量之间的关系，水 泵每天向水塔充水一次或两次，每次约两小时，试估计在任何 时刻，甚至包括水泵正在工作的时间内，水从水塔流出的流量 f(t)，并估计一天的总用水量

2、。实例（AMCM-91A题） 估计水塔的水流量时间 水位 时间 水位 时间 水位（秒） （0.01英尺） （秒） （0.01英尺） （秒） （0.01英尺）0 3137 35932 水泵工作 68535 28423316 3110 39332 水泵工作 71854 27676635 3054 39435 3550 75021 269710619 2994 43318 3445 79154 水泵工作13937 2947 46636 3350 82649 水泵工作17921 2892 49953 3260 85968 347521240 2850 53936 3167 89953 33972522

3、3 2797 57254 3087 93270 334028543 2752 60574 301232284 2697 64554 2927表1、某小镇某天的水塔水位表1给出了某个真实小镇的真实数据,水塔是一个 圆形柱体，高40英尺，直径57英尺，通常水塔的 水位降至约27英尺时水泵开始向水塔充水，而当 通常水塔的水位升至约35.5英尺时水泵停止向水 塔充水。二、基本假设1、影响水从水塔流出的流率的唯一因素是公众对水的传统要求。2、水塔中水的水位不影响水流量的大小。（因为物理学的定律指 出：水塔的最大水流量与水位高度的平方 根成正比，由表中数据有说明最高水位和最底水位的两个流量几乎相等）3、水

4、泵工作的起止时间由水塔的水位决定，水泵工作性能效率总 是一定的，没有工作时需维修、使用次数多影响使用效率问题， 水泵充水量远大于水塔水流量。一、问题的重述（略）对离散数据的处理，可以用数据逼近的方法来解决，本问 题要想到用数值逼近来建模，数值逼近的方法有很多，如 Lagrange插值、分段插值、样条插值、曲线拟合等.本问题分三步：1、先决定所给数据点处的水流量。（数据转换即可）2、找一个水从水塔流出的水流量光滑逼近函数3、处理水泵工作时的充水水流量及一天该镇的总用水量下面介绍样条插值理论4、表中的时间数据准确在一秒以内。5、水塔水流量与水泵状态独立，不因水泵工作而增加 或减少水流量的大小。6、

5、水塔水流量曲线可以用一条光滑的曲线来逼近。样条插值分段插值存在着一个缺点,就是会导致插值函数在子区间的端 点(衔接处)不光滑,即导数不连续,对于一些实际问题,不但要求 一阶导数连续,而且要求二阶导数连续。为了满足这些要求,人 们引入了样条插值的概念。所谓“样条”(SPLINE)是工程绘图中的一种工具,它是有弹 性的细长木条,绘图时,用细木条连接相近的几个结点,然后再 进行拼接,连接全部结点,使之成为一条光滑曲线,且在结点处 具有连续的曲率。样条函数就是对这样的曲线进行数学模拟 得到的。它除了要求给出各个结点处的函数值外,只需提供两 个边界点处导数信息,便可满足对光滑性的不同要求。1、样条函数的

6、定义设f(x)是区间a,b上的一个连续可微函数,在区间a,b上给 定一组基点:a=x0x1x2xn=b 设函数s(x)满足条件(1) s(x)在每个子区间xi , xi+1(i=0,1,2,n-1)上是次数不 超过m的多项式;(2) s(x)在区间a , b上有m-1阶连续导数;则称s(x)是定义在a ,b上的m次样条函数。x0，x1，x2, 称为样条结点,其中x1,xn-1称为内结点, x0 , xn 称为边界结点 。当m=3时,便成为最常用的三次样条函数。设y = f(x)在点 x0，x1，x2, xn的值为y0，y1， y2, yn，若函数S(x)满足下列条件S(xi)=f(xi) =y

7、i , i=0,1,2,n (1.1)则称S(x)为函数f(x)的三次样条插值函数, 简称 三次样条。2、三次样条插值函数构造三次样条插值函数的方法有很多，这里介绍一个常用 的方法：三弯矩插值法记Mi = S(xi), f(xi)= fi= yi ,考虑它在任一区间xi,xi+1上的 形式.根据三次样条的定义可知 ,S(x)的二阶导数S(x)在每一个 子区间xi,xi+1 ( i=0,1,2,n-1)上都是线性函数. 于是在xi,xi+1 上S(x)=Si(x)的二阶导数表示成(1.2) 其中 hi= xi+1xi .对S(x)连续积分两次,并利用插值条件S(xi)= yi ,得到三次样条函数

8、的构造x x i , x i+1S”(x) M i , M i+1因此，只要能求出所有的M i ，就能求出样条插值函数S(x).下面考虑Mi的求法则由连续性S(xi-)= S(xi+) ,(i=1,2,n-1) 得iMi-1+2Mi+iMi+1= di其中 上面的方程组有n-1个方程，但有n+1个变量Mi，故需 两个方程才能求唯一解，为此引入下列边界条件下面介绍几种常用的边界条件第一型边界条件：（可以用数值微分获得端点导数值）已知f(x)在两端点的导数f(a)和f(b) ，要求S(a) = f(a) , S(b) = f(b)第二型边界条件：已知f(x)在两端点的二阶导数f“(a)和f“(b)

9、 ，要求S(a)=M0 = f(a) , S(b)=Mn= f(b)特别当 S(a)= S(b) =0时，S(x)称为自然三次样条第三型边界条件：已知f(x)是以b -a为周期的周期函数 ，要求S(x)满 足周期条件S (a) = S(b) , S(a+)= S(b-) , S(a+)= S(b-)三次样条插值问题加上第i型边界条件称为第i型插值问题（ i，）可以证明第i型插值问题的解是存在且唯 一的。他们对应如下的三对角方程组：2 0 M0 d01 2 1 M1 d1. . . . . . . . = . (*). . . . .n-1 2 n-1 Mn-1 dn-1n 2 Mn dn对于第

10、一型插值问题，取 0=1，n=1,对于第二型插值问题，取0=0，n=0对于第三型插值问题，利用周期性，可导出其中 以上各组条件与方程组(*)联立，可以解出未知参数M0，M1 ,Mn，然后代入S(x) 表达式，即可求得样条函数 。上面构造方法中Mi相应于力学中细梁在xi处截面的弯矩，每 一个方程中又至多出现相邻的三个Mi，通常称为三弯矩法。总结以上论述，可得求三次样条的步骤为：（1）确定边界条件，判定是第几型插值问题；（2）根据所确定的条件计算各值，形成方程组(*)；（3）解三对角方程组(*)，求得M0， M1 ， M2, Mn ；（4）将求得的Mi值代回S(x)的表达式中，从而可求得函数y=f

11、(x)在任一点的近似值S(x)。估计水塔的水流量三、符号约定及说明h：水塔中水位的高度，是时间的函数，单位为英尺V ：水塔中水的体积，是时间的函数，单位为加仑t：时间，单位为小时f ：水塔水流量，是时间的函数，单位为加仑小时p：水泵工作时充水的水流量，是时间的函数，单位为加仑小时四、问题分析与建模采用三次样条插值来做曲线逼近，为 形象化，将表中数据描点画图时间水位1、补充充水的开始和截止数据由假设知水塔的水位降至约27英尺时水泵开始向水塔充水 ，水塔的水位升至约35.5英尺时水泵停止向水塔充水，水泵每天向水塔充水一次或两次，每次约两小时由表1 ，有第一次充水期的一些数据为32284（秒） 26

12、.97 27 （英尺）35932 水泵工作39332 水泵工作39435 （秒） 35.5（英尺）由 39435-32284=7048 1.958（小时）满足每次约两小时的条件可推断在32284（秒） 为充水的开始时间, 在39435（秒） 为充水的截 止时间.75021（秒） 26.97 27 （英尺）推断在75021(秒) 为充水的开始时间79154 水泵工作82649 水泵工作（补充数据35.5英尺） 85968 （秒） 34.75（英尺）（与35.5英尺相差太多）但 85968-75021=3.041（小时）而 82649-75021=7628 2.11889（小时）满足每次约两小时的

13、条件推断在82649（秒） 为充水的截止时间, 获得一个额外数据.第二次充水期的一些数据为时间 水体积 时间 水体积 时间 水体积（小时） （加仑） （小时） （加仑） （小时） （加仑）0 606125 9.9811 水泵工作 19.0375 5425540.9211 593716 10.9256 水泵工作 19.9594 5282361.8431 583026 10.9542 677715 20.8392 5148722.9497 571571 12.0328 657670 22.0150 水泵工作3.8714 562599 12.9544 639534 22.9581 6777154.9

14、781 552099 13.8758 622352 23.8800 6633975.9000 544081 14.9822 604598 24.9869 6485067.0064 533963 15.9039 589325 25.9083 6376257.9286 525372 16.8261 5750088.9678 514872 17.9317 5587812、数据转换表2(V=r2 h )用数学软件绘图如下3、由数据（ti,Vi)产生水流量f(t)的方法有：1.由对水体积的微商数据点直接获得水流量f(t)的近似函数值2.先获得水体积V(t)的近似函数，再对其求导我们利用第一种方法，为获得水体积的微商数据点，选用数 值微分公式用程序计算出流量函数点集f(tk), (给出计算的数据 表和散点图)。（略）用样条插值或数据拟合的数学软件，可以得到水流 量f(t)的近似函数，这里也记为f(t),这里，样条函数