插值_数据的统计描述.pptVIP

插值 一、机床加工问题 待加工零件的外形根据工艺要求由一组数据(x,y) 给出(平面情况下)，用程控机床加工时每一刀只能沿x方向和y方向移动非常小的一步。试根据下表数据给出横坐标改变0.1时的纵坐标(x介于0和15之间）。 0 3 5 7 9 11 12 13 14 15 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 二、插值 1.插值问题： 　已知n+1个节点(xi,yi), j=0,1,2,…,n,其中xj互异。不妨设a=x0x1…xn=b.求任一插值点x*处的插值y*. 构造相对简单的便于计算函数y=f(x)，使f(x)满足f(xi)=yi,然后利用求得的函数y=f(x)计算y*=f(x*). 2.插值方法 拉格朗日插值、分段线性插值、三次样条插值 三、Matlab命令及实现 yi = interp1(x,Y,xi,method) method: ‘linear’—分段线性插值,只能内插 　　　　‘spline’—三次样条插值 　　　　‘pchip’—分段三次Hermite插值 　　　　‘cubic’ 思考：1.对于n+1个节点，为什么用n次 多项式？ 2.若ff(x)=1,则sigma(0,n)li(x)=1 3.若f(x)为不超过n次的多项式，则ff(x)=f(x) 拉格朗日插值 function y=lagrint(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end 机床加工： x0=[0 3 5 7 9 11 12 13 14 15]; y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6]; x=0:0.1:15; y_lag=lagrint(x0,y0,x); y_plin=interp1(x0,y0,x); y_spline=interp1(x0,y0,x, spline); subplot(3,1,1) plot(x0,y0, k+,x,y_lag, r),grid Title(lagrange),subplot(3,1,2) plot(x0,y0, k+,x,y_plin, r),grid title(piecewise linear),subplot(3,1,3) plot(x0,y0, k+ ,x,y_spline, r ),grid title( spline ) 四、插值与拟合异同 　相同：根据一组数据构造一个函数 　　　　作为近似。 　不同： 　　1) 数学方法不同； 　　2) 插值构造的函数不通过数据点,而拟 合构造的函数不一定通过数据点，使 数据点与曲线总体来说尽量接近。 问题：针对实际问题，如何确定用插 　　　值还是拟合？ 五、练习： 1.内燃机转角与升程的关系 　内燃机的气门开启状态是由与发电机转轴连接 的凸轮来控制的，气门开启的大小决定了发动机 的性能，因此驱动气门的凸轮加工曲线是非常重 要的，下表给出了某凸轮转角与升程的几组实测 数据，为了机床加工的需要，试确定转角介于91 ~180度每隔1度的升程。 转角θ 升程x 91 105 110 115 120 124 128 0 1.0869 1.971 2.7555 3.3986 4.9073 8.3409 2.随高度变化的大气压 　大气压强p随高度x变化的理论公式为p=1.0332 exp{-(x+500)/7756}.为验证这一公式，测得某地大气 压强随高度变化的一组数据如下表，试用插值法和 拟合法进行计算并绘图，说明那种方法较合理。 转角θ 升程x 0 300 600 1000 1500 2000 0.9689 0.9322 0.8969 0.8519 0.7989 0.7491 六、多元函数的插值 　矩形温箱的温度 　　测量一个矩形温箱，已知其3×5个测试点的温度 如下表所示。求沿宽度x和深度y细分为x=1:0.2:3,y=1:0.2:5时交点上各点的温度。 x y 1 2 3 4 5 1 2 3 82 81 80 82 84 81 65 61 63 79 84 84 82