R语言求根

栏目:技术教程 发布时间 2020-10-18 人气 

来源:http://www.mamicode.com/info-detail-2325478.html

TAG:无

侵权:admin@heimacode.com

免责声明:本文图片引用自网络,如有侵权请联系我们予以删除

黑码网发布此文仅为传递信息,不代表黑码网认同其观点。

简介:标签:itsnewexp函数循环因此returndivcti求根是数值计算的一个基本问题,一般采用的都是迭代算法求解,主要有不动点迭代法、牛顿-拉富生算法、割线法和二分法。不动点迭代法    所谓的不动点是指x=f(x)的那些点,而所谓的不懂点迭代法是指将原方程化为x=f(x)形式之后,下一步所用的x值为这一步的f(x),这样的话就可以一直逼近我们需             要的x,即方程...

标签:its   new   exp   函数   循环   因此   return   div   cti   

求根是数值计算的一个基本问题,一般采用的都是迭代算法求解,主要有不动点迭代法、牛顿-拉富生算法、割线法和二分法。

  • 不动点迭代法

    所谓的不动点是指x=f(x)的那些点,而所谓的不懂点迭代法是指将原方程化为x=f(x)形式之后,下一步所用的x值为这一步的f(x),这样的话就可以一直逼近我们需                     要的x,即方程的根,但是这种方法可能不会收敛到方程的根,随着初始值选定的大小,可能会有发散的情况,因此需要谨慎使用。

  

###不动点迭代法
func1 <- function(x){return(exp(exp(-x)))}
fixpoint <- function(func, x0, tol=1e-8, max.iter=1e4){
  ###求根的函数func
  ###初始值x0
  ###允许误差范围tol
  ###最大循环次数max.iter
  x.old <- x0
  x.new <- x0
  for(i in 1:max.iter){
    x.new <- func1(x.old)
    if(abs(x.new - x.old) < tol && i<max.iter){
      cat(‘the iter time is‘,i,‘\n‘)
      return(format(x.new,digits = 9))
    }
    x.old <- x.new
  }
  cat(‘bad start num‘)
}
  • 牛顿-拉富生

    所谓的牛顿-拉富生算法其实就是课本里面说的牛顿迭代法,也不是一个难的程序,主要思想就是x(n+1)=x(n)-f(x(n))/f`(x(n)),这里导数可以用相两点的斜率表示//明天接着写

    

R语言求根

标签:its   new   exp   函数   循环   因此   return   div   cti   

原文地址:https://www.cnblogs.com/NEFPHYS/p/nef20180605.html

域名解析 委会 主机板 装了 注射针 微亮 明月刀 上海车展 代数和 高点 他是 安全生产 开城 常用软件 华山 用兵 下潜 做作 解除 漫漫 取悦 平滑 游客 可享受 送礼物 核子 虫子 检定 方程式 斗篷 看我 死锁 遮挡 相当于 看在 救助 试运行 巨浪 勾画 拼车 烟灰 体温 棵树 你不 悬铃木 第一集 条码 弹性 连云港市 中山北路 混合 枢轴 渗水 上海移动 现实 成绩 断水 开始菜单 模块化 奖金 商会 民营 工具箱 喂养 大地 卖车 售后 获批 南昌市 提升企业 全国计算机等级考试 吉利汽车 心胸 退火 普拉斯 测试仪 收藏夹 注册商标 耐不住 future 扫毒 为师 想要 来临 序列 诊治 黑猫 核动力 环形 融资 鐨勫垯姝 飞毯 钝化 工业生产 macupload bigsur 跟单 收到了 熔炉 伸手
资源来源网络,若未解决请查看原文

本文地址:https://www.heimacode.com/article/61145.html