Question

java 数值 分析 的程序

public double[] gauss(Numbers number) //消去法
{
int x = number.x;
int y = number.y;
double[] res = new double[x];
double[][] num = number.getNumbers();
for(int i = 1; i<= x-1 ; i++)
{
for(int j = i+1 ; j<= x ; j++ )
{
double a = (num[j-1][i-1])/(num[i-1][i-1]);
for(int m = 1; m <= y ; m++ )
{
num[j-1][m-1] -= a*num[i-1][m-1];
}
}
}
System.out.println("变换之后的矩阵为");
for(int i = 1 ; i<= x ; i++)
{
for(int j = 1 ; j <= y ; j++)
{
System.out.print(num[i-1][j-1]+" ");
}
System.out.println();
}
for(int i = x ; i >= 1 ; i-- )
{
double result = 0;
for(int j = 1 ; j <= x-i ; j++)
{
result += num[i-1][y-j-1]*res[x-j];
}
result = num[i-1][y-1] - result;
res[i-1] = result/num[i-1][i-1];
}
return res;
}
额...我写了一个程序来算消去法 就是一个矩阵解方程组的那种
试验了一些小的方程 都能解对
然后进行书上要求的10乘11的矩阵
X1 到X10 是10个解 有八个对 剩下两个是这样的：
方程组的解为
X3=-1.1398289719484943E-14
X6=-4.952844806906863E-16
其中X3 和 X6 让我很费解 不知道是算错了呢 还是计算机数据类型表示的问题

如果不是，就是程序算法错了
附上算法代码.....
比较闹心的代码 大侠们指教

Answer 1

可能是浮点数精度的问题，double、float直接运算是不能得到精确值的
要么改为BigDecimal，要么保留小数位数

Answer 2

在学校时用c 写过一个 多种方法解线性方程的程序。 其中用到了，迭代法，高斯消元法，及约当法。
可以参考。
确实有楼上说的精度问题，但那时处理了下，没有完全解决