位置: 编程技术 - 正文

Python编程实现二分法和牛顿迭代法求平方根代码(python erf)

编辑:rootadmin

推荐整理分享Python编程实现二分法和牛顿迭代法求平方根代码(python erf),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:python2//4,python2ide,python 2to3,python2ide,python2怎么用,python 2,python 2to3,python 2to3,内容如对您有帮助,希望把文章链接给更多的朋友!

求一个数的平方根函数sqrt(int num) ,在大多数语言中都提供实现。那么要求一个数的平方根,是怎么实现的呢?实际上求平方根的算法方法主要有两种:二分法(binary search)和牛顿迭代法(Newton iteration)

1:二分法

求根号5

a:折半: 5/2=2.5b:平方校验: 2.5*2.5=6.>5,并且得到当前上限2.5c:再次向下折半:2.5/2=1.d:平方校验:1.*1.=1.<5,得到当前下限1.e:再次折半:2.5-(2.5-1.)/2=1.f:平方校验:1.*1.=3.<5,得到当前下限1.

每次得到当前值和5进行比较,并且记下下下限和上限,依次迭代,逐渐逼近平方根:

运行结果:1 2. 1. 1. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2...[Finished in 0.1s]

经过次二分法迭代,得到的值和系统sqrt()差别在0.,精度在亿分之一,

0.需要迭代8次

因此,在对精度要求不高的情况下,二分法也算比较高效的算法。

2:牛顿迭代

仔细思考一下就能发现,我们需要解决的问题可以简单化理解。

从函数意义上理解:我们是要求函数f(x)=x²,使f(x)=num的近似解,即x²-num=0的近似解。

Python编程实现二分法和牛顿迭代法求平方根代码(python erf)

从几何意义上理解:我们是要求抛物线g(x)=x²-num与x轴交点(g(x)=0)最接近的点。

我们假设g(x0)=0,即x0是正解,那么我们要做的就是让近似解x不断逼近x0,这是函数导数的定义:

可以由此得到

从几何图形上看,因为导数是切线,通过不断迭代,导数与x轴的交点会不断逼近x0。

对于一般情况:

将m=2代入:

运行结果:1 2. 2. 2...

精确到亿分之一,牛顿法只迭代了3次,是二分法的十倍

3:利用牛顿法求开立方

微积分、概率、线代是高级算法的基础课。可是,这么多年,已经忘得差不多了..............................

总结

标签: python erf

本文链接地址:https://www.jiuchutong.com/biancheng/377056.html 转载请保留说明!

上一篇:Python编程给numpy矩阵添加一列方法示例(Python编程给定a=1,b=4,输出a+b的值)

下一篇:Python数据可视化正态分布简单分析及实现代码(python数据可视化课后题答案)

  • 个体户也要做实收资本吗
  • 制造业的加工费可以放制造费用的二级科目加工费吗
  • 一般纳税人租赁费税率5%是简易征收吗?
  • 兼职人员怎么申报工资
  • 小微企业免增值税2023年政策
  • 员工报销没有发票怎么入账
  • 公司茶水间有什么
  • 企业入统的利弊
  • 培训产生的住宿费算培训费吗
  • 个体工商户能否转让
  • 会计政策变更与会计估计的审计
  • 二手车交易怎么办理过户手续
  • 固定资产出租的收入计入什么科目
  • 当月未认证发票分录
  • 帮其他企业推广产品
  • 企业如何申请定增资金
  • 不含税单价怎么换成含税单价
  • 新注册的个体户怎么开发票
  • 对方给我们开发票
  • 银行对账单不平衡
  • 以前年度损益调整会计分录
  • 租房税收缴款书在哪里打印
  • win10开机无限werfault一直刷新
  • 财务原始凭证
  • php的编辑工具有哪些
  • php利用云片网实现短信验证码功能的示例代码
  • php string
  • 出口退税运费的账务处理
  • 房地产开发企业销售自行开发的房地产项目
  • 仙岩寺庙
  • php中的九大缓存技术
  • 小规模纳税人的企业所得税怎么算
  • 财务的几张报表
  • 成功解决用英语怎么说
  • 收到投资款的会计科目怎么入账
  • vue全家桶学多久能上手项目
  • chkconfig命令参数
  • linux symbolic link
  • 现金流量表的附表如何编制
  • 帐没做完怎么报增值税
  • 个体工商户税务登记需要哪些资料
  • 进料余料结转 报关单
  • webpack怎么优化
  • mysql慢查询日志查看
  • sqlyog提示
  • 社保工伤退回分录
  • 普票需要填写全额发票吗
  • 建筑企业结转成本依据准则
  • mysql 缓冲区
  • 应付账款的入账时间为
  • 代金券抵帐怎么退款
  • 当月无销售要计提工资吗
  • 建筑业异地预缴税款最新规定
  • 个体工商户可以变更法人吗
  • 发票要不回来怎么办
  • aspnet开发和运行环境
  • 注册表修改大全
  • linux的文件和目录的区别
  • win10预览版
  • service5.exe - service5是什么进程 有什么用
  • win 8系统怎么样
  • centos7安全
  • 安装win7旗舰版变了家庭版
  • 微软新版windows
  • 在linux系统中拥有
  • 常用的加锁方式
  • windows下用CMD调用COM口
  • node实战
  • unity引擎用什么模拟器
  • jquery自定义事件
  • js获取tbody
  • linux反弹shell方法
  • javascript判断题
  • 套接字java
  • 代码实现观察者模式
  • 车船税的纳税期限是
  • 建筑企业要交哪几项环保税
  • 企业所得税率2023年
  • 苏州工业园区国资委
  • 广东省电子税务局官网
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设