第4章 - 函数的数值逼近.ppt

第四章 函数的数值逼近 引言 代数多项式插值 分段线性插值与“保形”插值 三次样条函数插值 曲线拟合的最小二乘法 §1 引言 (1)复杂函数的计算; (2)函数表中非表格点计算 (3)光滑曲线的绘制; (4)提高照片分辩率算法 (5)定积分的离散化处理; (6)微分方程的离散化处理; (7)积分方程的离散化处理; 2、插值多项式的存在唯一性 一、线性插值与抛物线插值 二、Lagrange 多项式插值(n次) 定理（Lagrange）插值多项式 三、插值多项式的余项 例 §4 三次样条插值 二、样条函数定义 ?Chebyshev多项式 常用 正交多项式(续1) 是区间[-1,1]上关于权函数?(x)= 的正交多项式,且满足: a) b) 有递推公式 c) Tn(x)在[-1,1]上的n个零点为 ?Hermite多项式 常用 正交多项式(续2) 是区间[- ?,+ ?]上关于权函数?(x) = 的正交多项式,且满足: a) b) 有递推公式 3. 正交多项式数据拟合 已知点集 且 数表 步骤 a) 选取Hn中正交基 选取Hn中关于点集X以及权系数 的正交多项式组 即 b) 求解aj*得到Pn*(x) 接下页 步骤 c) 计算最小平方误差 接下页 于是，随n增加 优点：增加一个节点，只需计算an+1* ，计算量小 利用多项式{?0(x), ?1(x),…, ?m(x)} 的离散正交性易知，此时法方程组（41）的系数矩阵G为对角阵 。 如此，无需解线性方程组，可减少含入误差，避免病态情况出现。 利用离散正交多项式求所给数据表的二次拟合多项式 例 解：按三项递推式及αk , ?k的计算公式可得： 而由系数ak的计算公式有 ： ?i 1 2 3 4 5 6 7 8 9 ?xi -1 -0.75 -0.5 -0.25 0 0.25 0.5 0.75 1 ?yi -0.2209 0.3295 0.8826 1.4392 2.0003 2.5645 3.1334 3.7601 4.2836 无需解正规方程组，只需要计算内积，避免出现病态方程组的可能。并且当逼近次数增加一次时，只需在原有的多项式中增加一项，即在上例还可求三次、四次拟合多项式。 由此可得最小二乘二次拟合多项式为： (2) 求解方程组(26)(或(29)),求 。 存在唯一性 且 定理 （三次样条插值函数存在唯一） 三次样条插值函数f (x) ，且满足给定的边界条件。 计算步骤： 先计算(26)式中的 （若是第二类 或第三类边界条件，要计算 ） (3) 用(20)及(21)式进行插值计算 (先确定x所在区间) (1)如果f(x)是定义在[a, b]上的函数，且已知y=f (x)有函数表 (2)给定边界条件(a)(或(b)或(c))，则f(x)在[a, b]存在唯一的 因为Sj(x)是三次样条插值函数，所以 是一次函数。 2、二阶导数表示 由两点Lagrange插值得 参数 对上式积分,得 再积分,得 (30) (31) (32) 令 令 由条件 ，确定积分常数 即得 将上式代入(32)得到3次样条插值函数的表达式 (33) (34) 将(33)代入(31),得 由条件 满足的线性 方程组 两边同除以 (35) (36) 两个条件。 上式有n-1个方程,要确定n+1个未知量 需增加 ——三弯矩方程 (37) 则令j=0，由(35),得 令j=n，由(36),得 (1) 若 已知， n-1个方程 (38) (2) 若 已知，代入方程(37)，只 需解n-1个方程 (39) 满足方程 计算步骤 : 同三次样条插值函数的一阶导数表示的计算步骤。 说明： 此方程组(38)和(39)有唯一解{Mj} 。 （2）三弯矩方程(37 )与(25)比较， 与 交换了位置。 （1） 方程组(38)和(39)的系数矩阵都是严格对角占优矩 阵，因 （3） Mj 在力学上为细梁在 xj 处截面处的弯矩, 且弯矩与相邻 的两个弯矩有关, 故方程组(38)和(39)称为三弯矩方程。Mj在 数学上称为曲率。而方程(25)叫做三转角方程。Mj 在力学上 叫做细梁上 xj 处的转角，在数学上叫做变化率。 一、问题的提出 1．函数逼近 设f(x)为[a, b]上的连续函数，寻求一个近似函数 P(x) ，在[a, b]上均匀逼近 f(x) 。 2．数据逼近（实验数据） 已知 , 求n次多项式 Pn(x), nm, 使Pn(x) 能