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

《法和最速下降法的C语言程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、0．618法的c语言程序#include#includefloatm=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("jingduise:t");在t1，t2点处的函数值为f1，f2。sca

2、nf("%f",&e);while(b-a>e)每次都判断b-a与e的关系，如果b-a>e，进行以下循环{t1=a+(1-m)*(b-a);t2=a+m*(b-a);f1=fun(t1);f2=fun(t2);if(f1>f2)每次迭代后生成新的f1，f2判断他们的大小，如果f1>f2{a=t1;进行以下计算，得到新的、更小的定义区间。t1=t2;f1=f2;t2=a+m*(b-a);f2=fun(t2);}else如果f1

3、f1;t1=a+(1-m)*(b-a);f1=fun(t1);}}t=(a+b)/2;直到得到的b-a不再满足b-a>e，将(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为目标函数对

4、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

5、，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(m>e)判断，当m>e时进行以下循环{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

6、，fun2，并将值b=fun2(x1,x2);赋给a，bm=sqrt(a*a+b*b);计算新的m}y=fun3(x1,x2);当m不满足m>e的时候退出循环，并计算fun3，printf("Theminis%f",y);将值赋给y,并输出。}