实施具有明确限制的比率

计算科学 参考请求 稳定
2021-12-28 09:00:17

我正在实现一个函数由两个封闭形式方程的比率给出。 f(x1,x2,xn)

f(x1,x2,,xn)=g(x1,x2,,xn)h(x1,x2,,xn)

不幸的是,当时,都是但是有一个明确的限制 x1=0gh0

k(x2,x3,xn)=limx10f(x1,x2,xn)

实现这一点的明智方法是什么?我在想

if |x_1| > cutoff:
   return f(x_1, ..., x_n)/g(x_1, ..., x_n)
else:
   return k(x_2, ... x_n)

如何确定合适的截止值?

我试图研究这个问题,我遇到了海厄姆关于数值算法的准确性和稳定性的书。

我在第 1.14.1 节中找到了一个相关示例,他在其中讨论了 他对不稳定性的解决方案与我上面的幼稚方法不同;他使用坐标的变化来以某种方式减少不准确性。f(x)=(ex1)/xy=ex

我将不胜感激有关如何解决此问题的任何建议/参考。

1个回答

的幂展开函数并计算误差项。当错误低于您想要的阈值时使用扩展。fx1

线性化您将拥有: 作为最后一项错误选择作为从的绝对值最大的值f

f=f0+fx0+12f(ξ)x02
ϵξf(x0)f(ξ))=0

if epsilon<threshold
    f=f_0+f'x_0
else
    f=g/h