算法设计五.docVIP

算法设计五 ??? 一、用分段线性Lagrange插值多项式求函数值的近似值。 ??? （一）算法名称 ??? ??? （二）存贮方式 ??? ???????? 一维数组，存放节点 ??? ????????? 一维数组，存放节点上的函数值 ??? （三）分段线性Lagrange插值多项式算法 ??? 1. 输入 （节点个数减1）； ；插值点U ??? 2. 如果 ，则 车5。否则 ??? 3. 如果 则 转5。否则 ??? 4. 对于 ??? 如果 ，则 ??? 5. ?? ??????????????? ??? 6. 输出V ??? 二、非等距节点分段Newton插值 ??? 设给定函数 的函数表如下： ? ? ??? 并且 ，用分段m次Newton插值多项式计算 在 处的近似值，其中 的整数值）。 ??? （一）算法名称 ??? ??? （二）存贮方式 ??? ????????????? 一维数组，存节点 ??? 一维数组，存函数值 ??? ??二维数组，保存各阶均差，其中N是最后一个节点的下标值，M为计算的均差的最高阶数。 ??? 数组元素 表示以 为最后一个节点的 阶均差。即 ????????????? 由 ????????????? 有??? ??? 本算法主要组成部分及计算公式有 ??? （1）计算各阶均差。 ??? ?????????????????????????? ??? （2）确定表头节点下标 。 ??? 对于 ??? ?????????????? ??? （3）计算插值多项式的值P(U) ??? 据Newton插值公式，可得出以 为表头的计算公式如下： ??????????????????????????? 即???????????????????????? ??????????????????????????? （当 ） 及 ??????????????????????????? 即???????????????????????? ??????????????????????????? （当 ） ??? （三）非等距节点的分段M次Newton插值算法 ??? 1. 输入N，M， ??? 2. 对于 ????????????? ??? 3. 对于 ??? 对于 ??? ??? 4. 对于 ??? 4.1 如果 ??? 4.1.1如果 ，则 ??? 并转7，否则 ??? 4.1.2 ,转10 ??? 5. 如果 ，则 ??? 转10 ??? 6. 如果 则 ??? ，转10 ??? 7. ??? 8. 对于 ??? ??? 9. 转12 ??? 10. ??? 11. 对于 ??? ??? 12. 输出 ??? 13. 停机 ??? 四、算例 ??? 已知某直升飞机旋转机翼外形轮廓线部分坐标值如下表：??? xi 0.52 3.10 8.00 17.95 28.65 39.62 50.65 78 104.6 yi 5.28749 9.40 13.84 20.20 24.90 28.44 31.1 35 36.5 xi 156.6 208.6 260.7 312.5 364.4 416.3 468 494 507 520 yi 36.6 34.6 31.0 26.34 20.9 14.8 7.8 3.7 1.5 0.2 ??? 试用分段的Newton插值计算在 处的函数值。 ??? ??????????????????????????? 从而采用分段的三次Newton插值，输入 ， ， ， ，运行程序NEWTN（X，Y，u，n，m，A），输出各插值点上的数值结果如下： u 5 25 65 130 170 v 11.2480 23.45437 33.5074 36.8432 36.2650 u 220 320 380 440 500 v 33.9213 25.5844 19.1201 11.7780 2.6972 而该函数表的五阶均差满足 ?????????????????????????????????? 若采用四次Newton插值公式，即 输出各插值点上的数值结果如下： u 5 25 65 130 170 v 11.2571 23.4589 33.5340 36.8432 36.2724 u 220 320 380 440 500 v 33.9284 25.5991 19.0888 11.4585 2.6966 ??? 以上可以看出Newton插值法很容易在计算机上实现，而且要想将插值多项式的次数提高一次，只要重新输入插值多项式次数M即可。 ???