估计水塔水流量.ppt

一、问题美国某州的各用水管理机构要求各社区提供以每小时多少加仑计的用水率以及每天所用的总水量，但许多社区并没有测量流入或流出当地水塔的水量的设备，他们只能代之以每小时测量水塔中的水位，其精度在0.5%以内。更为重要的是，无论什么时候，只要水塔中的水位下降到某一最底水位L时，水泵就启动向水塔重新充水直至到某一最高水位H，但也无法得到水泵的供水量的测量数据。因此，在水泵正在工作时，人们不容易建立水塔中的水位与水泵工作时的用水量之间的关系，水泵每天向水塔充水一次或两次，每次约两小时，试估计在任何时刻，甚至包括水泵正在工作的时间内，水从水塔流出的流量f(t)，并估计一天的总用水量。实例（AMCM-91A题）估计水塔的水流量时间水位时间水位时间水位（秒）（0.01英尺）（秒）（0.01英尺）（秒）（0.01英尺）0313735932水泵工作6853528423316311039332水泵工作7185427676635305439435355075021269710619299443318344579154水泵工作13937294746636335082649水泵工作179212892499533260859683475212402850539363167899533397252232797572543087932703340285432752605743012322842697645542927表1、某小镇某天的水塔水位表1给出了某个真实小镇的真实数据,水塔是一个圆形柱体，高40英尺，直径57英尺，通常水塔的水位降至约27英尺时水泵开始向水塔充水，而当通常水塔的水位升至约35.5英尺时水泵停止向水塔充水。二、基本假设1、影响水从水塔流出的流率的唯一因素是公众对水的传统要求。2、水塔中水的水位不影响水流量的大小。（因为物理学的定律指出：水塔的最大水流量与水位高度的平方根成正比，由表中数据有说明最高水位和最底水位的两个流量几乎相等）3、水泵工作的起止时间由水塔的水位决定，水泵工作性能效率总是一定的，没有工作时需维修、使用次数多影响使用效率问题，水泵充水量远大于水塔水流量。1275.35一、问题的重述（略）对离散数据的处理，可以用数据逼近的方法来解决，本问题要想到用数值逼近来建模，数值逼近的方法有很多，如Lagrange插值、分段插值、样条插值、曲线拟合等.本问题分三步：1、先决定所给数据点处的水流量。（数据转换即可）2、找一个水从水塔流出的水流量光滑逼近函数3、处理水泵工作时的充水水流量及一天该镇的总用水量下面介绍样条插值理论4、表中的时间数据准确在一秒以内。5、水塔水流量与水泵状态独立，不因水泵工作而增加或减少水流量的大小。6、水塔水流量曲线可以用一条光滑的曲线来逼近。样条插值分段插值存在着一个缺点,就是会导致插值函数在子区间的端点(衔接处)不光滑,即导数不连续,对于一些实际问题,不但要求一阶导数连续,而且要求二阶导数连续。为了满足这些要求,人们引入了样条插值的概念。所谓“样条”(SPLINE)是工程绘图中的一种工具,它是有弹性的细长木条,绘图时,用细木条连接相近的几个结点,然后再进行拼接,连接全部结点,使之成为一条光滑曲线,且在结点处具有连续的曲率。样条函数就是对这样的曲线进行数学模拟得到的。它除了要求给出各个结点处的函数值外,只需提供两个边界点处导数信息,便可满足对光滑性的不同要求。1、样条函数的定义设f(x)是区间[a,b]上的一个连续可微函数,在区间[a,b]上给定一组基点:a=x0x1x2xn=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)=yi,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+1–xi.对S″(x)连续积分两次,并利用插值条件S(xi)=yi,得到三次样条函数的构造],[)(''111SiiiiiiiixxxhxxMhxxMxxxi,xi+1S”(x)Mi,Mi+1))(())(()()()(iiiiiiiiiiiiiiiixxhMhyxxhMhyhxxMhxxMx6666S1113131因此，只要能求出所有的{Mi}，就能求出样条插值函数S(x).下面考虑Mi的求法],[)()()(1112121622S'iiiiiiiiiiiiiixxxhMMhyyhxxMhxxMx则由连续性S'(xi-)=S'(xi+),(i=1,2,……,n-1)得μiMi-1+2Mi+λiMi+1=di其中1111116)((,iiiiiiiiiiiiihhhyyhyydhhhu上面的方程组有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)，要求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λ0M0d0μ12λ1M1d1.........=.(**).....μn-12λn-1Mn-1dn-1μn2Mndn对于第一型插值问题，取λ0=1，μn=1,对于第二型插值问题，取λ0=0，μn=0对于第三型插值问题，利用周期性，可导出其中'',''nnydyd2200nnnnnndMMMMM2110)'(),'(nnnnnnhyyyhdyhyyhd10101106611110116))((,nnnnnnnnnnhhhyyhyydhhhu以上各组条件与方程组(**)联立，可以解出未知参数M0，M1,……,Mn，然后代入S(x)表达式，即可求得样条函数。上面构造方法中Mi相应于力学中细梁在xi处截面的弯矩，每一个方程中又至多出现相邻的三个Mi，通常称为三弯矩法。总结以上论述，可得求三次样条的步骤为：（1）确定边界条件，判定是第几型插值问题；（2）根据所确定的条件计算各值，形成方程组(**)；（3）解三对角方程组(**)，求得M0，M1，M2,Mn；（4）将求得的Mi值代回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（秒）为充水的截止时间,获得一个额外数据.第二次充水期的一些数据为时间水体积时间水体积时间水体积（小时）（加仑）（小时）（加仑）（小时）（加仑）06061259.9811水泵工作19.03755425540.921159371610.9256水泵工作19.95945282361.843158302610.954267771520.83925148722.949757157112.032865767022.0150水泵工作3.871456259912.954463953422.95816777154.978155209913.875862235223.88006633975.900054408114.982260459824.98696485067.006453396315.903958932525.90836376257.928652537216.82615750088.967851487217.93175587812、数据转换表2(V=r2h)用数学软件绘图如下dttdVtf)()(水流量3、由数据（ti,Vi)产生水流量f(t)的方法有：1.由对水体积的微商数据点直接获得水流量f(t)的近似函数值2.先获得水体积V(t)的近似函数，再对其求导我们利用第一种方法，为获得水体积的微商数据点，选用数值微分公式)(243)()(243)()(,)(1288)(12112112112iiiiiiiiiiiiiiiiiiiiittVVVtfttVVVtftVVttVVVVtf计算右端点：计算左端点：用程序计算出流量函数点集{f(tk)},(给出计算的数据表和散点图)。（略）用样条插值或数据拟合的数学软件，可以得到水流量f(t)的近似函数，这里也记为f(t),这里，样条函数所需的边界条件可以由数值微分公式得出。于是我们得到了水流量的估计函数模型。注：使用样条插值时，得出的水流量f(t)不必给出具体的函数关系式，画出它的图形即可。使用拟合时，得出的水流量f(t)有具体的函数关系式，此时要选好拟合函数类本题可选为8次多项式。4、水泵充水期间的水流量处理水泵充水期间的水流量用平均水流量代替：第一次充水期间充满水的水量）小时加仑量第一次充水的平均用水充水期间的用水量小时）充水的时间加仑）充满水的体积(97576))((1)((9864.19678.89542.10(1628435148726777159542.109678.81119542.109678.811dttfVtpdttftV同理可以得出第二次充水期间的平均水流量p2=91910加仑/小时于是有充水期间的