我试图在矩形区域上集成解析函数(sqrt和三角函数的复合函数)。它在区域内没有奇点,似乎是使用dblquad
的完美候选者。我的问题是如何评估Matlab提供给我的数值的准确性。在不知道积分的确切值的情况下,我们如何证明有效数字的合理性?当您需要给出具有特定精度数字的值时,您应该能够进行校验。是否有可能实现这一给定的值是通过使用Matlab计算?
除非您另行设置,否则dblquad
对绝对正交误差使用默认容差阈值(最新版本中为10-6)。积分的近似值将在不大于规定公差的误差范围内。
你可以看看dblquad
的源代码,在某个地方它将使用一定数量的"步骤"。我猜你可以做一个新的m文件,用重要的位来计算积分,然后调整步数,直到它花了计算机很长时间,而不改变结果。就我个人而言,我使用自定义的辛普森规则进行数值积分,只是将N(步数)更改为一些较大的数字。