当前位置: 首页>編程日記>正文

数值计算之 插值法(3)多项式插值法的解,范德蒙矩阵,龙格现象

数值计算之 插值法(3)多项式插值法的解,范德蒙矩阵,龙格现象

数值计算之 插值法(3)多项式插值法的解,范德蒙矩阵,龙格现象

  • 前言
  • 多项式插值法的解与范德蒙矩阵
  • 龙格现象

前言

上两篇分别是拉格朗日插值法和牛顿插值法。

拉格朗日插值的思想是将一个多项式拆分为多个多项式,每个多项式完成一个节点插值。

牛顿插值法的思想是一种递推法,当出现新的抽样点时,增量更新一次插值函数。

从形式上来看,两种插值法得到的多项式不同,牛顿插值因为使用增量更新,因此计算量相对较小。

多项式插值法的解与范德蒙矩阵

未知表达式的函数f(x)f(x)f(x)满足y0=f(x0),y1=f(x1),y2=f(x2),…,yn=f(xn)y_0=f(x_0),y_1=f(x_1),y_2=f(x_2),\dots,y_n=f(x_n)y0=f(x0),y1=f(x1),y2=f(x2),,yn=f(xn),插值多项式的形式表示为:
P(x)=a0+a1x+a2x2+⋯+anxnP(x)=a_0+a_1x+a_2x^2+\dots+a_nx^n P(x)=a0+a1x+a2x2++anxn

如果直接使用方程组求解,则有以下方程组:
P(x0)=a0+a1x0+a2x02+⋯+anx0n=y0P(x1)=a0+a1x1+a2x12+⋯+anx1n=y1P(x2)=a0+a1x2+a2x22+⋯+anx2n=y2…P(xn)=a0+a1xn+a2xn2+⋯+anxnn=ynP(x_0)=a_0+a_1x_0+a_2x_0^2+\dots+a_nx_0^n=y_0 \\ P(x_1)=a_0+a_1x_1+a_2x_1^2+\dots+a_nx_1^n=y_1 \\ P(x_2)=a_0+a_1x_2+a_2x_2^2+\dots+a_nx_2^n=y_2 \\ \dots \\ P(x_n)=a_0+a_1x_n+a_2x_n^2+\dots+a_nx_n^n=y_n \\ P(x0)=a0+a1x0+a2x02++anx0n=y0P(x1)=a0+a1x1+a2x12++anx1n=y1P(x2)=a0+a1x2+a2x22++anx2n=y2P(xn)=a0+a1xn+a2xn2++anxnn=yn
用矩阵方程表示:
Xα=γ[1x0x02…x0n1x1x12…x1n1x2x22…x2n……………1xnxn2…xnn][a0a1a2…an]=[y0y1y2…yn]X\alpha = \gamma \\ \quad \\ \begin{bmatrix} 1 & x_0 & x_0^2 & \dots & x_0^n \\ 1 & x_1 & x_1^2 & \dots & x_1^n \\ 1 & x_2 & x_2^2 & \dots & x_2^n \\ \dots & \dots & \dots & \dots & \dots \\ 1 & x_n & x_n^2 & \dots & x_n^n \\ \end{bmatrix}\begin{bmatrix} a_0 \\ a_1 \\ a_2 \\ \dots \\ a_n \\ \end{bmatrix} =\begin{bmatrix} y_0 \\ y_1 \\ y_2 \\ \dots \\ y_n \\ \end{bmatrix} Xα=γ1111x0x1x2xnx02x12x22xn2x0nx1nx2nxnna0a1a2an=y0y1y2yn
其中,矩阵XXX是范德蒙矩阵,其行列式∣X∣=∏0≤j<i≤n(xi−xj)|X|=\prod_{0\le j<i\le n}(x_i-x_j)X=0j<in(xixj)。由于采样点横坐标都不相同,因此∣X∣≠0|X|\ne 0X=0XXX可逆,则α\alphaα必然是有唯一解的。这说明,同阶多项式插值的解实际上是唯一的,拉格朗日插值和牛顿插值的结果只是形式表达上不相同。

龙格现象

在我们的直觉上,使用多项式插值法时,多项式阶数越高,插值效果越好。但是,当多项式的阶数较大时,可能会出现龙格现象即在插值区间的边缘,插值结果与真实函数差值巨大。这个现象表明在使用多项式插值时,尽量避免使用阶数过高的多项式。
在这里插入图片描述
除了不使用高阶多项式插值外,也可以通过插值节点的选择来避免龙格现象。当均匀取节点时更有可能出现龙格现象,而选择切比雪夫或者高斯节点就可能避免龙格现象。


https://www.fengoutiyan.com/post/14587.html

相关文章:

  • 插值多项式的次数
  • 三次lagrange插值多项式
  • 二次lagrange插值多项式
  • lagrange插值多项式实验报告
  • 龙格函数多项式插值Lagrange
  • 埃尔米特插值多项式
  • 一次插值多项式
  • lagrange插值多项式例题
  • 鏡像模式如何設置在哪,圖片鏡像操作
  • 什么軟件可以把圖片鏡像翻轉,C#圖片處理 解決左右鏡像相反(旋轉圖片)
  • 手機照片鏡像翻轉,C#圖像鏡像
  • 視頻鏡像翻轉軟件,python圖片鏡像翻轉_python中鏡像實現方法
  • 什么軟件可以把圖片鏡像翻轉,利用PS實現圖片的鏡像處理
  • 照片鏡像翻轉app,java實現圖片鏡像翻轉
  • 什么軟件可以把圖片鏡像翻轉,python圖片鏡像翻轉_python圖像處理之鏡像實現方法
  • matlab下載,matlab如何鏡像處理圖片,matlab實現圖像鏡像
  • 圖片鏡像翻轉,MATLAB:鏡像圖片
  • 鏡像翻轉圖片的軟件,圖像處理:實現圖片鏡像(基于python)
  • canvas可畫,JavaScript - canvas - 鏡像圖片
  • 圖片鏡像翻轉,UGUI優化:使用鏡像圖片
  • Codeforces,CodeForces 1253C
  • MySQL下載安裝,Mysql ERROR: 1253 解決方法
  • 勝利大逃亡英雄逃亡方案,HDU - 1253 勝利大逃亡 BFS
  • 大一c語言期末考試試題及答案匯總,電大計算機C語言1253,1253《C語言程序設計》電大期末精彩試題及其問題詳解
  • lu求解線性方程組,P1253 [yLOI2018] 扶蘇的問題 (線段樹)
  • c語言程序設計基礎題庫,1253號C語言程序設計試題,2016年1月試卷號1253C語言程序設計A.pdf
  • 信奧賽一本通官網,【信奧賽一本通】1253:抓住那頭牛(詳細代碼)
  • c語言程序設計1253,1253c語言程序設計a(2010年1月)
  • 勝利大逃亡英雄逃亡方案,BFS——1253 勝利大逃亡
  • 直流電壓測量模塊,IM1253B交直流電能計量模塊(艾銳達光電)
  • c語言程序設計第三版課后答案,【渝粵題庫】國家開放大學2021春1253C語言程序設計答案
  • 18轉換為二進制,1253. 將數字轉換為16進制
  • light-emitting diode,LightOJ-1253 Misere Nim
  • masterroyale魔改版,1253 Dungeon Master
  • codeformer官網中文版,codeforces.1253 B
  • c語言程序設計考研真題及答案,2020C語言程序設計1253,1253計算機科學與技術專業C語言程序設計A科目2020年09月國家開 放大學(中央廣播電視大學)
  • c語言程序設計基礎題庫,1253本科2016c語言程序設計試題,1253電大《C語言程序設計A》試題和答案200901
  • 肇事逃逸車輛無法聯系到車主怎么辦,1253尋找肇事司機