选矿过程模拟与优化第三章插值.pptVIP

IF (JI) THEN IF((ABS(X[i]-X[j]))=EPS) THEN BEGIN Writeln (′数据点有误！′)； GOTO EXIT END ELSE A：=A*（x1-x[j]）/(x[i]-x[j]) END; ; EXIT; END; END; 例：给定列表函数 求 时的五位数近似插值P(23)。 解： 20 22 24 26 0.34202 0.37461 0.40674 0.43837 作业3： 某厂原煤浮沉资料如下，试分别用线性插值，拉格朗日一元三点插值和拉格朗日插值计算密度为1.35, 1.45和1.55时的浮物累积产率和灰分。 δ －1.3 1.3～1.4 1.4～1.5 1.5～1.6 1.6～1.8 +1.8 y% 23.74 47.56 9.39 3.28 2.10 13.93 Ag% 5.76 10.11 19.31 28.75 38.96 80.18 §3-4 样条插值的应用 1．为何要建立样条插值 分段插值的缺点是在连接处不能保证光滑地过渡，如果采用其插值结果绘制曲线，则很难保证曲线光滑，用高阶插值虽可弥补这个缺点，但因次数越高、计算越繁，误差往往也大，要保证连接处有一定的光滑性，一般采用样条插值。 2．什么是样条插值 样条插值是一种改进的分段插值，它在每个由相邻节点组成的小区间，都构造一个三次或两次函数，为了在连接处保持光滑，在节点上保持一阶连续导数。 3．三次样条函数的建立 设给定节点的序列为： 其中， 要求构造一个函数 ，使其满足三次样条函数的三个条件： (1)函数通过相应的节点，即 (2)在区间[ ]上， 具有一阶和二阶连续导数。 (3)在每个小区间[ ]上， 是x的三次多项式。 根据上述条件可以导出三次样条函数表达式： 设三次样条函数为 ，在[ ]上， 若要满足插值条件，则三次多项式可构造为： ① 根据第一个条件，显然 将①式对x求导： ② 若用 表示 在 处的一阶导数，则： ， 将 和 分别代入②式得 式中 解上面方程组可得： 的二阶导数为： 于是有： 因此同样有： 在区间 的右端点 的二阶左导数为： 而在子区间 左端点 的二阶右导数为： 根据 在 连接点处二阶导数连续, 则有： 即: 式中 该式对所有连接点 均成立，由此可以得到N－1个含 的方程式。 将a、b值代入①式得 ： 式中： ——插值节点及函数值 ; ——系数， ; ——表示 在 处的一阶导数。 只要求出了 ，即求出了三次样 条函数。 可以通过二阶导数连续来求。 注意：二个边界点的二阶导数取 ， 即按自然样条处理。可得三条角线方程组： 补充：“追赶法”解三对角线方程组： 简记 其中： 当 时， (a) (b) (c) 解：首先将系数矩阵A分解为二个矩阵，即 A=LU 其中L为下三角矩阵，U为单位上三角矩阵，则 有： 其中 为待定系数。 比较上面矩阵两边和已知条件，可得分解的二个 矩阵系数。 ∴求解 等价于解二个三角形方程组 （1） ，求y （2） ，求x 从而得到解三对角线方程组的追赶法公式： 1．计算 的递推公式 2．解 3．解 我们将计算系数 及 过程称为“追”的过程，将计算方程组