1、黃金分割法的C語言程序應(yīng)用#include #include #define e 0.001 （給定收斂精度e=0.001）float f(float x)float y=pow(x,2)+2*x; (我們需要求解的目標(biāo)函數(shù)：)return(y);main() （主函數(shù)）float a=-3,b=5; （給定初始單峰區(qū)間 a, b）float x1,x2,f1,f2,c,min;x1=a+0.382*(b-a); （在區(qū)間 a, b內(nèi)取兩個內(nèi)插點(diǎn)x1和x2并計(jì)算其函數(shù)值）x2=a+0.618*(b-a);f1=f(x1);f2=f(x2);printf(n %f %f,a,b); while(

2、b-a)e) （當(dāng)新的b-ae=0.001時，就執(zhí)行下列這些運(yùn)算）if (f1f2)（若 f1 f2 ，則取a,x2為新區(qū)間，而x1則作為新區(qū)間內(nèi)的第一個試算點(diǎn)，即令b=x2,x2=x1,f2=f1,而另一試算點(diǎn)可按x1=a+0.382*(b-a)算出）b=x2; x2=x1;x1=a+0.382*(b-a);f2=f1;f1=f(x1);else （若 f1 f2 ，則取x1,b為新區(qū)間，而x2則作為新區(qū)間內(nèi)的第一個試算點(diǎn)，即令a=x1,x1=x2,f1=f2,而另一試算點(diǎn)可按x2=a+0.618*(b-a)算出）a=x1;x1=x2;x2=a+0.618*(b-a);f1=f2;f2=f(

3、x2);printf(n %f %f,a,b); c=(a+b)/2; min=f(c);printf(n The area is %f to %f,a,b); printf(n The min is %f and the result is %f,c,min);二次插值法的C語言程序應(yīng)用#include#include double f(double x) return x*x*x-3*x+2;void main(void) double a1,a2,a3,a0; double f1,f2,f3,f0; double q=0.; （給定計(jì)算精度q=0.） int k=1; a1=0; (給定

4、初始搜索區(qū)間【a1,a3】) a3=3; a2= 0.5*(a1+a3); (采用等距原則取點(diǎn)) f1=f(a1); （計(jì)算三點(diǎn)的函數(shù)值） f2=f(a2); f3=f(a3); while(1) a0=0.5*(a2*a2-a3*a3)*f1+(a3*a3-a1*a1)*f2+(a1*a1-a2*a2)*f3)/(a2-a3)*f1+(a3-a1)*f2+(a1-a2)*f3); （計(jì)算插值函數(shù)的極小點(diǎn)a0） f0=f(a0); (計(jì)算極小值f0) if(fabs(a0-a2)a2) （用進(jìn)退法確定搜索區(qū)間） if(f0=f2) (進(jìn)一步縮小區(qū)間，以得到新的單峰區(qū)間) a1=a2; a2=a0; f1=f2; f2=f0; k+; else a3=a0; f3=f0; k+; else （這里的否定是表示：如果當(dāng)a0a2） if(f0=f2) a3=a2; a2=a0; f3=f2; f2=f0; k