【题4】1元3次方程求解.docVIP

【题4】一元三次方程求解 有形如：ax3+bx2+cx+d＝0这样的一个一元三次方程。给出该方程中各项的系数（a，b，c，d 均为实数），并约定该方程存在三个不同实根（根的范围在—100至100之间），且根与根之差的绝对值≥1。 要求由小到大依次在同一行输出这三个实根（根与根之间留有空格），并精确到小数点后2位。 提示：记方程f（x）=0，若存在2个数x1和x2，且x1x2，f（x1）*f（x2）＜0，则在（x1，x2）之间一定有一个根。 输入：a，b，c，d 输出：三个实根（根与根之间留有空格） 输入输出样例 输入： 1 -5 -4 20 输出： -2.00 2.00 5.00 题解 这是一道有趣的解方程题。为了便于求解，将原方程 f(x)=ax3+bx2+cx+d＝0 变换成 f’(x)=x3+b’*x2+c’*x+d’=0 的形式（其中b’=，c’=，d’=），f(x)和f’(x)的根不变。设x的值域（-100至100之间）中有x， 其左右两边相距0.0005的地方有x1和x2两个数，即 x1=x-0.0005 x2=x+0.0005 x1和x2间的距离（0.001）满足精度要求（精确到小数点后2位）。若出现如图10.1.2所示的两种情况之一，则确定x为f’(x)=0的根。 图10.1.2 有两种方法计算f’(x)=0的根x: 1．枚举法 根据根的值域和根与根之间的间距要求（≥1），我们不妨将根的值域扩大100倍（-10000≤x≤10000），依次枚举该区间的每一个整数值x，并在题目要求的精度内设定区间：x1=，x2=。若区间端点的函数值f’(x1)和f’(x2)异号或者在区间端点x1的函数值f’(x1)=0，则确定为f’(x)=0的一个根。由此得出算法： 输入方程中各项的系数a，b，c，d ； b←b/a； c←c/a； d←d/a； a←1； {将方程变换为x3+b’*x2+c’*x+d’=0的形式} for x←-10000 to 10000 do {枚举当前根*100的可能范围} begin x1←(x-0.05)/100； x2←(x+0.05)/100； {在题目要求的精度内设定区间} if (f(x1)*f(x2)0) or (f(x1)=0){若在区间两端的函数值异号或在x1处的函数值为0，则确定x/100为根} then write(x/100:0:2，’’)； end；{for} 其中函数f(x)计算x3+b*x2+c*x+d： function f(x:extended):extended； {计算x3+b*x2+c*x+d} begin f←x*x*x+b*x*x+c*x+d； end；{f} 分治法 枚举根的值域中的每一个整数x(-100≤x≤100)。由于根与根之差的绝对???≥1，因此设定搜索区间[x1，x2]，其中x1=x，x2=x+1。若 ⑴f’(x1)=0，则确定x1为f’(x)的根； ⑵f’(x1)*f’(x2)0，则确定根x不在区间[x1，x2]内，设定[x2，x2+1]为下一个搜索区间 ⑶f’(x1)*f’(x2)0，则确定根x在区间[x1，x2]内。 问题是，如果确定根x在区间[x1，x2]内的话（f’(x1)*f’(x2)0），如何在该区间找到根的确切位置。采用二分法，将区间[x1，x2]分成左右两个子区间：左子区间[x1，x]和右子区间[x，x2]（其中x=）： 如果f’(x1)*f’(x)≤0，则确定根在左区间[x1，x]内，将x设为该区间的右指针（x2=x），继续对左区间进行对分；如果f’(x1)*f’(x)0，则确定根在右区间[x，x2]内，将x设为该区间的左指针（x1=x），继续对右区间进行对分； 上述对分过程一直进行到区间的间距满足精度要求为止（x2-x10.001）。此时确定x1为f’(x)的根。由此得出算法： 输入方程中各项的系数a，b，c，d ； b←b/a； c←c/a； d←d/a；a←1； {将方程变换为x3+b’*x2+c’*x+d’=0的形式} for x←-100 to 100 do {枚举每一个可能的根} begin x1←x；x2←x+1； {确定根的可能区间} if f(x1)=0