匿名用户
1级
2012-03-21 回答
这是高次方程是没有解析解的,但可以用逼近法求解。
关于这类题最经典的是二分法:
问题二分法解决方程求解问题
利用二分法,求方程 的实根,精确到两位小数。
分析二分法是一种典型的迭代问题,前面已经介绍了二分法定义,这里为了便于计算函数值 编制函数float function(float x)。在主函数中首先给出了有根区间 ,在程序中用[x1,x2]表示。由于不确定函数需要执行的次数,因此使用do-while循环,循环条件为区间中点的函数值小于 ,当函数值小于 时,近似认为当前的值为方程根。
实现
#include "stdio.h"
#include "math.h"
float function(float x)
{
float f;
f= x*x*(x-A)-B;
return f;
}
void main()
{
float x1,x2,x0,fx1,fx2,fx0;
x1=10;x2=-10;
fx1=function(x1);
fx2=function(x2);
do
{
x0=(x1+x2)/2.0;/*计算中点*/
fx0=function(x0);/*计算中点处的函数值*/
if(fx0*fx1<0)/*计算新的区间*/
{/*区间中点的函数值与x1的函数值正负号相反*/
/*区间中点的y坐标与x1点的y坐标在不同y半轴上*/
x2=x0;/*新区间为[x1,x0]*/
fx2=fx0;
}
else
{/*区间中点的y坐标与x1点的y坐标在相同y半轴上*/
x1=x0;/*新区间为[x0,x2]*/
fx1=fx0;
}
}while(fabs(fx0)>=1e-6); /*常量1e – 6为精度,可根据需要修改*/
printf("The root is %f",x0);
}
——————————————————————————————————————————
使用时,你直接把A,B换做常数即可,比如你取A=1,B=2,
那么就会输出:The root is 2.000000
这是个精确解,很多这种高次方程求不出精确解的只能得到近似解,你可以试试。
追问:
不可以,要的结果是x,不是f,我改成了x,但是还不对,你再帮帮忙看看哪里是不是编写错了
追答:
程序没错。你仔细看看,这里面涉及到一个式子变形,求x²(x-A)=B的根x,就是求使得f= x*x*(x-A)-B=0的x的值,你把f改成了x干嘛,只替换A,B程序输出就是x的解。
另外需要注意的一点是,二分法需要给定有根区间(解区间)即程序中的(x1,x2)=(10,-10),这个要根据函数的不同有所变化,具体判断方法可用函数求导和特殊值代入法。