matlab数值分析插值法
1拉格朗日插值法functionyh=lagrange(x，y，xh ) n=length(x ) x； m=length(xh； yh=Zeros(1，m )； for j=1:m； fori=1: nxp=x ([ 1: I-1i 1: n ]； yh(j )=yh ) j ) y(I ) *prod ) ) xh(j )-xp )./(x ) I-XP ) )； %yh和y endend 调用程序
x=[ 11，12，13 ]； y=[2.3979、2.4849、2.5649]； xh=11.75； yh=lagrange(x，y，xh )2 gjdyl插值法
功能yh=Newton pol (x，y，xh ) n=length(x ) x； p 65:1 )=x； p (3360，2 )=y； forj=:n1p(1:n2-j，j )=diff ) 1:n3-j，j-1 )./(x(j-1:n )-x (1: N2-j ) ) ) %有求差商表(此处)符号，注意与差商表的不同之处) endq=p (1，2: n1 )； 求出gjdyl法系数---取第1行的yh=0； m=1； yh=q(1； forI=2:nm=q(I； forj=2:Im=m*(xh-x(j-1 )； gjdyl法中各多项式的值(xh-x0 )…(xh-xn-1 ) end yh=yh m； %合计end 调用程序
x=[ 11，12，13 ]； y=[2.3979、2.4849、2.5649]； xh=11.75； yh=Newtonpol(x，y，xh ) )。