数值分析——多项式插值之Lagrange插值

2023-11-20 16:10

本文主要是介绍数值分析——多项式插值之Lagrange插值,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、引言

  考虑这样一个实际例子,当我们按下计算器的正弦按钮时,会发生什么?我们都知道计算器有可以处理加法和乘法的硬件,但是,它是如何计算一个数的正弦值呢?多项式插值法就可以解决这样的问题。我们将在未来重新审视这个问题。目前,我们先来学什么是插值以及如何插值。

二、什么是插值

  如下图所示,假定我们收集了一组数据点$(x, y)$,譬如$(0, 1), (2, 2), (3, 4)$。有一条经过这三点的抛物线,我们把这条抛物线称为经过这3点的二次插值多项式

 

 

这样就引出了插值的数学定义,如下:

【插值的定义】  如果$P(x_i) = y_i (1 \leqslant i \leqslant n)$,那么函数$y = P(x)$插值了数据点

$(x_1, y_1), \cdot \cdot \cdot , (x_n, y_n)$

 简单来讲就是,如果一个函数通过了一组数据点,那么就称这个函数插值了这组数据点。

二、Lagrange插值

2.1 讨论

  现在我们知道了什么是插值,请大家考虑一个问题,如果我只知道一组n个数据点$(x_1, y_1), \cdot \cdot \cdot , (x_n, y_n)$,我们想要求出一个多项式,能够插值这一组所有的数据点。并且这个多项式的次数是$d = n - 1$次的,该怎么做?

  Lagrange插值公式给出了这个问题的解答方案。例如,假设给定点$(x_1, y_1), (x_2, y_2), (x_3, y_3)$,那么其2次插值多项式可以由Lagrange插值多项式给出,如下:

$P_2(x) = y_1 \frac{(x - x_2)(x - x_3)}{(x_1 - x_2)(x_1 - x_3)} + y_2 \frac{(x - x_1)(x - x_3)}{(x_2 - x_1)(x_2 - x_3)} + y_3 \frac{(x - x_1)(x - x_2)}{(x_3 - x_1)(x_3 - x_2)}$

   那么有人可能会问,这样的多项式一定是正确的吗,答案是:是的。我们可以验证一下:

1. 当$x = x1$时,$P_2(x1) = y1$;

2. 当$x = x2$时,$P_2(x2) = y2$;

3. 当$x = x3$时,$P_2(x3) = y3$;

   我们只考虑这三个点,因为我们只有这三个点,在这3个点上,这个多项式都成功的插值了,因此,这个多项式一定是正确的。(注意这个多项式关于变量$x$是2次的)

 

2.2 数学定义

一般地,假设给出$n$个点$(x_1, y_1), \cdot \cdot \cdot , (x_n, y_n)$,则对于1和n之间的每一个$k$可定义

$L_k(x) = \frac{(x - x_1) \cdot \cdot \cdot (x - x_{k -1})(x - x_{k + 1}) \cdot \cdot (x - x_n)}{(x_k - x_1) \cdot \cdot \cdot (x_k - x_{k -1})(x_k - x_{k + 1}) \cdot \cdot (x_k - x_n)}$

 $L_k$的一个有趣性质是:

1. $L_k(x_k) = 1$

2. $L_k(x_j) = 0 \qquad (j \neq k)$

  因此,定义$n - 1$次Lagrange多项式

$P_{n - 1}(x) = y_1 L_1(x) + \cdot \cdot \cdot + y_n L_n(x)$

2.3 存在性和及唯一性

  有人会问,对于给定的n个数据点,其插值多项式是唯一的吗?即只能是由一个多项式才能插值这n个点吗?答案是:不是

  大家想想就知道,对于二维平面的$n$个坐标点,我们肯定能画出无穷条线来穿过这些点,每一条线都对应这一个多项式。那么这个问题的意义何在?

  多项式是无穷的,但是,对于插值$n$个数据点的多项式,其最高次数是小于等于$n - 1$的,这样的多项式,只能是只有一个。用数学来描述这个问题如下:

【定理】

  设$(x_1, y_1), \cdot \cdot \cdot, (x_n, y_n)$是平面上$x_i$互不相同的$n$个点,那么存在一个而且仅存在一个次数小于等于$n - 1$次的多项式,满足

$P(x_i) = y_i,  \qquad i = 1,  \cdot \cdot \cdot, n$

【证明】

  (1) 存在性:存在性已由Lagrange插值的显式公式得出。

  (2) 唯一性:假定有存在两个这样公式,譬如$P(x)$及$Q(x)$,它们最多是$n - 1$次,而且都插值所有$n$个点,即有:

$P(x_1) = Q(x_1) = y_1, P(x_2) = Q(x_2) = y_2, \cdot \cdot \cdot, P(x_n) = Q(x_n) = y_n$。

  则有$H(x) = P(x) - Q(x)$,显然,$H$的次数最多也是$n - 1$,而且注意到

$H(x_1) = H(x_2) = \cdot \cdot \cdot = H(x_n)$

  即$H$有$n$个不同的零点。按照代数学基本定理,一个$d$次多项式,除了它恒等于零多项式,最多可能有$d$个零点。因此有

$H \equiv 0$

  于是,

$P(x) \equiv Q(x)$

因此,存在唯一的次数小于等于$n - 1$的多项式$P(x)$插值与$n$个点$(x_i, y_i)$。

 

2.4 范例

【题目】求插值于点$(0, 2), (1, 1), (2, 0), (3, -1)$的次数小于等于3的多项式。

【解】Lagrange形式如下:

$P(x) = 2 \frac{(x - 1)(x - 2)(x - 3)}{(0 - 1)(0 - 2)(0 - 3)} + 1 \frac{(x - 0)(x - 2)(x - 3)}{(1 - 0)(1 - 2)(1 - 3)} + 0 \frac{(x - 0)(x - 1)(x - 3)}{(2 - 0)(2 - 1)(2 - 3)} + (-1) \frac{(x - 0)(x - 1)(x - 2)}{(3 - 0)(3 - 1)(3 - 2)} = -x + 2$

转载于:https://www.cnblogs.com/Qling/p/9764941.html

这篇关于数值分析——多项式插值之Lagrange插值的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/395843

相关文章

性能分析之MySQL索引实战案例

文章目录 一、前言二、准备三、MySQL索引优化四、MySQL 索引知识回顾五、总结 一、前言 在上一讲性能工具之 JProfiler 简单登录案例分析实战中已经发现SQL没有建立索引问题,本文将一起从代码层去分析为什么没有建立索引? 开源ERP项目地址:https://gitee.com/jishenghua/JSH_ERP 二、准备 打开IDEA找到登录请求资源路径位置

SWAP作物生长模型安装教程、数据制备、敏感性分析、气候变化影响、R模型敏感性分析与贝叶斯优化、Fortran源代码分析、气候数据降尺度与变化影响分析

查看原文>>>全流程SWAP农业模型数据制备、敏感性分析及气候变化影响实践技术应用 SWAP模型是由荷兰瓦赫宁根大学开发的先进农作物模型,它综合考虑了土壤-水分-大气以及植被间的相互作用;是一种描述作物生长过程的一种机理性作物生长模型。它不但运用Richard方程,使其能够精确的模拟土壤中水分的运动,而且耦合了WOFOST作物模型使作物的生长描述更为科学。 本文让更多的科研人员和农业工作者

MOLE 2.5 分析分子通道和孔隙

软件介绍 生物大分子通道和孔隙在生物学中发挥着重要作用,例如在分子识别和酶底物特异性方面。 我们介绍了一种名为 MOLE 2.5 的高级软件工具,该工具旨在分析分子通道和孔隙。 与其他可用软件工具的基准测试表明,MOLE 2.5 相比更快、更强大、功能更丰富。作为一项新功能,MOLE 2.5 可以估算已识别通道的物理化学性质。 软件下载 https://pan.quark.cn/s/57

衡石分析平台使用手册-单机安装及启动

单机安装及启动​ 本文讲述如何在单机环境下进行 HENGSHI SENSE 安装的操作过程。 在安装前请确认网络环境,如果是隔离环境,无法连接互联网时,请先按照 离线环境安装依赖的指导进行依赖包的安装,然后按照本文的指导继续操作。如果网络环境可以连接互联网,请直接按照本文的指导进行安装。 准备工作​ 请参考安装环境文档准备安装环境。 配置用户与安装目录。 在操作前请检查您是否有 sud

线性因子模型 - 独立分量分析(ICA)篇

序言 线性因子模型是数据分析与机器学习中的一类重要模型,它们通过引入潜变量( latent variables \text{latent variables} latent variables)来更好地表征数据。其中,独立分量分析( ICA \text{ICA} ICA)作为线性因子模型的一种,以其独特的视角和广泛的应用领域而备受关注。 ICA \text{ICA} ICA旨在将观察到的复杂信号

【软考】希尔排序算法分析

目录 1. c代码2. 运行截图3. 运行解析 1. c代码 #include <stdio.h>#include <stdlib.h> void shellSort(int data[], int n){// 划分的数组,例如8个数则为[4, 2, 1]int *delta;int k;// i控制delta的轮次int i;// 临时变量,换值int temp;in

三相直流无刷电机(BLDC)控制算法实现:BLDC有感启动算法思路分析

一枚从事路径规划算法、运动控制算法、BLDC/FOC电机控制算法、工控、物联网工程师,爱吃土豆。如有需要技术交流或者需要方案帮助、需求:以下为联系方式—V 方案1:通过霍尔传感器IO中断触发换相 1.1 整体执行思路 霍尔传感器U、V、W三相通过IO+EXIT中断的方式进行霍尔传感器数据的读取。将IO口配置为上升沿+下降沿中断触发的方式。当霍尔传感器信号发生发生信号的变化就会触发中断在中断

kubelet组件的启动流程源码分析

概述 摘要: 本文将总结kubelet的作用以及原理,在有一定基础认识的前提下,通过阅读kubelet源码,对kubelet组件的启动流程进行分析。 正文 kubelet的作用 这里对kubelet的作用做一个简单总结。 节点管理 节点的注册 节点状态更新 容器管理(pod生命周期管理) 监听apiserver的容器事件 容器的创建、删除(CRI) 容器的网络的创建与删除

PostgreSQL核心功能特性与使用领域及场景分析

PostgreSQL有什么优点? 开源和免费 PostgreSQL是一个开源的数据库管理系统,可以免费使用和修改。这降低了企业的成本,并为开发者提供了一个活跃的社区和丰富的资源。 高度兼容 PostgreSQL支持多种操作系统(如Linux、Windows、macOS等)和编程语言(如C、C++、Java、Python、Ruby等),并提供了多种接口(如JDBC、ODBC、ADO.NET等

OpenCV结构分析与形状描述符(11)椭圆拟合函数fitEllipse()的使用

操作系统:ubuntu22.04 OpenCV版本:OpenCV4.9 IDE:Visual Studio Code 编程语言:C++11 算法描述 围绕一组2D点拟合一个椭圆。 该函数计算出一个椭圆,该椭圆在最小二乘意义上最好地拟合一组2D点。它返回一个内切椭圆的旋转矩形。使用了由[90]描述的第一个算法。开发者应该注意,由于数据点靠近包含的 Mat 元素的边界,返回的椭圆/旋转矩形数据