数值分析复习:Newton插值

文章目录

  • 牛顿(Newton)插值
    • 插值条件
    • 基函数
    • 插值多项式
    • 差商
      • 差商的基本性质
      • 差商估计
      • 差商的Leibniz公式
    • 余项估计

本篇文章适合个人复习翻阅,不建议新手入门使用

牛顿(Newton)插值

插值条件

n+1个插值节点 x 0 , x 1 , … , x n x_0,x_1,\dots,x_n x0,x1,,xn 处函数值相同

基函数

{ ω i ( x ) } i = 0 n \{\omega_i(x)\}_{i=0}^n {ωi(x)}i=0n,其中 ω i ( x ) = ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x i − 1 ) \omega_i(x)=(x-x_0)(x-x_1)\cdots(x-x_{i-1}) ωi(x)=(xx0)(xx1)(xxi1)称之为节点多项式

插值多项式

∏ n f ( x ) = ∑ i = 0 n f [ x 0 , x 1 , … , x k ] ω k ( x ) \prod_nf(x)=\sum\limits_{i=0}^nf[x_0,x_1,\dots,x_k]\omega_k(x) nf(x)=i=0nf[x0,x1,,xk]ωk(x)其中 f [ x 0 , x 1 , … , x k ] f[x_0,x_1,\dots,x_k] f[x0,x1,,xk]称为 f f f 关于点 x 0 , x 1 , … , x k x_0,x_1,\dots,x_k x0,x1,,xk的k阶牛顿差商

差商

差商的基本性质

  • n阶差商为n次插值多项式的首项系数
  • 差商值与节点排列顺序无关
  • f [ x 0 , x 1 , … , x n ] = f ( x n ) − ∏ n − 1 f ( x n ) ω n ( x n ) = ∏ n f ( x n ) − ∏ n − 1 f ( x n ) ω n ( x n ) = ∑ i = 0 n f ( x i ) ω n + 1 ′ ( x i ) \begin{split} f[x_0,x_1,\dots,x_n]&=\frac{f(x_n)-\prod_{n-1}f(x_n)}{\omega_n(x_n)}\\ &=\frac{\prod_nf(x_n)-\prod_{n-1}f(x_n)}{\omega_n(x_n)}\\ &=\sum\limits_{i=0}^n\frac{f(x_i)}{\omega'_{n+1}(x_i)}\\ \end{split} f[x0,x1,,xn]=ωn(xn)f(xn)n1f(xn)=ωn(xn)nf(xn)n1f(xn)=i=0nωn+1(xi)f(xi)
  • f [ x 0 , x 1 , … , x n ] = f [ x 0 , … , x n − 1 ] − f [ x 1 , … , x n ] x 0 − x n f[x_0,x_1,\dots,x_n]=\frac{f[x_0,\dots,x_{n-1}]-f[x_1,\dots,x_n]}{x_0-x_n} f[x0,x1,,xn]=x0xnf[x0,,xn1]f[x1,,xn]

证明思路:

第二条性质:
前两个等号容易得到;第三个等号:只需注意到

  • n阶差商是n次Newton插值的首项系数
  • 等式右端是Lagrange插值多项式的首项系数
  • Newton插值、Lagrange插值是同一插值多项式的不同表达
  • 多项式插值的唯一性(由Vandermonde行列式的性质易证)

第三条性质:归纳法可证

差商估计

f [ x 0 , x 1 , … , x n ] = f ( m ) ( ξ ) m ! f[x_0,x_1,\dots,x_n]=\frac{f^{(m)}(\xi)}{m!} f[x0,x1,,xn]=m!f(m)(ξ)其中 ξ ∈ ( min ⁡ { x i } , max ⁡ { x i } ) \xi\in(\min\{x_i\},\max\{x_i\}) ξ(min{xi},max{xi})

证明思路:构造辅助函数 f ( x ) − ∏ n f ( x ) f(x)-\prod_nf(x) f(x)nf(x),使用 n n n次Rolle中值定理

差商的Leibniz公式

f ( x ) = ϕ ( x ) ψ ( x ) f(x)=\phi(x)\psi(x) f(x)=ϕ(x)ψ(x),则
f [ x 0 , x 1 , … , x n ] = ∑ i = 0 n ϕ ( x 0 , … , x i ) ψ ( x i , … , x n ) f[x_0,x_1,\dots,x_n]=\sum\limits_{i=0}^n\phi(x_0,\dots,x_i)\psi(x_i,\dots,x_n) f[x0,x1,,xn]=i=0nϕ(x0,,xi)ψ(xi,,xn)

证明思路:对 f , ϕ , ψ f,\phi,\psi f,ϕ,ψ 分别进行Newton插值即可

余项估计

R n ( x ) = f ( x ) − ∏ n f ( x ) = f [ x 0 , x 1 , … , x n , x ] ∏ i = 0 n ( x − x i ) = f [ x 0 , x 1 , … , x n , x ] ω n + 1 ( x ) \begin{split} R_n(x)&=f(x)-\prod_nf(x)\\ &=f[x_0,x_1,\dots,x_n,x]\prod\limits_{i=0}^n(x-x_i)\\ &=f[x_0,x_1,\dots,x_n,x]\omega_{n+1}(x)\\ \end{split} Rn(x)=f(x)nf(x)=f[x0,x1,,xn,x]i=0n(xxi)=f[x0,x1,,xn,x]ωn+1(x)

参考书籍:《数值分析》李庆扬 王能超 易大义 编

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/751872.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

使用 pnpm 搭建 monorepo 项目

引言 在我之前的开发经历中,并没有实际使用过 Monorepo 管理项目,尽管之前对此有所了解,但并未深入探究。然而,如今许多开源项目都采纳了 Monorepo 方式,对于不熟悉它的开发者来说,阅读和理解这些项目的源…

【HarmonyOS】ArkUI - 向左/向右滑动删除

核心知识点:List容器 -> ListItem -> swipeAction 先看效果图: 代码实现: // 任务类 class Task {static id: number 1// 任务名称name: string 任务${Task.id}// 任务状态finished: boolean false }// 统一的卡片样式 Styles func…

C语言快速入门之内存函数的使用和模拟实现

1.memcpy 它可以理解为memory copy的组合,memory有记忆的意思,这里指的是内存,copy是拷贝,这个函数是针对内存块进行拷贝的 函数原型 void* memcpy(void* destination,const void* source, size_t num); 从source位置开始&am…

ChatGPT国内镜像站大全

#今天在知乎看到一个问题:“平民不参与内测的话没有账号还有机会使用ChatGPT吗?” 从去年GPT大火到现在,关于GPT的消息铺天盖地,真要有心想要去用,途径很多,别的不说,国内GPT的镜像站到处都是&…

基于sortablejs实现拖拽element-ui el-table表格行进行排序

可以用原生的dragstart、drag、dragend、dragover、drop、dragleave实现这个效果&#xff0c;但是有现成的轮子就不要重复造了&#xff0c;看效果&#xff1a; <template><el-table :class"$options.name" :data"tableData" ref"table"…

Docker进阶教程 - 1 Dockerfile

更好的阅读体验&#xff1a;点这里 &#xff08; www.doubibiji.com &#xff09; 1 Dockerfile Dockerfile 是做什么的&#xff1f; 我们前面说到&#xff0c;制作镜像的方法主要有两种方式&#xff1a; 使用 docker commit 命令&#xff1b;使用 Dockerfile 文件。 但是…

leetcode每日一题310.最小高度树

目录 一.题目原型 二.题目思路 三.代码实现 一.题目原型 二.题目思路 首先&#xff0c;我们看了样例&#xff0c;发现这个树并不是二叉树&#xff0c;是多叉树。 然后&#xff0c;我们可能想到的解法是&#xff1a;根据题目的意思&#xff0c;就挨个节点遍历bfs&#xff0c;…

瑞_Redis_短信登录_Redis代替session的业务流程

文章目录 项目介绍1 短信登录1.1 项目准备1.2 基于Session实现登录流程1.3 Redis代替session的业务流程1.3.1 设计key的结构1.3.2 设计key的具体细节1.3.3 整体访问流程1.3.4 代码实现 &#x1f64a; 前言&#xff1a;本文章为瑞_系列专栏之《Redis》的实战篇的短信登录章节的R…

springboot项目读取excel表格内容到数据库,excel表格字段为整数的读取方法

在我昨天的项目中&#xff0c;我需要把excel表格中字段为整数的字段读取到数据库中进行保存&#xff0c;但是在内置方法中并没有读取整数的方法&#xff08;也有可能是我没发现&#xff0c;太菜了~~&#xff09;&#xff0c;那接下来我就提供给大家一个简单地方法来读取excel表…

Apache-Doris基础概念

OLAP数据库Doris 一、Doris架构二、基本概念1. Row & Column2. Partition & Tablet3. 建表示例&#xff08;1&#xff09;列的定义&#xff08;2&#xff09;分区分桶&#xff08;3&#xff09;多列分区&#xff08;4&#xff09;PROPERTIES&#xff08;5&#xff09;E…

【LabVIEW FPGA入门】单周期定时循环

单周期定时循环详解 单周期定时环路是FPGA编程中最强大的结构之一。单周期定时循环中的代码更加优化&#xff0c;在FPGA上占用更少的空间&#xff0c;并且比标准While循环中的相同代码执行得更快。单周期定时环路将使能链从环路中移除&#xff0c;以节省FPGA上的空间。…

windows下使用tree指定层数生成项目结构

windows自带的tree太辣鸡了&#xff0c;我们需要找东西代替 工具链接 Tree for Windows (sourceforge.net) 点击这里下载 置入Git 解压下载的压缩包&#xff0c;将bin目录下的exe复制下来 进入你的Git目录&#xff0c;将其放入Git目录下的usr/bin目录下 打开Git Bash 输入…

如何在Mac中删除照片?这里有详细步骤

前言 本文介绍如何从Mac中删除照片,以释放硬盘空间或更好地组织文件和文件夹。 如何使用废纸篓删除Mac上的图片 在Mac上删除图片的最简单方法之一是使用废纸篓功能。学习只需几秒钟。下面是如何删除单个图片以及如何在Mac上删除多个图片,以及一些关键和有用的提示,以使该…

Selenium-webdriver_manager判断是否已经下载过驱动(复用缓存驱动)

1,谷歌浏览器默认位置 2&#xff0c;ChromeDriverManager 下载的驱动位置 其中admin为机器的用户名 def installDriver(self):"""判断是否需要下载driver""""""找到本机谷歌浏览器版本""""""C:\P…

【学习心得】Python数据分析的基本思路

一、什么是数据分析&#xff1f; 数据分析是指通过一些方法&#xff0c;对一些数据进行分析&#xff0c;从中提取出有价值的信息并形成结论进行展示。 &#xff08;1&#xff09;一些方法 数学和统计学方法&#xff1a;例如回归分析、聚类分析、主成分分析、时间序列分析等&a…

【代码随想录算法训练营第二十五天 | LeetCode216.组合总和III、17.电话号码的字母组合】

代码随想录算法训练营第二十五天 | LeetCode216.组合总和III、17.电话号码的字母组合 一、77. 组合 解题代码C&#xff1a; class Solution { private:vector<vector<int>> result; // 存放结果集vector<int> path; // 符合条件的结果void backtracking(in…

Shell脚本——实用工具类

删除用户脚本.sh #!/bin/bash # # Delete_User - Automates the 4 steps to remove an account # ################################################################# # # Define Functions # ################################################################# function…

机器人路径规划:基于改进型A*算法的机器人路径规划(提供Python代码)

一、A*算法介绍 A*算法最早可追溯到1968年&#xff0c;在IEEE Transactions on Systems Science and Cybernetics中的论文A Formal Basis for the Heuristic Determination of Minimum Cost Paths中首次提出。 https://blog.csdn.net/weixin_46204734/article/details/1367905…

C++/CLI学习笔记5(快速打通c++与c#相互调用的桥梁)

3.6 句柄与指针 在标准C中&#xff0c;指针容纳的是另一个变量或函数的内存地址。这意味着可用指针间接引用变量。 但在C/CLI中&#xff0c;是“运行时”帮你管理内存&#xff0c;所以它保留了将内存里的东西移来移去以最大化可用内存的权利。这意味着对象可能不老是呆在一个位…

Windows 网络质量测试

Windows 网络质量测试 References 保证网络稳定性&#xff0c;建议最大数据包延迟 200ms 以内&#xff0c;数据包最大和最小延迟差 100ms 以内&#xff0c;丢包率最好不丢包或 5% 以内。 ping www.baidu.com -t 调出 运行 (快捷键 Win R)&#xff0c;输入 cmd&#xff0c;pi…