我有一个表示概率分布的数据点列表。我需要在这个发行版上进行集成。但是,由于我没有函数,并且只有一组数据点,所以我想出了以下代码来表示概率分布:
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];
);其中数据点的集合在文本文件中。
执行以下命令:
Plot[dDist[x, True], {x, -1, 1}]提供了以下内容:

然而,执行此操作时:
NIntegrate[dDist[x, True], {x, -1, 1}]计算结果为零,同时显示以下警告:

我已经尝试增加MinRecursion,但无济于事。我不确定该怎么做,并且愿意接受任何建议,包括修改dDist函数。
发布于 2012-12-13 10:47:40
如果没有数据可用,我建议使用表格来创建在不连续的地方分开的分段函数。NIntegrate应该可以处理那些没有问题的问题。
https://stackoverflow.com/questions/13831385
复制相似问题