请阅读奇异值分解维基百科文章。下面的代码是第二部分示例的代表。
import Jama.Matrix;
import Jama.SingularValueDecomposition;
public class JAMATest {
static public void printMatrix(Matrix m){
double[][] d = m.getArray();
for(int row = 0; row < d.length; row++){
for(int col = 0; col < d[row].length; col++){
System.out.printf("%6.4f\t", m.get(row, col));
}
System.out.println();
}
System.out.println();
}
public static void main(String[] args) {
double[][] vals = { {1., 0., 0., 0., 2.},
{0., 0., 3., 0., 0.},
{0., 0., 0., 0., 0.},
{0., 4., 0., 0., 0.}
};
Matrix A = new Matrix(vals);
SingularValueDecomposition svd = new SingularValueDecomposition(A);
System.out.println("A = ");
printMatrix(A);
System.out.println("U = ");
printMatrix(svd.getU());
System.out.println("Sigma = ");
printMatrix(svd.getS());
System.out.println("V = ");
printMatrix(svd.getV());
}
}
并产生输出:
A =
1.0000 0.0000 0.0000 0.0000 2.0000
0.0000 0.0000 3.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 4.0000 0.0000 0.0000 0.0000
U =
0.0000 0.0000 -1.0000 0.0000
0.0000 1.0000 -0.0000 0.0000
0.0000 0.0000 -0.0000 1.0000
1.0000 0.0000 -0.0000 0.0000
Sigma =
4.0000 0.0000 0.0000 0.0000 0.0000
0.0000 3.0000 0.0000 0.0000 0.0000
0.0000 0.0000 2.2361 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
V =
0.0000 -0.0000 -0.4472 -0.8944 -0.0000
0.0000 -0.0000 -0.0000 -0.0000 -0.0000
0.0000 1.0000 -0.0000 -0.0000 -0.0000
0.0000 -0.0000 -0.0000 -0.0000 1.0000
1.0000 -0.0000 -0.8944 0.4472 -0.0000
希望这可以帮到你。另外,就此问题,以下是Matlab的输出结果:
>> A = [1.0000, 0.0000, 0.0000, 0.0000, 2.0000
>> A
A =
1 0 0 0 2
0 0 3 0 0
0 0 0 0 0
0 4 0 0 0
>> [U, S, V] = svd(A)
>> U
U =
0 0 1 0
0 1 0 0
0 0 0 -1
1 0 0 0
>> S
S =
4.0000 0 0 0 0
0 3.0000 0 0 0
0 0 2.2361 0 0
0 0 0 0 0
>> V
V =
0 0 0.4472 0 -0.8944
1.0000 0 0 0 0
0 1.0000 0 0 0
0 0 0 1.0000 0
0 0 0.8944 0 0.4472
关于您的第一个问题,以下代码不会产生错误:
import Jama.Matrix;
public class JAMATest {
public static void main(String[] args) {
double[][] vals = {{1.,1.,0},{1.,0.,1.},{1.,3.,4.},{6.,4.,8.}};
Matrix A = new Matrix(vals);
}
}
所以你正在做的其他事情一定会导致它出现异常。尝试使用我的printMatrix方法替换你正在使用的任何东西,看看是否有帮助。