Question

假设我有

 fu3 = g.^5.*qQd./(exp(g.^2/T(i)));
    I_value = (1/T^3)*trapz(g,fu3);

所以，这取决于T和s（包含在qQd中）。这个计算需要几分钟。所以我想把这个表达式放到一个简单的多项式表达式中。我认为名称适合或插值。我想找到一个满足I_value积分的多项式表达式。有什么帮助吗？

Answer 1

您的整合范围是什么？让我假设它非常大。然后你的被积函数变化了很多个数量级（你是否绘制了它？你可能需要用对数标度来理解问题的全部范围）。这是梯形积分需要几分钟的原因。用多项式代替被积函数不是一个选项。

解决方案1：通过分析解决问题（部分整合，替换g ^ 2-> u）。

解决方案2：如果您想以数字方式进行，请使用比等距梯形更好的集成方案。您需要一个积分器，根据被积函数的幅度变化调整步长。