注册 登录
编程论坛 C++教室

用二分法求根:陷入了死循环(已解决)

yqiong 发布于 2008-03-20 17:18, 2188 次点击
用二分法求根:陷入了死循环
#include<iostream>
using namespace std;
int main()
{
    double fabs(double  ,double);
    double f(double );
       double x2,x1,x;
    do{
             cout<<"请输入x1和x2的值:";
            cin>>x1>>x2;
            cout<<"x1="<<x1<<"  x2="<<x2;
    }while(f(x1)*f(x2)>0);

    do{
        x=(x1+x2)/2;
        if(f(x)*f(x1)>0) x2=x;
        else x1=x;
    }while(fabs(f(x),0)>=1e-5);

    cout<<"结果为:"<<x;
    return 0;
}

double f(double x)
{      return(((x+-5)*x+16)*x+-80);}

double fabs(double a,double b)
{    return (a-b>0?a-b:b-a);}

[[it] 本帖最后由 yqiong 于 2008-3-21 15:41 编辑 [/it]]
7 回复
#2
PcrazyC2008-03-20 18:09
二分法有一个严格的限制就是要在一个单调区间上,如果你输入的数不在一个单调区间就会出现死循环,你这个题中明显有一个答案是5,所以我输入4,6就行了,其它的,自己再找
#3
yqiong2008-03-20 22:50
[bo]以下是引用 [un]PcrazyC[/un] 在 2008-3-20 18:09 的发言:[/bo]

二分法有一个严格的限制就是要在一个单调区间上,如果你输入的数不在一个单调区间就会出现死循环,你这个题中明显有一个答案是5,所以我输入4,6就行了,其它的,自己再找 ...


答案是5,但我输入2 6时怎么也陷入死循环了呀?
#4
eagleboycn2008-03-21 14:43
LZ,输入2,6是f(2)<0,f(6)>0,而且f((2+6)/2)<0,则f(2)*f(4) 负负得正
我改了两个while判断语句
#include<iostream>
using namespace std;
int main()
{
    double fabs(double  ,double);
    double f(double );
    double x2,x1,x;
    do{
       cout<<"请输入x1和x2的值:";
       cin>>x1>>x2;
       cout<<"x1="<<x1<<"  x2="<<x2;
    }while(f(x1)>0&&f(x2)>0);/////////////////////////////////////////

    do{
        x=(x1+x2)/2;
        if(f(x)>0&&f(x1)>0) ////////////////////////////////////
            x2=x;
        else x1=x;
    }while(fabs(f(x),0)>=1e-5);

    cout<<"结果为:"<<x<<endl;
    return 0;
}

double f(double x)
{      
    cout<<"x:"<<x<<"  ";
    cout<<"f(x)"<<(((x-5)*x+16)*x-80)<<endl;
    return(((x-5)*x+16)*x-80);
}

double fabs(double a,double b)
{    return (a>b?(a-b):(b-a));}
#5
yqiong2008-03-21 15:33
do{
       cout<<"请输入x1和x2的值:";
       cin>>x1>>x2;
       cout<<"x1="<<x1<<"  x2="<<x2;
    }while(f(x1)>0&&f(x2)>0)//不明白。只有当异号时才能确保在x1和x2间有一值使函数值为0呀!
    do{
        x=(x1+x2)/2;
        if(f(x)>0&&f(x1)>0)  x2=x;
        else x1=x;
    }while(fabs(f(x),0)>=1e-5);

运行时输入2 6确实能正确输出答案5,但改输入3 6又出问题了呢!
#6
yqiong2008-03-21 15:40
do{
        x=(x1+x2)/2;
        if(f(x)*f(x1)>0) x1=x;     
       else x2=x;  //找到原因了,就错在这!
  }while(fabs(f(x),0)>=1e-5);
#7
newyj2008-03-21 19:55
double fabs(double  ,double);
    double f(double );
放在主函数外和 放在里面有什么区别
#8
yqiong2008-03-24 12:10
没区别吧,放在里面的是旧版的。两种形式C都支持
1