计算方法引论实验报告模板 2.docx
- 文档编号:30449041
- 上传时间:2023-08-15
- 格式:DOCX
- 页数:12
- 大小:92KB
计算方法引论实验报告模板 2.docx
《计算方法引论实验报告模板 2.docx》由会员分享,可在线阅读,更多相关《计算方法引论实验报告模板 2.docx(12页珍藏版)》请在冰豆网上搜索。
计算方法引论实验报告模板2
《水塔的流量估算》
实验报告
系别资源环境与旅游系
课程计算方法引论
学期2013-2014学期
班级11测绘2班
学号110714238
姓名张树义
池州学院资源环境与旅游系
【实验目的】
1.了解有关数据处理的基本概念和原理。
2.初步了解处理数据插值与拟合的基本方法,如样条插值、分段插值等。
3.学习掌握用MATLAB命令处理数据插值与拟合问题。
【实验内容】
某居民区有一供居民用水的圆形水塔,一般可以通过测量其水位来估计水的流量。
但面临的困难是,当水塔水位下降到设定的最低水位时,水泵自动启动向水塔供水,到设定的最高水位时停止供水,这段时间是无法测量水塔的水位和水泵的供水量。
通常水泵每天供水一两次,每次约两小时。
水塔是一个高12.2米、直径17.4米的正圆柱。
按照设计,水塔水位降到约8.2米时,水泵自动启动,水位升到约10.8米时水泵停止工作。
某一天的水位测量记录如表1所示,试估计任何时刻(包括水泵正供水时)从水塔流出的水流量,及一天的总用水量。
时刻(h)
水位(cm)
0
968
0.92
948
1.84
931
2.95
913
3.87
898
4.98
881
5.90
869
7.01
852
7.93
839
8.97
822
时刻(h)
水位(cm)
9.98
//
10.92
//
10.95
1082
12.03
1050
12.95
1021
13.88
994
14.98
965
15.90
941
16.83
918
17.93
892
时刻(h)
水位(cm)
19.04
866
19.96
843
20.84
822
22.01
//
22.96
//
23.88
1059
24.99
1035
25.91
1018
曲线拟合的基本方法
曲线拟合问题是指:
已知平面上
个点(
,
),
=0,1,…,
,
互不相同,寻求函数
=
,使
在某种准则下与所有数据点最为接近,即曲线拟合得最好。
线性最小二乘法是解决曲线拟合最常用的方法,其基本思路是,令
=
+
+…+
(5)
其中
是事先选定的一组函数,系数
(
=0,1,…,
<
)待定。
寻求
,使得残差平方和
=
(6)
达到最小。
这里的建模原理实质上与实验七中的回归分析是一致的。
【实验方法与步骤】
1.引例问题的分析
流量是单位时间流出的水的体积,由于水塔是圆柱形,横截面积是常数,在水泵不工作时段,流量很容易从水位对时间的变化率算出,问题是如何估计水泵供水时段的流量。
水泵供水时段的流量只能靠供水时段前后的流量拟合得到,作为拟合的原始数据,我们希望水泵不工作时段的流量越准确越好。
我们可以考虑先用表中数据拟合水位~时间函数然后对之求导即可得到各时段的流量。
流量是时间的连续函数,只取决于水位差,与水位本身无关,与水泵是否工作无关,计算中将流量定义为单位时间流出的水的高度,即水位对时间变化率的绝对值(水位是下降的)。
水泵第1次供水时段为
=9.0到
=11.0(小时),第2次供水时段为
=20.8到
=23.0(小时)。
这是根据最高和最低水位分别为10.8米和8.2米,及表1的水位测量记录作出的假设,其中前3个时刻直接取自实测数据(精确到0.1小时),最后1个时刻来自每次供水约两小时的已知条件(从记录看,第2次供水时段应在记录的22.96小时之后不久结束)。
水泵工作时单位时间的供水量大致为常数,这个常数应该大于单位时间的平均流量。
首先考虑拟合水位~时间函数,从表1测量记录看,一天有两个供水时段(以下称第1供水时段和第2供水时段),和三个水泵不工作时段(简称第1时段
=0到
=8.97,第2时段
=10.95到
=20.84,第3时段
=23以后)。
对第1、2时段的测量数据可直接分别作多项式拟合,得到水位函数。
为使拟合曲线比较光滑,多项式次数不要太高,一般在3~6次。
由于第3时段只有3个测量记录,无法对这一时段的水位作出较好的拟合。
接着确定流量~时间函数,对于第1、2时段只需将水位函数求导数即可,对于两个供水时段的流量,则用供水时段前后(水泵不工作时段)的流量拟合得到,并将拟合得到的第2供水时段流量外推,将第3时段流量包含在第2供水时段内。
最后一天总用水量等于两个水泵不工作时段和两个供水时段(将第3时段包含在第2供水时段内)用水量之和,它们都可以由流量对时间的积分再乘以水塔截面积得到。
1.拟合水位--时间函数
1.第一时段(0~8.97)
x=[00.921.842.953.874.985.907.0417.938.97];
y=[968948931913898881869852839822];
a=polyfit(x,y,3)
fval=polyval(a,x);plot(x,fval),holdon
plot(x,y),holdoff
2第二时段(10.95~20.84)
x=[10.9512.0312.9513.8814.9815.9016.8317.9319.0419.9620.84];
y=[108210501021994965941918892866843822];
a=polyfit(x,y,3)
fval=polyval(a,x);
plot(x,fval),holdon
3第三时段(23.88~25.91)
x=[23.8824.9925.91];
y=[105910351018];
a=polyfit(x,y,3)
fval=polyval(a,x);
plot(x,fval),holdon
2.拟合流量--时间函数
对于1,2时段只需将1,2时段的水位--时间函数求导即可,对于供水时段的流量,则用供水时段前后(水泵不工作时段)的流量拟合的到,并且将拟合得到的第二供水时段流量外推,将第3时段流量包含在低2供水时段内。
1.第一时段(0~8.97)
对第一时段求导得得到第一时段流量--时间函数y=0.2451*x.^2-2.7984*x+22.2227
x=linspace(0,8.98);
y=0.2451*x.^2-2.7984*x+22.2227;
plot(x,y)
2第二时段(10.95~20.84)
对第二时段求导得得到第一时段流量--时间函数y=0.137*x.^2-5.2517*x+74.1016
x=linspace(10.95,20.84);
y=0.137*x.^2-5.2517*x+74.1016;
plot(x,y)
2第三时段(23.88~25.91)
对第三时段求导得得到第一时段流量--时间函数y=-0.4848*x.^2+21.0788*x-203.8489
x=linspace(23.88,25.91);
y=-0.4848*x.^2+21.0788*x-203.8489;
plot(x,y)
4将三时段函数在一个图中做出
x=linspace(0,8.98);
y1=0.2451*x.^2-2.7984*x+22.2227;
plot(x,y1);
holdon;
x=linspace(10.95,20.84);
y2=0.137*x.^2-5.2517*x+74.1016;
plot(x,y2);
holdon;
x=linspace(23.88,25.91);
y3=-0.4848*x.^2+21.0788*x-203.8489;
plot(x,y3);
【实验心得】
本次试验,我学会了matlab的一些基本的操作,了解了一些简单的命令和用法。
并且用这些简单算法去处理老师布置的作业。
实验的结果还是不太理想,中间有些处理不好。
对于matlab的数据处理方法运用上还是一知半解,大多都是依靠XX去查询遇到的问题然后根据自己的理解进行实验的处理其中有一些处理很不理想。
最后将三段函数进行拟合连接时更是一头雾水,中间那段空白不知该如何连接起来尝试了多次无果。
同时反映出对于matlab的认知还差很多。
对其中的数据进行检核是发现有些数据误差还是很大的。
经过这次实验还是有许多的收获,首先对于matlab的认知和简单程序的操作有所了解,对将来的数据处理带来了方便。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算方法引论实验报告模板 计算方法 引论 实验 报告 模板
![提示](https://static.bdocx.com/images/bang_tan.gif)