程序问答   发布时间:2022-06-02  发布网站:大佬教程  code.js-code.com
大佬教程收集整理的这篇文章主要介绍了为多项式函数查找 x 的最大值大佬教程大佬觉得挺不错的,现在分享给大家,也给大家做个参考。

如何解决为多项式函数查找 x 的最大值?

开发过程中遇到为多项式函数查找 x 的最大值的问题如何解决?下面主要结合日常开发的经验,给出你关于为多项式函数查找 x 的最大值的解决方法建议,希望对你解决为多项式函数查找 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,请注明来意。