牛顿迭代法

简介: 牛顿迭代法(Newton's method)又称为牛顿-拉夫逊(拉弗森)方法(Newton-Raphson method),它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。牛顿迭代公式    设r是\(f(x)=0\)的根,选取\(x_0\)作为r的初始近似值,过点\((x_0,f(x_0))\),做曲线\(y=f(x)\)的切线L,L的方程为\(y=f(x_0)+f’(x_0)(x-x_0)\),求出L与x轴交点的横坐标\[x_1=x_0-\frac{f(x_0)}{f’(x_0)}\]    称\(x_1\)为r的一次近似值。

     牛顿迭代法Newton's method)又称为牛顿-拉夫逊(拉弗森)方法(Newton-Raphson method),它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。

牛顿迭代公式

    设r是\(f(x)=0\)的根,选取\(x_0\)作为r的初始近似值,过点\((x_0,f(x_0))\) ,做曲线 \(y=f(x)\)的切线L,L的方程为\(y=f(x_0)+f’(x_0)(x-x_0)\) ,求出L与x轴交点的横坐标

\[x_1=x_0-\frac{f(x_0)}{f’(x_0)}\]

    称\(x_1\)为r的一次近似值。过点\((x_1,f(x_1))\) 做曲线 \(y=f(x)\)的切线,并求该切线与x轴交点的横坐标

\[x_2=x_1-\frac{f(x_1)}{f’(x_1)}\]

    称\(x_2\)为r的二次近似值。重复以上过程,得r的近似值序列,其中,

\[x_{n+1}=x_n-\frac{f(x_n)}{f’(x_n)}\]

    称为r的\(n+1\)次近似值,上式称为牛顿迭代公式

    用牛顿迭代法解非线性方程,是把非线性方程\(f(x)=0\)线性化的一种近似方法。把 \(f(x)\)在点 \(x_0\)的某邻域内展开成泰勒级数

\[f(x)=f(x_0)+f’(x_0)(x-x_0)+\frac{f’’(x_0)(x-x_0)^2}{2!}+…+\frac{f^{(n)}(x_0)(x-x_0)^n}{n!}+R_n(x)\]

     取其线性部分(即泰勒展开的前两项),并令其等于0,即 \(f(x_0)+f’(x_0)(x-x_0)=0\),以此作为非线性方程 \(f(x)=0\)的近似方程,若\(f’(x_0)\neq0\),则其解为

\[x_1=x_0-\frac{f(x_0)}{f’(x_0)}\]

     这样,得到牛顿迭代法的一个迭代关系式:

\[x_{n+1}=x_n-\frac{f(x_n)}{f’(x_n)}\]

    从下面的图中,我们可以看到牛顿迭代的几何意义,每次迭代,都会更加逼近\(f(x)=0\)的解。

img_a3eb93fc70dd437946ab3784c0fcf33a.jpg

     下面用牛顿迭代法在matlib解出方程\(x^2+2xe^x+e^{2x}=0\)的根,首先画出函数的图像,猜测根的大致位置。

函数图像如下图所示:

image

近似结果:

image

x=-1:0.01:1;
y= x.^2+2*x.*exp(x)+exp(2*x)
plot(x,y);
grid on;
clc;clear;
%syms x;
%diff(x^2+2*x*exp(x)+exp(2*x),x,1)
%clear;
x=0.0
for i=1:100
    x1=x-(x^2+2*x*exp(x)+exp(2*x))/(2*x + 2*exp(2*x) + 2*exp(x) + 2*x*exp(x))
    if(abs(x1-x)>0.0001)
        x=x1;
    else
        break;
    end
end
i
View Code
相关文章
|
6月前
|
存储 自然语言处理 算法
109_噪声鲁棒微调:对抗训练
在当今大语言模型(LLM)的广泛应用中,模型的鲁棒性问题日益凸显。对抗性攻击通过在输入中添加微小但精心设计的扰动,能够误导模型产生错误输出,这对依赖LLM的关键系统构成了严重威胁。噪声鲁棒微调作为提升模型抵抗对抗攻击能力的重要技术,正成为大模型安全性研究的核心方向之一。
735 2
|
机器学习/深度学习 算法 数据挖掘
最优化--梯度下降法--牛顿法(详解)
最优化--梯度下降法--牛顿法(详解)
2409 1
|
Java Maven
idea中maven项目pom文件Could not acquire lock(s)
idea中maven项目pom文件Could not acquire lock(s)
8745 2
|
存储 监控 安全
数据泄露后的应对措施:个人和企业的行动指南
数据泄露后的应对措施:个人和企业的行动指南
1859 2
|
Ubuntu
一分钟在Ubuntu 20.04安装QEMU-KVM + Virt-Manage
一分钟在Ubuntu 20.04安装QEMU-KVM + Virt-Manage
|
机器学习/深度学习 算法 数据库
基于深度学习的多人步态识别系统(目前数据集大小124人,准确率96.5%)
基于深度学习的多人步态识别系统(目前数据集大小124人,准确率96.5%)
1961 0
基于深度学习的多人步态识别系统(目前数据集大小124人,准确率96.5%)
|
JSON 小程序 JavaScript
超详细微信小程序开发学习笔记,看完你也可以动手做微信小程序项目
这篇文章是一份全面的微信小程序开发学习笔记,涵盖了从小程序介绍、环境搭建、项目创建、开发者工具使用、文件结构、配置文件、模板语法、事件绑定、样式规范、组件使用、自定义组件开发到小程序生命周期管理等多个方面的详细教程和指南。
|
算法 Java
KMP算法详解及其在字符串匹配中的应用
KMP算法详解及其在字符串匹配中的应用
|
数据库
1NF | 2NF | 3NF的区分以及什么是函数依赖、部分函数依赖、值传递依赖(最详细的讲解1NF、2NF、3NF的关系)
这篇文章详细讲解了数据库范式中的1NF、2NF和3NF,包括它们的定义、区分方法和如何判断部分函数依赖和传递函数依赖,以及如何将数据表规范化到相应的范式。
1NF | 2NF | 3NF的区分以及什么是函数依赖、部分函数依赖、值传递依赖(最详细的讲解1NF、2NF、3NF的关系)
|
数据采集 算法 编译器
倚天710规模化应用 - 性能优化 -自动反馈优化分析与实践
编译器优化分成静态优化与动态优化,静态优化指传统编译器gcc/llvm时,增加的优化等级,如O1,O2,O3,Ofast,此时,编译器会依据编译优化等级增加一些优化算法,如函数inline、循环展开以及分支静态预测等等。一般情况下,优化等级越高,编译器做的优化越多,性能会更会好。在阿里生产环境中,单纯依赖于静态优化,并不能达到程序运行流畅目的,通过分析CPU硬件取指令、执行指令,往往会出现一些分支预测失败导致iCacheMiss率高的场景,限制了程序的性能进一步提升。基于此,业务引入了动态反馈优化工具,依据生产环境的实际运行数据,反哺指导编译器对程序代码进一步调整编译优化策略,提高分支预准确率