水塔水流量的估计

《水塔水流量的估计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数学实验MATLAB水塔水流量的估计一．实验问题某居民区有一供居民用水的圆柱形水塔，一般可以通过测量其水位来估计水的流量。但面临的困难是，当水塔水位下降到设定的最低水位时，水泵自动启动向水塔供水，到设定的最高水位时停止供水，这段时间无法测量水塔的水位和水泵的供水量。通常水泵每天供水一次，每次约2h。水塔是一个高为12.2m，直径为17.4m的正圆柱。按照设计，水塔水位降至约8.2m时，水泵自动启动，水位升到约10.8m时水泵停止工作。表1是某一天的水位测量纪录（符号“//”表示水泵启动），试估计任何时刻（包括水泵正供水时）从水塔流出的水流量，及一天的总用水量。   表1水

2、位测量纪录时刻（h）00.921.842.953.874.985.907.017.93水位（cm）968948931913898881869852839时刻（h）8.979.9810.9210.9512.0312.9513.8814.9815.90水位（cm）822////108210501021994965941时刻（h）16.8317.9319.0419.9620.8422.0122.9623.8824.99水位（cm）918892866843822////10591035时刻（h）25.91        水位（cm）1018        -11-数学实验MATL

3、AB二．问题分析根据以上数据的形式和以往经验，适合采用线性拟合的方式进行数据处理。对第1、2、3未供水时段可直接进行用五次多项式进行拟合。对第1、2供水时段分别在两端各取两个点用前后时刻的流速拟合得到。结果可以用分段函数表示分为5段，分别是第一未供水时段，第一供水时段，第二未供水时段，第二供水时段，第三未供水时段。得出流速之后再乘以水塔横截面积即得任何时刻与水塔流出水流量的关系，即流速与时间的关系。对流速进行分段积分并求和，即得一天的总水流量。三．程序的设计与求解方法1.数据的单位转换时间（h）水位（m）时间（h）水位（m）09.6812.9510.210.929.481

4、3.889.931.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.90水泵开启10.93水泵开启23.8810.5910.9510.8224.9910.3512.0310.5025.9110.18-11-数学实验MATLAB水塔的横截面积为A=（17.4）^2*pi/4=237.0661(平方米)。2.拟合水位——时间函数（1）对第

5、1未供水时段的数据进行拟合。t=[00.921.842.903.874.985.907.007.938.9710.9512.0312.9513.8814.9815.9016.8317.9319.0419.9620.8423.8824.9925.91]h=[9.689.489.319.138.988.818.698.528.398.2210.8210.5010.219.949.659.419.188.928.668.438.2210.5910.3510.18]f1=polyfit(t(1:10),h(1:10),5);tm1=0:0.1:9.0;y1=polyval(f1,

6、tm1);plot(tm1,y1)（2）对第2未供水时段的数据进行拟合。f2=polyfit(t(11:21),h(11:21),5);tm2=10.9:0.1:20.9;-11-数学实验MATLABy2=polyval(f2,tm2);plot(tm2,y2)3．确定流量——时间函数（1）对第1，2未供水时段的水位求导可得流量，用5次多项式拟合第1，2未供水时段的流速与时间关系曲线。f1=polyfit(t(1:10),h(1:10),5);f2=polyfit(t(11:21),h(11:21),5);b1=polyder(f1);b2=polyder(f2);tm1

7、=0:0.01:9.00;tm2=10.95:0.01:20.84;g1=-polyval(b1,tm1);s1=trapz(tm1,g1);%计算第1未供水时段的总用水量-11-数学实验MATLABg2=-polyval(b2,tm2);s2=trapz(tm2,g2);subplot(1,2,1)plot(tm1,g1)subplot(1,2,2)plot(tm2,g2)（2）为使流速函数在水泵开启时连续，取4个点，用5次多项式拟合得第1供水时段流速与时间关系曲线。q1=-polyval(b1,[7.938.97]);q2