数学实验_第四章.docVIP

云南大学数学与统计学实验教学中心 实 验 报 告 课程名称：数学实验 学期： 成绩： 指导教师： 学生姓名： 学生学号： 实验名称： 插值问题 实验要求： 实验学时： 实验编号：04 实验日期： 完成日期： 学院： 数学与统计学院 专业 ：数学与应用数学 年级： 一、实验目的 1、学会用Matlab求简单微分方程的解析解. 2、学会用Matlab的dslove 语句。 3、学会应用插值思想以及插值方法解决实际问题，进行建模、求解与分析。 二、实验环境 Matlab7.0版本 三、实验内容 P88 4.7.1 1.火车行驶的路程，速度数据如下表，计算从静止开始20分钟内走过的路程。 t(分) 2 4 6 8 10 12 14 16 18 20 v（公里/小时） 10 18 25 29 32 20 11 5 2 0 P89 4.7.1 4.确定地球与金星之间的距离 天文学家在1914年8月份的7次观察中，测得地球与金星之间的距离（单位：米），并取其常用对数值与日期的一组历史数据如下表，试推测何时金星与地球的距离（米）的对数为9.9351799 日期 18 20 22 24 距离对数 9.9617724 9.9543645 9.9468069 9.9390950 日期 26 28 30 距离对数 9.9312245 9.9231915 9.9149925 P89 4.7.2 3. 山区地貌图 在某山区（平面区域（0，2800）(（0，2400）内，单位：米）测得一些地点的高程（单位：米）如表，试作出该山区的地貌图和等高线图。 表5.5 2400 2000 1600 1200 800 400 0 1430 1450 1470 1320 1280 1200 1080 940 1450 1480 1500 1550 1510 1430 1300 1200 1460 1500 1550 1600 1550 1600 1600 1600 1370 1500 1200 1100 1550 1600 1550 1380 1270 1500 1200 1100 1350 1450 1200 1150 1230 1390 1500 1500 1400 900 1100 1060 1180 1320 1450 1420 1400 1300 700 900 Y/X 0 400 800 1200 1600 2000 2400 2800 四、实验过程 程序如下： 利用一维三次样条插值方法可以得到时间与速度的曲线. 输入： t=[2/60,4/60,6/60,8/60,10/60,12/60,14/60,16/60,18/60,20/60]; v=[10,18,25,29,32,20,11,5,2,0]; ti=0:0.01:20/60; vi=interp1(t,v,ti,spline); plot(t,v,o,ti,vi);xlabel(t(分));ylabel(v(公里/小时)); title(速度曲线) 图形所围成面积即火车从静止开始20分钟内走过的路程，用梯形公式进行数值积分 输入： sum=0; for i=1:1:9; sum=sum+1/2*(v(i)+v(i+1))*(1/30); end sum 输出： sum = 4.9000 利用一维三次样条插值可以得到距离与时间的曲线图形，程序如下 输入： （1）画散点图 t=[18 20 22 24 26 28 30]; L=[9.9617724 9.9543645 9.9468069 9.9390950 9.9312245 9.9231915 9.9149925]; plot(t,L,*) 图如下 （2）计算 t=18:2:30; L=[9.9617724 9.9543645 9.9468069 9.9390950 9.9312245 9.9231915 9.9149925]; L1=9.9351799; t1=interp1(L,t,L1,spline) 输出 t1 = 25.0000 P4.7.2 山区地貌图问题下面采用双三次插值方法来得到山区地貌图和等高线图，程序如下 输入： x=0:400:2800; y=0:400:2400; [X,Y]=meshgrid(x,y); z=[1180 1320 1450 1420 1400 1300 700 900; 1230 1390 1500 1500 1400 900 1100 1060; 1270