c语言坐标轮换法_优化设计-鲍威尔法程序(c语言)

#include

#include

#define m 10 /*数组长度m >= 维数n */

float f(float x[]);

void mjtf(int n,float x0[],float h,float s[],float a[],float b[]);

void mhjfgf(int n,float a[],float b[],float flag,float x[]);

void mbwef(int n,float x0[],float h,float flag,float a[],float b[],float x[]);

float f(float x[])

{

float result;

result=60-10*x[0]-4*x[1]+x[0]*x[0]+x[1]*x[1]-x[0]*x[1];

return result;

}

/*多维进退法子程序*/

void mjtf(int n,float x0[],float h,float s[],float a[],float b[])

{

int i;

float x1[m],x2[m],x3[m],f1,f2,f3;

for(i=0;i

{

x1[i]=x0[i];

x2[i]=x0[i]+h*s[i];

}

f1=f(x1);

f2=f(x2);

if(f2>=f1) /*判断搜索方向*/

{ /*搜索方向为反向，转身*/

h=(-1)*h;

for(i=0;i

x3[i]=x1[i];

f3=f1;

for(i=0;i

x1[i]=x2[i];

f1=f2;

for(i=0;i

x2[i]=x3[i];

f2=f3;

} /*搜索方向为正向*/

for(i=0;i

x3[i]=x2[i]+h*s[i];

f3=f(x3);