0.618法和最速下降法的C语言程序

0．618法的c语言程序#includestdio.h#includemath.hfloatm=0.618;floatfun(floatt)定义目标函数{floaty;y=3*t*t-21.6*t-1;returny;}main(){floata=0,b=10,e;定义域从a=0开始，到b=10；精度为e；floatt1,t2,t,f1,f2,min;定义变量：按照算法要求定义变量t1为t2，并且目标printf(\njingduise:\t);在t1，t2点处的函数值为f1，f2。scanf(%f,&e);while(b-ae)每次都判断b-a与e的关系，如果b-ae，进行以下循环{t1=a+(1-m)*(b-a);t2=a+m*(b-a);f1=fun(t1);f2=fun(t2);if(f1f2)每次迭代后生成新的f1，f2判断他们的大小，如果f1f2{a=t1;进行以下计算，得到新的、更小的定义区间。t1=t2;f1=f2;t2=a+m*(b-a);f2=fun(t2);}else如果f1f2时，进行一下计算，得到新的、更小的定义区间{b=t2;t2=t1;f2=f1;t1=a+(1-m)*(b-a);f1=fun(t1);}}t=(a+b)/2;直到得到的b-a不再满足b-ae，将(a+b)/2赋给t，即为min=fun(t);所求的最小点的因变量。并求出在这点的函数值minprintf(t=%f.,t);将t和min输出。printf(Theminimumis%f.,min);}最速下降法的C语言程序#include(stdio.h)#include(math.h)floatfun1(floatx1,floatx2)定义函数fun1为目标函数对x1求偏导{floata;a=6*x1+2*x2+2;returna;}floatfun2(floatx1,flaotx2)定义函数fun2为目标函数对x2求偏导{floatb;b=2*x1+2*x2-2;returnb;}floatfun3(floatx1,flaotx2)定义函数fun3为目标函数{flaoty;y=3*x1*x1+2*x1*x2+x2*x2+2*x1-2*x2+1;returny;}main(){flaott,e=0.01,x1=0,x2=1,a,b,y,m;定义起始点为x1=0，x2=1，并定义精度为e=0.01a=fun1(x1,x2);把fun1的值赋给ab=fun2(x1,x2);把fun2的值赋给bm=sqrt(a*a+b*b);对a*a+b*b求开方，将值赋给mwhile(me)判断，当me时进行以下循环{t=(a*a+b*b)/(8*a*a+2*b*b+4*a*b);根据梯度法（最速下降法），利用梯度和海赛矩阵x1=x1-a*t;求步长t。根据步长和梯度方向求出新的x1，x2x2=x2-b*t;a=fun1(x1,x2);根据得到的x1，x2求出新的fun1，fun2，并将值b=fun2(x1,x2);赋给a，bm=sqrt(a*a+b*b);计算新的m}y=fun3(x1,x2);当m不满足me的时候退出循环，并计算fun3，printf(Theminis%f,y);将值赋给y,并输出。}