大佬教程收集整理的这篇文章主要介绍了为多项式函数查找 x 的最大值,大佬教程大佬觉得挺不错的,现在分享给大家,也给大家做个参考。
我正在使用一个简单的多项式来拟合曲线。
poly <- function(a,b,c,X) a * x^2 + b * x + c
我想找到导致曲线最大值的 x 值。目前,我创建了一个 x 范围从 20000 到 50000 的网格,为每一行运行该函数,然后对结果使用 @H_516_5@max()。它有效,但我有很多组,每次我这样做时它都会创建一个大数据框。它非常笨重,我觉得必须有更好的方法。
一些典型的系数是:
a <- -0.000000179
b <- 0.011153167
c <- 9.896420781
如果你重新排列你的函数,那么你想要最大化的变量是第一个,你像这样设置默认值:
poly <- function(x,a,b,C) a * x^2 + b * x + c
formals(poly)$a <- -0.000000179
formals(poly)$b <- 0.011153167
formals(poly)$c <- 9.896420781
然后您可以使用 optimize
函数在您的时间间隔内最大化:
optimize(poly,c(20000,50000),maximum = T)
$`maximum`
[1] 31154.1
$objective
[1] 183.6298
其中 $maximum
是出现最大值的 x
值,$objective
是高度。
您可以使用优化。我认为该线程中回答的其他解决方案更有吸引力,但为了完整性,我会写下来:
a <- -0.000000179
b <- 0.011153167
c <- 9.896420781
o <- optim(
par=list(x=0),fn=function(X){ -poly(a,c,x=X) },method="Brent",lower=-50e3,upper=50e3
)
输出:
> o
$par
[1] 31154.1
$value
[1] -183.6298
$counts
function gradient
NA NA
$convergence
[1] 0
$message
NULL
,
如果 a
为负,则在 a * x^2 + b * x + c
处达到抛物线 -b/(2*a)
的最大值:
a<0
#[1] TRUE
-b/(2*a)
#[1] 31154.1
以上是大佬教程为你收集整理的为多项式函数查找 x 的最大值全部内容,希望文章能够帮你解决为多项式函数查找 x 的最大值所遇到的程序开发问题。
如果觉得大佬教程网站内容还不错,欢迎将大佬教程推荐给程序员好友。
本图文内容来源于网友网络收集整理提供,作为学习参考使用,版权属于原作者。
如您有任何意见或建议可联系处理。小编QQ:384754419,请注明来意。