博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU.2899.Strange fuction(牛顿迭代)
阅读量:4591 次
发布时间:2019-06-09

本文共 1305 字,大约阅读时间需要 4 分钟。

\(Description\)

  求函数\(F(x)=6\times x^7+8\times x^6+7\times x^3+5\times x^2-y\times x\)\(x\in \left[0,100\right]\)时的最小值。

\(Solution\)

  \(x\geq 0\)\(F(x)\)为单峰凹函数,三分即可。

  而且由此可知\(F(x)\)的导数应是单增的。函数最值可以转化为求导数零点问题,于是也可以二分求\(F'(x)\)的零点,或者用牛顿迭代求。
  峰值函数最值也可以用模拟退火求。
  练习下牛顿迭代。其它代码可以。

  牛顿迭代:\[x=x_0-\frac{F(x_0)}{F'(x_0)}\]

  对\(F(x)\)泰勒展开,\(F(x)=F(x_0)+F'(x_0)(x-x_0)+\frac{F''(x_0)}{2!}(x-x_0)^2+\ldots+\frac{F^{(n)}(x_0)}{n!}(x-x_0)^n+R_n(x)\)
  为方便计算?只保留线性部分\(F(x)=F(x_0)+F'(x_0)(x-x_0)\),令其等于\(0\)
  就可以得到\(x=x_0-\frac{F(x_0)}{F'(x_0)}\)
  多次迭代、多次选取\(x_0\)即可。

//0MS   1628K#include 
#include
#include
#define eps (1e-7)double y;inline double f(double x){ return 6*pow(x,7)+8*pow(x,6)+7*pow(x,3)+5*x*x-y*x;}inline double fd(double x){ return 42*pow(x,6)+48*pow(x,5)+21*x*x+10*x-y;}inline double fdd(double x){ return 252*pow(x,5)+240*pow(x,4)+42*x+10;}double Get_zero(double x)//求导函数零点 { double las=x+1; while(fabs(las-x)>eps) las=x, x=x-fd(x)/fdd(x); return x;}int main(){ int T; scanf("%d",&T); while(T--) { scanf("%lf",&y); double ans=1e15; for(int i=0; i<=100; i+=10) ans=std::min(ans,f(Get_zero(i))); printf("%.4lf\n",ans); } return 0;}

转载于:https://www.cnblogs.com/SovietPower/p/9162128.html

你可能感兴趣的文章
checkbox和文字对齐
查看>>
JConsole远程连接配置 服务器监控工具
查看>>
了解HTTP协议栈(实践篇)
查看>>
loj10035. 「一本通 2.1 练习 1」Power Strings
查看>>
%s的用法
查看>>
调用底层不能直接访问的类和方法
查看>>
清理缓存的方法 #DF
查看>>
JAVA array,map 转 json 字符串
查看>>
2017-12-27练习
查看>>
NET设计规范(二) 命名规范
查看>>
VMware 9.0.1安装Mac OS X Mountain Lion 10.8.2
查看>>
SSL延迟
查看>>
android新手关于左右滑动的问题,布局把<android.support.v4.view.ViewPager/><ImageView/> 放在上面就不行了。...
查看>>
深入理解DIP、IoC、DI以及IoC容器
查看>>
赋值文件
查看>>
Vue 数组 字典 template v-for 的使用
查看>>
蓝牙模块选择经验谈
查看>>
java中==和equals
查看>>
CCActionPageTurn3D
查看>>
python random
查看>>