机械优化设计C语言程序.doc

PAGE / NUMPAGES 目 录 黄金分割法 二次插值法 最速下降法（阶梯法） 变尺度法 鲍威尔法 一、黄金分割法 #includestdio.h #includemath.h #define r 0.618 #define f(x) x*x+2*x golden(float,float,float); main() { float a,b,e; printf(\n请输入区间和收敛精度:a,b,e\n); scanf(%f,%f,%f,a,b,e); golden(a,b,e); } golden(float a,float b,float e) { float y1,y2,a1,a2,A,Y; int n=0; a1=b-r*(b-a); a2=a+r*(b-a); y1=f(a1); y2=f(a2); printf(黄金分割法的搜索过程:); do { printf(\n %d a=%f,b=%f,a1=%f,a2=%f,y1=%f,y2=%f,n,a,b,a1,a2,y1,y2); if(y1=y2) { a=a1; a1=a2; y1=y2; a2=a+r*(b-a); y2=f(a2); } else { b=a2; a2=a1; y2=y1; a1=b-r*(b-a); y1=f(a1); } n++; }while(fabs((b-a)/b)=e||fabs((y2-y1)/y2)=e); A=(a+b)/2; Y=f(A); printf(\n %d a=%f,b=%f,a1=%f,a2=%f,n,a,b,a1,a2); printf(\n结果:\n极值点及其函数值:A=%f,Y=%f\n, A,Y); } 二、二次插值法 #includestdio.h #includemath.h #includeconio.h void main() { float*area(float a1,float p,float a[3]); float f(float x); float ar,fr; float a1=10,p=0.01,e=0.000001; float pa[3]; area(a1,p,pa); a1=pa[0]; float a2=pa[1]; float a3=pa[2]; float f1=f(a1); float f2=f(a2); float f3=f(a3); do { ar=((a3*a3-a2*a2)*f1+(a1*a1-a3*a3)*f2+(a2*a2-a1*a1)*f3); ar=ar/2/((a3-a2)*f1+(a1-a3)*f2+(a2-a1)*f3); fr=f(ar); if(ara2) { if(frf2) { a3=ar; f3=fr; } else if(frf2) { a1=a2; f1=f2; a2=ar; f2=fr; } else { a3=ar; a1=a2; a2=(a1+a3)/2; f1=f2; f3=fr; f2=f(a2); } } else if(ara2) { if(frf2) { a1=ar; f1=fr; } else if(frf2) { a3=a2; f3=f2; a2=ar; f2=fr; } else { a1=ar; a3=a2; a2=(a1+a3)/2; f1=fr; f3=f2; f2=f(a2); } } if (fabs(a1-a3)=e) break; } while(1); if(f2fr) { ar=a2; fr=f2; } printf(\nx=%f\nf(x*)=%f,ar,fr); } float*area(float a1,float p,float a[3]) { float f(float x); float a2,f2,a3,f3,temp; float acc=0.00001; float f1=f(a1); float storep=p; while(1) { a2=a1+p; f2=f(a2); if(f2=f1) { if(fa