黄金分割算法求函数的极值C++实现
生活随笔
收集整理的這篇文章主要介紹了
黄金分割算法求函数的极值C++实现
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
黃金分割算法求極值迭代的最終判斷是兩次函數(shù)值的結(jié)果差是否很小,如果小,則說明該處函數(shù)的斜率很小,就是一個極值點(diǎn)。
過程如下:
?
以下是代碼:
//黃金分割法! /* ================================================================ 題目:用黃金分割法求解3*x*x*x-2*x*x-16=0的根。 ================================================================ */ #include<stdio.h> #include <math.h> #define E 0.5 //精度 //方程的解,可以通過修改方程求解不同的問題 double hs(double x) {return x*x-x+2; } double gen(double a, double b) //黃金分割算法求解函數(shù),a為左區(qū)間,b為右區(qū)間 {double g,a1,b1;g = (sqrt(5) - 1.0) / 2; //黃金分割的值a1 = a + (1 - g)*(b - a);b1 = a + g * (b - a);while (fabs(b-a) > E){if (hs(a1)>hs(b1)){a = a1;a1 = b1;b1= a + g * (b - a);}else {b = b1;b1 = a1;a1= a + (1 - g)*(b - a);}}return (a+b)/2; } int main() {double a, b, t, fa, fb;printf("輸入估計(jì)區(qū)間: a="); //左區(qū)間scanf("%lf", &a);printf("輸入估計(jì)區(qū)間: b="); //有區(qū)間scanf("%lf", &b);t = gen(a, b);printf("解得:(%.4lf,%0.4lf) ", t,hs(t));return 0; }結(jié)果:
?
?
?
?
《新程序員》:云原生和全面數(shù)字化實(shí)踐50位技術(shù)專家共同創(chuàng)作,文字、視頻、音頻交互閱讀總結(jié)
以上是生活随笔為你收集整理的黄金分割算法求函数的极值C++实现的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: Yen 的k_shortest pat
- 下一篇: css,css,javascript实现