黄金分割法优化子程序

风花雪月

本程序包含4个C文件
mhjfgf.c
funct.c
jtf.c
hjfgf.c


本程序由heizi友情提供

风花雪月

mhjfgf.c代码如下：

1. #include "hjfgf.c"
   
2. main()
   
3. {double xx[1],a[1],b[1],ff;
   
4. double s[]={1},x0[]={0};
   
5. jtf(x0,0.1,s,1,a,b);
   
6. ff=gold(a,b,0.00001,1,xx);
   
7. printf("\nx[0]=%f,,ff=%f",xx[0],ff);
   
8. getch();
   
9. }

复制代码

风花雪月

funct.c代码如下：

1. 
   
2. #include "stdio.h"
   
3. #include "stdlib.h"
   
4. #include "math.h"
   
5. double objf(double x[])
   
6. {double ff;
   
7. ff=8*pow(x[0],3)-2*x[0]*x[0]-7*x[0]+3;
   
8. return(ff);
   
9. }
   

复制代码

jtf.c代码如下：

1. 
   
2. #include "funct.c"
   
3. void jtf(double x0[],double h0,double s[],int n,double a[],double b[])
   
4. {int i;
   
5. double *x[3],h,f1,f2,f3;
   
6. for(i=0;i<3;i++)
   
7.   x[i]=(double *)malloc(n*sizeof(double));
   
8.   h=h0;
   
9. for(i=0;i<n;i++)
   
10.   *(x[0]+i)=x0[i];
    
11. f1=objf(x[0]);
    
12. for(i=0;i<n;i++)
    
13.   *(x[1]+i)=*(x[0]+i)+h*s[i];
    
14.   f2=objf(x[1]);
    
15. if(f2>=f1)
    
16.   {h=-h0;
    
17.     for(i=0;i<n;i++)
    
18.     *(x[2]+i)=*(x[0]+i);
    
19.    f3=f1;
    
20.     for(i=0;i<n;i++)
    
21.     {*(x[0]+i)=*(x[1]+i);
    
22.      *(x[1]+i)=*(x[2]+i);
    
23.     }
    
24.    f1=f2;
    
25.    f2=f3;
    
26.    }
    
27.    for(;;)
    
28.    {h=2*h;
    
29.      for(i=0;i<n;i++)
    
30.      *(x[2]+i)=*(x[1]+i)+h*s[i];
    
31.    f3=objf(x[2]);
    
32.    if(f2<f3) break;
    
33.    else
    
34.     { for(i=0;i<n;i++)
    
35.        {*(x[0]+i)=*(x[1]+i);
    
36.         *(x[1]+i)=*(x[2]+i);
    
37.        }
    
38.       f1=f2;
    
39.       f2=f3;
    
40.     }
    
41.    }
    
42.    if(h<0)
    
43.     for(i=0;i<n;i++)
    
44.     {a[i]=*(x[2]+i);
    
45.      b[i]=*(x[0]+i);
    
46.     }
    
47.    else
    
48.     for(i=0;i<n;i++)
    
49.     {a[i]=*(x[0]+i);
    
50.      b[i]=*(x[2]+i);
    
51.      }
    
52.    for(i=0;i<3;i++)
    
53.    free(x[i]);
    
54. }
    

复制代码

hjfgf.c代码如下：

1. 
   
2. #include "jtf.c"
   
3. double gold(double a[],double b[],double eps,int n,double xx[])
   
4. {int i;
   
5. double f1,f2,*x[2],ff,q,w;
   
6. for(i=0;i<2;i++)
   
7.   x[i]=(double *)malloc(n*sizeof(double));
   
8. for(i=0;i<n;i++)
   
9.   {*(x[0]+i)=a[i]+0.618*(b[i]-a[i]);
   
10.    *(x[1]+i)=a[i]+0.382*(b[i]-a[i]);
    
11.   }
    
12.   f1=objf(x[0]);
    
13.   f2=objf(x[1]);
    
14.   do
    
15.    {if(f1>f2)
    
16.      {for(i=0;i<n;i++)
    
17.       {b[i]=*(x[0]+i);
    
18.        *(x[0]+i)=*(x[1]+i);
    
19.        }
    
20.      f1=f2;
    
21.      for(i=0;i<n;i++)
    
22.       *(x[1]+i)=a[i]+0.382*(b[i]-a[i]);
    
23.      f2=objf(x[1]);
    
24.      }
    
25.     else
    
26.      { for(i=0;i<n;i++)
    
27.        {a[i]=*(x[1]+i);
    
28.        *(x[1]+i)=*(x[0]+i);}
    
29.      f2=f1;
    
30.     for(i=0;i<n;i++)
    
31.      *(x[0]+i)=a[i]+0.618*(b[i]-a[i]);
    
32.     f1=objf(x[0]);
    
33.      }
    
34.   q=0;
    
35.   for(i=0;i<n;i++)
    
36.    q=q+(b[i]-a[i])*(b[i]-a[i]);
    
37.   w=sqrt(q);
    
38.   }while(w>eps);
    
39.   for(i=0;i<n;i++)
    
40.    xx[i]=0.5*(a[i]+b[i]);
    
41.   ff=objf(xx);
    
42.   for(i=0;i<2;i++)
    
43.   free(x[i]);
    
44.   return(ff);
    
45. }
    

复制代码

gmda191

如题：好帖，感谢，顶

sunshinexj

这也行 不错 学习了 谢过

5526232

这个怎么组合啊
我怎出现: fatal error C1083: Cannot open include file: 'hjfgf.c': No such file or directory