数学建模——水塔流量问题.docVIP

实验十四水塔流量问题

【实验目得】

１．了解有关数据处理得基本概念与原理。

2。初步了解处理数据插值与拟合得基本方法，如样条插值、分段插值等。

３。学习掌握用MATLAB命令处理数据插值与拟合问题。

【实验内容】

某居民区有一供居民用水得圆形水塔，一般可以通过测量其水位来估计水得流量.但面临得困难就是,当水塔水位下降到设定得最低水位时,水泵自动启动向水塔供水,到设定得最高水位时停止供水,这段时间就是无法测量水塔得水位与水泵得供水量。通常水泵每天供水一两次,每次约两小时。水塔就是一个高1２、2米、直径17、4米得正圆柱。按照设计,水塔水位降到约8、2米时,水泵自动启动，水位升到约1０、8米时水泵停止工作.

某一天得水位测量记录如表１所示,试估计任何时刻（包括水泵正供水时）从水塔流出得水流量，及一天得总用水量。

表１水位测量启示录(//表示水泵启动）

时刻（h）

水位（ｃm）

0

968

０、92

９48

1、８4

931

2、95

9１3

3、８7

89８

4、98

８８1

5、90

869

7、01

852

7、９3

8３9

8、97

8２2

时刻（h)

水位（ｃｍ)

9、９8

//

１０、9２

//

10、９5

1082

12、03

10５0

12、95

１０21

13、88

994

14、９8

965

1５、90

941

1６、８３

9１8

17、9３

8９2

时刻(h）

水位(cｍ)

１9、０4

866

１9、9６

84３

20、８４

8２2

22、01

//

２2、９6

/／

23、8８

10５9

２４、99

103５

2５、91

1018

【实验准备】

在生产实践与科学研究中，常常遇到这样得问题:由实验或测量得到得一批离散样点,需要确定满足特定要求得曲线或曲面（即变量之间得函数关系或预测样点之外得数据)。如果要求曲线(面）通过所给得所有数据点(即确定一个初等函数通过已知各数据，一般用多项式或分段多项式），这就就是数据插值。在数据较少得情况下，这样做能够取得好得效果.但就是,如果数据较多，那么插值函数就是一个次数很高得函数，比较复杂。如果不要求曲线（面)通过所有得数据点，而就是要求它反映对象整体得变化趋势,可得到更简单实用得近似函数,这就就是数据拟合.函数插值与曲线拟合都就是要根据一组数据构造一个函数作为近似，由于近似得要求不同，二者在数学方法上就是完全不同得.

1。数据插值得基本方法

拉格朗日插值

若知道函数＝在互异得两个点与处得函数值与,而想估计该函数在另一点处得函数值，最自然得想法就是作过点（，）与点（,)得直线=，用作为准确值得近似值,如果得到得结果误差太大,还可增加一点得函数值,即已知＝在互异得三个点，与处得函数值，与，可以构造过这三点得二次曲线=，用作为准确值得近似值。

一般得，若已知=在互异得＋１个点,,…，处得函数值，,…，，则可以考虑构造一个过这+１个点得次数不超过得多项式

＝＋+…+＋(1）

通过所有+１个点,即满足

=,=0,1，…，（2）

然后用作为准确值得近似值。这样构造出来得多项式称为得次拉格朗日插值多项式或插值函数。

分段插值

多项式历来都被认为就是最好得逼近工具之一，它插值光滑，但不具有收敛性,会随着节点数目增多而次数升高，一般不宜采用高次多项式（如＞７)插值，否则逼近得效果往往就是不理想得,甚至发生龙格振荡（当节点数目不断增大时,在区间中部趋于,但对于区间两端得,并不趋于,也称龙格现象）。

在插值范围较小,用低次插值往往就能奏效。最直观得办法就就是将各数据点用折线连接起来，这种增加节点，用分段低次多项式插值得化整为零得处理方法称作分段插值法,即不去寻求整个插值区间上得一个高次多项式，而就是把区间划分为若干个小区间。如果

＝＜…=(3）

那么分段线性插值公式为

＝＋，〈≤,＝0，１,…，（4）

分段线性插值通常有较好得收敛性与稳定性,算法简单，克服了龙格现象,其缺点就是不如拉格朗日插值多项式光滑。

样条插值

分段线性插值函数在节点得一阶导数一般不存在，且不光滑,这就导致了样条插值函数得提出。在机械制造、航海、航空工业中，经常需要解决下列问题:已知一些数据点（,），（,),（,）,如何全部通过这些数据点作一条比较光滑得曲线呢？

绘图员解决了这一问题，首先把数据描绘在平面上，再把一根富有弹性得细直条（称为样条)弯曲，使其一边通过这些数据点,用压铁固定其形状，沿样条边绘出一条光滑得曲线，往往要用几根样条，分段完成上述工作,同时也应让连接点处保持光滑.对绘图员用样条画出得曲线,进行数学模拟，就导出了样条函数得概念。如今已经