机械优化设计实验报告 - 图文 

}

*(xx[0]+i)=*(xx[3]+i);

else

{for (i=0;i

{*(ss+(i+1)*(n+1))=x[i]-(*(xx[1]+i)); *(s+i)=*(ss+(i+1)*(n+1)); }

f=oneoptim(xx[0],s,h0,epsg,n,x);

for(i=0;i

*(xx[0]+i)=x[i];

for (j=m+1;j<=n;j++) for (i=0;i

*(ss+i*(n+1)+j-1)=*(ss+i*(n+1)+j); } } } main() { }

double p[]={1,1}; double ff,x[2];

ff=powell(p,0.3,0.001,0.0001,2,x); printf(\输出最优点及其目标函数值：\\n\printf(\

五、计算结果及分析

由图可知程序求出来的最优点为x*=(2.500861 2.499306)T,目标函数值为0.000003.可见已经足够接近极值点x*=(2.5 2.5)T及极小值法f(x*) = 0。程序满足要求。此方法是一种有效的共轭方向法。