函数的数值逼近培训资料.ppt

三、数据拟合的Matlab实现 1、polyfit() p=polyfit(x,y,n) 其中x,y是数据点，n是拟合多项式次数，p是多项式系数（由高到低） 2、polyval() y=polyval(p,x)其中p是拟合多项式系数，y是拟合多项式在x点的值。 注意：当已知n个点处的函数值，而拟合多项式的次数为n-1时，则拟合函数polyfit()就是计算n-1阶插值多项式的值。 Least.m 应用实例：人口预测 概念：内插、外推 censusgui.m 例题 4 考虑在区间[-1,1]上的函数 构造等距节点上的Lagrange插值多项式。 目的 观察当插值节点逐渐增多时，插值多项式对原函数的拟合变化情况。 tryrunge.m §4.3 分段插值与保形插值 1、分段插值 引例 fenduanex2.m 给定n+1 个结点a=x0x1…xn =b上的函数值y0, y1,…,yn ，求一折线函数Ih(x)满足： 则称Ih(x)为分段线性插值函数。 由定义可知Ih(x)在每一个小区间[xk,xk+1]上可表示为 若用插值基函数表示，则在整个区间[a,b]上Ih(x)为 2、 保形插值 给定[a, b] 上的一种分划 及f(x) 在分划点上的函数值f(xk)=yk, 及f’(xk)=dk (k=0,1,…,n),则作一个分段三次Hermite 插值多项式 P(x) , 满足条件 ① P(x)∈C1[a, b]； ② P(xk)=yk, P’(xk)=dk, k=0,1,…,n； ③ P(x) 在每个小区间上是三次多项式。 根据前面的结果，不难作出各个点上的插值基函数。 dk未知？ dk的取法 记 第k个小区间的区间长度 第k个小区间分段线性插值斜率 如果 与 符号相同，则令 如果 与 符号相反，则令 极值点 保持原来曲线的走向 12 15 17 21 18 16 6 5 4 3 2 1 考虑函数表 生成它的保形插值多项式。 例题 5 Chap4_3.m §4.4 样条插值 对于给定函数表 x x0 x1 … xn y=f(x) y0 y1 … yn 表4.2 （其中 ），若函数S(x) ∈C2[a,b]满足条件： ① S(x) 在每个子区间[xi-1, xi] (i=1,2,…,n) 上都是不高于三次的多项式； ② S(x), S’(x), S’’(x) 在区间[a, b] 上连续； ③ S(xi)=yi (i=0,1,…,n). （A－1） 则称S(x) 为函数f(x) 关于结点x0, x1,…,xn 的三次样条插值函数。 条件①表明S(x) 是一个分段三次多项式，用Si(x) 表示S(x) 在第i 个子区间上的表达式，则Si(x) 形如 其ai, bi, ci, di 待定，子区间共有n 个，需确定4n 个系数。 另外，根据条件②③，有 共有4n-2个方程，显然要确定S(x) 还差两个条件。 ◆ 边界条件的取法 Ⅰ. Ⅱ. 特别地，当 时，称为自然边界条件。 Ⅲ . 当原函数f(x)以b-a 为周期的周期函数，则要求S(x)也是周期函数。 ◆ 样条插值函数的建立 构造满足插值条件（A-1）及相应边界条件的三次样条插值函数S(x)的表达式可以有多种方法。例如，可以直接利用分段三次Hermite插值，只要假定S’(xj)=mj (j=0,1,…,n),于是 其中 是插值基函数。 利用插值条件及相应的边界条件，则可得关于mj(j=0,1,…,n)的三对角方程组，求出mj则得到所求的三次样条函数S(x). 利用S(x)的二阶导数值S’’(x)=Mj(j=0,1,…,n)表示S(x) 交互实验：插值方法比较 目的：通过比较几种不同插值方法，使大家对四种插值效果在感性上有所认识，便于以后使用。 tryinterp.m 实例：考察某种纤维的强度与其拉伸倍数的关系,下表 是实际测定的24个纤维样品的强度与相应的拉伸倍数 是记录: §4.5 曲线拟合的最小二乘方法 纤维强度随拉伸 倍数增加而增加 并且24个点大致分 布在一条直线附近 ---------