1、 任意t时刻流出水箱的流量水泵的泵水速度四、模型建立与求解4.1模型分析问题要求是分析水箱流量与时间的关系,因此我们需要得到具体时间点所对应的流量数据,由于原始数据中只有一个时间段所对应的水量变化值,于是我们用一个时间段的平均流量作为该时间段时间中点所对应的流量值,然后再通过曲线插值拟合得出时间和流量的函数关系式。4.2数据处理首先我们要将表中数据换算为标准单位制,其中:时间用小时(h)、水箱水量用加仑(G)换算公式有:24m,, 用EXCEL公式进行换算,结果如表一:表一:时间与水量表时间(h)水量()606.112.95639.510.92593.6913.88622.321.845831
2、4.98604.572.95571.5515.9589.33.87562.5716.83574.984.98552.0717.93558.765.9544.0619.04542.537.01533.5619.96528.217.93525.3520.84514.858.9722.02泵水9.9822.9610.9323.88663.3710.95677.6824.99648.4812.03657.6425.91637.6平均流量用EXCEL公式进行计算,计算结果为表二:表二:时间中点与平均流量表时间中点(h)平均流量(0.4613.491.3811.622.410.323.419.764.43
3、9.465.448.716.467.478.928.4510.19.48#VALUE!10.4610.9411.4918.5612.4919.7113.4218.4814.4316.1415.4416.616.3715.417.3814.7518.4914.6219.515.5720.415.1821.4322.4923.4224.4413.4125.4511.834.3数据插值在数据中水泵工作时的流量数据并没有给出,为了模型拟合的准确性,我们采用三次样条插值作出水泵泵水时水从水箱流出的流量值(即时间点9.48h、10.46h、10.94h、21.43h、22.49h、23.42h对应的流量值
4、),并作出时间流量的散点图(程序见附录一),散点如图一;表三:插值所得空缺的流量值流量(G)127511571117117147141449014190图一:时间中点平均流量散点图4.4曲线拟合从图一中可以看出数据分布不均匀,局部紧密,因此不能采用插值多项式进行拟合,而应用曲线拟合的最小二乘法。我们选择作为基函数拟合得到时间流量的函数关系式(程序见附录二),拟合图像如图二:图二:函数拟合图像得到的时间流量函数式为:其中 (h)。该模型拟合得到的是一天内流量随时间的函数式,根据经验如果该天没发生特殊原因每天的用水量和用水时段应该大致一样,因此我们可以将此模型推广到任意一天的用水量随时间的变化关系
5、:若,则。4.5误差估计用MATLAB软件做出拟合结果的残差图(图三):从残差图可以看出,所有数据的残差离零点均较近,且残差的置信区间均包含零点,这说明回归模型能较好的符合原始数据。4.6验证泵水速度为常数如果该模型确实准确地模拟了这些数据,那么在不同的泵水周期中,按此模型计算出的水泵泵水速度应近似为常数;而在水泵工作时水的体积平均变化速度应近似等于水泵我泵水速度减去此段时间水从水箱流出的平均速度。即:此处为在泵水区间两端的定积分。于是通过检验两个泵水期间的泵水速度是否为同一常数来检验该模型的误差(定积分求解过程见附录三)。第一段泵水的始停时间及水量:故第二段泵水的始停时间及水量:相对误差:所
6、以相差仅为7.9%,此结果足以支持该模型。六、模型评价与推广6.1模型的优点1、该模型验证了泵水速度为常数的这一假设。2、如果所给数据准确无误,那么该模型可以计算任何一天任何时刻的用水量。3、该模型应用的数学知识简单易懂,操作方便。4、该模型具有很好的通用性。6.2模型的缺点 1、该模型提供的数据太少,只有一天的观测数据,在实际建模中最好应有不同条件下很多天所采集得数据。 2、由于水泵工作时间不确定导致模型有一定的误差,如果能知道水泵的具体工作时间将会很多地提高模型的精确度。6.3模型的推广此模型根据一些水量和时间的数据得到流量随时间的变化关系,可以预测一天中任意时刻的用水量,这对居民避免用水
7、高峰、供水公司合理分配供水量有很大的帮助,同时该模型还可以推广到用电供给和天然气供给等领域。七、参考文献【1】张磊,毕靖,郭莲英,MATLAB实用教程,北京,人民邮电出版社,2008年【2】亨塞尔曼.美,精通MATLAB,北京,清华大学出版社,2006年八、附录附录一插值clearx=0.46,1.38,2.4,3.41,4.43,5.44,6.46,7.47,8.45,11.49,12.22,13.14,14.43,15.44,16.37,17.38,18.49,19.5,20.4,24.44,25.45;%时间中点y=13490,11620,10320,9760,9460,8710,946
8、0,8920,10100,18560,19710,18480,16140,16600,15400,14750,14620,15570,15180,13410,11830; %流量y0=interp1(x,y,9.48,10.46,10.94,21.43,22.49,23.42,spline) %三次样条插值散点图x1=0.46,1.38,2.4,3.41,4.43,5.44,6.46,7.47,8.45,9.48,10.46,10.94,11.49,12.22,13.14,14.43,15.44,16.37,17.38,18.49,19.5,20.4,21.43,22.49,23.42,24.
9、44,25.45; %插值过后时间y1=13490,11620,10320,9760,9460,8710,9460,8920,10100,12751,15712,17117,18560,19710,18480,16140,16600,15400,14750,14620,14729,14264,14714,14490,14190 ,13410,11830;%插值过后流量plot(x1,y1,*xlabel(时间中点ylabel(平均流量附录二非线性拟合建立m-文件volum.m如下functionyhat=volum(beta,x1) yhat=beta(1)-beta(2)*x1.3+beta
10、(3)*x1.5-beta(4)*cos(0.1*x1)+beta(5)*sin(0.1*x1)输入程序:clcbeta0=112000,200,1,10000,30000;beta,r ,J=nlinfit(x1,y1,volum,beta0)YY,delta=nlpredci(,x1,beta,r ,J);k+,x1,YY,r) %拟合图z=YY-delta, -YY+deltarcoplot(r,z) % 残差图附录三定积分求解syms xy=97566-16.8*x.3+0.013*x.5-83143*cos(0.1*x)-27478*sin(0.1*x);p=int(y,8.97,10.93);p1=vpa(p)q=int(y,20.84,22.96);p2=vpa(q)
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1