[自然科学]水流量的估计.ppt

图12.7 混合时段时间—流速示意图 图12.8是用分段线性及三次样条插值方法得到的整个过程的时间—流速函数示意图。 表12-4 各时段及一天的总用水量（用水高度） 539.4450 81.7699 53.3334 258.6547 145.6870 三次样条插值法 532.1866 76.4688 49.6051 258.9697 147.1430 分段线性插值法 550.6921 92.1337 54.0689 258.8664 145.6231 拉格朗日插值法 全天 混合时段 第3 供水段 第2 未供水段 第1 未供水段 表12-5是对一天中任取的4个时刻分别用3种方法得到的水塔水流量近似值。 29.41733175863551 25.54908920557423 33.70895535953259 15.05278581965820 ③ 25.4715789473684 25.44655913978495 32.99762626262627 14.82724137931034 ② 34.70996210551694 25.56622418180477 33.74260090853463 15.98266712348514 ① 22.88 15.88 10.88 6.88 时间 注：①拉格朗日插值法②分段线性插值法③三次样条插值法 2）拟合法 （1）拟合水位—时间函数 从表12-2中的测量记录看，一天有两次供水时段和三次未供水时段，分别对第1，2未供水时段的测量数据直接作多项式拟合，可得到水位函数（注意，根据多项式拟合的特点，此处拟合多项式的次数不宜过高，一般以3～6次为宜）。对第3未供水时段来说，数据过少不能得到很好的拟合。 设t，h分别为已输入的时刻和水位测量记录（由表12.2提供，水泵启动的4个时刻不输入），这样第1未供水时段各时刻的水位可由如下MATLAB程序完成： * 水流量的估计 12.1 实验目的 12.2 实验内容 美国某州的用水管理机构要求各社区提供以每小时多少加仑计的用水量以及每天所用的总水量。许多社区没有测量流入或流出水塔水量的装置，只能代之以每小时测量水塔中的水位，其误差不超过5%。需要注意的是，当水塔中的水位下降到最低水位L时，水泵就自动向水塔输水直到最高水位H，此期间不能测量水泵的供水量，因此，当水泵正在输水时不容易建立水塔中水位和用水量之间的关系。水泵每天输水一次或两次，每次约2小时。 本实验的主要目的是使学生学会用MATLAB软件进行插值和曲线拟合计算并解决一些具体的实际问题。通过实际问题的解决，使学生了解如何利用曲线插值及曲线拟合解决实际问题的全过程。 2.1 实验问题 试估计任何时刻（包括水泵正在输水时间）从水塔流出的水流量f(t)，并估计一天的总用水量。已知该水塔是一个高为40ft（英尺），直径为57ft（英尺）的正圆柱，表5-1给出了某个小镇一天水塔水位的真实数据，水位降至约27.00ft水泵开始工作，水位升到35.50ft停止工作。（注：1ft（英尺）=0.3048m（米）） 表12.1 某小镇某天水塔水位 3340 93270 3445 43318 3397 89953 3550 39435 3475 85968 水泵开动 39332 水泵开动 82649 水泵开动 35932 水泵开动 79254 2697 32284 2697 75021 2752 28543 2767 71854 2795 25223 2842 68535 2850 21240 2927 64554 2892 17921 3012 60574 2947 13937 3087 57254 2994 10619 3167 53936 3054 6635 3260 49953 3110 3316 3350 46636 3175 0 水位/0.01ft 时间/s 水位/0.01ft 时间/s 2.2 问题分析 本实验所指流量可视为单位时间内流出水的体积。由于水塔是正圆柱形，横截面积是常数，所以在水泵不工作时段，流量很容易根据水位相对时间的变化算出。问题的难点在于如何估计水泵供水时段的流量。 水泵供水时段的流量只能靠供水时段前后的流量经插值或拟合得到。作为用于插值或拟合的原始数据，我们希望水泵不工作时段的流量越准确越好。这此流量大体上可由两种方法计算，一是直接对表12-1中的水量用数值微分算出各时段的流量，用它们拟合其它时刻或连续时间的流量；二是先用表中数据拟合水位一时间函数，求导数即可得到连续时间的流量。 有了任何时刻的流量，就不难计算一天的总用水量。其实，水泵不工作时段的用水量可以由测量记录直接得到，