水箱流量数学建模论文.docx

问题重述： 许多供水单位由于没有测量流入或流出水箱流量的设备，而只能测量水箱中的水位．试通过测得的某时刻水箱中的水位的数据，估计在任意时刻（包括水泵灌水期间）t流出水箱的流量f(t). 模型假设： 除题目中的假设，另加如下假设： 影响水箱流量的唯一因素是该区民众对水的普通需求； 水塔中的水位不影响水流量的大小。可根据物理学的Torricelli定律得出； 水泵起止工作时间由水塔的水位决定； 水塔的水流量与水泵状态独立,并不因水泵工作而增加或减少水流量的大小； 符号说明： t:时间 Y：水量 V:水流量 I：总用水量 模型求解： 首先，将时间换算成小时，水位高度变换成水的体积，如表1所示： 表1：某一时刻水的体积 因为水量达到514×103g时开始泵水；达到667.6×103g时停止泵水。根据表1中的数据，可以近似认为在8.968h时开始第一次泵水，10.954h时第一次泵水停止；在20.839h时开始第二次泵水，22.880h时第二次泵水停止。 根据表1中的数据，计算相邻时间区间中点以及相应的时间段内平均水流量： 记水量为y，时间为t，流量为v： （ti，ti+1）区间内平均流量vi=（yi+1-yi）/（ti+1-ti） 区间首尾点流量为：v0= (3y0-4y1+y2)/(t2-t0) vn= (-3yn+4yn-1-yn-2)/(tn-tn-2) 计算结果如表2所示： 时间中点（h） 水流量（10^3g/h） 0 14.21 0.46 13.49 1.38 11，62 2.4 10.32 3.41 9.76 4.43 9.46 5.44 8.71 6.46 9.46 7.47 8.92 8.45 10.1 8.97 11.99 10.96 20.75 11.49 18.56 12.49 19.71 13.42 18.48 14.43 16.41 15.44 16.6 16.37 15.4 17.38 14.75 18.49 14.62 19.5 15.57 20.4 15.18 20.84 14.32 23.88 16，67 24.44 13.41 25.45 11.83 25.91 8.78 表2:某一时刻的对应流量 将表二中的数据写入名为“data2”的txt文件中，以供计算使用。 方法1：插值与拟合法： 在数据中水泵工作时的流量数据并没有给出，为了模拟拟合的准确性，采用插值拟合方法作出水箱中水的流速图，如图1所示： 对流速积分，计算一天中的总流量，部分计算结果如图2所示： 计算出24小时之内总用水量： I=363.4×103g=364×3.78=1373.65m3 图1：插值与拟合计算出的每个时间水流量 图2：插值与拟合运行结果部分显示 计算的matlab程序如下： clc,clear a=load(data2.txt); tt=a(:,1) dv2=a(:,2); plot(tt,dv2,*)%画出流速的散点图 pp=csape(tt,dv2);%对流量进行插值 tt0=0:0.1:tt(end);%给出插值点 fdv=ppval(pp,tt0);%计算各插值点的流量值 hold on,plot(tt0,fdv)%画出插值曲线 I=trapz(tt0(1:260),fdv(1:260))%计算24h内总流量的数值积分 方法2：灰色理论： 对表2中计算出来的流量，因为是非单调的摆动发展序列，故采用灰色理论中的GM(2,1)模型，预测流量随时间的变化关系： 计算得到的时间影响序列为： x =482.495*exp(0.0239263*t) - 0.0463653*exp(0.276532*t) - 468.238 计算出来的预测数据、残差以及相对误差如表3所示： 表3：灰色理论预测出的流量值 画出原始数据以及预测值的图像如图3所示： 图3：流量原始值以及灰色理论预测值图 由图3和表3可以看出用灰色理论预测得到的流量值与原始值差别较大，主要是由于模型建立不完善引起的，仍需进一步修改与提高。 计算的matlab程序如下： clc,clear a=load(data2.txt);%提取数据 tt=a(:,1); x0=a(:,2) x0=x0 n=length(x0); x1=cumsum(x0) a_x0=diff(x0) z=0.5*(x1(2:end)+x1(1:end-1)); B=[-x0(2:end),-z,ones(n-1,1)]; u=B\a_x0 x=dsolve(D2x+a1*Dx+a2*x=b,x(0)=c1,x(26)=c2); x=subs(x,{a1,a2,b,c1,c2},{u(