原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
水流量的估计(插值)
* 图12.7 混合时段时间—流速示意图 * 图12.8是用分段线性及三次样条插值方法得到的整个过程的时间—流速函数示意图。 * 表12-4 各时段及一天的总用水量（用水高度） 第1 未供水段 第2 未供水段 第3 供水段 混合时段 全天 拉格朗日插值法 145.6231 258.8664 54.0689 92.1337 550.6921 分段线性插值法 147.1430 258.9697 49.6051 76.4688 532.1866 三次样条插值法 145.6870 258.6547 53.3334 81.7699 539.4450 * 表12-5是对一天中任取的4个时刻分别用3种方法得到的水塔水流量近似值。 时间 6.88 10.88 15.88 22.88 ① 15.98266712348514 33.74260090853463 25.56622418180477 34.70996210551694 ② 14.82724137931034 32.99762626262627 25.44655913978495 25.4715789473684 ③ 15.05278581965820 33.70895535953259 25.54908920557423 29.41733175863551 注：①拉格朗日插值法②分段线性插值法③三次样条插值法 * 2）拟合法 （1）拟合水位—时间函数 从表12-2中的测量记录看，一天有两次供水时段和三次未供水时段，分别对第1，2未供水时段的测量数据直接作多项式拟合，可得到水位函数（注意，根据多项式拟合的特点，此处拟合多项式的次数不宜过高，一般以3～6次为宜）。对第3未供水时段来说，数据过少不能得到很好的拟合。 设t，h分别为已输入的时刻和水位测量记录（由表12.2提供，水泵启动的4个时刻不输入），这样第1未供水时段各时刻的水位可由如下MATLAB程序完成： * t=[0,0.92,1.84,2.95,3.87,4.98,5.90,7.00,7.93,8.97,10.95,12.03, 12.95,13.88,14.98,15.90,16.83,17.93,19.04,19.96,20.84,23.88 24.99,25.66] h=[9.68,9.48,9.31,9.13,8.98,8.81,8.69,8.52,8.39,8.22,10.82,10.50, 10.21,9.94,9.65,9.41,9.18,8.92,8.66,8.43,8.22,10.59,10.35, 10.18]； c1=polyfit(t(1:10),h(1:10),3)； tp1=0:0.1:8.9； x1=polyval(c1,tp1)； plot(tp1,x1) * 图12.9给出的是第1未供水时段的时间—水位拟合函数图形。 * 变量x1中存放了以0.1为步长算出的各个时刻的水位高度。同样地，第2未供水时段时间—水位图可由如下MATLAB程序完成，读者可自己上机运行查看。 c2=polyfit(t(11:21),h(11:21),3)； tp2=10.9:0.1:20.9； x2=-polyval(c2,tp2)； plot(tp2,x2) （2）确定流量—时间函数 对于第1，2未供水时段的流量可直接对水位函数求导 ，程序如下： * c1=polyfit(t(1:10),h(1:10),3)； c2=polyfit(t(11:21),h(11:21),3)； a1=polyder(c1)； a2=polyder(c2)； tp1=0:0.01:8.97； tp2=10.95:0.01:20.84； x13=-polyval(a1,tp1)； x113=-polyval(a1,[0:0.01:8.97])； wgsysl1=100*trapz(tp1,x113)； */计算第1未供水时段的总用水量/* x14=-polyval(a1,[7.93,8.97])； */为下面的程序准备数据/* x23=-polyval(a2,tp2)； x114=-polyval(a2,[10.95:0.01:20.84]) wgsysl2=100*trapz(tp2,x114)； */计算第2未供水时段的总用水量/* x24=-polyval(a2,[10.95,12.03])； */为下面的程序准备数据/* x25=-polyval(a2,[19.96,20.84])； */为下面的程序准备数据/* subplot(1,2,1) plot(tp1,x13*100) */与图12.10单位保持一致/* sub
文档评论（0）