实验二、函数的插值法.doc

实验二、函数的插值法 一、插值算法的基本思想： 函数的变化规律往往是通过一组实验数据给出，为了研究此变化规律往往需要求出不在表上的一些函数值。因此我们希望根据给定的函数表做一个既能反函数的特性，又便于计算的简单函数，用近似。通常选一类较简单的函数作为，并使这样确定的就是我们希望得到的插值函数，主要有①Lagrange 插值；②Newton插值；③Hermit插值；④分段线性插值；⑤三次样条插值。 二、实验要求： 1、利用Lagrange插值公式 编写插值多项式程序。 2、给出分段二次和三次多项式的一般表达式，并根据节点选取规则对上述两种插值法进行比较分析。 3、对此插值问题若用牛顿插值多项式结果如何？牛顿插值多项式公式如下： 三、目的和意义： 1、学会常用的插值计算方法以求函数的近似表达式， 以解决其它科学实验的计算问题； 2、通过比较加深理解插值多项式与分段插值多项式的优劣问题； 3、熟悉各种插值方法的程序设计； 4、通过画出的各种插值曲线比较起插值的光滑度问题。 四、算法步骤： 1、Lagrange 插值 输入N各节点数组； 定义初始值和； 3)利用公式 求 N 次插值基函数； 4)利用多项式求 插值函数; 流程图： 开始 、 否 是 输出 五、MATLAB程序源代码 1、Lagrange 插值 数值实验例题 数值插入：x=0:3:10 　y=sin(x) 运行： 数值例题 % p60 clear x=[.25 .3 .39 .45 .53]; y=[.5 .5477 .6245 .6708 .728]; [F,a]=scyt(x,y) x=[.25 .3 .39 .45 .53]; y=[.5 .5477 .6245 .6708 .728]; [F,a]=scyt(x,y) x=[27.7 28 29 30]; y=[4.1 4.3 4.1 3]; [F,a]=scyt(x,y) n=7; x=linspace(-pi,pi,n) y=sin(x); [F,a]=scyt(x,y) 或采用： function [S,y]=bkj(x0,y0,p,x,s1,s2) (注释: 第一种边界条件) X=[27.75 28.85 29.95] Y= [4.2222 4.0877 3.1888] 得到的S(x)= CPU 时间: 0.0250s 结果分析：