《水塔水流量估计问题》由会员分享，可在线阅读，更多相关《水塔水流量估计问题（4页珍藏版）》请在金锄头文库上搜索。
1、-水 塔 水 流 量 估 计 问 题一问题描述*居民区有一供居民用水的园柱形水塔，一般可以通过测量其水位来估计水的流量，但面临的困难是，当水塔水位下降到设定的最低水位时，水泵自动启动向水塔供水，到设定的最高水位时停顿供水，这段时间无法测量水塔的水位和水泵的供水量 通常水泵每天供水一两次，每次约两小时.水塔是一个高12.2米，直径17.4米的正园柱按照设计，水塔水位降至约8.2米时，水泵自动启动，水位升到约10.8米时水泵停顿工作 表1 是*一天的水位测量记录，试估计任何时刻包括水泵正供水时从水塔流出的水流量，及一天的总用水量拟合水位时间函数二 流量估计的解题思路拟合水位时间函数拟合水位时间函数
2、1.拟合水位时间函数测量记录看，一天有两个供水时段以下称第1供水时段和第2供水时段，和3个水泵不工作时段以下称第1时段t=0到t=8.97，第2次时段t=10.95到t=20.84和第3时段t=23以后。对第1、2时段的测量数据直接分别作多项式拟合，得到水位函数为使拟合曲线比较光滑，多项式次数不要太高，一般在36由于第3时段只有3个测量记录，无法对这一时段的水位作出较好的拟合。2. 确定流量时间函数对于第1、2时段只需将水位函数求导数即可，对于两个供水时段的流量，则用供水时段前后水泵不工作时段的流量拟合得到，并且将拟合得到的第2供水时段流量外推，将第3时段流量包含在第2供水时段3. 一天总用水
3、量的估计总用水量等于两个水泵不工作时段和两个供水时段用水量之和，它们都可以由流量对时间的积分得到。三算法设计与编程1、拟合第1时段的水位，并导出流量设t，h为已输入的时刻和水位测量记录水泵启动的4个时刻不输入，第1时段各时刻的流量可如下得：1 c1=polyfitt1：10，h1：10，3； %用3次多项式拟合第1时段水位，c1输出3次多项式的系数2a1=polyderc1； % a1输出多项式系数为c1导数的系数3tp1=0：0.1：9；*1=-polyvala1，tp1；% *1输出多项式系数为a1在tp1点的函数值取负后边为正值，即tp1时刻的流量4） 流量函数为：2、 拟合第2时段的水
4、位，并导出流量设t，h为已输入的时刻和水位测量记录水泵启动的4个时刻不输入，第2时段各时刻的流量可如下得：1 c2=polyfit(t(10.9:21),h(10.9:21),3)； %用3次多项式拟合第2时段水位，c2输出3次多项式的系数2 a2=polyder(c2)； % a2输出多项式系数为c2导数的系数3tp2=10.9:0.1:21; *2=-polyval(a2,tp2); % *2输出多项式系数为a2在tp2点的函数值取负后边为正值，即tp2时刻的流量4流量函数为：3、拟合供水时段的流量在第1供水时段t=911之前即第1时段和之后即第2时段各取几点，其流量已经得到，用它们拟合第
5、1供水时段的流量为使流量函数在t=9和t=11连续，我们简单地只取4个点，拟合3次多项式即曲线必过这4个点，实现如下： *1=-polyvala1，8 9； %取第1时段在t=8,9的流量 *2=-polyvala2，11 12； %取第2时段在t=11,12的流量 *12=*1 *2； c12=polyfit8 9 11 12，*12，3； %拟合3次多项式 tp12=9：0.1：11； *12=polyvalc12，tp12； % *12输出第1供水时段各时刻的流量拟合的流量函数为：在第2供水时段之前取t=20，20.8两点的流水量，在该时刻之后第3时段仅有3个水位记录，我们用差分得到流量
6、，然后用这4个数值拟合第2供水时段的流量如下： dt3=difft(22：24)； %最后3个时刻的两两之差 dh3=diffh(22：24)； %最后3个水位的两两之差 dht3=-dh3./dt3； %t(22)和t(23)的流量 t3=20 20.8 t(22) t(23)； *3=-polyval(a2，t3(1：2)，dht3； %取t3各时刻的流量 c3=polyfitt3，*3，3； %拟合3次多项式 tp3=20.8：0.1：24； *3=polyvalc3，tp3；% *3输出第2供水时段外推至t=24各时刻的流量拟合的流量函数为：4、一天总用水量的估计第1、2时段和第1、2
7、供水时段流量的积分之和，就是一天总用水量虽然诸时段的流量已表为多项式函数，积分可以解析地算出，这里仍用数值积分计算如下： y1=0.1*trapz(*1)； %第1时段用水量仍按高度计，0.1为积分步长 y2=0.1*trapz(*2)； %第2时段用水量 y12=0.1*trapz(*12)； %第1供水时段用水量 y3=0.1*trapz(*3)； %第2供水时段用水量 y=y1+y2+y12+y3*237.8*0.01； %一天总用水量 计算结果：y1=146.1815, y2=266.4409, y12=48.5004, y3=74.8064，y=1274.45、流量及总用水量的检验计
8、算出的各时刻的流量可用水位记录的数值微分来检验用水量y1可用第1时段水位测量记录中下降高度968-822=146来检验，类似地，y2用1082-822=260检验供水时段流量的一种检验方法如下：供水时段的用水量加上水位上升值260是该时段泵入的水量，除以时段长度得到水泵的功率单位时间泵入的水量，而两个供水时段水泵的功率应大致相等第1、2时段水泵的功率可计算如下： p1=(y12+260)/2； %第1供水时段水泵的功率水量仍以高度计 tp4=20.8：0.1：23； *p2=polyvalc3，tp4； % *p2输出第2供水时段各时刻的流量 p2=(0.1*trapz(*p2)+260)/2.2；%第2供水时段水泵的功率水量仍以高度计计算结果：p1=154.2502 ，p2=142.3670四.计算结果n1,n2y1y2y12y3yp1p2(3,4)146.1815266.440948.500474.80641274.4154.2502142.3670(5,6)146.5150265.541746.131772.60571262.2153.0659141.4479(3,4)流量函数为：(56)流量函数为：画图n=5,6n=3,4. z.
