优化设计-梯度法

《优化设计-梯度法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、优化设计作业梯度法用梯度法求法min（x12+2x22-4x1-2x1x2），s.t.x1+x2≤6。程序如下：#include"stdio.h"#include"math.h"voidmain(){floatx1,x2,x3,x4,a,b,c,d,y,y1,y2,y3,y4;x1=1,x2=1;do{y1=2*x1-4-2*x2;y2=4*x2-2*x1;x3=x1-a*y1;x4=x2-a*y2;b=y1*y1+2*y2*y2-2*y1*y2;c=-y1*x1-x1*y1-y2*x2*4+4*y1+2*x1*y2+2*y1*x2;a

2、=-c/(2*b);d=x3+x4;y3=x1*x1+2*x2*x2-4*x1-2*x1*x2;y4=x3*x3+2*x4*x4-4*x3-2*x3*x4;if(d<6&&y3>y4){a=a*2;x3=x1-a*y1;x4=x2-a*y2;}else{a=0.5*a;x3=x1-a*y1;x4=x2-a*y2;}printf("x1=%f,x2=%f,a=%f",x1,x2,a);x1=x3;x2=x4;}while(sqrt(y1*y1+y2*y2)>0.05);y=x1*x1+2*x2*x2-4*x1-2*x1*x2;prin

3、tf("minf(x)=%f",y);运行结果：