前续“问题六十二:怎么求一元十次方程在区间内的所有不相等的实根”和“问题五十九:怎么求一元六次方程在区间内的所有不相等的实根(3)——修正一个问题”
不管是求解一元六次方程还是一元十次方程,都存在这样一个区间端值问题:
当区间端值为0时,在求该区间上不等实根总个数时出错。
下方以一元六次方程为例说明问题:
如上图显示,该函数在区间[-2.5,2.5]上有六个不等实根,分别是-1/2, -1/3, -1/4, 1/4, 1/3, 1/2
所以,该函数在区间[-2.5,0]上应该有三个不等实根,分别是-1/2, -1/3, -1/4
但是,程序输出的结果是这样的:
显示不等实根个数还是6个,这明显是错误的。
我们将程序修改如下:(将端点零值用1e-6表示)
bool roots_num_equation_6th(float ee6[7], double a, double b, int &num) {
double temp1, temp2
, aaa, bbb;
double ee77[7][7] = {0.0};
double fun_a[7] = {0.0};
double fun_b[7] = {0.0};
int change_num_a = 0;
int change_num_b = 0;
if (fabs(a)<1e-6) {aaa = 1e-6;
}
else {
aaa = a;
}
if (fabs(b)<1e-6) {
bbb = 1e-6;
}
else {
bbb = b;
}
for (int i=0; i<6; i++) {
//determine f0, f1
ee77[0][i] = ee6[i];
ee77[1][i] = (6-i)*ee6[i];
}
ee77[0][6] = ee6[6];
for (int i=2; i<6; i++) {
//determine f2, f3, f4, f5
temp1 = ee77[i-2][0] / ee77[i-1][0];
temp2 = (ee77[i-2][1] - temp1*ee77[i-1][1]) / ee77[i-1][0];
for (int j=0; j<(6-i); j++) {
ee77[i][j] = temp1*ee77[i-1][j+2] + temp2*ee77[i-1][j+1] - ee77[i-2][j+2];
if (fabs(ee77[i][j])<1e-5) {
ee77[i][j] = 0.0;
}
}
ee77[i][6-i] = temp2*ee77[i-1][6-i+1] - ee77[i-2][6-i+2];
if (fabs(ee77[i][6-i])<1e-5) {
ee77[i][6-i] = 0.0;
}
if ((ee77[i][0] == 0.0) && (ee77[i][1] == 0.0) && (ee77[i][2] == 0.0) &&
(ee77[i][3] == 0.0) && (ee77[i][4] == 0.0) && (ee77[i][5] == 0.0) &&
(ee77[i][6] == 0.0)) {
break;
}
}
if (fabs(ee77[5][0])>1e-5) {
//determine f6
temp1 = ee77[4][0] / ee77[5][0];
temp2 = (ee77[4][1] - temp1*ee77[5][1]) / ee77[5][0];
ee77[6][0] = temp2*ee77[5][1] - ee77[4][2];
if (fabs(ee77[6][0])<1e-5) {
ee77[6][0] = 0.0;
}
}
else {
ee77[6][0] = 0.0;
}
for (int j=0; j<7; j++) {
fun_a[0] = fun_a[0] + ee77[0][j]*pow(
aaa, (6-j));
fun_b[0] = fun_b[0] + ee77[0][j]*pow( bbb, (6-j));
}
for (int i=1; i<7; i++) {
for (int j=0; j<(7-i); j++) {
fun_a[i] = fun_a[i] + ee77[i][j]*pow(
aaa, ((6-i)-j));
fun_b[i] = fun_b[i] + ee77[i][j]*pow( bbb, ((6-i)-j));
}
if ((fun_a[i]*fun_a[i-1]) < 0) {
change_num_a ++;
}
if ((fun_b[i]*fun_b[i-1]) < 0) {
change_num_b ++;
}
}
num = abs(change_num_a - change_num_b);
return true;
}
如上修改之后,程序运行结果如下: