1水箱水流量问题.docx

1水箱水流量问题.docx

《1水箱水流量问题.docx》由会员分享，可在线阅读，更多相关《1水箱水流量问题.docx（14页珍藏版）》请在冰豆网上搜索。

1水箱水流量问题

例题

估计水箱的水流量模型

（1991年美国大学生数学建模竞赛的A题.问题中使用的长度单位为E（英尺,1E=30.24cm）,容积单位是G（加仑,1G=3.785L））.

问题的提出：

某些州的用水管理机构需估计公众的用水速度（单位:

G/h）和每天的总用水量.许多供水单位由于没有测量流入或流出量的设备,而只能测量水箱中的水位（误差不超过5%）.当水箱水位低于水位L时,水泵开始工作将水灌入水箱,直至水位达到最高水位H为止.但是依然无法测量水泵灌水流量,因此,在水泵工作时无法立即将水箱中的水位和水量联系起来.水泵一天灌水1~2次,每次约2h.

试估计在任一时刻（包括水泵灌水期间）t流出水箱的流量

并估计一天的总用水量.

表1给出了某镇某一天的真实用水数据.水箱是直径为57E,高为40E的正圆柱体.当水位落到27E以下,水泵自动启动把水灌入水箱;当水位回升至35.5E时,水泵停止工作。

表1

时间/s

水位

E

时间/s

水位

E

0

3316

6635

10619

13937

17921

21240

25223

28543

32284

35932

39332

39435

43318

3175

3110

3054

2994

2947

2892

2850

2795

2752

2697

泵水

泵水

3550

3445

46636

49953

53936

57254

60574

64554

68535

71854

75021

79254

82649

85968

89953

93270

3350

3260

3167

3087

3012

2927

2842

2767

2697

泵水

泵水

3475

3397

3340

模型假设

（1）影响水从水塔流出的流率的唯一因素是公众对水的传统要求．因为附表只给出了某一天（实际是近26小时）水塔的水位数据，并没有对这些数据的产生有影响的因素作出具体说明，我们只能假定所给数据反映了有代表性的一天，而不包括任何特殊情况，如自然灾害、火灾、水塔溢水、水塔漏水等对水的特殊要求．

（2）水塔中水的水位不影响水流量的大小．据物理学的Torricelli（托里查里）定律，水塔最大水流量是与水位的高度的平方根成正比的．针对表8－1所给的数据，最大高度是35.50英尺,最小高度是27.00英尺,所以两个高度的最大水流量之比是

接近于1,所以我们假定水位不影响水流量,类似地,我们假定气候条件、条件变化等也不直接影响水流量.

（3）水泵工作起止时间由水塔的水位决定.我们总是假定水位大约27.00英尺时,水泵就开始工作,直到水位升至大约35.50英尺时停止工作,每次充水时间约为两小时.水泵工作性能、效率总是一定的,不因使用次数多少而变化,水泵工作时不需要维修,也不中途停止工作.当然,水泵充水的水流量远大于水塔的水流量,以保证人们对水的需求.

（4）表1中水位数据取得的时间准确在1秒以内.

（5）水塔的水流量与水泵状态独立,并不因水泵工作而增加或减少水流量的大小.

（6）水塔的水流量曲线可以用一条光滑的曲线来逼近.这时,在每一个数据点,水流量的两阶导数是连续的,因为水的消耗是基于社区公众一天的活动,如洗澡、做饭、洗衣服等,每一个使用者的要求与整个社会的要求相比是微不足道的,而整个社会的需求是不可能同时增加或减少的,由于水的消耗的自然性,可以设想水流量曲线是一条连续光滑的曲线.

问题分析与模型建立

1、为方便起见,记

V----表示水的容积;

----表示时刻

（单位:

h）时水的容积;

------表示流出水箱的水的流速（单位;G/h）,它是时间的函数;

p---------表示水泵的灌水速度（G/h）.

先将表1中数据作变换,时间单位用小时（h）,水位高转换成水的体积，具体数据如表2所示（

单位:

）.

表2

时间/h

水量/

G

时间/h

水量/

G

0.

0.921111

1.84306

2.94972

3.87139

4.97806

5.9

7.00639

7.92861

8.96778

9.98111

10.9256

10.9542

12.0328

606.098

593.69

583.

571.546

562.574

552.074

544.057

533.557

525.349

514.849

/

/

677.685

657.64

12.9544

13.87558

14.9822

15.9039

16.8261

17.9317

19.0375

19.9594

20.8392

22.015

22.9581

23.88

24.9869

25.9083

639.505

622.324

604.571

598.299

574.982

558.756

542.529

528.212

514.849

/

/

663.367

648.477

637.593

根据假定，当水塔水位降至约27.00英尺时，开始充水，而水位升至约35.50英尺时停止充水，从所给数据自然无法知道水泵开始和停止工作的准确时间，但我们发现第一次充水前的最后一个测量数据是32284秒时水位为26.97英尺，可见水泵在32284位秒后不久开始充水。

39332秒时水泵仍在工作，而39435秒时水为35.50英尺，水泵在这两个时刻之间停止了工作，这两个时刻的差距是103秒≈0.028小时，很短的时间，所以我们可以假设水泵开始工作的时间为32284秒，结束工作的时间为39332秒，充水时间约为7048秒≈1.95778小时，符合每次抽水约两小时的假设。

现在再来分析第二次充水期间的数据，充水前最后一个测量数据是75021秒时水位26.97英尺，与第一次充水前水位相同，所以可以假设75012秒后水泵即开始工作，但充水后的第一个数据显示了一个与水泵应该停止工作时的不同的水位，说明水泵停止工作已有一段时间，水泵停止工作的准确时间在82649秒和85968秒之间，但这两个时刻相差近一个小时，更靠近哪能一个时刻？

我们发现82649秒时水泵已工作约82649-75021=7628秒，与第一次充水时间相比，第二次充水时间也约为2小时，由此即知82649秒后水泵即停止了工作，所以我们可选取82649秒≈22.95806小时为第二次充水停止的时间，并可假定此时刻的水位也为约35.50英尺，这就解决了水泵起止工作时间问题。

即有

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

t始=8.968（h），v始=514.8×103（G）

t末=10.926（h），v末=677.6×103（G）

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

t始=20.839（h），v始=514.8×103（G）

t末=22.958（h），v末=677.6×103（G）

2、由于要求的是水箱流量与时间的关系，因此须由上表的数据计算出相邻时间区间的中点及在时间区间内水箱中流出的水的平均速度：

平均流速=（区间左端点的水量―区间右端点的水量）/区间的长度

得下表：

表3

时间区间的中点值/h

平均水流量/

G/h

时间区间的中点值/h

平均水流量/

G/h

0.460556

1.38208

2.39639

3.41056

4.42472

5.43903

6.45319

7.4675

8.44819

9.47444

10.4533

10.9399

11.4935

12.4936

13.471

11.5953

10.3498

9.73471

9.48735

8.69649

9.48974

8.90086

10.1036

/

/

/

18.5833

19.6766

13.4151

14.429

15.4431

16.365

17.3789

18.4846

19.4985

20.3993

21.4271

22.4865

23.419

24.4335

25.4476

18.6466

16.0463

16.5697

15.5248

14.677

14.6733

15.5294

15.1898

/

/

/

13.4514

11.8095

3、建立模型

1）步骤

①输入数据

；

②作出散点图；

③进行曲线拟合；

④作出拟合曲线图，输出拟合函数；

⑤进行误差估计。

2）具体实现（以多项式拟合为例）

由散点图可知，可采用多项式拟合曲线（如图）：

以8次多项式拟合（其中拟合优度可达96%），可得：

3）误差估计

误差估算时，由于已假定水泵的灌水速度为一常数，同时知道在水泵抽水时，水箱中水的体积的平均变化速度

应近似等于水泵的灌水速度P减去此段时间从水箱中流出的平均速度。

即

此处f（t）在Δt区间的两端点间进行积分。

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

下面通过水泵开始和停止工作的两段区间，即t∈[8.968，10.926]及t∈[20.839，22.958]来进行检验。

第一段：

对应于t始=8.968（h），t末=10.926（h），

水量分别为v始=514.8×103（G），v末=677.6×103（G）

故

（G）

（h）

（G/h）

第二段：

对应于t始=20.839（h），t末=22.958（h），

水量分别为v始=514.8×103（G），v末=677.6×103（G）

故

（G）

（h）

（G/h）

（G/h）

（G/h）

模型求解

将

h和

h代入到水的流速拟合函数

我们得到这两时刻的流速分别近似为13688.8G/h和13335.6G/h,相差仅2.58%,从而可以认为

能近似表达一天的用水流量.于是,一天里的用水总量近似地等于函数

在24小时周期内的积分.有

（G）

（G）

（G）

可见误差均控制在能接受的范围之内，故该社区一天的总用水量约为：

（G）

若按常规每1000人的用水量为105000G/d,因此估计出这个地区大约有3200人.同时可得流出水箱的平均流速为：

（G/h）

水泵灌水的平均速度为：

模型评价

从建模的设想及实施过程我们可以看出有如下优点及不足之处，优点：

①这模型很灵活,能被有一个正水箱的任何小镇使用,输入的数据可为任何规则逼近时间区间的水位,即时间分布可以是随机的.

②模型用到的数学知识是简单易懂的，其计算过程完全可由拟合工具箱来完成，是很容易实现的。

③只要有一台计算机甚至计算器，输入数据，这模型就容易完成。

④模型不仅提供了水流量及一天用水量的较为准确的估计。

还可以估计任何时刻的水流量，包括水泵工作时的水流量。

⑤该模型可推广到用电分布的情况。

缺点：

1模型最大的缺点是无法准确估计结果的误差。

对不同的输入数据，误差则不同。

2用光滑曲线拟合的方法无法模拟真实水流量曲线的微小变化，除非时间区间取得充够小。

③数据太少，精确度稍有欠缺。

其中具体的计算过程如下：

曲线拟合过程（还可用样条插值求解splinetool）

>>t=[0.4605561.382082.396393.410564.424725.43903…6.453197.46758.4481911.493512.493613.415114.4290…15.443116.365017.378918.484619.498520.399324.4335…25.4476];

>>f=[13.471011.595310.34989.734719.487358.69649…9.489748.9008610.103618.583319.676618.646616.0463…16.569715.524814.67714.673315.529415.189813.451411.8095]*10^3;

>>scatter（t,f）

>>plot（t,f）

>>f=0.00024547-0.0248441.0109-21.114240.11-1468.84690.8-7839.916281

>>formatlongg

>>a=polyfit（t,f,8）;

>>cftool

LinearmodelPoly8:

fittedmodel1（x）=p1*x^8+p2*x^7+p3*x^6+p4*x^5+

p5*x^4+p6*x^3+p7*x^2+p8*x+p9

Coefficients（with95%confidencebounds）:

p1=0.0002455（0.0001133,0.0003777）

p2=-0.02484（-0.03856,-0.01113）

p3=1.011（0.4274,1.594）

p4=-21.11（-34.24,-7.988）

p5=240.1（72.52,407.7）

p6=-1469（-2678,-259.2）

p7=4691（72.96,9309）

p8=-7840（-1.58e+004,117.9）

p9=1.628e+004（1.203e+004,2.054e+004）

定积分计算过程：

>>q=int（'0.00024547*x^8-0.024844*x^7+1.0109*x^6-21.114*x^5+240.11*x^4-1468.8*x^3+4690.8*x^2-7839.9*x+16281,8.968,10.926）；

>>q=q/1.958；

>>p=int（'0.00024547*x^8-0.024844*x^7+1.0109*x^6-21.114*x^5+240.11*x^4-1468.8*x^3+4690.8*x^2-7839.9*x+16281',20.839,22.958）；

>>p=p/2.119；

>>q=83150+q；

>>p=76830+p；

>>r=（p-q）/p

r=0.0268

模型求解过程

上述所有过程可直接在cftool工具箱中一次性完成。

水箱水流量问题练习题

试将水箱水流量问题的建模方法推广到闭路电视的普及预测模型，下表列出了美国自1952年至1978年闭路电视的统计数据。

年

家庭拥有电视数

家庭拥有闭路电视数

闭路电视百分比

电视台数目

闭路电视系统数目

1952

15300

14

108

70

1954

26000

65

402

300

1956

34900

300

0.9

496

450

1958

41424

450

1.1

556

525

1960

45750

850

1.4

579

640

1962

48850

1085

1.7

571

800

1964

51600

1085

2.1

582

1200

1966

53850

1575

2.9

613

1570

1968

56670

2800

4.4

642

2000

1970

59550

4500

7.7

686

2490

1972

—

6000

9.7

690

2841

1974

—

8700

12.7

694

3158

1976

71460

10800

14.8

701

3681

1978

74700

13000

17.1

708

3997

根据上表中的数据可以绘制美国家庭采用闭路电视的增长曲线。

请利用已有历史资料来预测未来闭路电视在家庭采用的百分比。