数值分析复习(二)拉格朗日插值法、插值余项与误差估计

拉格朗日插值法

数值分析复习(一)线性插值、抛物线插值中我们讨论过线性插值与二次插值,其实都是接下来要讲的拉格朗日插值的特殊情况,接下来我们一一分析:

定义插值基函数:

n次多项式l_j(x)(j=0,1,\cdots,n)在n+1个节点x_0<x_1<\cdots<x_n上满足条件:

l_j(x_k)= \left\{ \begin{array}{l} 1,k=j, \\ 0,k\neq j,\\ \end{array} \right. j,k=0,1,\cdots,n

就称这n+1个n次多项式l_0(x),l_1(x),\cdots l_n(x)为节点x_0,x_1,\cdots,x_n上的n次插值基函数。

l_k(x)=\frac{(x-x_0) \cdots (x-x_{k-1})(x-x_{k+1})\cdots (x-x_n)}{(x_k-x_0) \cdots (x_k-x_{k-1}) (x_k-x_{k+1} \cdots (x_k-x_n)},k=0,1,\cdots,n

引入记号:

拉格朗日插值多项式可变换为:

L_n(x)=\sum_{k=0}^{n}y_k l_k(x)=\sum_{k=0}^{n}y_k \frac{ \omega_{n+1}(x)}{(x-x_k)\omega_{n+1}^{'}(x_k)}

当n=1时,L_1(x)=y_0\frac{x-x_1}{x_0-x_1}+y_1\frac{x-x_0}{x_1-x_0},为线性插值

当n=2时,L_2(x)=y_0\frac{(x-x_1)(x-x_2)}{(x_0-x_1)(x_0-x_2)}+y_1\frac{(x-x_0)(x-x_2)}{(x1-x_0)(x_1-x_2)}+y_2\frac{(x-x_0)(x-x_1)}{(x_2-x_0)(x_2-x_1)},展开后可得抛物线插值

 

注:n次插值多项式L_n(x)通常是次数为n的多项式,特殊情况下次数可能小于n,如当二次插值多项式插值的三点共线时L_2(x)将退化为一次多项式

插值余项与误差估计

R_n(x)=f(x)-L_n(x),x\in [a,b]为插值多项式的截断误差,也称余项

有如下定理:

通过余项表达式我们可以知道,若插值函数f_(x) \in H_nH_n代表次数小于等于n的多项式集合),由于f^{(n+1)}(x)=0,故R_n(x)=f(x)-L_n(x)=0,即它的插值多项式为其本身。 

实验一 误差分析 一、实验目的及要求 1.了解误差分析对数值计算的重要性。 2.掌握避免或减小误差的基本方法。 、实验设备 安装有C、C++或MATLAB的计算机。 三、实验原理 误差是指观测值真值之差,偏差是指观测值平均值之差。根据不同的算法,得到的结果的精度是不一样的。 四、实验内容及步骤 求方程ax2+bx+c=0的根,其中a=1,b= -(5×108+1),c=5×108 采用如下两种计算方案,在计算机上编程计算,将计算结果记录下来,并分析产生误差的原因。 ////////////////////////////// 实验 Lagrange插值 一、实验目的及要求 1.掌握利用Lagrange插值法及Newton插值法求函数值并编程实现。 2.程序具有一定的通用性,程序运行时先输入节点的个数n,然后输入各节点的值( ),最后输入要求的自变量x的值,输出对应的函数值。 、实验设备和实验环境 安装有C、C++或MATLAB的计算机。 三、算法描述 1. 插值的基本原理(求解插值问题的基本思路) 构造一个函数y=f(x)通过全部节点,即 (i=0、1、… n) 再用f(x)计算插值,即 2. 拉格朗日(Lagrange)多项式插值 Lagrange插值多项式: 3.牛顿(Newton)插值公式 //////////////////////////////////// 实验三 高斯消去法解方程组 一、实验目的及要求 1.掌握求解线性方程组的高斯消去法---列选主元在计算机上的算法实现。 2.程序具有一定的通用性,程序运行时先输入一个数n表示方程含有的未知数个数,然后输入每个线性方程的系数和常数,求出线性方程组的解。 、实验设备和实验环境 安装有C、C++或MATLAB的计算机。 三、算法描述 1.高斯消去法基本思路 设有方程组 ,设 是可逆矩阵。高斯消去法的基本思想就是将矩阵的初等行变换作用于方程组的增广矩阵 ,将其中的 变换成一个上三角矩阵,然后求解这个三角形方程组。 2. 利用列选主元高斯消去法求解线性方程组
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值