多项式中的牛顿迭代法

news/2025/3/29 3:06:28/文章来源:https://www.cnblogs.com/EmilyDavid/p/18792888

牛顿迭代法(\(Newton's Method\)

快速求零点

对于一条曲线,我们要求它的零点。

先猜测一个零点,再求出该点处的切线,然后取切线的零点重复以上过程,可以快速逼近零点。

快速求函数

我们已知\(G(F(x))=0\),求\(F(x)\)

先猜测一个函数\(F_0(x)\),套用求零点的过程,可以得到\(F_1(x)=F_0(x)-\frac{G(F_0(x))}{G'(F_0(x))}\),再不断重复提高精度。

正确性证明:

最开始手动求出\([x^0]F(x)\),记作\(F_0(x)\)

假设现在已知\(F_k(x)\equiv F(x)\pmod{x^n}\)

我们把\(G(F(x))\)\(F_k(x)\)处泰勒展开得到:

\[G(F(x))=\sum\limits_{i\ge 0}\frac{G^{(i)}(F_k(x))}{i!}(F(x)-F_k(x))^i=0 \]

注意到\(F(x)-F_k(x)\)的次数至少为\(n\),所以:

\[\forall_{i\ge 2}(F(x)-F_k(x))^i\equiv 0\pmod{x^{2n}} \]

所以\(\sum\limits_{i=0}^1\frac{G^{(i)}(F_k(x))}{i!}(F(x)-F_k(x))^i\equiv 0\pmod{x^{2n}}\)

即:

\[G(F_k(x))+G'(F_k(x))(F(x)-F_k(x))\equiv 0\pmod{x^{2n}} \]

整理后得到:

\[F(x)\equiv F_k(x)-\frac{G(F_k(x))}{G'(F_k(x))}\pmod{x^{2n}} \]

我们令\(F_{k+1}(x)\equiv F_k(x)-\frac{G(F_k(x))}{G'(F_k(x))}\pmod{x^{2n}}\),就可以开始倍增求\(F(x)\)了。

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

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

相关文章

王炸!Kafka 4.0 重磅发布,Java 8 和 Zookeeper 彻底被抛弃了,一个时代结束了!

大家好,我是R哥。 Kafka 4.0 终于来了!这次更新可不只是常规的版本优化,而是一次重大架构调整,尤其是不再支持 Java 8,彻底移除了 Zookeeper,正式启用 KRaft 模式,让 Kafka 变得更加独立、高效。 除了这个重磅更新外,Kafka 4.0 还砍掉了一些历史包袱,还带来几个重磅新…

会计学-开篇(一)

介绍 是不是觉得会计学极其枯燥,死板,无趣?因为会计准则本身就是人为设定的游戏规则,整个学习过程感觉就是一直在背诵各种条条框框,完全没有创造性。相比之下,不管是数学,编程,物理学,生物学,投资学,都显得很有意思,这些学科都是帮助我们通向创造。而会计学似乎只会…

会计学-开篇

介绍 是不是觉得会计学极其枯燥,死板,无趣?因为会计准则本身就是人为设定的游戏规则,整个学习过程感觉就是一直在背诵各种条条框框,完全没有创造性。相比之下,不管是数学,编程,物理学,生物学,投资学,都显得很有意思,这些学科都是帮助我们通向创造。而会计学似乎只会…

2 小时,我搭好了“一物一码” 的设备巡检管理系统!

说实话,以前每次看到设备巡检表上那些手写的记录,我就头疼——字迹潦草、容易丢、查起来还费劲。直到昨天下午,我实在忍不了了,决定自己动手搞个"一物一码"的巡检系统。 从零到落地,现在所有设备贴个二维码,手机一扫就能查记录、报故障,还能自动生成报表! 老…

2022-PTA正式赛-L1-8 静静的推荐(思路)

未达到分数线passAcCode: #include<bits/stdc++.h> using namespace std; int vis[100010]; int main(){int N, K, S, ans = 0;cin >> N >> K >> S;while(N--){int sc1, sc2;cin >> sc1 >> sc2;if(sc1 < 175) continue;if(sc1 >= 17…

GlusterFS 三节点集群部署指南(Heketi 管理模式)

前言:随着信息技术的飞速发展,数据量呈爆炸式增长,传统的集中式存储系统已经难以满足现代应用对存储容量、性能和可靠性的苛刻需求。分布式存储技术应运而生,而 GlusterFS 作为其中的杰出代表,凭借其简单却强大的架构设计,成为了存储领域的一颗新星。 GlusterFS 的魅力不…

性价比拉满!最新的ISP图像处理方案来啦,基于瑞芯微RK3562J全国产工业平台!

今天为大家带来基于瑞芯微RK3562J工业平台的ISP图像处理方案,不仅低成本、低功耗,更能让图像清晰呈现!下面,一起来深入探究看它如何实现!ISP的作用与优势 何为ISP?全称为Image Signal Processor(图像信号处理器),其主要作用是处理前端图像传感器输出的信号,主要功能有…

什么是受管文件传输(MFT)?对企业有何价值和作用?

一、什么是受管文件传输(MFT)? 受管文件传输(Managed File Transfer, MFT),也可以称为“托管文件传输”,是一种专门用于安全、可靠地传输文件的解决方案。与传统的文件传输方式(如电子邮件、FTP)相比,MFT提供了更高的安全性、可管理性和自动化能力,能够满足企业对数…

FreeSWITCH 异常重启后会话快速恢复

FreeSWITCH 异常重启后会话快速恢复 配置 会话追踪sip_profiles/*.xml<param name="track-calls" value="true"/>(官方文档描述会有轻微性能影响) 存储恢复数据、修改数据库位置 重要 fs默认使用sqlite进行恢复数据的存储,也可以正常使用,但此处…

SpringSecurity5(10-动态权限管理)

Spring Security动态权限管理通过实时更新权限配置,支持按需调整用户权限,实现灵活的访问控制。结合注解、表达式等方式,能够动态加载和检查用户权限,无需重启系统,确保应用在复杂业务场景下的安全性与可扩展性,提升用户体验与管理效率。授权流程SpringSecurity 的授权流…

【2025年企业必备】这款跨网跨域传输软件 支持多达6种传输方式!

在全球化经济的推动下,企业的业务范围不断扩展,跨网跨域的数据传输需求日益增长。无论是跨国企业的内部协作,还是与合作伙伴之间的数据交换,高效、安全、可靠的跨网跨域传输软件已成为企业运营不可或缺的一部分。 传统的文件传输方式,如电子邮件附件、FTP服务器、甚至是物…