j6/example/matrixEx.js

Javascript scientific library (like R, NumPy, Matlab)

var R = require("../lib/matrix");
var M = R.M, V = R.V;
var m = [[1,2,3], [4,5,6]];
print('m=\n', m.strM());
print('m.rowSum()=\n', m.rowSum().strM());
print('m.colSum()=\n', m.colSum().strM());
print('m.tr()=\n', m.tr().strM());
print('m.madd(m)=\n', m.madd(m).strM());

print('m*mt=', m.mdot(m.tr()));

var v = [1,2,3];
print('vadd(v,v)=\n', R.vadd(v,v).str());
// print("v.sin()=", v.sin().str());
// print("v.norm()=%s", v.norm().str());
print("v.diag()=\n%s\n", v.diag().strM());
var A = [[1,2,3],[4,5,6],[7,3,9]];
// var AiA = M.mdot(M.minv(A), A);
var AiA = A.inv().mdot(A);
print("dim(A)=", V.dim(A));
print("det(A)=", M.det(A)); // det(A)=-30
print("AiA=\n", AiA.strM());
print("AiA.tr()=\n", M.tr(AiA).strM());
print("A=\n", A.strM());
// print("A.mul(0.1)=\n", A.mul(0.1).strM()); 
// 目前只允許矩陣乘矩陣，要加上矩陣 *+-/ 數量 (R, C, ....)
print("A.row(1)=", A.row(1));
print("A.col(1)=", A.col(1));
print("A.sum()=", A.sum());
// print("A.rowSum()=", A.rowSum());
// print("A.colSum()=", A.colSum());
// print("A.mean(row)=", A.rowMean().str());
// print("A.mean(col)=", A.colMean().str());

var D = v.diag();
print("D=", D);

print('I([2])=', M.identity([2]));


print("===========eigen================");
var Eλ = A.eig();
print("E=", Eλ.E.strM(), "λ=", Eλ.lambda.strM());
var E = Eλ.E, λ=Eλ.lambda;
print("E*[λ]*E-1=", E.mdot(λ.diag()).mdot(E.inv()).strM());

print("===========LU================");
var lu = M.lu([[0,0,1,1,1,2,2],[0,1,0,1,2,1,2],[2,-1,-1,2,-1,-1,2]]);
print('lu:', M.strM(lu));
var luSolve = M.luSolve(lu,[5,-8,13]);
print('luSolve:', luSolve.strM());

/*
print("===========Sparse================");
var S = [
[0,0,0,0,0,0],
[0,3,0,0,0,0],
[0,0,0,6,0,0],
[0,0,9,0,0,0],
[0,0,0,0,12,0],
[0,0,0,0,0,5],
[0,0,1,1,0,0],
[0,0,0,0,0,0],
];
print("sparse(S)=", M.sparse(S)); 
// The relation between A and its sparse representation SA is:
//  A[i][SA[1][k]] = SA[2][k] with SA[0][i] ≤ k < SA[0][i+1]
// 雖然有點不同，但可參考： http://openhome.cc/Gossip/AlgorithmGossip/SparseMatrix.htm
// sparse(S)=[[0,  0, 1,     3,    5, 6,  7 ], // (行)
//                 [ 1,  3, 6, 2, 6, 4,  5 ],  // (列)
//                 [ 3,  9, 1, 6, 1, 12, 5 ] ] // (值)
*/