水塔用水量的估计-插值

1、重庆大学学生实验报告实验课程名称数学实验开课实验室DS1421学院年级专业班学生姓名学 号开课时 间 2013 至 2014 学年第 2 学期总成绩教师签名数学与统计学院制开课学院、实验室：数统学院 DS1421实验时间：2014年5月28日课程名称数学实验实验项目名称水塔用水量的估计插值实验项目类型验证演示综合设计其他指导教师肖剑成 绩实验目的1234了解插值的基本原理了解拉格朗日插值、线性插值、样条插值的基本思想； 了解三种网格节点数据的插值方法的基本思想； 掌握用MATLAB十算三种一维插值和两种二维插值的方法;通过范例展现求解实际问题的初步建模过程；基础实验一、实验内容1 编写拉格朗日

2、插值方法的函数M文件2用三种插值方法对已知函数进行插值计算，通过数值和图形输出，比较它们的效果;3.针对实际问题，试建立数学模型，并求解。二、实验过程(一般应包括实验原理或问题分析，算法设计、程序、计算、图表等,验结果及分析)1. 一维插值，利用以下一些具体函数，考察分段线性插值、三次样条插值和拉格朗日多项式插值等三 种插值方法的差异。1 101) 2 , xW-5,5 ;2 ) sin x, x“0,2 二;3 ) cos x, x“0,2 二.1 x2注意：适当选取节点及插值点的个数；比较时可以采用插值点的函数值与真实函数值的差异，或采用两个函数之间的某种距离。(1) 程序:21 ) 1/

3、(1+x ) , x 三-5,5x=-5:0.1:5;y=1./(1+(x.A2);x0=-5:1:5;y0=1./(1+(x0.A2);y1=i nterp1(x0,y0,x,'li near'y2=i nterp1(x0,y0,x,'spli ne'y3=lagr1(x0,y0,x);plot(x0,y0,'k*',x,y1,'b'gtext( 'li near');gtext( 'spli ne');gtext( 'Lagr');,x,y2,'g',x,y3,&

4、#39;r'););2 ) sin x, x 三0,2 二x=0:0.1:pi;y=si n( x); x0=0:0.5:pi; y0=s in( x0);y1=i nterp1(x0,y0,x,'li near');y2=i nterp1(x0,y0,x,'spli ne');y3=lagr1(x0,y0,x); plot(x0,y0,'k*',x,y1,'b',x,y2,'g' ,x,y3,gtext( 'li near');gtext( 'spli ne');gtext

5、( 'Lagr');'r')3) cos10x, x 0,2 二'lin ear');x=0:0.1:2*pi; y=(cos(x)"0; x0=0:0.5:2*pi; y0=(cos(x0).A10; y1=i nterp1(x0,y0,x,y2=i nterp1(x0,y0,x,'spli ne' );y3=lagr1(x0,y0,x);plot(x0,y0,'k*',x,y1,'b',x,y2,'g',x,y3,'r')gtext( 'li n

6、ear');gtext( 'spli ne');gtext( 'Lagr');(2) 结果:File0splineQlinear-2口0linear0.423.60 8y=sin xFigure 1Edit Qew Insert Tools desktop window HelpQS J2y=1/(1+x )Edit Viei<v Insert Tools desktop window HelpFigure 110y=cos x(3) 分析：由图可以看出，函数y=1心+x 2)使用三次样条插值效果最好，函数y=sinx使用拉格朗日插值效果最好，y=

7、cos 10x使用分段线性插值效果最好，可见，三种插值方法各有各自最适用的函数。2. 轮船的甲板成近似半椭圆面形，为了得到甲板的面积。首先测量得到横 向最大相间8.534米；然后等间距地测得纵向高度，自左向右分别为：0.914, 5.060, 7.772, 8.717, 9.083, 9.144, 9.083, 8.992, 8.687, 7.376, 2.073，计算甲板的面积。(1) 程序：x=li nspace(0,8.534,13);y=0 0.914 5.060 7.772 8.717 9.083 9.144 9.083 8.992 8.687 7.376 2.073 0;x0=0:

8、0.001:8.534;y1=i nterp1(x,y,x0);figure,plot(x,y,'k*',x0,y1,'-r')S=trapz(y1)*0.001(2) 结果：S = 54.6894Figure 1File Edit View Insert Tools Desktop Window Help理tod A也fe -v C ®* 巳百 0X(3) 分析：甲板横向最大相间为 8.534米，然后等间距地测得纵向高度 ，共有11个值，所以应该是吧 8.534米分 成12分，对应的值为纵向高度；以左边零点位坐标原点，建立坐标系。线性插值得到图形，再

9、用数值积 分可求面积。3.火车行驶的路程、速度数据如表 7.2，计算从静止开始 20分钟内走过的路程。表7.2t(分)2468101214161820v(km/h)10182529322011520(1) 程序：x=0:2:20;y=0 10 18 25 29 32 20 11 5 2 0; x0=0:0.001:20;y1= in terp1(x,y,xO,'spli ne' );plot(x,y, 'k*',xO,y1,'r')S=trapz(y1)*0.001(2) 结果：S = 304（3）分析：用线性插值的方法作出火车行驶的v-t关系图

10、，则火车行驶的路程为图形的面积，用数值积分的方法可以求出。（单4 .确定地球与金星之间的距离，天文学家在1914年8月份的7次观测中，测得地球与金星之间距离位：米），并取其常用对数值，与日期的一组历史数据如表7.3。表7.3日期（号）18202224262830距离对数9.96177249.95436459.94680699.93909509.93122459.92319159.9149925由此推断何时金星与地球的距离(米)的对数值为9.9351799 ?(1)程序：x=18:2:30;y=9.9617724 9.9543645 9.9468069 9.9390950 9.9312245 9

11、.9231915 9.9149925; x0=18:0.1:30;y1= in terp1(x,y,xO,'spli ne' );plot(x,y, 'k*',xO,y1,'r')gridgtext( 'y=9.9351799');(2)结果:从图像中可以看出在25号时金星与地球的距离的对数值为9.9351799 。(单5.日照时间分布表7.4的气象资料是某一地区1985-1998年间不同月份的平均日照时间的观测数据位：小时/月)，试分析日照时间的变化规律。表7.4月份123456789101112日照80.967.267.150

12、.532.033.636.646.852.362.064.171.2(1)程序：mo nth=1:12;temps=80.9 67.2 67.1 50.5 32.0 33.6 36.6 46.8 52.3 62.0 64.1 71.2; m=1:0.1:12;t=in terp1(m on th,temps,m,'spli ne');plot(m on th,temps,'k*',m,t,'r')xlabel('mon th'),ylabel( 'degrees clsius')(2)结果:从图像中可以看出，日照时

13、间从一月到五月逐渐递减，从六月到十二月逐渐增加，近似于二次抛物线。6.山区地貌图 位：米）如表 7.5，在某山区（平面区域 （0, 2800）（ 0, 2400）内，单位：米）测得一些地点的高程（单试作出该山区的地貌图和等高线图。表7.52400200016001200800400143014501470132012801200108094014501480150015501510143013001200146015001550160015501600160016001370150012001100155016001550138012701500120011001350145012001150

14、1230139015001500140090011001060118013201450142014001300700900Y/X0 400 800 1200 1600 2000 2400 2800（1）程序：x=0:400:2800;y=0:400:2400;1320 1280 1200 1080 940z=1430 1450 1470 1450 1480 1500 1550 1510 1430 1300 1200 1460 1500 1550 1600 1550 1600 1600 16001370 1500 1200110015501600 155013801270 1500 120011

15、0013501450 120011501230 1390 150015001400900 110010601180 1320 1450142014001300 700900;figure(1);xi,yi=meshgrid(0:50:2800,0:50:2400);zi=i nterp2(x,y,z,xi,yi, meshz(xi,yi,zi) xlabel( X ),ylabel( figure(2);C,H=c on tour(xi,yi,zi); clabel(C,H)xlabel( 'x' ),ylabel('spli ne');'Y' ),zlabel('Z&