科学计算0445.pptVIP

一、Lagrange插值多项式 思考: 如何编程计算Lagrange插值多项式? 算法 (类程序的形式) Lagrange (n, xn , yn , x ) 1. L=0.0 2. for i=1, ···,n 2.1 l =1.0 2.2 for j=1,2, ···,n if j ≠ i then l =l* (x-xj)/(xi - xj) 2.3 L=L+l yk 3. 输出”L(x)=L” 二、拟合多项式 * 下页 上页 可知, 除xi 点外, 其余节点xj处都是0,故可设 考虑下面最简单` 最基本的插值问题。求n-1 次多项式 l i(x) (i=1, ···, n)，使其满足条件 其中A为常数, 由li (xi)=1可得 , 所以有 称之为n-1次拉格朗日基函数. 设 可知其满足 注:Lagrange插值多项式优点是形式整齐对称,便于记忆,易于编程,缺点是若补充或删去一个节点时,全盘皆动,原来的插值基函数都不可用. 1. Lagrange插值函数m文件为lag.m function Lx=L