首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >NIntegrate会将一个明显非零的积分计算为零。被积函数是一组变形为分布的数据点

NIntegrate会将一个明显非零的积分计算为零。被积函数是一组变形为分布的数据点
EN

Stack Overflow用户
提问于 2012-12-12 09:31:44
回答 1查看 289关注 0票数 1

我有一个表示概率分布的数据点列表。我需要在这个发行版上进行集成。但是,由于我没有函数,并且只有一组数据点,所以我想出了以下代码来表示概率分布:

代码语言:javascript
复制
dList1 = Import["Z-1.txt", "Table"];
dList2 = Import["Z_over-1.txt", "Table"];

dDist[X_,sym_] := (

  dList = 0;

  If[sym,
    dList = dList1;
  ,
    dList = dList2;
  ];

  val = 0;

  If[Abs[X] < Pi,
    i = 2;
    While[dList[[i]][[1]] < X, i++];

    width = dList[[i]][[1]] - dList[[i-1]][[1]];

    difX = dList[[i]][[1]] - X;
    difY = dList[[i]][[2]] - dList[[i-1]][[2]];

    val = dList[[i-1]][[2]] + (1-(difX/width)) difY;
  ];

  Return[val];
);

其中数据点的集合在文本文件中。

执行以下命令:

代码语言:javascript
复制
Plot[dDist[x, True], {x, -1, 1}]

提供了以下内容:

然而,执行此操作时:

代码语言:javascript
复制
NIntegrate[dDist[x, True], {x, -1, 1}]

计算结果为零,同时显示以下警告:

我已经尝试增加MinRecursion,但无济于事。我不确定该怎么做,并且愿意接受任何建议,包括修改dDist函数。

EN

回答 1

Stack Overflow用户

发布于 2012-12-13 10:47:40

如果没有数据可用,我建议使用表格来创建在不连续的地方分开的分段函数。NIntegrate应该可以处理那些没有问题的问题。

票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/13831385

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档