用0.618法求函数f(x)=x^2+2*x的最小值,-3<=x<=5。是算法有错吗?求解决!!
测试数据与正确数据不符,什么原因呢?测试数据:输入 0.176 输出 -1.11139 -0.940987,但正确输出为-1.112 -0.936,麻烦高手看下,是算法有错吗?
源代码如下:
//用0.618法求minf(x)=x^2+2*x,-3<=x<=5
# include <iostream>
using namespace std;
float fx(float x)
{
float f;
f=x*x+2*x;
return f;
}
void main()
{
float a=-3,b=5,f1,f2,e,p,v,u;
p=0.618;
cout<<"用0.618法求minf(x)=x^2+2*x,-3<=x<=5,请输入精度e:";
cin>>e;
v=a+(1-p)*(b-a);
u=a+p*(b-a);
for(a=-3,b=5;b-a>e;)
{
f1=fx(v);
f2=fx(u);
if(f1>=f2)
{
a=v;
v=u;
f1=f2;
u=a+p*(b-a);
f2=fx(u);
}
else
{
b=u;
u=v;
f2=f1;
v=a+(1-p)*(b-a);
f1=fx(v);
}
}
cout<<a<<" "<<b<<endl;
}