Mathematica中同参数下N[CDF[NormalDistribution[m,s],x]]与CDF结果差异疑问
问题原因分析
核心触发点
这是Mathematica新手高频遇到的两类场景,你可以按以下逻辑逐一排查:
- 变量残留导致两次计算参数不一致
Mathematica的全局变量会默认保留最后一次赋值结果,除非手动清除。大概率你在直接计算CDF[NormalDistribution[m,s],x]之前,m/s/x三个变量中至少有一个已经被其他运行过的代码覆盖为其他值,而你运行带N的表达式之前又重新给三个变量赋了正确的m=12/s=25/x=10,两次计算的实际参数不同,结果自然有差异。
你可以运行以下代码验证:
正常情况下ClearAll[m, s, x]; m = 12; s = 25; x = 10; res1 = N[CDF[NormalDistribution[m, s], x]] res2 = CDF[NormalDistribution[m, s], x] N[res2]res1和N[res2]结果完全一致,均为0.468119。 - 符号表达式读取错误
不带N的CDF会返回符号形式的结果1/2 Erfc[1/(25 Sqrt[2])],如果你没有调用数值转换函数,而是手动估算、或者用了非官方的数值提取工具读取该符号表达式的结果,很可能得到错误的0.324195。
规避方案
每次运行独立计算前先用ClearAll[变量名]清除对应变量的残留值,需要数值结果时统一用N[]包裹表达式,避免手动计算符号表达式的数值。
内容的提问来源于stack exchange,提问作者T123
相关产品推荐
相关产品推荐

