数值分析 一元插值.ppt

0－阶差商1－阶差商2－阶差商k－阶差商第30页，共42页，星期日，2025年，2月5日数值分析数值分析一元插值第1页，共42页，星期日，2025年，2月5日第十二讲

插值与逼近

第五章插值与逼近第2页，共42页，星期日，2025年，2月5日插值和逼近是解决什么问题的？用某个简单函数在满足一定条件下，在某个范围内近似代替另一个较为复杂或者解析表达式未给出的函数，以便于简化对后者的各种计算或揭示后者的某些性质。第3页，共42页，星期日，2025年，2月5日(1)复杂函数的计算;(2)函数表中非表格点计算;(3)光滑曲线的绘制;(4)提高照片分辩率算法;(5)定积分的离散化处理;(6)微分方程的离散化处理;(7)积分方程的离散化处理;插值方法的应用:3/18第4页，共42页，星期日，2025年，2月5日例如：在工程实践和科学实验中，经常需要建立函数关系，即y=f(x)。但通常只能观测到它的部分信息，这些值构成了观测数据：而不知道函数在其他点x上的取值，这时只能用一个经验函数y=g(x)对真实函数y=f(x)作近似。xix1x2…xnf(xi)f(x1)f(x2)…f(xn)第5页，共42页，星期日，2025年，2月5日插值的任务就是由已知的观测点(xi,yi)为物理量(未知量),建立一个简单的、连续的解析模型g(x)，以便能根据该模型推测该物理量在非观测点处的特性。oxy●●●●●y0x1x2xny1y2ynx0y=f(x)g(x)第6页，共42页，星期日，2025年，2月5日机翼断面的下轮廓线如图所示，下表给出了下轮廓线上的部分数据。用程控铣床加工时，每一刀只能沿x方向和y方向走非常小的一步。例如，如果工艺要求铣床沿x方向每次只能移动0.1单位，这时就需要求出当x坐标每改变0.1单位时的y坐标。试完成加工所需的数据，画出曲线。例：机床加工xi035791112131415yi01.21.72.02.12.01.81.21.01.6第7页，共42页，星期日，2025年，2月5日第8页，共42页，星期日，2025年，2月5日插值法:由实验或测量的方法得到所求函数y=f(x)在互异点x0,x1,...,xn处的值y0,y1,…,yn构造一个简单函数F(x)作为函数y=f(x)的近似表达式 y=f(x)?F(x)使F(x0)=y0,F(x1)=y1,?,F(xn)=yn(a)这类问题称为插值问题。f(x)称为被插值函数，F(x)称为插值函数，x0,x1,...,xn称为插值节点。(a)式称为插值条件。函数插值的基本问题第9页，共42页，星期日，2025年，2月5日插值函数的类型第10页，共42页，星期日，2025年，2月5日第11页，共42页，星期日，2025年，2月5日代数插值一元函数插值（一元Lagrange插值）二元函数插值（二元Lagrange插值）Hermite插值样条插值第12页，共42页，星期日，2025年，2月5日解法：设所求n次多项式为pn(x)=a0+a1x+…+anxn,由pn(xk)=yk我们有问题:给定n+1对互异的点对(xi,yi)i=0,1,2,…,n，求一个次数不高于n的多项式pn(x)，使得pn(xi)=yi,i=0,1,2,…,n一元函数插值第13页，共42页，星期日，2025年，2月5日ai(i=0,1,2,…,n)的系数行列式是Vandermonde行列式由于xi互异，所以上式右端不为零，从而方程组的解a0,a1,…,an存在且唯一。但遗憾的是方程组是病态方程组时,阶数n越高，病态越严重。为此我们从另一途径寻求获得Pn(x)的方法----Lagrange插值和Newton插值。（这两种方法称为基函数法）。第14页，共42页，星期日，2025年，2月5日函数组在点集上线性无关的概念：设函数组{φk(x),k=1,2,…}是次数不高于n的多项式组，称在点集{x0,x1,x2,…,xm}(n=m)上线性无关，如果向量组线性无关。第15页，共42页，星期日，2025年，2月5日问题的严格数学提法:给定