給定一個數字向量，是否可以在 R 中創建一個用于連分數估計的函式？
公式為：
uj5u.com熱心網友回復：
我們可以定義一個遞回函式f來計算連分數
f <- function(v) {
ifelse(length(v) == 1, v, v[1] 1 / f(v[-1]))
}
然而，缺點是遞回深度有限制，例如，
> f(1:1e5)
Error: node stack overflow
因此，如果您有一個大陣列v，更好的選擇可能是使用for回圈，例如，
f <- function(v) {
if (length(v) == 1) {
return(v)
}
s <- tail(v, 1)
for (k in (length(v) - 1):1) {
s <- v[k] 1 / s
}
s
}
你會看到
> f(1:1e5)
[1] 1.433127
uj5u.com熱心網友回復：
使用遞回一個選項來實作您想要的結果可能如下所示：
cfrac <- function(x) {
n <- length(x)
if (n > 1) {
res <- x[1] 1 / cfrac(x[2:n])
} else {
res <- x
}
return(res)
}
identical(cfrac(2), 2)
#> [1] TRUE
identical(cfrac(c(2, 3)), 2 1 / 3)
#> [1] TRUE
identical(cfrac(c(-3, 2, 18)), -3 1 / (2 1 / 18))
#> [1] TRUE
轉載請註明出處，本文鏈接：https://www.uj5u.com/qiye/388512.html
上一篇：計數事故的發生，直到下一次事故