法和最速下降法的C语言程序(共2页).doc

《法和最速下降法的C语言程序(共2页).doc》由会员分享，可在线阅读，更多相关《法和最速下降法的C语言程序(共2页).doc（2页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上0618法的c语言程序#include #include float m=0.618; float fun(float t) 定义目标函数float y; y=3*t*t-21.6*t-1; return y;main() float a=0,b=10,e; 定义域从a=0开始，到b=10；精度为e； float t1,t2,t,f1,f2,min; 定义变量：按照算法要求定义变量t1为t2，并且目标 printf(njingdu is e:t); 在t1，t2点处的函数值为f1，f2。 scanf(%f,&e); while(b-ae) 每次都判断b-a与e的关系，

2、如果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 如果f1e，将(a+b)/2赋给t，即为 min=fun(t); 所求的最小点的因变量。并求出在这点的函数值min printf(t=%f.,t); 将t和min输出。 printf(The minimum is %f.,min);最速下降

3、法的C语言程序#include (stdio.h)#include (math.h)float fun1(float x1,float x2) 定义函数fun1为目标函数对x1求偏导float a; a=6*x1+2*x2+2; return a;float fun2(float x1,flaot x2) 定义函数fun2为目标函数对x2求偏导float b; b=2*x1+2*x2-2; return b;float fun3(float x1,flaot x2) 定义函数fun3为目标函数flaot y; y=3*x1*x1+2*x1*x2+x2*x2+2*x1-2*x2+1; retur

4、n y;main() flaot t, e=0.01, x1=0, x2=1, a, b, y, m; 定义起始点为x1=0，x2=1，并定义精度为e=0.01 a=fun1(x1,x2); 把fun1的值赋给a b=fun2(x1,x2); 把fun2的值赋给b m=sqrt(a*a+b*b); 对a*a+b*b求开方，将值赋给m while(me) 判断，当me时进行以下循环 t=(a*a+b*b)/(8*a*a+2*b*b+4*a*b); 根据梯度法（最速下降法），利用梯度和海赛矩阵 x1=x1-a*t; 求步长t。 根据步长和梯度方向求出新的x1，x2 x2=x2-b*t; a=fun1(x1,x2); 根据得到的x1，x2求出新的fun1，fun2，并将值 b=fun2(x1,x2); 赋给a，b m=sqrt(a*a+b*b); 计算新的m y=fun3(x1,x2); 当m不满足me的时候退出循环，并计算fun3， printf(The min is %f,y); 将值赋给y,并输出。专心-专注-专业