第1章常用数值分析方法§3插值法与曲线拟合.pptVIP

§3 插值法与曲线拟合;3.1 实验数据统计处理;3.1.2 数据的统计分析 ; 函数常被用来描述客观事物变化的内在规律——数量关系，如宇宙中天体的运行，地球上某地区平均气温的变化等等，但在生产和科研实践中碰到的大量的函数中，不仅仅是用解析表达式表示的函数，还经常用数表和图形来表示函数，其中函数的数表形式在实际问题中应用广泛，主要原因是有相当一部分函数是通过实验或观测得到的一些数据，这些数据只是某些离散点 xi 上的值（包括函数值f (xi)，导数值 f?(xi)等,i = 1,2,…,n)，虽然其函数关系是客观存在的，但却不知道具体的解析表达式，因此不便于分析研究这类数表函数的性质，也不能直接得出其它未列出点的函数值，我们希望能对这样的函数用比较简单的表达式近似地给出整体的描述。; 另一方面，有些函数，虽然有解析表达式，但因其过于复杂，不便于计算和分析，同样希望构造一个既能反映函数的特性又便于计算的简单函数，近似代替原来的函数。 如在积分 中，当f (x)很复杂，要计算积分 I 是很困难的，构造近似函数使积分容易计算，并且使之离散化能上机计算求出积分I，都要用到插值逼近。; 解决上述问题的方法有两类：一类是对于一组离散点(xi,f (xi)) (i = 1, 2, …,n)，选定一个便于计算的函数形式?(x)，如多项式，分段线性函数，有理式，三角函数等，要求?(x)通过点?(xi)=f (xi) (i = 1, 2,…,n),由此确定函数?(x)作为f (x)的近似。这就是插值法。这里的 g(x) 称为f(x) 的插值函数。最常用的插值函数是 …? ; 已知：一系列离散的（互不相同的）点xi ， yi（i = 1，2，…n） 求：给定点 x 对应的函数值 y 或近似函数表达式。;3.2.2 线性插值 ;3.2.3 抛物线插值 ;3.2.4 Lagrange插值的一般形式 ;3.2.5 分段插值（分段抛物线插值） ;3.3 曲线拟合?;两类问题;3.3.2 求经验公式—多项式拟合;3.3.2.2 最小二乘技术 ;skj;3.3.2.3 多项式的次数（m）的确定 ;3.3.2.4 多项式拟合举例; 显然，上述方程组中只有两个未知参数，任取两个方程即可解出a, b。其几何意义就是任选两点，连成一线，故该法又称“选点法”。然而，由于试验误差的存在，这些点并不会在同一条直线上，因此随着选点不同，得到的a, b的值也有差异，如：;仍设试验数据满足直线关系：;（3）进而得到线性方程组：;两类问题; ;3.3.3.1 线性模型的参数确定 ;利用最小二乘技术求出模型中的各待定参数aj（j = 0, 1, 2, …, m） ;m+1元线性方程组;3.3.3.2 简单非线性模型的参数确定 ;A;例2. 已知一活化极化体系的阳极极化过程符合如下规律： 现通过实验测得其n组阳极极化数据为：ΔEi、iai（i = 1, 2, …, n），求该体系腐蚀速度icor和Tafel斜率ba。 ;3.3.3.3 一般非线性模型的参数确定 ;迭代的最小二乘技术（以充电曲线方程式的拟合为例 ）;b) 近似处理变超越方程为线性方程;c) 用最小二乘法拟合线性模型，求出其中的参数u、v、z ;迭代的最小二乘技术 高斯—牛顿曲线拟合法。 ;作 业;2、编制最小二乘多元线性回归的通用程序，上机调试，并对以下试验数据进行回归分析（误差限自定）。;3、分别编制Lagrange分段抛物线插值和Lagrange一般多点插值的通用程序。 已知CO2热容试验数据如表3.4所示，试用上述程序求出在1708K时之热容。