电子科大MATLAB第10节 分段插值解析.pptVIP

第四章 函数的数值逼近 第二节 分段插值 分段插值 问题： 随着插值结点数增加，插值多项式的次数也相应增加，而对于高次插值容易带来剧烈振荡，带来数值不稳定。为了既要增加插值节点，减小插值区间，以便更好的逼近被插值函数，又要不增加插值多项式的次数以减少误差，可采用分段插值的办法。 线性分段插值 线性分段插值收敛性 线性分段插值收敛性 分段Hermite插值 分段Hermite插值 或者直接分段求解线性方程组计算系数。 样条插值 三次样条函数是由分段三次曲线连接而成，且二阶导数连续的函数。 样条插值 样条插值 样条插值 样条插值 该矩阵严格对角占有，故非奇异，所以样条插值存在且唯一。 样条插值 样条插值的收敛性 定理给出了f (x)与S (x)及其1、2、3阶导数的差， 0、1、2阶误差：步长h越小，误差越小； 3阶导数误差：步长相等误差最小。 样条插值 分段线性插值 样条插值 一维插值：yi = interp1(x, y, xi, ‘method ’) 二维插值：zi= interp2(x, y, z, xi, yi, ‘method’) 三维插值：vi = interp3(x,y,z,v,xi,yi,zi, ‘method’) 插值比较 clear all close all clc % 插值比较 x=-5:1.25:5; y=1./(1+x.*x); xi=-5:0.2:5; Yi(:,1)=1./(1+xi.*xi); Yi(:,2)=interp1(x,y,xi, linear); Yi(:,3)=interp1(x,y,xi, spline); Yi(:,4)=interp1(x,y,xi, cubic); for iii = 1 : 3 figure hold on plot(x,y) plot(xi,Yi(:, 1)) plot(xi,Yi(:, iii + 1)) end xxxx = -5 : 0.01 : 5; yyy = interp1(x,y,xxxx, spline); 应用举例-自由绘制工具 插值节点 约束条件