水箱的水流量估计Word文档格式.docx

水箱的水流量估计Word文档格式.docx

《水箱的水流量估计Word文档格式.docx》由会员分享，可在线阅读，更多相关《水箱的水流量估计Word文档格式.docx（11页珍藏版）》请在冰豆网上搜索。

任意t时刻流出水箱的流量

水泵的泵水速度

四、模型建立与求解

4.1模型分析

问题要求是分析水箱流量与时间的关系，因此我们需要得到具体时间点所对应的流量数据，由于原始数据中只有一个时间段所对应的水量变化值，于是我们用一个时间段的平均流量作为该时间段时间中点所对应的流量值，然后再通过曲线插值拟合得出时间和流量的函数关系式。

4.2数据处理

首先我们要将表中数据换算为标准单位制，其中：

时间用小时（h）、水箱水量用加仑（G）

换算公式有：

24m，

用EXCEL公式进行换算，结果如表一：

表一：

时间与水量表

时间（h）

水量（

）

606.1

12.95

639.51

0.92

593.69

13.88

622.32

1.84

583

14.98

604.57

2.95

571.55

15.9

589.3

3.87

562.57

16.83

574.98

4.98

552.07

17.93

558.76

5.9

544.06

19.04

542.53

7.01

533.56

19.96

528.21

7.93

525.35

20.84

514.85

8.97

22.02

泵水

9.98

22.96

10.93

23.88

663.37

10.95

677.68

24.99

648.48

12.03

657.64

25.91

637.6

平均流量

用EXCEL公式进行计算，计算结果为表二：

表二：

时间中点与平均流量表

时间中点（h）

平均流量（

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

7.47

8.92

8.45

10.1

9.48

#VALUE!

10.46

10.94

11.49

18.56

12.49

19.71

13.42

18.48

14.43

16.14

15.44

16.6

16.37

15.4

17.38

14.75

18.49

14.62

19.5

15.57

20.4

15.18

21.43

22.49

23.42

24.44

13.41

25.45

11.83

4.3数据插值

在数据中水泵工作时的流量数据并没有给出，为了模型拟合的准确性，我们采用三次样条插值作出水泵泵水时水从水箱流出的流量值（即时间点9.48h、10.46h、10.94h、21.43h、22.49h、23.42h对应的流量值），并作出时间—流量的散点图（程序见附录一），散点如图一;

表三：

插值所得空缺的流量值

流量（G）

12751

15711

17117

14714

14490

14190

图一：

时间中点—平均流量散点图

4.4曲线拟合

从图一中可以看出数据分布不均匀，局部紧密，因此不能采用插值多项式进行拟合，而应用曲线拟合的最小二乘法。

我们选择

作为基函数拟合得到时间流量的函数关系式（程序见附录二），拟合图像如图二：

图二：

函数拟合图像

得到的时间流量函数式为：

其中

（h）。

该模型拟合得到的是一天内流量随时间的函数式，根据经验如果该天没发生特殊原因每天的用水量和用水时段应该大致一样，因此我们可以将此模型推广到任意一天的用水量随时间的变化关系：

若

则

。

4.5误差估计

用MATLAB软件做出拟合结果的残差图（图三）：

从残差图可以看出，所有数据的残差离零点均较近，且残差的置信区间均包含零点，这说明回归模型

能较好的符合原始数据。

4.6验证泵水速度为常数

如果该模型确实准确地模拟了这些数据，那么在不同的泵水周期中，按此模型计算出的水泵泵水速度应近似为常数；

而在水泵工作时水的体积平均变化速度

应近似等于水泵我泵水速度

减去此段时间水从水箱流出的平均速度。

即：

此处

为

在泵水区间

两端的定积分。

于是通过检验两个泵水期间的泵水速度是否为同一常数来检验该模型的误差（定积分求解过程见附录三）。

第一段泵水的始停时间及水量：

故

第二段泵水的始停时间及水量：

相对误差：

所以

相差仅为7.9%，此结果足以支持该模型。

六、模型评价与推广

6.1模型的优点

1、该模型验证了泵水速度为常数的这一假设。

2、如果所给数据准确无误，那么该模型可以计算任何一天任何时刻的用水量。

3、该模型应用的数学知识简单易懂，操作方便。

4、该模型具有很好的通用性。

6.2模型的缺点

1、该模型提供的数据太少，只有一天的观测数据，在实际建模中最好应有不同条件下很多天所采集得数据。

2、由于水泵工作时间不确定导致模型有一定的误差，如果能知道水泵的具体工作时间将会很多地提高模型的精确度。

6.3模型的推广

此模型根据一些水量和时间的数据得到流量随时间的变化关系，可以预测一天中任意时刻的用水量，这对居民避免用水高峰、供水公司合理分配供水量有很大的帮助，同时该模型还可以推广到用电供给和天然气供给等领域。

七、参考文献

【1】张磊，毕靖，郭莲英，MATLAB实用教程，北京，人民邮电出版社，2008年

【2】亨塞尔曼.美，精通MATLAB，北京，清华大学出版社，2006年

八、附录

附录一

插值

clear

x=[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,9460,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.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（3）*x1.^5-beta（4）*cos（0.1*x1）+beta（5）*sin（0.1*x1）

输入程序：

clc

beta0=[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+delta]

rcoplot（r,z）%残差图

附录三

定积分求解

symsx

y=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）