正在加载图片...
实验10样条插值 参考答案 1.各边界条件的插值样条生成函数见代码文件: (1)自然边界样条:functionS=-cnsfit(X,Y) (2)端点曲率调整样条:function S=cepcasfit(X,Y,dx20,dx2n) (3)紧压样条:functionS=csfit(X,Y,dx0,dxn) (4)外推样条:function S=cesfit((X,Y) (5)抛物线终结样条:functionS=cptsfit(X,Y) 利用上述各函数计算及绘图,得到如下结果: (1)自然边界样条计算及绘图代码段,输出如图1。 x05 y=x.*exp(-x). nS-=cnsfit(x,.y) x1=0:.01:1:y1=polyval(nS(I.).x1-x(1)): x2=1:.01:2y2=polyval(nS(2,),x2-x2 x3=2:.01:3y3=polyval(nS(3,),x3-x(3) x4=3:.01:4y4=polyval(nS(4,),x4-x(4) x5-4:.01:5:y5-polyval(nS(5.).x5-x(5)). figure,plot(x1,y1,x2.y2,x3.y3.x4.y4,x5.y5,x.y,'ro),title('natural spline) (2)端点曲率调整样条代码段(与(1)类似),输出如图2。 (3)紧压样条代码段(与(1)类似),输出如图3。 (4)外推样条代码段(与(1)类似),输出如图4。 (5)抛物线终结样条代码段(与(1)类似),输出如图5。 实验 10 样条插值 参考答案 1. 各边界条件的插值样条生成函数见代码文件: (1)自然边界样条:function S=cnsfit(X,Y) (2)端点曲率调整样条:function S=cepcasfit(X,Y,dx20,dx2n) (3)紧压样条:function S=csfit(X,Y,dx0,dxn) (4)外推样条:function S=cesfit(X,Y) (5)抛物线终结样条:function S=cptsfit(X,Y) 利用上述各函数计算及绘图,得到如下结果: (1) 自然边界样条计算及绘图代码段,输出如图 1。 x=0:5; y=x.*exp(-x); nS=cnsfit(x,y); x1=0:.01:1;y1=polyval(nS(1,:),x1-x(1)); x2=1:.01:2;y2=polyval(nS(2,:),x2-x(2)); x3=2:.01:3;y3=polyval(nS(3,:),x3-x(3)); x4=3:.01:4;y4=polyval(nS(4,:),x4-x(4)); x5=4:.01:5;y5=polyval(nS(5,:),x5-x(5)); figure,plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5,x,y,'ro'),title('natural spline') (2) 端点曲率调整样条代码段(与(1)类似),输出如图 2。 (3) 紧压样条代码段(与(1)类似),输出如图 3。 (4) 外推样条代码段(与(1)类似),输出如图 4。 (5) 抛物线终结样条代码段(与(1)类似),输出如图 5
向下翻页>>
©2008-现在 cucdc.com 高等教育资讯网 版权所有