请各位大侠帮看看小弟的程序.........这个程序有错吗?除了结构之外.... - C语言论坛

| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛 |

编程论坛 → 开发语言 → 『 C语言论坛』 → 请各位大侠帮看看小弟的程序.........这个程序有错吗?除了结构之外....

我的收件箱(0)

共有 611 人关注过本帖

标题：请各位大侠帮看看小弟的程序.........这个程序有错吗?除了结构之外....

只看楼主加入收藏

w527705090

等　级：贵宾
威　望：11
帖　子：441
专家分：1882
注　册：2011-6-28
结帖率：100%

楼主

问题点数：0 回复次数：6

请各位大侠帮看看小弟的程序.........这个程序有错吗?除了结构之外....

程序代码：

#include<stdio.h>
#include<math.h>
#define f(x) x*x+2*x+1 //在这里定义函数
main()
{
    float a,b,x1,x2,f1,f2,e;
    int i;
    printf("intput a b e:");
    scanf("%f%f%f",&a,&b,&e);//输入区间和精度
    for(i=0;;i++)
    {
        if(fabs(a-b)<=e)
       
        {
            printf("%f\n%d\n",f(0.5*(a+b)),i);
            break;//符合这个条件后就结束循环
        }
        else
        {
            x1=a+0.382*(b-a);
            x2=a+0.618*(b-a);
            f1=f(x1);
            f2=f(x2);
            if(f(x1)<f(x2))
      
            {
                b=x2;
           
            }
            else
            {
                a=x1;
            }
       
            continue;//这个用来再次循环
        }
    }
}

小弟是想用黄金分割法求极值.....不知这样编程有没有不对的地方....希望各位给点意见.....谢谢~~~
调试了没有提示错误.....但是感觉不对.....

[ 本帖最后由 w527705090 于 2011-9-22 22:36 编辑 ]

搜索更多相关主题的帖子: 成都　

2011-09-22 22:33

w527705090

等　级：贵宾
威　望：11
帖　子：441
专家分：1882
注　册：2011-6-28

第 2 楼

得分:0

先自己顶起~~~

有心者，千方百计；无心者，千难万难。

2011-09-22 22:34

w527705090

等　级：贵宾
威　望：11
帖　子：441
专家分：1882
注　册：2011-6-28

第 3 楼

得分:0

小弟自己找出原因了....

程序代码：

#include<stdio.h>
#include<math.h>
#define f(x) x*x+2*x+1 //在这里定义函数
main()
{
    double a,b,x1,x2,f1,f2,e; //把这些字符类型改成双精度
    int i;
    printf("intput a b e:");
    scanf("%lf%lf%lf",&a,&b,&e);//输入区间和精度 这里同样要改一下
    for(i=0;;i++)
    {
        if(fabs(a-b)<=e)
      
        {
            printf("%f\n%d\n",f((a+b)/2),i);
            break;//符合这个条件后就结束循环
        }
        else
        {
            x1=a+0.382*(b-a);
            x2=a+0.618*(b-a);
            f1=f(x1);
            f2=f(x2);
            if(f(x1)<f(x2))
     
            {
                b=x2;
          
            }
            else
            {
                a=x1;
            }
      
            continue;//这个用来再次循环
        }
    }
}