在R中,我们可以将停止准则作为一个输入参数来添加到函数中。例如,假设我们有一个找到方程解的函数,我们想让它在满足某个条件时停止,我们可以通过添加一个停止准则来实现。
以下是一个示例函数,该函数将一个方程的解作为输入,并返回它的根,直到满足停止准则:
find_root <- function(f, lower, upper, stopping_criterion=function(x) abs(f(x)) < 1e-06) {
if (stopping_criterion(lower) || stopping_criterion(upper)) {
return(NULL) # no root in the range
}
mid <- (lower + upper) / 2
while (!stopping_criterion(mid) && (upper - lower) > 1e-06) {
if (f(mid) * f(lower) < 0) {
upper <- mid
} else {
lower <- mid
}
mid <- (lower + upper) / 2
}
return(mid)
}
在这个函数中,我们使用了一个匿名函数来定义停止准则,如果没有指定任何停止准则,则默认使用 $|f(x)|<1e^{-6}$ 作为停止准则。您可以通过向函数传递另一个函数来更改默认停止准则。
下面是一个使用该函数的示例:
f <- function(x) x^2 - 2
find_root(f, 0, 2)
# output: 1.414214
stopping_criterion <- function(x) abs(f(x)) < 1e-03
find_root(f, 0, 2, stopping_criterion)
# output: 1.414176