数学建模插值法.ppt

插值法 问题 一水库上游河段降暴雨，根据预报测算上游流入水库的流量为Q(t) (102立方米/秒) ： t (时) 8 12 16 24 30 44 48 56 60 Q（t） 36 54 78 92 101 35 25 16 13 通过这个预报值，估计14 和 20 时上游流入水库的流量. 一、问题的提出 在生产和实验中，关于函数f(x)，经常存在两种情况： （1）其表达式不便于计算； （2）无表达式， 而只有函数在给定点的函数值， 怎样计算其它点的函数值？ 二、问题的解决 两种方法： （1）函数插值； （2）曲线拟合． 三、插值法 定义：当精确函数 y = f(x) 非常复杂或未知时，在一系列节点 x0 … xn 处测得函数值 y0 = f(x0), … ，yn = f(xn)， 由此构造一个简单易算的近似函数 p(x) ? f(x)，满足条件 p(xi) = f(xi) (i = 0, … n)------ 插值条件 这里的 p(x) 称为f(x) 的插值函数； 构造插值函数的方法为插值法。 多项式插值法包含： 拉格朗日插值 牛顿插值 三次埃尔米特插值法 分段线性插值 分段三次埃尔米特插值法 三次样条插值 1、 拉格朗日插值公式 （１）定义 对给定的n+1个节点x0 , x1,x2,…,xn及对应的函数值y0 , y1,y2,…,yn， 构造一个n次插值多项式： 即为拉格朗日插值公式，其中 （2）拉格朗日插值的matlab实现 function y=lagrange(x0,y0,x) % x0插值节点， y0插值节点处的函数值，x要计算函数值的点； n=length(x0); %计算x0的长度 m=length(x); %计算x的长度 for i=1:m s=0; z=x(i); for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); %计算插值基函数 end end s=p*y0(k)+s; end y(i)=s; %计算在x(i)处的函数值（拉格朗日） end （3） 函数的调用 给定数据表： 用lagrange插值公式计算 clear x0=0.4:0.1:0.8； y0=log(x0)； x=0.54； lagrange(x0,y0,x) ln(0.54)的函数值 ans = -0.6161 （4）练习 2、牛顿插值法 牛顿插值公式： Nn(x)=f(x0)+f[x0,x1](x-x0)+…+f[x0,x1,x2,…,xn](x-x0)(x-x1)…(x-xn) 其中： f[x0,x1] 一阶差商 f[x0,x1,x2,…,xn] n阶差商 注：牛顿插值法与拉格朗日插值法，同一个多项式，不同的表达方式，但是计算量不一样，牛顿插值法的计算量小。 这两种方法在MATLAB 中没有现成的命令和函数，因为它们存在龙格现象。 牛顿插值法的matlab实现 function Y=newpoly(x,y,X) %牛顿插值 n=length(x); D=zeros(n); D(:,1)=y; for j=2:n %计算差商矩阵 for k=j:n D(k,j)=(D(k,j-1)-D(k-1,j-1))/(x(k)-x(k-j+1)); end end C=D(n,n); for k=(n-1):-1:1 %构造插值多项式 C=conv(C,poly(x(k))); m=length(C); C(m)=C(m)+D(k,k); end Y=polyval(C,X); 龙格现象 Runge在上个世纪初发现： 在[-5,5]上用n+1个等距节点作n次插值多项式Pn(x)，