当前位置: 首页 > news >正文

牛顿均差知识

牛顿均差


牛顿均差公式

f ( x ) = ∑ i = 0 n [ f [ x 0 , ⋯ , x i ] ⋅ ∏ k = 0 i − 1 ( x − x k ) ] + f [ x 0 , ⋯ , x n , x ] ⋅ ∏ k = 0 n ( x − x k ) f(x) = \sum_{i=0}^{n} \left[ f[x_0, \cdots, x_i] \cdot \prod_{k=0}^{i-1} (x - x_k) \right] + f[x_0, \cdots, x_n, x] \cdot \prod_{k=0}^{n} (x - x_k) f(x)=i=0n[f[x0,,xi]k=0i1(xxk)]+f[x0,,xn,x]k=0n(xxk)

性质证明

均差( D i v i d e d D i f f e r e n c e s \mathcal{Divided\quad Differences} DividedDifferences)的定义为:
f [ x i , x i + 1 , … , x i + k ] = f [ x i + 1 , … , x i + k ] − f [ x i , … , x i + k − 1 ] x i + k − x i f[x_i, x_{i+1}, \ldots, x_{i+k}] = \frac{f[x_{i+1}, \ldots, x_{i+k}] - f[x_i, \ldots, x_{i+k-1}]}{x_{i+k} - x_i} f[xi,xi+1,,xi+k]=xi+kxif[xi+1,,xi+k]f[xi,,xi+k1]

性质1:均差与导数的关系

若函数 f ( x ) f(x) f(x) 在包含节点 x 0 , x 1 , … , x n x_0, x_1, \ldots, x_n x0,x1,,xn 的区间内存在 n n n 阶导数,则存在一点 ξ \xi ξ 使得:
f [ x 0 , x 1 , … , x n ] = f ( n ) ( ξ ) n ! , ξ ∈ 区间 ( x 0 , x 1 , … , x n ) f[x_0, x_1, \ldots, x_n] = \frac{f^{(n)}(\xi)}{n!}, \quad \xi \in \text{区间}(x_0, x_1, \ldots, x_n) f[x0,x1,,xn]=n!f(n)(ξ),ξ区间(x0,x1,,xn)

证明

通过牛顿插值多项式的误差项推导。已知 n n n 次插值多项式 N n ( x ) N_n(x) Nn(x) 的误差为:
f ( x ) − N n ( x ) = f [ x 0 , x 1 , … , x n , x ] ⋅ ∏ i = 0 n ( x − x i ) f(x) - N_n(x) = f[x_0, x_1, \ldots, x_n, x] \cdot \prod_{i=0}^{n} (x - x_i) f(x)Nn(x)=f[x0,x1,,xn,x]i=0n(xxi)
x = x n + 1 x = x_{n+1} x=xn+1 时( x n + 1 x_{n+1} xn+1 为区间内另一点):
f ( x n + 1 ) − N n ( x n + 1 ) = f [ x 0 , x 1 , … , x n , x n + 1 ] ⋅ ∏ i = 0 n ( x n + 1 − x i ) f(x_{n+1}) - N_n(x_{n+1}) = f[x_0, x_1, \ldots, x_n, x_{n+1}] \cdot \prod_{i=0}^{n} (x_{n+1} - x_i) f(xn+1)Nn(xn+1)=f[x0,x1,,xn,xn+1]i=0n(xn+1xi)
根据罗尔定理,存在 ξ \xi ξ 使得:
f [ x 0 , x 1 , … , x n , x n + 1 ] = f ( n + 1 ) ( ξ ) ( n + 1 ) ! f[x_0, x_1, \ldots, x_n, x_{n+1}] = \frac{f^{(n+1)}(\xi)}{(n+1)!} f[x0,x1,,xn,xn+1]=(n+1)!f(n+1)(ξ)
特别地,当 x n + 1 → x n x_{n+1} \to x_n xn+1xn 时,均差退化为导数形式。

性质2:均差的对称性

均差的值与节点的排列顺序无关,即:
f [ x 0 , x 1 , … , x n ] = f [ x π ( 0 ) , x π ( 1 ) , … , x π ( n ) ] f[x_0, x_1, \ldots, x_n] = f[x_{\pi(0)}, x_{\pi(1)}, \ldots, x_{\pi(n)}] f[x0,x1,,xn]=f[xπ(0),xπ(1),,xπ(n)]
其中 π \pi π 是任意排列函数。

证明

通过数学归纳法:

  • 基础情况:

n = 1 n=1 n=1 时, f [ x 0 , x 1 ] = f ( x 1 ) − f ( x 0 ) x 1 − x 0 f[x_0, x_1] = \frac{f(x_1) - f(x_0)}{x_1 - x_0} f[x0,x1]=x1x0f(x1)f(x0),显然对称。

  • 归纳假设:

假设对于 n = k n=k n=k 时成立。

  • 归纳步骤:对于 n = k + 1 n=k+1 n=k+1,利用均差递推公式:
    f [ x 0 , x 1 , … , x k + 1 ] = f [ x 1 , … , x k + 1 ] − f [ x 0 , … , x k ] x k + 1 − x 0 f[x_0, x_1, \ldots, x_{k+1}] = \frac{f[x_1, \ldots, x_{k+1}] - f[x_0, \ldots, x_k]}{x_{k+1} - x_0} f[x0,x1,,xk+1]=xk+1x0f[x1,,xk+1]f[x0,,xk]
    根据归纳假设,分子中的两项均对称,因此整体对称。

性质3:均差的线性性

f ( x ) = a ⋅ g ( x ) + b ⋅ h ( x ) f(x) = a \cdot g(x) + b \cdot h(x) f(x)=ag(x)+bh(x),则:
f [ x 0 , x 1 , … , x n ] = a ⋅ g [ x 0 , x 1 , … , x n ] + b ⋅ h [ x 0 , x 1 , … , x n ] f[x_0, x_1, \ldots, x_n] = a \cdot g[x_0, x_1, \ldots, x_n] + b \cdot h[x_0, x_1, \ldots, x_n] f[x0,x1,,xn]=ag[x0,x1,,xn]+bh[x0,x1,,xn]

证明

通过均差的定义直接推导:
f [ x 0 , x 1 ] = f ( x 1 ) − f ( x 0 ) x 1 − x 0 = a ⋅ g ( x 1 ) + b ⋅ h ( x 1 ) − a ⋅ g ( x 0 ) − b ⋅ h ( x 0 ) x 1 − x 0 = a ⋅ g ( x 1 ) − g ( x 0 ) x 1 − x 0 + b ⋅ h ( x 1 ) − h ( x 0 ) x 1 − x 0 = a ⋅ g [ x 0 , x 1 ] + b ⋅ h [ x 0 , x 1 ] \begin{aligned} f[x_0, x_1] &= \frac{f(x_1) - f(x_0)}{x_1 - x_0} \\ &= \frac{a \cdot g(x_1) + b \cdot h(x_1) - a \cdot g(x_0) - b \cdot h(x_0)}{x_1 - x_0} \\ &= a \cdot \frac{g(x_1) - g(x_0)}{x_1 - x_0} + b \cdot \frac{h(x_1) - h(x_0)}{x_1 - x_0} \\ &= a \cdot g[x_0, x_1] + b \cdot h[x_0, x_1] \end{aligned} f[x0,x1]=x1x0f(x1)f(x0)=x1x0ag(x1)+bh(x1)ag(x0)bh(x0)=ax1x0g(x1)g(x0)+bx1x0h(x1)h(x0)=ag[x0,x1]+bh[x0,x1]
对高阶均差同理可证。
\end{proof}

性质4:均差与多项式的关系

f ( x ) f(x) f(x) m m m 次多项式,则:

  • n ≤ m n \leq m nm 时, n n n 阶均差 f [ x 0 , x 1 , … , x n ] f[x_0, x_1, \ldots, x_n] f[x0,x1,,xn] m − n m-n mn 次多项式;
  • n > m n > m n>m 时, n n n 阶均差 f [ x 0 , x 1 , … , x n ] = 0 f[x_0, x_1, \ldots, x_n] = 0 f[x0,x1,,xn]=0

证明

通过导数性质:

  • m m m 次多项式的 m m m 阶导数为常数, m + 1 m+1 m+1 阶导数为零。
  • 根据均差与导数的关系:
    f [ x 0 , x 1 , … , x n ] = f ( n ) ( ξ ) n ! f[x_0, x_1, \ldots, x_n] = \frac{f^{(n)}(\xi)}{n!} f[x0,x1,,xn]=n!f(n)(ξ)
    n > m n > m n>m 时, f ( n ) ( ξ ) = 0 f^{(n)}(\xi) = 0 f(n)(ξ)=0,故均差为零。

性质5:均差的递推性质

均差满足递推公式:
f [ x 0 , x 1 , … , x n ] = f [ x 1 , … , x n ] − f [ x 0 , … , x n − 1 ] x n − x 0 f[x_0, x_1, \ldots, x_n] = \frac{f[x_1, \ldots, x_n] - f[x_0, \ldots, x_{n-1}]}{x_n - x_0} f[x0,x1,,xn]=xnx0f[x1,,xn]f[x0,,xn1]

证明
通过牛顿插值多项式的构造:

  • n n n 次插值多项式 N n ( x ) N_n(x) Nn(x) 可表示为:
    N n ( x ) = N n − 1 ( x ) + f [ x 0 , x 1 , … , x n ] ⋅ ∏ i = 0 n − 1 ( x − x i ) N_n(x) = N_{n-1}(x) + f[x_0, x_1, \ldots, x_n] \cdot \prod_{i=0}^{n-1} (x - x_i) Nn(x)=Nn1(x)+f[x0,x1,,xn]i=0n1(xxi)
  • 代入 x = x n x = x_n x=xn 并整理,可得递推公式。

差分与位移算子


例1

Δ ( f k g k ) = f k Δ g k + g k + 1 Δ f k \Delta(f_kg_k)=f_k\Delta g_k + g_{k+1}\Delta f_k Δ(fkgk)=fkΔgk+gk+1Δfk

证明
Δ ( f k g k ) = f k + 1 g k + 1 − f k g k = f k + 1 ( g k + 1 − g k ) + g k + 1 ( f k + 1 − f k ) = f k + 1 Δ g k + g k + 1 Δ f k \begin{align*} \Delta(f_kg_k) &= f_{k+1}g_{k+1} - f_kg_k \\ &= f_{k+1}(g_{k+1} - g_k) + g_{k+1}(f_{k+1} - f_k) \\ &= f_{k+1}\Delta g_k + g_{k+1}\Delta f_k \end{align*} Δ(fkgk)=fk+1gk+1fkgk=fk+1(gk+1gk)+gk+1(fk+1fk)=fk+1Δgk+gk+1Δfk

例2

∑ k = 0 n − 1 f k Δ g k = f n g n − f 0 g 0 − ∑ k = 0 n − 1 g k + 1 Δ f k \sum_{k=0}^{n-1}f_{k}\Delta g_{k}=f_ng_n-f_0g_0-\sum_{k=0}^{n-1}g_{k+1}\Delta f_{k} k=0n1fkΔgk=fngnf0g0k=0n1gk+1Δfk

这个是著名的阿贝尔变换

∑ k = m n f k Δ g k = f n + 1 g n + 1 − f m g m − ∑ k = m n − 1 g k + 1 Δ f k \sum_{k=m}^{n}f_{k}\Delta g_{k}=f_{n+1}g_{n+1}-f_mg_m-\sum_{k=m}^{n-1}g_{k+1}\Delta f_{k} k=mnfkΔgk=fn+1gn+1fmgmk=mn1gk+1Δfk

证明

∑ k = m n f k Δ g k = ∑ k = m n f k g k + 1 − ∑ k = m n f k g k = ∑ k = m n f k g k + 1 − ∑ k = m + 1 n + 1 f k g k + f n + 1 g n + 1 − f m g m = f n + 1 g n + 1 − f m g m − ∑ k = m n − 1 g k + 1 ( f k + 1 − f k ) ) = f n + 1 g n + 1 − f m g m − ∑ k = m n − 1 g k + 1 Δ f k \begin{align*} \sum_{k=m}^{n}f_{k}\Delta g_{k} &=\sum_{k=m}^{n}f_kg_{k+1}-\sum_{k=m}^{n}f_kg_k \\ &=\sum_{k=m}^{n}f_kg_{k+1}-\sum_{k=m+1}^{n+1}f_kg_k+f_{n+1}g_{n+1}-f_{m}g_{m} \\ &=f_{n+1}g_{n+1}-f_mg_m-\sum_{k=m}^{n-1}g_{k+1}(f_{k+1}-f_{k})) \\ &=f_{n+1}g_{n+1}-f_mg_m-\sum_{k=m}^{n-1}g_{k+1}\Delta f_{k} \end{align*} k=mnfkΔgk=k=mnfkgk+1k=mnfkgk=k=mnfkgk+1k=m+1n+1fkgk+fn+1gn+1fmgm=fn+1gn+1fmgmk=mn1gk+1(fk+1fk))=fn+1gn+1fmgmk=mn1gk+1Δfk

赋值: m = 0 , n = n − 1 m=0,n=n-1 m=0,n=n1即可得到:

∑ k = m n f k Δ g k = f n + 1 g n + 1 − f m g m − ∑ k = m n − 1 g k + 1 Δ f k \sum_{k=m}^{n}f_{k}\Delta g_{k}=f_{n+1}g_{n+1}-f_mg_m-\sum_{k=m}^{n-1}g_{k+1}\Delta f_{k} k=mnfkΔgk=fn+1gn+1fmgmk=mn1gk+1Δfk

相关文章:

  • 负进制转换
  • STC32G12K128实战:串口通信
  • 计算机网络-MPLS LDP基础实验配置
  • 易学探索助手-个人记录(十)
  • Kuka AI音乐AI音乐开发「人声伴奏分离」 —— 「Kuka Api系列|中文咬字清晰|AI音乐API」第6篇
  • python打卡day25
  • 数字高程模型(DEM)公开数据集介绍与下载指南
  • vscode extention踩坑记
  • 七部门:设立“国家创业投资引导基金”,优先支持取得关键核心技术突破的科技型企业上市融资
  • 学会使用ai作图
  • c++在头文件中声明全局的变量
  • TNNLS-2020《Autoencoder Constrained Clustering With Adaptive Neighbors》
  • Deep Learning(手写字识别 - CNN)
  • 阿里的库存秒杀实现与Inventory Hint技术解析
  • Python 接入DeepSeek
  • 高等数学第七章---微分方程(§7.1-§7.3微分方程概念、一阶微分方程、一阶微分线性方程)
  • 生成式人工智能认证(GAI认证)官网 - 全国统一认证中文服务平台上线
  • 一种资源有限单片机处理cJSON数据的方法
  • WordPress 文章和页面:它们的区别是什么?
  • Qt 的.pro文件配置详解
  • 新任美国驻华大使庞德伟抵京履职,外交部回应
  • 创同期历史新高!1至4月全国铁路发送旅客14.6亿人次
  • 赖清德为“临阵脱逃”作准备,国台办:绝不会任“台独”祸首逍遥法外
  • 联合国秘书长欢迎中美经贸高层会谈成果
  • 云南威信麟凤镇通报“有人穿‘警察’字样雨衣参与丧事”:已立案查处
  • 山西省委常委李金科添新职