共轭梯度法

《共轭梯度法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、用C语言编写的共轭梯度法计算程序如下：#include“math.h”#include“stdio.h”floatx[10],y[10],p[10],f,h;intn;voidfun(){inti;for(I=1;I

2、(“%d,%d”,&n,&e”);for(I=1;I

3、;h=h;fun();f1=f;p2:u=0;t=t0;h2=h1+t;h=h2;fun();f2=f;if(f1>f2){t=t+t;u=u+1;}else{t=-t;h3=h1;f3=f1;h1=h2;f1=f2;h2=h3;f2=f3;}p1:h3=h2+t;=h3;fun();f3=f;if(f2>f3){t=t+t;u=u+1;h1=h2;f1=f2;h2=h3;f2=f3;gotop1;}else{if(u>0){h4=0.5*(h2+h3);h=h4;fun();f4=f;if(f4>f2){h3=h4;f3=f4;}else{h1=h2;f1=f2;h2=h4;f2=

4、f4;}}c1=(f3-f1)/(h3-h1);c2=((f2-f1)/(h2-h1)-c1)/(h2-h3);if((f2<1)f5=1;elsef5=f2;if((fabs(f4-f2)/f5)f2){h1=h2;f1=f2;}else{h1=h4;f1=f4;}t0=v*t0;gotop2;}}}}}gotop4;p3:h=0;fun();printf(“OBJ.FUNCF=%f”,f);for(I=1;I

5、f(“=%f”,x[I]);}}