问题标签 [integral]

坚持对这个积分进行数值评估，并且只有做简单正交方案的经验，所以如果这看起来很业余，请多多包涵。下面的嵌套积分使用 Gauss-Hermite(-inf,+inf) 和 Gauss Laguerre 方案的变体 (0,+inf) - 在高斯和 gamma(v/2,2) 密度上。我到了几乎完成积分的地步，中间步骤看起来不错，但我坚持如何组合权重来评估整体积分。对于修改代码/其他想法以编写更好的正交方案来解决问题的建议，我将非常感激。\begin{方程} \int^{\infty} {-\infty}\int^{\infty} {0} \prod_{i=1}^n\Phi \left(\frac{\sqrt{w/v }\,C{i}-a{i}Z}{\sqrt{1a {i^2}}\right)f{z}(Z)f{w}(W)dwdz \end{方程}

% 脚本定义节点、权重、参数然后调用一个主函数和一个子函数

ρ=0.3；nfirms=10；h=repmat(0.1,[1,nfirms]); T=1；R=0.4；v=8；阿尔法=v/2；高斯点=15；

我需要计算积分公式给出的两个变量的退化超几何函数：

我使用 Matlab 进行数值积分：

我得到的结果是 54.9085，我知道这个值是错误的！所以请帮助我计算上述积分的真值，奇点为 0。

我想在 Matlab 中计算以下类型的积分。

它是函数 e^-(u)*u 的积分，边界为零和无穷大。这个积分应该返回 1。

我怎样才能在 Matlab 中做到这一点？

我知道 Erlang 通常不用于这类事情，但是否有已经编写了积分/导数的数值函数（如 GNU Octave 中的四边形，取函数和两个数字）？

有人可以告诉我如何在matlab中计算以下内容：

在哪里

对于某些整数类型，即使浮点值远远超出整数的可表示范围，我如何才能找到最接近浮点类型的某个值的值。

或者更准确地说：

设F为浮点类型（可能是float、double或long double）。设为I整数类型。

假设两者F和I都有有效的专业化std::numeric_limits<>。

给定一个可表示的值F，并且只使用 C++03，我怎样才能找到最接近的可表示值I？

我追求的是一种纯粹、高效和线程安全的解决方案，除了 C++03 所保证的之外，它对平台没有任何假设。

如果不存在这样的解决方案，是否可以使用 C99/C++11 的新特性找到一个解决方案？

由于报告域错误的方式不平凡，使用lround()C99 似乎存在问题。能否以可移植和线程安全的方式捕获这些域错误？

注意：我知道 Boost 可能通过它的boost::numerics::converter<>模板提供了一个解决方案，但是由于它的高度复杂性和冗长，我无法从中提取要领，因此我无法检查他们的解决方案是否使超出 C++03 的假设。

以下天真的方法失败了，因为I(f)当 的整数部分f不是 的可表示值时，C++03 未定义 的结果I。

然后考虑以下方法：

F(std::numeric_limits<I>::min())这也失败了，因为和的组成部分F(std::numeric_limits<I>::max())可能仍然无法在 中表示I。

最后考虑这第三种方法，它也失败了：

这个时间I(f)总是会有一个明确定义的结果，但是，因为F(std::numeric_limits<I>::max())可能比 小得多std::numeric_limits<I>::max()，所以我们可能会返回std::numeric_limits<I>::max()一个浮点值，它是低于 的多个整数值std::numeric_limits<I>::max()。

请注意，之所以会出现所有麻烦，是因为未定义转换是F(i)向上舍入还是向下舍入到最接近的可表示浮点值。

以下是 C++03（4.9 浮点整数转换）的相关部分：

整数类型或枚举类型的右值可以转换为浮点类型的右值。如果可能，结果是准确的。否则，它是下一个较低或较高可表示值的实现定义选择。

我有很多数据点，我想计算滑动窗口曲线下的面积。但它应该很快。我用谷歌搜索了一下，在 Java 中找到了一个 NewtonCotes 实现，但我不知道是否有更快的方法。

有任何想法吗？

我刚刚定义了一个函数，它是给定阶跃函数的积分，现在我想对获得的三角形函数进行积分。 jerk:1,-1,1-1 步进函数对应于给定的连续 jerk 列表。为了获得三角形函数，

我使用了这个公式：a(t)=kt+ack 是 Jerk，ac 是一个常数，它确保了每一步的连续性和变化。但是要知道，如果我使用公式v(t)=(1/2)kt^2+ac*t+vc，我会得到不正确的速度负值。

我找到了一些 scipy 函数，它返回一个整数的数字，但我需要一条曲线，你知道如何整合它吗？谢谢你

伙计们，我需要有人解决这个问题吗？当我编译该代码时，出现此错误：

我在这部分有问题：

那我该如何解决

我正在计算电容器电压和电流。现在我也想确定能量。能量只是力量的积分，但是我无法整合我的幂函数：

我收到错误：被积函数必须返回与输入向量长度相同的输出向量。

任何人都可以帮忙吗？