第四章_多项式与插值--席.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 分段线性插值在MATLAB中的函数调用： 1）横纵坐标都是向量； 2）x是n维向量，y为矩阵时，行数必须是x的维数。 例：x0 =-6:6; y0 =sin(x0); xi = -6:.25:6; %插值点 yi = interp1(x0,y0,xi); x=-6:0.001:6; y=sin(x);plot(x0,y0,o,xi,yi,x,y,:), legend(节点(x_{i},y_{i}),分段线性插值函数,被插值函数y=sinx) title(y=sinx及其分段线性插值函数和节点的图形, Color, w) 分段Hermite插值：即在每个小区间上用分段的 三次Hermite插值。 在matlab里有专门的调用函数： y1=interp1(x,y,x1,’pchip’) 例: x0 =-6:6; y0 =sin(x0); xi = -6:.25:6; yi = interp1(x0,y0,xi,pchip); x=-6:0.001:6; y=sin(x); plot(x0,y0,o,xi,yi,x,y,:), legend(节点(xi,yi),分段埃尔米特插值函 数,被插值函数y=sinx) title(y=sinx及其分段埃尔米特插值函数和 节点的图形, Color, w) x=-5:5; y=1./(1+x.^2); plot(x,y,x,y,o) x=-5:5; y=1./(1+x.^2); xi=-5:.05:5; yi=spline(x,y,xi); plot(xi,yi,b,x,y,ro) 被插值函数: -5≤ x ≤ 5 七、样条插值 分段线性插值只能保证曲线的连续性，不能保证曲线在连接点处的光滑性。 Hermite 插值虽然在连接点处一阶光滑，但整体插值由于节点多次数高而有可能发生Runge现象。 若用分段三次Hermite插值，节点处一阶导数连续光滑，但在曲线的凹凸性变化比较大的地方，误差较大。 希望： 既想分段插值，又想在节点处保持光滑，甚至二阶光滑--三次样条。 定义：设有 n +1个点的点列 若函数 满足： (3) 在整个区间 内具有二阶连续导数。 (2) 在每个小区间 上是三次多项式； ---此时 叫插值函数； 则称 为点列的三次样条插值函数或三次样条多 项式，简称三次样条。 定义中的一阶导数连续意味着曲线没有急弯，二 阶导数连续意味着曲线每一点的曲率半径有意义。 在MATLAB中，三次样条插值运算实现如下： yi=interp1( x, y, xi, ’spline’) 或 yi=spline( x, y, xi ) 其中 x，y 都是向量形式的点，xi是进行插值的点的 横坐标向量，yi 为插值函数值。 例 求三次样条插值并作图比较。 %以s为参数，分别做x，y的样条函数 xx = [-1 -0.866 -0.5 0 0.5 0.866 1 1.0402 1.1500 1.3... 1.54 1.8280 2.1736 2.5883 3.0860]; yy = [0 -0.25 -0.433 -0.5 -0.433 -0.25 0 0.15 0.2598 0.3... 0.3 0.3 0.3 0.3 0.3]; s=1:length(xx);sp=1:length(xx)/100:length(xx); xp=spline(s,xx,sp); yp=spline(s,yy,sp); %以s为参数，分别对已知的两组数据 %点进行三次样条插值 plot(xp,yp); hold on plot(xx,yy, ro);xlabel(x); ylabel(y); axis([-1 3.5 -1 1]) %若不取s