R：求多项式方程的根

Question

我在 R 中有这个方程。

f <- function(x) {first +second*x +third*x^2  +fourth*filter_factor - log(myBITRATE)}

在哪里

first= -5.219078
second = 0.7613156
third = -0.01298033
fourth = -0.05218249
filter_factor = 1
myBITRATE = 184.47

有没有办法找到这个方程的根？

我需要一个 newton-raphson 方法的起点。

Answer 1

您的函数是具有一个变量的二阶多项式，因此您可以进行简单的计算：

delta = second^2 - 4*third*(first + fourth*filter_factor - log(myBITRATE))

由于您的增量是严格正数：

x1 = (- second - sqrt(delta))/(2*third)
x2 = (- second + sqrt(delta))/(2*third)

#> x1
#[1] 36.53336
#> f(x1)
#[1] 0
#> x2
#[1] 22.11812
#> f(x2)
#[1] 8.881784e-16