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

【学习笔记】拉格朗日插值

EZ、什么是拉格朗日插值?

众所周知,\(n+1\) 个点可以唯一确定一个 \(n\) 次多项式。

拉格朗日插值法要解决的就是给定 \(n+1\) 个点确定一个多项式 \(f(x)\),求出在自变量 \(x=k\) 时多项式的取值。

拉格朗日插值法的思想和 CRT 非常像——把每一个维度独立拆开来。

考虑对一个点 \((x_i,y_i)\) 构造一个子函数 \(f_i(x)\),使得 \(f_i(x_i)=y_i\),且对于每一个 \(j \neq i\),使得 \(f_i(x_j)=0\)

如果你上过奶猫老师的数论课,就会意识到这些子函数和 CRT 的 \(\omega_i\) 的思想是完全一样的!

然后就是怎么算子函数了。不需要 exgcd,不需要神秘科技,就是朴素的 \(O(n^2)\) 做法。

\(f_i(x_j)=0\) 是容易的,只要让 \(f_i(x)=\displaystyle{\prod_{j=1}^{n}[j \neq i](x-x_j)}\) 就行了。

那怎么修改使得满足另外一个条件呢?其实很简单,加个系数。

最终的结果就是 \(f_i(x)=\displaystyle{\dfrac{y_i\prod_{j=1}^{n}[j \neq i](x-x_j)}{(x_i-x_j)}}\)

全部加起来就做完了。多简单。

其实复杂度可以更低,但是太超标了。我不会。哭哭。

http://www.wxhsa.cn/company.asp?id=6903

相关文章:

  • 一种将离散化状态方程映射为并行多处理器计算机的方法
  • 基本数据类型题目
  • 一种基于动作指令交互的动态活体检测技术,提升人脸识别安全性
  • [系统] Windows 已有office版本和visio不兼容的解决方案
  • CF 2127F Hamed and AghaBalaSar
  • AT_agc055_b [AGC055B] ABC Supremacy
  • “Sequential Thinking MCP Server 和codex等AI工具本身任务拆解功能对比
  • 基于错误xsleak 悬空标记 运用css利用帧计数 -- Pure leak ASIS CTF 2025
  • 网易伏羲:当算法遇见社交,解码游戏世界的连接密码
  • 在 CentOS 7 上安装Nginx和配置http代理
  • 题解:P2624 [HNOI2008] 明明的烦恼
  • 在AI技术快速实现创想的时代,挖掘新需求成为核心竞争力——某知名DevOps学习平台需求洞察
  • Windows Powershell 获取版本version
  • XXL-JOB (1)
  • 记录---Vue3对接UE,通过MQTT完成通讯
  • 《Real-Time Rendering》第一章 介绍
  • C语言基础
  • 公益站Agent Router注册送200刀额度竟然是真的
  • 数据集中valid的作用
  • 深入 RocketMQ 核心源码:从环境搭建到高可用设计的全方位解析
  • 单例模式
  • apache修改默认位置
  • 实用指南:YOLOv11的旋转目标检测改进-(扩展检测头支持旋转框预测,适配遥感场景)
  • 从零到顶会:NLP科研实战手册 - 实践
  • 肝不好能喝酒吗
  • ROS中如何将日志格式设置为行号的形式
  • USB相关的sysfs文件(重要的)【转】
  • 25上第一周
  • 深入解析:RxJava在Android中的应用
  • 模型选择与配置说明