为何差距这么大

一个数用两种类型表示算出的结果为何差距这么大：
N[Integrate[
Exp[(10^(-4))*t]*ArcTan[((10^(-4)) - t)^2], {t, 0, 10^(-4)}]/
Integrate[2*t*((10^(-4)) - t), {t, 0, 10^(-4)}], 10]
1.0000000025+ 0.*10^-11i
N[Integrate[Exp[0.0001*t]*ArcTan[(0.0001 -t)^2], {t, 0, 0.0001}]/
Integrate[2*t*(0.0001 - t), {t, 0, 0.0001}], 10]
-32.7418
请教各位这是怎么回事啊

在mathematica语言中，相同小数的工作精度是不如分数的，因此下面一种写法的误差会大于上面一种。
另外，对于这种数值较小的定积分你可以尝试使用NIntegrate来代替N[Integrate]，在你的例子中这两个式子在NIntegrate下给出一样的结果1.

N这个函数不能提升精度，只能把高精度数舍入为低精度数，而默认的MachinePrecision数字（不知道我在说啥的执行下 Precision[0.0001] ）舍入优先级是最低的。
这个问题不用NIntegrate也是可以的，只要你精度设置得当：
With[{coef = 0.0001`32},
Integrate[Exp[coef*t]*ArcTan[(coef - t)^2], {t, 0, coef}]/
Integrate[2*t*(coef - t), {t, 0, coef}]]
(* 1.000000002500 + 0.*10^-12 I *)

另参百度知道问题《Mathematica实数保留n位有效数字和保留n位小数的最简单的方法分别是什么？》（编号588921103153505885）

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频！

7回复贴，共1页

<<返回mathematica吧

分享到:

日	一	二	三	四	五	六