我目前正在进行一些数值传热计算,这需要一个很好的互补误差函数(erfc)的近似。其中一个传热方程由exp(s^2)和erfc(s)的乘积组成,这意味着(对于较大的s值)一个非常大的数乘以一个非常小的数。因此,电流变流控制的精度需要相当高。否则,对大量s的计算将是不稳定的。在

由于处理not的数学函数(例如sqrt()、exp()和pow())不能取双倍数,所以我使用float,这给了我对于大值s的稳定性问题。对于大值,我的意思是大约4到5。例如,4.5的erfc应该是1.97*10^(-10)。在

…现在来回答问题。在

1)是否可以在处理的数学函数(例如平方根、exp和幂)中使用双精度来获得互补误差函数值的更精确表示?在

2)据我所知,Python也在使用float,而不是double,这意味着(至少在我看来)Processing和Python对于互补错误函数应该得到相同的结果。然而,事实并非如此。Python对错误函数的逼近比处理更精确。怎么会这样?在

下面显示了Python和互补误差函数近似的处理版本。在

Python:from math import sqrt

from math import pi

from math import exp

n=1000 #number of intervals

b=float(1) #integration end

h=b/n #Interval width

leftRiemann=0

for s in range(0, n):

leftRiemann=leftRiemann+h*2/sqrt(pi)*exp(-((s*h)**2))

rightRiemann=0

for s in range(1, n+1):

rightRiemann=rightRiemann+h*2/sqrt(pi)*exp(-((s*h)**2))

trapez=0.5*leftRiemann+0.5*rightRiemann

centrePoint=0

for s in range(1, n+1):

centrePoint = centrePoint+h*2/sqrt(pi)*exp(-((h*(s-0.5))**2))

erfc=1 - (1.0/3*trapez+2.0/3*centrePoint)

print erfc

正在处理:

^{pr2}$

Python:erfc(1)=0.15729920705

处理:erfc(1)=0.15729982

表值:erfc(1)=0.157299

python计算和差积商_Python与处理计算的区别相关推荐

  1. python计算和差积商_Python计算股票投资组合的风险价值(VaR)

    原文链接:http://tecdat.cn/?p=17758​tecdat.cn 什么是风险价值(VaR)? 风险价值(VaR)用于尝试量化指定时间范围内公司或投资组合中的财务风险水平.VaR提供了一 ...

  2. python 编辑距离 2组匹配_Python文本相似性计算之编辑距离详解

    编辑距离 编辑距离(Edit Distance),又称Levenshtein距离,是指两个字串之间,由一个转成另一个所需的最少编辑操作次数.编辑操作包括将一个字符替换成另一个字符,插入一个字符,删除一 ...

  3. python日历下拉框_python日历来计算月份倒退

    我们正在尝试在python中创建日历函数.我们创建了一个小型的内容管理系统,要求是,网站的右上角会有一个下拉列表,它会给出选项 – 月 – 1个月,2个月,3个月等等. ..,如果用户选择8个月,那么 ...

  4. python计算auc的代码_python是怎么计算auc指标的?

    1.安装scikit-learn 1.1Scikit-learn 依赖Python (>= 2.6 or >= 3.3), NumPy (>= 1.6.1), SciPy (> ...

  5. python一个数的阶乘_python整数阶乘计算

    广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 问题本身很简单,主要是通过这个小问题来演示python的一些用法,例如测试代码运 ...

  6. python爬虫实例电商_Python实现爬取并分析电商评论

    现如今各种APP.微信订阅号.微博.购物网站等网站都允许用户发表一些个人看法.意见.态度.评价.立场等信息.针对这些数据,我们可以利用情感分析技术对其进行分析,总结出大量的有价值信息.例如对商品评论的 ...

  7. python根据生日求年龄_python根据出生日期计算年龄的代码

    python根据出生日期计算年龄的代码,运行后会提醒用户输出出生的年月日,然后输出年龄,可以改写为一个通用函数 from time import * #a function to find your ...

  8. python爬虫实例电商_Python案例:先爬虫爬取,再分析2000款避孕套,最后得出这个结论...

    在学过了Python爬虫.数据分析之后,曾几何时,这样一个问题一直缠绕着我. Python爬虫能用来干点什么?Python数据分析能不能让这剪不断理还乱的数据变得更有意义? 也许,想来想去,只有用来干 ...

  9. python计算时间装饰器_python 写一个计算运行时间的装饰器

    面试题之中的一个. 写一个装饰器wraps,它接收一个參数t.假设函数运行时间超过1秒,输出"bad",否则输出"goods". 首先,计算函数的运行时间: i ...

最新文章

  1. 中文 查询_查询商标,商标注册通过分析的几个小技巧
  2. 业务逻辑组件化android,AppJoint 极简 Android 组件化方案
  3. Unity中制作游戏的快照游戏支持玩家拍快照
  4. 微软Hololens设备 浅分析
  5. mysql主主mycat_MySQL数据库主主复制并基于MyCAT实现高可用
  6. 【框架】[Spring]XML配置实现AOP拦截-切点:JdkRegexpMethodPointcut
  7. 关于 AfxSocketInit() 的使用
  8. Java语言编程基础:数组
  9. nyoj--120--校园网络(scc+缩点)
  10. jzoj3085-图的计数【组合数,数论】
  11. 深入理解脚本化CSS系列第二篇——查询计算样式
  12. 输出质数(java)
  13. yum php5.6源码,5.2 YUM升级PHP5.6
  14. 开课吧:全栈工程师的关键开发技能(硬实力)
  15. ​使用 QuickBI 搭建酷炫可视化分析
  16. LINGO 11.0安装教程
  17. [ffmpeg][goav][issues] goav Dictionary crash double free
  18. mysql练习题库一
  19. 搭讪达人这样和女生发短信
  20. 女朋友让我深夜十二点催她睡觉,我用Python轻松搞定

热门文章

  1. QT的QWhatsThis类的使用
  2. 南邮微型计算机实验,南邮 微机原理 微型计算机与接口技术 微机原理实验.doc...
  3. php与tcp哪个快,PHP程序员如何理解TCP协议
  4. MySQL 5.7.18的安装与主从复制(转自:https://www.baidu.com/home/news/data/newspage?nid=9485770887287731252n_typ)
  5. 02_创建Git仓库,克隆仓库,git add,git commit,git push,git pull,同行冲突,不同行冲突的结局方案,git mergetool的使用
  6. 1.信号处理之:kill(),alarm(),pause()函数
  7. 数据库元数据数据字典查询_4_列出给定表的索引信息
  8. sudo apt install镜像_将Docker镜像安全扫描步骤添加到CI/CD管道
  9. pytorch之torch.backends.cudnn.benchmark=True——使用 GPU 来加速网络的训练
  10. C++11中range-based for loops中与的区别