重庆大学-数学实验-6插值.docVIP

开课学院、实验室：数学与统计DS1407 实验时间 ： 课程 名称 数学实验 实验项目 名 称 水塔用水量的估计 ——插值 实验项目类型 验证 演示 综合 设计 其他 指导 教师 肖剑 成 绩 实验目的 [1] 了解插值的基本原理 [2] 了解拉格朗日插值、线性插值、样条插值的基本思想； [3] [4] 掌握用MATLAB计算三种一维插值和两种 [5] 通过范例展现求解实际问题的初步建模过程； 1．编写M文件； 2． 通过数值和图形输出，比较它们的效果； 针对实际问题，试建立数学模型，并求解。 1. 一维插值 利用以下一些具体函数，考察分段线性插值、三次样条插值和拉格朗日多项式插值等三种插值方法的差异。 1）x([-5,5]； x=-5:10/(m-1):5; y=1./(1+x.^2);z=0*x; plot(x,z,r,x,y,LineWidth,1.5), gtext(y=1/(1+x^2)),pause n=7; x0=-5:10/(n-1):5; y0=1./(1+x0.^2); y1=lagr1(x0,y0,x);hold on, plot(x,y1,r),gtext(lagr),pause, y2=interp1(x0,y0,x,spline); plot(x,y2,k);gtext(scyt),pause, y3=interp1(x0,y0,x); plot(x,y3,g);gtext(fenduan),pause, hold off 运行结果： 图1：插值比较图（红色：拉格朗日插值 蓝色：函数值 黑色：三次样条插值 绿色：分断插值） 由图直观可以看出：分断插值和三次样条与函数值较接近，拉格朗日插值与函数值相差较大 分析：在三种插值对比运用的过程中，进一步强化了hold on，hold off 的运用。 （二）应用实验 2．火车行驶的路程、速度数据如表7.，计算从静止开始20 表7.1 t(分) 2 4 6 8 10 12 14 16 18 20 v(km/h) 10 18 25 29 32 20 11 5 2 0 v0=[0 10 18 25 29 32 20 11 5 2 0]; t=0:.1:20; v=interp1(t0,v0,t,spline); plot(t0,v0,r,t,v) xlabel(t(min)),ylabel(v(km/h)) grid trapz(t,v) 运行结果： ans = 306.3918 图形： 图2：火车速度与时间关系的三次样条图 分析：三次样条插值和分断插值的运用相对普遍，求解积分的trapz函数还不够熟悉。 确定地球与金星之间的距离 1914年8月份的7次观测中，测得地球与金星之间距离（单位：米），并取其常用对数值，与日期的一组历史数据如表7.。 7.2 日期（号） 18 20 22 24 26 28 30 距离对数 9.9617724 9.9543645 9.9468069 9.9390950 9.9312245 9.9231915 9.9149925 由此推断何时金星与地球的距离（米）的对数值为9.9351799？ s0=[9.9617724 9.9543645 9.9468069 9.9390950 9.9312245 9.9231915 9.9149925]; day=18:0.001:30; sequ=9.9351799; s=interp1(day0,s0,day); plot(day0,s0,r,day,s), hold on plot(day,sequ,g), title(地球与金星距离随时间变化图); xlabel(时间day),ylabel(距离s); grid 运行结果： 图3：地球与金星距离随时间变化图 分析：进一步巩固了插值的运用，回顾了利用放大法求解数值的过程。 4．山区地貌图 0，2800）(0，2400）内，单位：米）测得一些地点的高程（单位：米）如表7.，试作出该山区的地貌图和等高线图。 7.3 2400 2000 1600 1200 800 400 0 1430 1450 1470 1320 1280 1200 1080 940 1450 1480