插值与拟合应用举例.docxVIP

插值与拟合 1.插值与拟合的基本概念 插值与插值函数：已知由g（Q （可能未知或非常复朵）产生的一批离散数据 （兀i，X）,j = O,,且 n+1 个互异插值节点a = xQxl--xn_[xn=bf 在插值区 间内寻找一个相对简单的函数/（兀），使其满足下列插值条件： 再利用已求得的/（劝计算任一非插值节点/的近似值y = ，这就是插值。其屮 /（X）称为插值函数，g（x）称为被插函数。 最小二乘拟合：已知一批离散的数据（兀?,必）丿=0丄???,/7,兀互不相同，寻求一个拟合 函数/（X）,使/（不）与X的误差平方和在最小二乘意义下最小。在最小二乘意义下确定 的/（无）称为最小二乘拟合函数。 温度问题 在12小H寸内，每隔1小时测量一次温度。温度依次为：5, 8, 9, 15, 25, 29, 31, 30, 22, 25, 27, 24o （单位:°C） （1） 试估计在3.2h, 6.5h, 7」h, 11.7h的温度值,并画出其图形。 （2） 每隔l/10h估计一次温度值，并画出其图形。 请你找出跟上述12个数据拟合的最好的一条曲线，请分别用分段线性插值、三次样条插值方 法（至少用两条不同的曲线，并比较它们拟合好坏的程度） hours=l: 12; temps二[5,9,15,25,29,31,30,22,25,27,24]; t=interp 1 （hours,temps,[3.2,6.5,7.1,11.7]） %线性插值 T=interp 1 （hours,temps, [3.2,6.5,7.1,11.7],spline） % 三次样条插值 计算结果为 t = 10.2000 30.0000 30.9000 24.9000 T = 9.6734 30.0427 31.1755 25.3820每隔l/10h估计一次温度值并画出其图形: hours=l:12; temps 二[5,8,9,15,25,29,31,30,22,25,27,24]; h 二 1:0.1:12; 匸 interp 1 (hours,temps,h,spline); plot(hours,temps,+,h,t,hours,temps,r:) xlabel(时 I可),ylabel(温度? 三次多项式拟合：hours=l:12; 三次多项式拟合： hours=l:12; temps 二[5,9,15,25,29,31,30,22,25,27,24]; a=polyfit(hours5temps,3) temps 1 =polyval(a,hours); plot(hours,temps/ro\hours4emps 1 /b/) 得到 y = -0.0065%3 一0.3283兀$ +7.128lx-4.4343 ,图形如下: 35 30 25 20 15 10 5 °0 2 48 10 12 °0 2 4 8 10 12 四次多项式拟合：得到 y = 0.0273无4 一0.7158? + 5.力07” 一 12.2251兀 +12.5884 比较拟合的好坏： 设5；?为拟合函数的值，牙为测量值，则残差幺=工（必-刃）2 。 I 汽车刹车距离 美国的某些司机培训课程中有这样的规则：正常驾驶条件下车速没增加10英里/小吋， 后面与前面一辆车的距离应增加一个车身的长度。又云，实现这个规则的一种简便办法是所 谓2秒准则。即后车司机从前车经过某一标志开始默数2秒钟后到达同一标志，而不管车速 如何。 试判籬2秒准则与上述规则是一样的么？这个规则的合理性如何，是否有更好的规则。 车速和刹车距离的一组数据 车速(km/h) 20 40 60 80 100 120 140 刹车距离 6.5 17.8 33.6 57.1 83.4 118.0 153.5 实际距离的拟合多项式，d=kiv2+ k2v其M-file程序如下: v=[20:20:140J/3.6; v2=v.A2; x=[v;v2]*; d=[6.5,17.33?6,57丄83?4,118?0,153?5『； a=x\d dd=x*a ddd=[6.5,17.33?6,57.1,83?4,118.0,153.5]; b=polyfit(v,ddd,2) y=polyval(b,v) plot(v,ddd/ro*,v,dd/b*) 得到 k、= 0.0853 , k2 = 0.6522 图3计算刹车距离曲线 作岀实际距离的拟合2次多项式，d=k)v2+ k2v + k3其M?file程序如下: v=[20:20:1401/3.6; d=[6.5,l 7.8,33.6,57.1,83.4,118.0,153.5]; k二polyfit(v,d,2) dl=polyval(k,v) plot(v,