水箱水流量问题第二十章建立数学建模案例分析.docx
《水箱水流量问题第二十章建立数学建模案例分析.docx》由会员分享，可在线阅读，更多相关《水箱水流量问题第二十章建立数学建模案例分析.docx（10页珍藏版）》请在冰豆网上搜索。
水箱水流量问题第二十章建立数学建模案例分析
第一十五章
综合实验
[学习目的]
1．学习对数学知识的综合运用；
2．学习数学建模——数学应用的全过程；
3．培养实际应用所需要的双向翻译能力。
工科数学而言，学习数学的最终目的应落实在数学的实际应用上，尽管数学也应将训练学生的抽象思维能力为目的，但这也许作为课堂教学的重要内容更为实际可行些，数学实验应注重学生对数学的应用能力——数学建模能力的培养、注意科学研究方法上的培养。
§15.1水箱水流量问题
[学习目标]
1.能表述水箱水流量问题的分析过程；
2.能表述模型的建立方法；
3.会利用曲线拟合计算水箱的水流量；
4.会利用Mathematica进行数据拟合、作图和进行误差估计。
一、问题
许多供水单位由于没有测量流入或流出水箱流量的设备，而只能测量水箱中的水位。
试通过测得的某时刻水箱中水位的数据，估计在任意时刻（包括水泵灌水期间）t流出水箱的流量f（t）。
给出下面原始数据表，其中长度单位为E（1E=）。
水箱为圆柱体，其直径为57E。
时间（s）
水位（10–2E）
时间（s）
水位（10-2E）
0
3175
44636
3350
3316
3110
49953
3260
6635
3054
53936
3167
10619
2994
57254
3087
13937
2947
60574
3012
17921
2892
64554
2927
21240
2850
68535
2842
25223
2795
71854
2767
28543
2752
75021
2697
32284
2697
79254
泵水
35932
泵水
82649
泵水
39332
泵水
85968
3475
39435
3550
89953
3397
43318
3445
93270
3340
假设：
1．影响水箱水流量的唯一因素是该区公众对水的普通需求；
2．水泵的灌水速度为常数；
3．从水箱中流出水的最大流速小于水泵的灌水速度；
4．每天的用水量分布都是相似的；
5．水箱的流水速度可用光滑曲线来近似；
6．当水箱的水容量达到×103g时，开始泵水；达到×103g时，便停止泵水。
二、问题分析与建立模型
1.引入如下记号：
V——水的容积；
Vi——时刻ti（h）水的容积（单位G，1G=（升）；
f（t）——时刻ti流出水箱的水的流速，它是时间的函数（G/h）；
p——水泵的灌水速度（G/h）。
根据要求先将上表中的数据做变换，时间单位用小时（h），水位高转换成水的体积（V=πR2h），得下表。
时间（h）
水量（103G）
时间（h）
水量（103G）
0
/
/
/
/
注：
第一段泵水的始停时间及水量为
t始=（h），v始=χ103（G）
t末=（h），v末=χ103（G）
第二段泵水的始停时间及水量为
t始=（h），v始=χ103（G）
t末=（h），v末=χ103（G）
2.由于要求的是水箱流量与时间的关系，因此须由上表的数据计算出相邻时间区间的中点及在时间区间内水箱中流出的水的平均速度：
平均流速=（区间左端点的水量―区间右端点的水量）/区间中点值
得下表：
时间区间的中点值（h）
平均流量（×103G/h）
/
/
/
/
/
/
做出散点图如图15-1。
图15-1散点图
从图中可以看出数据分布不均匀，局部紧密，因此不能采用插值多项式处理数据，而用曲线拟合的最小二乘法。
三、计算过程
1．算法：
第1步输入数据{xi，yi}；
第2步进行拟合；
第3步作出散点图；
第4步作出拟合函数图；
第5步进行误差估算。
2．实现
在算法步2中使用Fit[]函数，步3、步4使用Plot[]，步5选用Integrate[]函数。
3．误差估计：
误差估算时，由于水泵的灌水速度为一常数，水箱中水的体积的平均变化速度
应近似等于水泵的灌水速度P减去此段时间从水箱中流出的平均速度。
即
此处f（t）在Δt区间的两端点间进行积分。
如果此模型确实准确地模拟了这些数据，那么在不同的灌水周期中，按此模型计算出的水泵灌水速度应近似为常数。
下面通过水泵开始和停止工作的两段区间，即t∈[，]及t∈[，]来进行检验。
第一段：
对应于t始=（h），t末=（h）
水量分别为v始=514800（G），v末=677600（G）
故
ΔV1=677600-514800=162800（G）
Δt1=（h）
=83150（G/h）
第二段：
对应于t始=（h），t末=（h）
水量分别为v始=514800（G），v末=677600（G）
故
ΔV2=677600-514800=162800（G）
Δt2=（h）
=76830（G/h）
P1=83150+
P2=76830+
四、结果分析
通过水泵开始和停止工作的两段时间检验水泵灌水速度应近似为常数；其中由{1，x，x2，x3，…，x8}拟合的函数f（t）所产生的误差为%，由{1，x3，x5，sin，cos}拟合达到%。
由此可见如选择不同的基函数，将得不同的误差。
但是只要基函数选择恰当，所产生的误差也可以保持为相对稳定最小常数来支持该模型。
同时，一旦确定了最佳f（t），我们便可通过Integrate[]函数估算出一天的用水总量，从而根据常规每1000人用水量来推测出该地区的人口数，另外，还可求得水箱的平均流速。
评价
1．优点：
（1）任意时刻从水箱中流出的水速都可通过该模型计算出来；
（2）可推测几天的流速；
（3）可以将该建模过程推广到用电及用气的估算上。
2．缺点：
（1）如能知道水泵的抽水速度，就能更准确地估算水泵灌水期间水的流速；
（2）通过考虑体积测量的差异建模，该作法包含着某种不准确性。
源程序：
L={{，}，{，}，{，}，{，}，
{，}，{，}，{，}，{，}，
{，}，{，}，{，}，{，}，
{，}，{，}，{，}，{，}，
{，}，{，}，{，}，{，}，
{，}}
fx=Fit[L，{1，x^3，x^5，Sin[]，Cos[]}，x]
graph1=ListPlot[L，DisplayFunction→Identity]
graph2=ListPlot[fx，{x，，}，DisplayFunction→Identity]；
Show[graph1，graph2，DisplayFunction→$DisplayFunction，
PlotRange→All]
图15-2水箱水流量拟合图
v1=677600-514800；
t2=；
m1=v1/t1；
v2=677600-514800；
t1=；
m2=v2/t2；
p1=m1+Integrate[fx，{x，，}]/t1
p2=m2+Integrate[fx，{x，，}]/t2
%=（p1-p2）p2
运行结果为：
习题
1．试将水箱水流量问题的建模方法推广到闭路电视的普及预测模型，下表列出了美国自1952年至1978年闭路电视的统计数据。
年
家庭拥有电视数
家庭拥有闭路电视数
闭路电视百分比
电视台数目
闭路电视系统数目
1952
15300
14
108
70
1954
26000
65
402
300
1956
34900
300
496
450
1958
41424
450
556
525
1960
45750
850
579
640
1962
48850
1085
571
800
1964
51600
1085
582
1200
1966
53850
1575
613
1570
1968
56670
2800
642
2000
1970
59550
4500
686
2490
1972
—
6000
690
2841
1974
—
8700
694
3158
1976
71460
10800
701
3681
1978
74700
13000
708
3997
根据上表中的数据可以绘制美国家庭采用闭路电视的增长曲线。
请利用已有历史资料来预测未来闭路电视在家庭采用的百分比。
2．某工地用两台正向铲挖土机挖土，用容量8m3的汽车运土，汽车到达的时间间隔相互独立并服从负指数分布，平均每小时到24辆。
装车时间服从正态分布，平均每小时装车15辆。
设每台挖土机每天的使用费为250元，汽车的使用费每天150元。
试模拟该工地的生产过程，并计算：
（1）挖土机的利用率；
（2）每小时的平均挖土量；
（3）挖土机空闲和汽车等待时每天平均的损失费。
3．某杂货店只有一个收款台，顾客到达收款台的间隔是服从均值为的负指数分布，每个顾客的服务时间服从均值为标准差的正态分布。
这里时间的单位是分钟，且服务时间不取负值。
对100个顾客去收款台缴款的排队过程进行仿真。
升级会员 