java 高斯消去法_Java中的高斯消元

我试图实现一个Matrix.class来学习一些Java。

现在，对于一个应该在高斯消元之后返回矩阵的方法，我将遇到一些困难，这些方法将在稍后用于查找矩阵的逆矩阵。

这是我到目前为止所提出的：

public Matrix gaussianElimination() {

Matrix inv = this.clone();

int i = 0;

int j = 0;

while (i

int pivot = i;

for (int k=i+1; k

if (Math.abs(inv.getArray()[k][j]) > Math.abs(inv.getArray()[pivot][j])) {

pivot = k;

}

}

if (inv.getArray()[pivot][j] != 0) {

inv = inv.swapRow(i, pivot);

double div = inv.getArray()[i][j];

for (double value : inv.getArray()[i]) {

value = value/div;

}

for (int u=i+1; u < inv.getHeight(); u++) {

double mult = inv.getArray()[u][j];

for (int l=0; l

inv.getArray()[u][l] = mult * inv.getArray()[i][l];

}

}

}

j++;

i++;

}

return inv;

}getArray()函数返回矩阵的double [] []，而getHeight()和getWidth()分别返回inv.length和inv [0] .length。

我跟着这个wikipedia页面的伪代码来实现这个算法

该方法返回顶部第一个元素的行，但不正确计算下部行的矩阵。

例如：

一个结果

0.2635522849474877 0.10001114673002853 0.442971040143471

0.2986277338922876 0.7517642579959294 0.09150190333830721

0.8913610667753092 0.8898546572478708 0.25592546060133237

INV结果

0.8913610667753092 0.8898546572478708 0.25592546060133237

0.26618513545092265 0.26573527978742995 0.07642644034471581

0.062426597261833985 0.06232109565941264 0.017923775508624545

我非常感谢任何帮助，因为我找不到解决方案。我可能在某处混合了一个指针，或者实现了错误的算法。