第四章 函数的数值逼近..pptVIP

三、数据拟合的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.6 函数的最佳平方逼近 设有函数 ，总可以在一定意义下，确定的空间中求已知函数的近似，这类问题称为函数逼近。 例如 是n+1维的线性空间。在Pn上定义内积 其中 称为权函数。 若存在 ，使得 则称 为函数 的最佳平方逼近多项式。 考虑一般的最佳平方逼近问题 是一线性空间。 若存在 ，使得 则称 为函数 在Φ上的最佳平方逼近多项式。 即求解函数 的极值问题。该函数取得极值的必要条件是 例题 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 样条插值 对于给定函数表 yn … y1 y0 y=f(x) xn … x1 x0 x 表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)的表达式可以有多种方法。例如，可以直接利用分段三次Her