估计水塔水流量(精品)

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,估计水塔的水流量,一、问题,美国某州的各用水管理机构要求各社区提供以每小时多少加仑计的用水率以及每天所用的总水量，但许多社区并没有测量流入或流出当地水塔的水量的设备，他们只能代之以每小时测量水塔中的水位，其精度在0.5%以内。更为重要的是，无论什么时候，只要水塔中的水位下降到某一最底水位,L,时，水泵就启动向水塔重新充水直至到某一最高水位,H，,但也无法得到水泵的供水量的测量数据。因此，在水泵正在工作时，人们不容易建立水塔中的水位与水泵工作时的用水量之间的关系，水泵每天向水塔充水一次或两次，每次约两小时，试估计在任何时刻，甚至包括水泵正在工作的时间内，水从水塔流出的流量,f(t),，,并估计一天的总用水量。,实例（,AMCM-91A,题）估计水塔的水流量,时间 水位 时间 水位 时间 水位,（秒）（0.01英尺）（秒）（0.01英尺）（秒）（0.01英尺）,0 3137 35932 水泵工作 68535 2842,3316 3110,39332,水泵工作 71854 2767,6635 3054,39435,3550,75021,2697,10619 2994 43318 3445,79154,水泵工作,13937 2947 46636 3350,82649,水泵工作,17921 2892 49953 3260,85968,3475,21240 2850 53936 3167 89953 3397,25223 2797 57254 3087 93270 3340,28543 2752 60574 3012,32284 2697 64554 2927,表,1、某小镇某天的水塔水位,表1给出了某个真实小镇的真实数据,水塔是一个圆形柱体，高40英尺，直径57英尺，通常水塔的水位降至约27英尺时水泵开始向水塔充水，而当通常水塔的水位升至约35.5英尺时水泵停止向水塔充水。,二、基本假设,1,、影响水从水塔流出的流率的唯一因素是公众对水的传统要求。,2、水塔中水的水位不影响水流量的大小。,（因为物理学的定律指 出：水塔的最大水流量与水位高度的平方根成正比，由表中数据有,说明最高水位和最底水位的两个流量几乎相等）,3、水泵工作的起止时间由水塔的水位决定，水泵工作性能效率总是一定的，没有工作时需维修、使用次数多影响使用效率问题，水泵充水量远大于水塔水流量。,一、问题的重述（略）,对离散数据的处理，可以用数据逼近的方法来解决，本问题要想到用数值逼近来建模，数值逼近的方法有很多，如,Lagrange,插值、分段插值、样条插值、曲线拟合等.,本问题分三步：,1、先决定所给数据点处的水流量。（数据转换即可）,2、找一个水从水塔流出的水流量光滑逼近函数,3、处理水泵工作时的充水水流量及一天该镇的总用水量,下面介绍样条插值理论,4、表中的时间数据准确在一秒以内。,5、水塔水流量与水泵状态独立，不因水泵工作而增加或减少水流量的大小。,6、水塔水流量曲线可以用一条光滑的曲线来逼近。,样条插值,分段插值存在着一个缺点,就是会导致插值函数在子区间的端点(衔接处)不光滑,即导数不连续,对于一些实际问题,不但要求一阶导数连续,而且要求二阶导数连续。为了满足这些要求,人们引入了,样条插值,的概念。,所谓“样条”(,SPLINE),是工程绘图中的一种工具,它是有弹性的细长木条,绘图时,用细木条连接相近的几个结点,然后再进行拼接,连接全部结点,使之成为一条光滑曲线,且在,结点处具有连续的曲率,。样条函数就是对这样的曲线进行数学模拟得到的。它除了要求给出各个结点处的函数值外,只需提供两个边界点处导数信息,便可满足对光滑性的不同要求。,1、样条函数的定义,设,f(x),是区间,a,b,上的一个连续可微函数,在区间,a,b,上给定一组基点:,a=x,0,x,1,x,2,x,n,=b,设函数,s(x),满足条件,(1),s(x),在每个子区间,x,i,x,i+1,(i=0,1,2,n-1),上是次数不超过,m,的多项式;,(2),s(x),在区间,a,b,上有,m-1,阶连续导数;,则称,s(x),是定义在,a,b,上的,m,次样条函数,。,x,0,，x,1,，x,2,称为,样条结点,其中,x,1,x,n-1,称为,内结点,x,0,x,n,称为,边界结点,。当,m=3,时,便成为最常用的,三次样条函数,。,设,y=f(x),在点,x,0,，x,1,，x,2,x,n,的值为,y,0,，y,1,，y,2,y,n,，,若函数,S(x),满足下列条件,S(x,i,)=,f(x,i,)=,y,i,i=0,1,2,n (1.1),则称,S(x),为函数,f(x),的,三次样条插值函数,简称,三次样条,。,2、,三次样条插值函数,构造三次样条插值函数的方法有很多，这里介绍一个常用的方法：,三弯矩插值法,记,M,i,=,S,(,x,i,),f(,x,i,)=f,i,=,y,i,考虑它在任一,区间,x,i,x,i+1,上的形式.根据三次样条的定义可知,S(x),的二阶导数,S,(,x,),在每一个子区间,x,i,x,i+1,(i=0,1,2,n-1),上都是线性函数.,于是在,x,i,x,i+1,上,S(x)=,S,i,(x,),的二阶导数表示成,(1.2),其中,h,i,=,x,i+1,x,i,.,对,S,(,x,),连续积分两次,并利用插值条件,S(,x,i,)=,y,i,得到,三次样条函数的构造,x x,i,x,i+1,S”(x)M,i ,M,i+1,因此，只要能求出所有的,M,i,，就能求出样条插值函数,S(x).,下面考虑,M,i,的求法,则由连续性,S,(x,i,-,)=,S,(x,i,+,),(i=1,2,n-1),得,i,M,i-1,+2M,i,+,i,M,i+1,=,d,i,其中,上面的方程组有,n-1,个方程，但有,n+1,个变量,M,i,，,故需两个方程才能求唯一解，为此引入下列边界条件,下面介绍几种常用的边界条件,第一型边界条件,：（可以用数值微分获得端点导数值）,已知,f(x),在两端点的导数,f(a),和,f(b),，要求,S,(a)=f,(a),S,(b)=f,(b),第二型边界条件,：,已知,f(x),在两端点的二阶导数,f(a),和,f(b),，要求,S,(a)=M,0,=f,(a),S,(b)=,M,n,=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,M,0,d,0,1,2,1,M,1,d,1,.,.=.(*),.,n-1,2,n-1,M,n-1,d,n-1,n,2,M,n,d,n,对于第一型插值问题，取,0,=1，,n,=1,对于第二型插值问题，取,0,=0，,n,=0,对于第三型插值问题，利用周期性，可导出,其中,以上各组条件与方程组,(*),联立，可以解出未知参数,M,0,，M,1,M,n,，,然后代入,S(x),表达式，即可求得样条函数。,上面构造方法中,Mi,相应于力学中细梁在,x,i,处截面的弯矩，每一个方程中又至多出现相邻的三个,M,i,，,通常称为三弯矩法。,总结以上论述，可得求三次样条的步骤为：,（1）确定边界条件，判定是第几型插值问题；,（2）根据所确定的条件计算各值，形成方程组,(*),；,（3）解三对角方程组,(*),，求得,M,0,，M,1,，M,2,M,n,；,（4）,将求得的,M,i,值代回,S(x),的表达式中，,从而可求得函数,y=f(x),在任一点的近似值,S(x),。,估计水塔的水流量,三、符号约定及说明,h,：,水塔中水位的高度，是时间的函数，单位为英尺,V：,水塔中水的体积，是时间的函数，单位为加仑,t：,时间，单位为小时,f：,水塔水流量，是时间的函数，单位为加仑,小时,p：,水泵工作时充水的水流量，是时间的函数，单位为加仑,小时,四、问题分析与建模,采用三次样条插值来做曲线逼近，为形象化，将表中数据描点画图,时间,水位,1、补充充水的开始和截止数据,由假设知水塔的水位降至约27英尺时水泵开始向水塔充水，,水塔的水位升至约35.5英尺时水泵停止向水塔充水，,水泵每天向水塔充水一次或两次，每次约两小时,由表1，有第一次充水期的一些数据为,32284,（秒）26.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（小时）满足,每次约两小时的条件,推断在,82649（秒）为充水的截止时间,获得一个额外数据.,第二次充水期的一些数据为,时间 水体积 时间 水体积 时间 水体积,（小时）（加仑）（小时）（加仑）（小时）（加仑）,0 606125 9.9811 水泵工作 19.0375 542554,0.9211 593716,10.9256,水泵工作 19.9594 528236,1.8431 583026,10.9542,677715,20.8392,514872,2.9497 571571 12.0328 657670,22.0150,水泵工作,3.8714 562599 12.9544 639534,22.9581,677715,4.9781 552099 13.8758 622352,23.8800,663397,5.9000 544081 14.9822 604598 24.9869 648506,7.0064 533963 15.9039 589325 25.9083 637625,7.9286 525372 16.8261 575008,8.9678 514872 17.9317 558781,2、数据转换,表2,(,V,=,r,2,h),用数学软件绘图如下,3、,由数据（,t,i,V,i,),产生水流量,f(t),的,方法有,：,1.由对水体积的微商数据点直接获得水流量,f(t),的,近似函数值,2.先获得水体积,V(t),的近似函数，再对其求导,我们利用第一种方法，为获得水体积的微商数据点，选用数值微分公式,用程序计算出流量函数点集,f(t,k,),(,给出计算的数据表和散点图)。（略）,用样条插值或数据拟合的数学软件，可以得到水流量,f(t),的,近似函数，这里也记为,f(t),这里，样条函数所需的边界条件可以由数值微分公式得出。于是我们得到了水流量的估计函数模型。,注：使用样条插值时，得出的水流量,f(t),不必给出具体的函数关系式，画出它的图形即可。使用拟合时，得出的水流量,f(t),有具体的函数关系式，此时要选好拟合函数类本题可选为8次多项式,。,4、水泵充水期间的水流量处理,水泵充水期间的水流量用平均水流量代替：,第一次充水期间充满水的水量,同理可以得出第二次充水期间的平均水流量,p,2,=91910,加仑/小时,于是有充水期间的平均水流量：,p=(p,1,+p,2,)/2=94743,加仑/小时,五、,一天的总用水量,为确定模型的准确性，做如下检验：,检验2：利用给定的数据检验（对非充水期间的用