曲线曲面的插值与拟合方法次课 PPT.pptxVIP

曲线曲面的插值与拟合方法次课函数查表与地图边界线绘制?(2.3457)=？引例1 函数查表问题:已知标准正态分布函数表,求表中没有的值 ?(2、34)=0、99036 ?(2、35)=0、99061求~?(2、3457) ≈(2、35-2、3457)/(2、35-2、34)* ?(2、34)+ (2、3457-2、34)/(2、35-2、34)* ?(2、35)引例2 地图绘制问题:假如我们在地图边界获取了一些边界点的坐标,连接这些边界点形成闭合曲线,能够用来近似表示真实边界线,如何更准确地逼近真实边界线？如何更准确地逼近真实边界线？插值在数码图像放大中的应用引例3 图像插值放大:数码相机运用插值的方法能够创造出比传感器实际像素更多的图像,这种处理称为“数码变焦”。106*40原始图像:左边:最近邻插值放大450%右边:双三次插值放大450%插值在图像三维重建中的应用Surface recostruction from scattered points cloud 分段线性插值与拉格朗日插值分段线性插值:用直线(线性)连接数据点列上相邻的两点。比如~在两点[xi-1,xi]上线性插值函数为~拉格朗日插值:用n次拉格朗日插值多项式连接数据点列上相邻的n+1个点。Pszjs71拉格朗日插值基函数的构造比如 在三个点[x0,x1,x2]上lagrange插值函数为(线性插值是拉格朗日插值最简单的情形) 分段三次埃尔米特插值条件数分段三次埃尔米特插值:线性插值在每一小段上(两点之间),用到2个条件q(xi)=yi,因此确定了一个线性插值函数;三次埃尔米特插值在每一小段上,用到4个条件q(xi)=yi, q(xi)=yi,因此确定一个3次多项式插值函数。 分段插值主要是为了幸免高次插值估计出现的大幅度振荡现象,在实际应用中通常采纳分段低次插值来提高近似程度,比如可用分段线性插值或分段三次埃尔米特插值来逼近已知函数,但它们的总体光滑性较差,为了克服这一缺点,三次样条插值成为比较理想的工具。三次样条(spline)插值的概念 样条的概念出自工程设计与机械加工(飞机、船舶外形曲线设计)中的绘图工具(曲线尺),简单说就是具有连续二阶导数的三次插值多项式函数。 三次样条(spline)插值的条件数x2x1q2q1x0 首先从段数n=2分析:我们明白在每一小段的三次多项式有4个系数,因此如下图,总共需要有4*2=8个方程来确定; 由q(xi)=yi能够确定2*2=4个方程,又由内部节点q1(xi)= q2(xi)与q1(xi)= q2(xi)能够确定2*(2-1)=2个方程,看来剩下的8-(4+2)=2个方程只有靠外部给定(边界条件)了一维曲线等距插值函数interpinterps syntaxOne-dimensional r times longer data interpolation y = interp(y,r)题例 在原始数据点中增倍插值x=0:0、001:1; y=sin(2*pi*30*x)+sin(2*pi*60*x); yi=interp(y,4);subplot(1,2,1); stem(y(1:30)); title(Original Points);subplot(1,2,2); stem(yi(1:120)); title(Interpolated Points);一维曲线等距插值函数interp1interp1s syntaxOne-dimensional data interpolation yi = interp1(x,y,xi,method)nearest Nearest neighbor interpolationlinear Linear interpolation (default)spline Cubic spline interpolationcubic Piecewise cubic Hermite interpolation题例 在一天24小时内,从零点开始每间隔2小时测得的环境温度,推测在15点6分的的温度x=0:2:24; y=[12,9,9,10,18,24,28,27,25,20,18,15,13];plot(x,y,-ro); hold on; xi=15、1; yi=interp1(x,y,xi,spline),xi=0:1/3600:24; yi=interp1(x,y,xi,spline); plot(xi,yi,b-);二维曲面等距插值函数interp2interp2s syntaxTwo-dimensional data interpolation ZI = interp2(X,Y,Z,XI,Y