Я пытаюсь умножить 2 матрицы. Я пробовал с программным обеспечением R: проверьте эту тему: Умножение двух матриц в R
Теперь я пытаюсь сделать то же самое на Java.
Я использую библиотеку Jama для матрицы.
моя силовая функция
public Matrix power(Matrix M, double p) {
EigenvalueDecomposition evd = M.eig();
Matrix D = evd.getD();
for (int i = 0; i < D.getColumnDimension(); i++) {
D.set(i, i, Math.pow(D.get(i, i), p));
}
Matrix V = evd.getV();
return V.times(D.times(V.transpose()));
}
double[][] matDouble1 = {{0.25, 0.25, 0.25, 0.25}};
double[][] matDouble2 = {{0, 0, 1, 0},
{0.333, 0, 0.333, 0.333},
{1, 0, 0, 0},
{0, 1, 0, 0}};
Matrix mat1 = new Matrix(matDouble1);
Matrix mat2 = new Matrix(matDouble2);
mat2 = power(mat2, 10000);
mat1 = mat1.times(mat2.transpose());
Когда отображается мат1, я получаю:
0.25018740608813655 0.2498123125257854 0.2501874060881363 0.24981231252578548
вместо
0.5 0 0.5 0
с R я делаю
mpower = function(M,p) {
A = as.matrix(M)
if (dim(A)[1] != dim(A)[2]) stop("not a square matrix")
# M^{-1} = the matrix inverse of M
if (p==-1) return(solve(A))
# M^0 = I
if (p==0) return(diag(1,dim(A)[1],dim(A)[2]))
# M^1 = M
if (p==1) return(A)
if (p < -1) stop("only powers >= -1 allowed")
if (p != as.integer(p)) stop("only integer powers allowed")
R = A
for (i in 2:p) {
R = R %*% A
}
return(R)
}
mat1<-matrix(c(0.25,0.25,0.25,0.25),nrow=1)
mat2<-matrix(c(0,0,1,0,0.3,0,0.3,0.3,1,0,0,0,0,1,0,0),nrow=4)
mat1%*%t(mpower(mat2,10000))
я получил
[1,] 0.4996252 0 0.4996252 0