计算物理第1讲插值和拟合.pptVIP

一、样条插值 三弯矩算法 (利用M(x) 的分段插值函数) 试验题1 线性模型和最小二乘拟合 多项式拟合 试验题2 Xu Zhongfeng, Xi’an Jiaotong University, 2008 Computational Physics School of Science, Xi’an Jiaotong University Xu Zhongfeng 第1讲 离散数据的插值和拟合 问题的提出 设对于区间[a,b]，给定一个划分：a= x0 x1 … xn=b，给定节点上的函数值f(xi) = yi ( i=0,1,…,n) ，如果函数S(x)满足条件： (1) 在[a,b]上具有二阶连续导数，即 (2) 在每个子区间[xi , xi+1] ( i=0,1,…,n-1)上是三次多项式 (3) 则称S(x)是y = f(x) 的三次样条插值函数 确定S(x)需要4n个条件, 而我们只给出了4n-2个条件 端点函数值：2个 内节点函数值及连续条件：2(n-1) 个 内节点一、二阶导数连续条件：2(n-1) 个 需补充2个边界条件，通常有以下几种： 1.转角条件： 即 2.弯矩条件： 3.周期边界条件： 且常采用自然边界条件 记： 则M(x)可表示为 利用分段函数 积分，得 S(x)可进一步写成 利用 为了计算S (x), 需要进一步求出Mi 。 S (x)的一阶导数为 由S (x)在内节点处一阶导数的连续条件 记 转角条件： 弯矩条件： ——三弯矩方程 周期边界条件： 记： 总结 ——三弯矩方程 转角条件： 弯矩条件： 周期边界条件： 对于三次样条插值的三弯矩方法，编制用于第一种和第二种边界条件的程序。设已知数据 0.7280 0.53 0.6708 0.45 0.6245 0.39 0.5477 0.30 0.5000 0.25 yi xi 第一种边界条件: s(0.25)= 1.0000, s(0.53)=0.6868 第二种边界条件: s(0.25)= s(0.53)=0 分别用所编程序求解，输出各插值节点的弯矩值{Mi}和插值中点的样条函数值，并作点列{xi , yi}和样条函数 y=s(x) 的图形。 问题提出 设给定一组离散数据： ，连续函数 y=f(x)显然是存在的，但不可能得到其确切关系。于是，寻找一近似表达式y=?(x)。为此，我们需解决两个问题： (1) 选择?(x)为何类型的函数？即，建立什么样的数学模型？如 当参数确定后， 称为均方误差。 点处与测得的函数值之差的平方和最小，即 (2) 对于选择的数学模型，如何确定其中的待定参数？ 线性模型：令?(x)= ?0+ ?1x+ ?2 x2，其中?0 ,?1和?2 是待定参数, ?(x) 是它们的线性函数； 非线性模型：令?(x)= ?e?/x, 其中?和?是待定参数, ?(x) 是它们的非线性函数。 通常采用最小二乘法，即选择?(x)中的参数，使得其在离散 二 、离散数据的拟合 对于已知的n+1对离散数据 和权函数 , 记 在连续函数空间 C[a,b] 中选定m+1 (mn) 个线性无关的基函数 记由它们张成的子空间为 若有 使得 则称 为离散数据 在子空间?中带权 的最小二乘拟合。 ?(x)在离散点处的值 记 这样，求极值问题 的解 等价于求多元二次函数I(?0, ?1,…, ?m )极小值点 使得 对于已知的n+1对离散数据 和权函数 , 选定m+1维连续函数空间?，如果其有一组基 在点列 处的值向量组 线性无关，则最小二乘问题 存在唯一解 , 其中 的解，且有平方差 为正则方程 其中 最常选用的数学模型是多项式 即在多项式空间 作曲线拟合。 选择子空间?的基函数为 其在自变量序列 处的值向量为 即可得到最小二乘问题的的解 由此，求出正则方程G? =d 的解 平方误差为 deta = 0.1351 arfa = -0.0024 0.2073 0.2001 deta = 0.0069 arfa = 0.0000 -0.0052 0.2594 0.0362 4.66 50 4.60 45 4.51 40 4.37 35 4.15 30 3.87 25 3.44 20 2.86 15 2.16 10 1.27 5 9