问题:

牛顿法公式:
xk+1=xk−f(xk)f′(xk)x_{k+1}=x_k-\frac{f(x_k)}{f'(x_k)}xk+1​=xk​−f′(xk​)f(xk​)​

从而求出来如下:

xk+1=xk−xk2−3xk−exk+22xk−3−exkx_{k+1}=x_k-\frac{x_k^2-3x_k-e^{x_k}+2}{2x_k-3-e^{x_k}}xk+1​=xk​−2xk​−3−exk​xk2​−3xk​−exk​+2​

初值为x0=1x_0=1x0​=1

设置与真实根的误差不超过10−510^{-5}10−5。

x0=1
x=x0-(x0*x0-3*x0-np.power(np.e,x0)+2)/(2*x0-3-np.power(np.e,x0))
dx=x-x0
max_dx=0.00001
while(dx>max_dx or dx<-1*max_dx):x0=xx=x0-(x0*x0-3*x0-np.power(np.e,x0)+2)/(2*x0-3-np.power(np.e,x0))dx=x-x0print(x)


我们发现,超级快就收敛了,而且从第二次和第三次之间前10位小数都是一样的,也就是说精度其实已经超过了10−510^{-5}10−5,而是10−1010^{-10}10−10,无心插柳柳成荫。这说明牛顿法是很好的方法,平方收敛牛逼。

牛顿法求解非线性方程的根相关推荐

  1. 牛顿法求解非线性方程

    基础知识 [步骤1] 准备:选定初始近似值x0x_0x0​,计算f0=f(x0),f0′=f′(x0)f_0=f(x_0),f_0'=f'(x_0)f0​=f(x0​),f0′​=f′(x0​) [步 ...

  2. 数值计算大作业:非线性方程求根(二分法、牛顿法、弦截法在Matlab实现)

    作为研究生的入门课,数值计算的大作业算是所有研究生开学的重要编程作业. 我把二分法.牛顿法.弦截法求解非线性方程求根的数值计算作业在MATLAB中编程实现.具体的程序详细标注后放在文章附录了,算法数学 ...

  3. 2021-01-07 matlab数值分析 非线性方程求根 牛顿法

    matlab数值分析 非线性方程求根 牛顿法 %牛顿法求非线性方程的根: % 输入:fun--非线性函数:dfun--非线性函数导数:x0--初始值:tol--精度: % 输出:x--非线性方程数值根 ...

  4. matlab 牛顿法 初值,非线性方程的数值解法:牛顿法及牛顿下山法(含Matlab程序)...

    牛顿法及牛顿下山法 简介:牛顿迭代法是求解单变量非线性方程f(x)=0中最实用的方法,该方法在单根附近二阶收敛.但应用时要选用较好的初值x0近似才能保证迭代收敛.为克服这一缺点,可使用牛顿下山法.下面 ...

  5. 非线性方程求根的牛顿法

    牛顿迭代法的推导: 线性方程容易求解,但对于非线性方程,若能用某个线性方程来近似,求出该线性方程的解,即可得到原非线性方程的一个近似解. 设已知非线性函数的一个近似零点是,用在该点的Taylor展开式 ...

  6. 求解非线性方程组的牛顿法c语言,牛顿下山法求解非线性方程(组)(C实现)...

    1.算法描述 (1)符号说明与基本假设 对于非线性方程组:                                                        (1) 引入向量: 可将(1) ...

  7. c语言编程非线性方程求解,c语言计算机编程三种方法求解非线性方程

    c语言计算机编程三种方法求解非线性方程 本 科 专 业 学 年 论 文题 目:非线性方程求解比较姓 名: 何 娟 专 业: 计算机科学技术系 班 级: 08 级本科(2)班 指 导 老 师: 刘 晓 ...

  8. 计算物理学(数值分析)上机实验答案6、非线性方程求根

    实验六.非线性方程求根 ​ 在科学研究与工程技术中常会遇到求解非线性方程的问题.二分法简单易行, 但收敛较慢,仅有线性收敛速度.而且该方法不能用于求偶数重根或复根,但可 以用来确定迭代法的初始值.牛顿 ...

  9. c语言编程非线性方程求解,c语言计算机编程三种方法求解非线性方程.doc

    c语言计算机编程三种方法求解非线性方程.doc 本 科 专 业 学 年 论 文题 目非线性方程求解比较姓 名 何 娟 专 业 计算机科学技术系 班 级 08 级本科(2)班 指 导 老 师 刘 晓 娜 ...

最新文章

  1. LeetCode: 107. Binary Tree Level Order Traversal II
  2. Cisco设备静态NAT基本配置步骤
  3. SQL查询 的一些原则
  4. Tensorflow入门——训练结果的保存与加载
  5. Serverless 解惑——函数计算如何访问 PostgreSQL 数据库
  6. 纠结的velocity log那些事(出现Permission denied)
  7. AliOS Things网络适配框架 - SAL
  8. 关于Chrome字体模糊解决方案
  9. numpy 创建数组
  10. 开源要安全!Linux 基金会、GitHub、GitLab、Google 联合起来了!
  11. javascript匿名函数及闭包深入理解及应用
  12. html5---资料查询
  13. 艾迈斯半导体推出新款高性能读取IC,推动医疗和工业数字化X射线设备制造商降本增效
  14. 上传Android library到JitPack
  15. Learn Git Branching 笔记
  16. mvp的全称_库里常规赛mvp是哪一年,分别于2014-15和2015-16赛季荣膺
  17. Java常见的加密解密
  18. 浅谈5G网络及其应用
  19. 内网渗透 关于GPO
  20. Windows下的Socket编程——基本步骤

热门文章

  1. 乘风破浪的马里奥!这个AI带你一口气通29关,你猜连AI都过不去的是哪3关?
  2. Google开源新TensorFlow运行时TFRT,将取代现有进行时
  3. 独家 | 在树莓派+Movidius NCS上运用YOLO和Tiny-YOLO完成目标检测(附代码下载)
  4. 最新计算机技能需求排名出炉:Python仅排第三,第一你猜得到吗?(附链接)...
  5. 只花5-10分钟评审,还不提供拒稿理由,IJCAI就“枪毙”42%论文,网友:一脸懵逼...
  6. 【秋招必备】LeetCode神器,算法刷题宝典.pdf
  7. 如何评价英伟达9月2日凌晨发布的最强消费级显卡RTX 3090?
  8. 推荐10款Windows系统必备的高效软件!!!
  9. 最新BERT相关论文汇总
  10. GPT Timeline