函数逼近的插值法Hermite插值多项式教学幻灯片讲义.ppt

函数逼近的插值法——Hermite插值多项式;Lagrange插值公式所求得L(x)保证了节点处的函数值相等，也就是保证了函数的连续性。 但不少实际问题还需要插值得光滑度，也就是还要求它在节点处的导数值也相等，导数的阶数越高则光滑度越高。 现代的仿生学就是一个典型的例子。在设计交通具的外形，就是参照海豚的标本上已知点及已知点的导数，做插值在计算机上模拟海豚的外形制成飞机、汽车等外形。;Hermite插值问题的提法;Hermite插值多项式的求法—Lagrange方法;;;;;;;;;;;;Hermite插值余项;三次Hermite插值多项式(n=1);;;;;;x=[0,pi/6];y=sin(x);z=cos(x);u=pi/12; d=hermitchazhi(x,y,z,u) d = 0.25876861681747 e=abs(sin(u)-d) e = 5.042828505269492e-005 ;Hermit插值问题的一般提法;;这样的插值问题称为Hermite 插值问题;Hermite插值多项式的求法——Newton方法;;;;;;;分段多项式插值;Runge现象;;[a,z,z1]=divided1(5);plot(a,z,r,a,z1,g);[a,z,z1]=divided1(10);plot(a,z,r,a,z1,g);[a,z,z1]=divided1(15);plot(a,z,r,a,z1,g);[a,z,z1]=divided1(20);plot(a,z,r,a,z1,g);;分段线性插值;;;;j=find(x-u0,1), 判断xj-1uxj;x=[1:0.1:7];y=sin(x);u=5.5;v=fenduanlinear(x,y,u) v = -0.70554032557039 e=abs(sin(u)-v) e = 0;分段三次Hermite插值;;;分段三次Hermite插值算法;x=[1:4];y=sin(x);z=cos(x);u=3.3; v=fenduanhermit(x,y,z,u) v = -0.15718060155343 e=abs(sin(u)-v) e = 5.650925898138537e-004; x=[1:0.1:4];y=sin(x);z=cos(x);u=3.3; v=fenduanhermit(x,y,z,u) v = -0.15774569414325 e=abs(sin(u)-v) e = 0;样条函数;;;;;;;;;;;;;;;;;;;;;;;;;;;;;t=linspace(0,2.25,10); y=sqrt(t); a=linspace(0,2.25,40); d = sspline1(t,y,a); z=sqrt(a); plot(a,d,o,a,z);;zz=z-d; ;plot(a,zz);三次样条插值;;;;;;;;;;;;;例题;;;;;[z,d]=sspline2(t,y,0.2,e,f) z = -2.0445 -1.7776 -1.1303 -0.4371 0.0841 d =0.9615;上机实习题