由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Computation版 - 为什么两个等效表达式数值上不相同,如何减少roundoff error?
相关主题
mathematica 提问 "integrand not numerical"请问tecplot可以画函数曲线吗?
如何在mathematica中将符号消去具体值[转载] matlab问题请教
[转载] 会数值积分的请帮忙看一下unidentified_title
Re: mathematica 提问 "integrand not numeri[转载] How to minimize this variance?
how to remove a singularity in an integratanh函数
怎么求这个转角?一个循环的问题
请教Ansys高手一个问题非线性优化的方法
[转载] 高维空间中,点到曲线最短距离有显性表达式吗?Mathematica下面做function fit
相关话题的讨论汇总
话题: 数值话题: roundoff话题: 函数话题: 表达式话题: 两个
进入Computation版参与讨论
1 (共1页)
b****e
发帖数: 906
1
解析上log(z/(z-1))等价于Integrate[1/(z-x),{x,0,1}],其中z为任意复数, 但如果
把具体的数值放进去似乎情况有点复杂。log(z/(z-1))作为一个基本函数直接调用
built-in函数即可,积分在Mathematica中调用NIntegrate[]来数值计算到收敛精度10^
(-17)。然后比较两个数值,结果发现本来应该完全相等的两个数可以差上10^(-15),
远远超出计算机的机器精度。请谁能告诉我是怎么回事?是通过基本函数计算出来的结
果更可靠还是应该相信通过数值收敛出来的结果?谢谢
m*****n
发帖数: 3575
2
定积分如果用的是柱图法就比较烂,梯形法好很多,当然还有更好的
一般来讲计算次数越少的结果是越准确的,解析的好
m*******s
发帖数: 3142
3
Compute log(1+x) accurately for small values of x
https://www.mathworks.com/help/matlab/ref/log1p.html
b****e
发帖数: 906
4
解析上log(z/(z-1))等价于Integrate[1/(z-x),{x,0,1}],其中z为任意复数, 但如果
把具体的数值放进去似乎情况有点复杂。log(z/(z-1))作为一个基本函数直接调用
built-in函数即可,积分在Mathematica中调用NIntegrate[]来数值计算到收敛精度10^
(-17)。然后比较两个数值,结果发现本来应该完全相等的两个数可以差上10^(-15),
远远超出计算机的机器精度。请谁能告诉我是怎么回事?是通过基本函数计算出来的结
果更可靠还是应该相信通过数值收敛出来的结果?谢谢
m*****n
发帖数: 3575
5
定积分如果用的是柱图法就比较烂,梯形法好很多,当然还有更好的
一般来讲计算次数越少的结果是越准确的,解析的好
m*******s
发帖数: 3142
6
Compute log(1+x) accurately for small values of x
https://www.mathworks.com/help/matlab/ref/log1p.html
v****r
发帖数: 11
7
原因很简单,在有奇点的区域上,integrate函数不准,需要用特定网格可以变准,但
是我怀疑这个函数有这种能力。
你觉得应该完全一样就是搞笑, 能手算的当然手算表达式最准,数值积分永远是逼近
算法。除了多项式和一些特定函数,绝大多数都不会和真实值完全一样,还有这些软件
里得看怎么写收敛条件的。
1 (共1页)
进入Computation版参与讨论
相关主题
请教一个模型表达式how to remove a singularity in an integr
附件中的数据拟合怎么做?怎么求这个转角?
如何求取最优值?请教Ansys高手一个问题
请教一个曲线拟合的问题[转载] 高维空间中,点到曲线最短距离有显性表达式吗?
mathematica 提问 "integrand not numerical"请问tecplot可以画函数曲线吗?
如何在mathematica中将符号消去具体值[转载] matlab问题请教
[转载] 会数值积分的请帮忙看一下unidentified_title
Re: mathematica 提问 "integrand not numeri[转载] How to minimize this variance?
相关话题的讨论汇总
话题: 数值话题: roundoff话题: 函数话题: 表达式话题: 两个