多项式插值中的一些定理证明
多项式插值中的一些定理证明
期末复习中整理的一些证明
- 多项式插值定理
- 多项式差值误差定理
- Hermite插值误差估计定理
- Hermite插值存在唯一性定理
多项式插值定理
若x0,x1,⋯,xnx_0,x_1,\cdots,x_n是不同实数,则对任意的数值y0,y1,⋯,yny_0,y_1,\cdots,y_n,存在唯一的次数最多是nn次的多项式满足**
p(xi)=yi,i=0,1,⋯,np(x_i)=y_i,i=0,1,\cdots,n
证明:设p(x)=a0+a1x+⋯+anxnp(x)=a_0+a_1x+\cdots+a_nx^n,则有方程:
a_0+a_1x_i+\cdots+a_nx_i^n=y_i,i=0,1,\cdots,n
将它们写成矩阵的形式
\left( \begin{matrix} 1&x_1&\cdots&x_1^n\\ 1&x_2&\cdots&x_2^n\\ \vdots&\vdots&\ddots&\vdots\\ 1&x_n&\cdots&x_n^n\\ \end{matrix} \right) \left( \begin{matrix} a_0\\ a_1\\ \vdots\\ a_n\\ \end{matrix} \right)= \left( \begin{matrix} y_0\\ y_1\\ \vdots\\ y_n\\ \end{matrix} \right)
又
\left| \begin{matrix} 1&x_1&\cdots&x_n\\ 1&x_2&\cdots&x_n\\ \vdots&\vdots&\ddots&\vdots\\ 1&x_n&\cdots&x_n\\ \end{matrix} \right|= \prod_{i\neq j}(x_i-x_j) \neq0
由线性代数的知识可知,存在唯一一组 a0,a1,⋯,ana_0,a_1,\cdots,a_n满足以上条件,证毕。
多项式插值误差定理
f(x)∈Cn+1[a,b]f(x)\in C^{n+1}[a,b],p(x)p(x)是ff在区间[a,b][a,b]上n+1n+1个不同点x0,x1,⋯,xnx_0,x_1,\cdots,x_n上的次数不超过nn插值多项式,则对每个x∈[a,b]x \in [a,b]都存在ξx∈[a,b]\xi _x \in [a,b],满足:
f(x)−p(x)=1(n+1)!f(n+1)(ξx)∏i=0n(x−xi)f(x)-p(x)=\frac{1}{(n+1)!}f^{(n+1)}(\xi _x)\prod_{i=0}^n(x-x_i)
证明:若xx是插值节点中的一个,上式显然成立,下设x≠xix \neq x_i。设ϕ(t)=f(t)−p(t)−λω\phi (t)= f(t) -p(t)-\lambda \omega,其中ω=∏ni=0(t−xi),λ=f(x)−p(x)ω(x)\omega = \prod_{i=0}^n(t-x_i),\lambda = \frac{f(x)-p(x)}{\omega(x)},则ϕ(t)\phi(t)至少有x,x0,x1,⋯,xnx,x_0,x_1,\cdots,x_n这n+2n+2个不同的根,由罗尔定理,ϕ(n+1)\phi^{(n+1)}至少有一个零点ξx\xi _x,即:
\phi^{(n+1)}(\xi _x)=f^{(n+1)}(\xi _x)-\lambda(n+1)!=0
可以得到:
f(x)-p(x)=\frac{1}{(n+1)!}f^{(n+1)}(\xi _x)\prod_{i=0}^n(x-x_i)
Hermite插值误差估计定理
f(x)∈C2n+2[a,b]f(x)\in C^{2n+2}[a,b],x0,x1,⋯,xnx_0,x_1,\cdots,x_n上是区间[a,b][a,b]上n+1n+1个不同点,次数不超过2n+12n+1的多项式p(x)p(x)满足:
p(xi)=f(xi),p′(xi)=f′(xi)p(x_i)=f(x_i),p'(x_i)=f'(x_i)
则对每个x∈[a,b]x \in [a,b]都存在ξx∈[a,b]\xi _x \in [a,b],满足:f(x)−p(x)=1(2n+2)!f(2n+2)(ξx)∏i=0n(x−xi)2f(x)-p(x)=\frac{1}{(2n+2)!}f^{(2n+2)}(\xi _x)\prod_{i=0}^n(x-x_i)^2
证明:若xx是插值节点中的一个,上式显然成立,下设x≠xix \neq x_i。设ϕ(t)=f(t)−p(t)−λω\phi (t)= f(t) -p(t)-\lambda \omega,其中ω=∏ni=0(t−xi)2,λ=f(x)−p(x)ω(x)\omega = \prod_{i=0}^n(t-x_i)^2,\lambda = \frac{f(x)-p(x)}{\omega(x)},则ϕ(t)\phi(t)至少有x,x0,x1,⋯,xnx,x_0,x_1,\cdots,x_n这n+2n+2个不同于的根,由罗尔定理,ϕ′\phi'至少有n+1n+1个不同于x,x0,x1,⋯,xnx,x_0,x_1,\cdots,x_n的根,又插值节点也都是ϕ′\phi'的根,则ϕ′\phi'至少有2n+22n+2个不同的根,由罗尔定理ϕ(2n+2)\phi^{(2n+2)}至少有一个零点ξx\xi _x,即:
\phi^{(2n+2)}(\xi _x)=f^{(2n+2)}(\xi _x)-\lambda(2n+2)!=0
可以得到:
f(x)-p(x)=\frac{1}{(2n+2)!}f^{(2n+2)}(\xi _x)\prod_{i=0}^n(x-x_i)^2
Hermite插值的存在唯一定理
存在唯一的次数至多是mm的多项式p(x)p(x)满足:
p(j)(xi)=cij,(j=0,1,⋯,ki−1 i=0,1,⋯,n)p^{(j)}(x_i)=c_{ij},(j=0,1,\cdots,k_i-1 \ i=0,1,\cdots,n)
m+1=k0+k1+⋯+knm+1=k_0+k_1+\cdots+k_n
证明:根据已知条件,可以列出一个由m+1m+1个未知数(即多形式的系数)和m+1m+1个方程组成的方程组,由线性代数的知识可知,方程组存在唯一的解等价于对应的齐次方程只有零解,即:
p^{(j)}(x_i)=0,(j=0,1,\cdots,k_i-1 \ i=0,1,\cdots,n)
只有零解。则 xix_i至少是 kik_i次根,设:
p(x)=q\prod_{i=0}^n(x-x_i)^{k_i}
次数至少是 m+1m+1,又它是 mm次多项式,则p=0p=0,证毕。
参考
1.《Numerical Analysis》——David Kincaid&Ward Cheney
多项式插值中的一些定理证明相关推荐
- 【Paper】论文中定义、定理、引理、证明分别的含义
论文中定义.定理.引理.证明分别的含义 1. 定义 Definition 2. 定理 Theorem 3. 证明 Proof 4. 引理 Lemma 5. 谈论/评论 Remark 6. 注解/说明 ...
- Latex中定义、定理、引理、证明 设置方法总结
Latex中定义.定理.引理.证明 设置方法总结 在LaTex中需要有关定理.公理.命题.引理.定义等时,常用如下命令 \newtheorem{定理环境名}{标题}[主计数器名] \newtheore ...
- 能否用计算机证明数学定理大全,那些教材中未给出证明的定理,到底能否使用?...
那些教材中未给出证明的定理,到底能否使用? 初中几何知识体系,是基于<几何原本>中的公理体系建立的,可以看作是公理体系的子集,在教材编写中,为了教学方便,并没有直接照搬,而是有选择性地&q ...
- latex怎么打区间_Latex中定义、定理、引理、证明 设置方法总结
Latex中定义.定理.引理.证明 设置方法总结 在LaTex中需要有关定理.公理.命题.引理.定义等时,常用如下命令 \newtheorem{定理环境名}{标题}[主计数器名] \newtheore ...
- 吉布斯现象与插值优化(上)Matlab实现多项式插值
前言 傅里叶变换是贯穿一系列课程的应用,从起初的不知其所以然到学习电路理论时的轮廓逐渐清晰. 在电路理论的结课复习阶段,我在图书馆找到国外教材翻阅,关注到行文中提到的吉布斯现象(Gibbs Pheno ...
- matlab全域基函数,多项式函数插值:全域多项式插值(一)单项式基插值、拉格朗日插值、牛顿插值 [MATLAB]...
全域多项式插值指的是在整个插值区域内形成一个多项式函数作为插值函数.关于多项式插值的基本知识,见"计算基本理论". 在单项式基插值和牛顿插值形成的表达式中,求该表达式在某一点处的值 ...
- 数值计算之 插值法(1)多项式插值——拉格朗日插值法
数值计算之 插值法(1)多项式插值--拉格朗日插值法 前言 什么是插值 多项式插值法 拉格朗日插值法 总结 前言 移动机器人有一个非常重要的任务,轨迹规划.轨迹规划需要满足运动学原理,即在路径规划给出 ...
- 数学建模准备 插值(拉格朗日多项式插值,牛顿多项式插值,分段线性插值,分段三次样条插值,分段三次Hermite插值)
文章目录 摘要(必看) 0 基础概念 什么是插值 插值用途 什么是拟合 插值和拟合的相同点 插值和拟合的不同点 1 常用的基本插值方法 1.1 多项式插值法 1.1.1 拉格朗日多项式插值法 多项式插 ...
- 插值与拟合 (一) : 拉格朗日多项式插值 、Newton插值 、分段线性插值、Hermite插值 、样条插值、 B 样条函数插值、二维插值
插值:求过已知有限个数据点的近似函数. 拟合:已知有限个数据点,求近似函数,不要求过已知数据点,只要求在某种意义下它在这些点上的总偏差最小. 插值和拟合都是要根据一组数据构造一个函数作为近似,由于近似 ...
最新文章
- Mapped Statements collection does not contain value for TaskMapper.selectByPrimaryKey
- Linux系统下不同机器之间拷贝文件的方法
- 功能与命令式编程。 Java 8中的斐波那契,素数和阶乘
- C语言extern用法
- hdu 4775 Infinite Go(暴力)
- springboot踩过的坑
- Concrete Mathematics A Foundation for Computer Science
- win10计算机怎么拨号上网,win10拨号连接怎么创建 win10宽带拨号连接如何设置
- 脱口秀在尝试处理一件难事
- 第004讲 浮动窗口 表单及表单控件
- java 毫秒转分钟和秒_毫秒转换为天、小时、分、秒
- 微型计算机主板usb电源损坏,自已动手彻底解决主板USB供电不足的问题
- background的用法
- 【bzoj1150】【CSTC2007】【数据备份】【贪心】
- 03弱引用内存泄露和hash冲突-ThreadLocal详解-并发编程(Java)
- 【堡垒机】2022年云堡垒机品牌排名大比拼
- 不同vlan间的通信--三层交换技术
- 在matlab中find是什么意思,Matlab中find函数的使用说明
- 千寻位置 开发demo_千寻位置发布全新操作系统,助力卫星导航产业化
- STM32学习笔记1-软件安装
热门文章
- Google Cloud 发起“Data Cloud Alliance”新联盟
- Android系统设置选项的包名
- Java---设计模式app小软件汇总应用
- 打印机服务器端口无显示器,什么是打印机端口?无法找到打印机端口是怎么回事?...
- win7 找不到 本地用户和组
- bezier.CSS_SVG_canvas画_贝兹曲线
- 利用fiddler和低版本的iTunes实现iOS抓包
- 面试宝典(一) - 让你不再错过“金九银十“的求职浪潮之简历包装篇
- python 坐标轴单位标注_matlab坐标轴如何带单位标注?
- 通过网页控制嵌入式设备