资源描述：
《水塔水流量估计2009.docx》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验水塔水流量的估计实验目的本次实验的主要目的是让学生会用数学软件进行插值计算并解决一些具体的实际问题。介绍一些经典的插值方法，包括拉格朗日插值法、埃尔米特插值法、分段插值法、三次样条插值法等等。实验内容1实验问题美国某州的各用水管理机构要求各社区提供以每小时多少加仑计的用水率以及每天所用的总水量。许多社区没有测量流入或流出水塔的水量装置,他们只能代之以每小时测量水塔中的水位,其误差不超过5%。更重要的是,当水塔中的水位下降到最低水位L时水泵就启动向水塔输水直到最高水位H,期间不能测量水泵的供水量。因此,当水泵正在输水时不容易建立水塔中水位和用水量之间的关系。水泵每天输水一次
2、或两次,每次约二小时。试估计任何时刻(包括水泵正在输水时间)从水塔流出的水流量f(t),并估计一天的总用水量。已知该水塔是一个高为40英尺（ft），直径为57英尺(ft)的正圆柱，表12.1给出了某个小镇一天水塔水位的真实数据，水位降至约27.00ft水泵开始工作,水位升到35.50ft停止工作。(注：1英尺（ft）=0.3024米（m）)表12-1某小镇某天水塔水位时间/s水位/0.01ft时间/s水位/0.01ft03175466363350331631104995332606635305453936316710619299457254308713937294760574
3、301217921289264554292721240285068535284225223279571854276728543275275021269732284269779254水泵开动35932水泵开动82649水泵开动39332水泵开动8596834753943535508995333974331834459327033402问题分析流量是单位时间内流出水的体积,由于水塔是正圆柱形,横截面积是常数,所以在水泵不工作时段,流量很容易根据水位相对时间的变化率算出。问题的难点在于如何估计水泵供水时段的流量。水泵供水时段的流量只能靠供水时段前后的流量经插值或拟合得到。作为用于插
4、值或拟合的原始数据,我们希望水泵不工作时段的流量越准确越好。这些流量大体上可由两种方法计算,一是直接对表12.1中的水量用数值微分算出各时段的流量,用它们拟合其它时刻或连续时间的流量;二是先用表中数据拟合水位—时间函数,求导数即可得到连续时间的流量。有了任何时刻的流量,就不难计算一天的总用水量。其实,水泵不工作时段的用水量可以由测量记录直接得到,由表12.1中下降水位乘以水塔的截面积就是这一时段的量这个数值可以用来检验数据插值或拟合的结果。模型建立为了表示方便,我们将问题中所给表12-1中的数据全部化为国际标准单位（表12-5）,时间用小时(h),高度用米(m):表12-5一
5、天内水塔水位记录时间(h)水位(m)时间(h)水位(m)09.6812.9510.210.929.4513.889.941.849.3114.989.652.959.1315.909.413.878.9816.839.184.988.8117.938.925.908.6919.048.667.008.5219.968.437.938.3920.848.228.978.2222.02水泵开动9.98水泵开动22.96水泵开动10.93水泵开动23.8810.5910.9510.8224.9910.3512.0310.5025.9110.18模型假设(1)流量只取决于水位差,与水
6、位本身无关，故由物理学中Torricelli定律：从小孔流出的液体的流速正比于水面高度的平方根。题目给出水塔的最低和最高水位分别是8.1648m（27×0.3024）和10.7352m（35.50×0.3024）(设出口的水位为零)，因为sqrt(10.7352/8.1648)≈1.1467，约为1,所以可忽略水位对流速的影响。（2）将流量看作是时间的连续函数，为计算简单,不妨将流量定义成单位时间流出水的高度,即水位对时间变化率的绝对值(水位是下降的),水塔截面积为S=(57*0.3048)^2*π/4=237.8(m^2).得到结果后乘以S即可。流量估计方法首先依照表12-
7、5所给数据,用MATLAB作出时间—水位散点图12.7。Fig.12.7下面我们来计算水箱流量与时间的关系。根据表12-5中数据散点图12.7，一种简单的处理方法为先将表12-5中的数据分为三段，然后对每一段的数据做如下处理：设某段数据为，相邻数据中点的平均流速用下面的公式(流速=（左端点的水位-右端点的水位）/区间长度)：，每段数据首尾点的流速用下面的公式计算：，，用以上公式求得时间与流速之间的数据表12-6如下:表12-6时间(h)流速(cm/h)时间(h)流速(cm/h)029.8912.4931
