本文整理汇总了Java中org.ujmp.core.Matrix.svd方法的典型用法代码示例。如果您正苦于以下问题:Java Matrix.svd方法的具体用法?Java Matrix.svd怎么用?Java Matrix.svd使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类org.ujmp.core.Matrix
的用法示例。
在下文中一共展示了Matrix.svd方法的11个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Java代码示例。
示例1: testSVDWikipedia
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDWikipedia() throws Exception {
Matrix a = createMatrixWithAnnotation(4, 5);
if (!isSupported(a, MatrixLibraries.SVD, MatrixLayout.FAT, Size.SMALL, null)) {
return;
}
a.setAsDouble(1, 0, 0);
a.setAsDouble(2, 0, 4);
a.setAsDouble(3, 1, 2);
a.setAsDouble(4, 3, 1);
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例2: testSVDSquareSmall
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDSquareSmall() throws Exception {
Matrix a = createMatrixWithAnnotation(5, 5);
if (!isSupported(a, MatrixLibraries.SVD, MatrixLayout.SQUARE, Size.SMALL, null)) {
return;
}
for (int r = 0, v = 1; r < a.getRowCount(); r++) {
for (int c = 0; c < a.getColumnCount(); c++) {
a.setAsDouble(v++, r, c);
}
}
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例3: testSVDSquareLarge
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDSquareLarge() throws Exception {
if (!isTestLarge()) {
return;
}
Matrix a = createMatrixWithAnnotation(107, 107);
for (int r = 0, v = 1; r < a.getRowCount(); r++) {
for (int c = 0; c < a.getColumnCount(); c++) {
a.setAsDouble(v++, r, c);
}
}
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例4: testSVDSquareRandSmall
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDSquareRandSmall() throws Exception {
Matrix a = createMatrixWithAnnotation(10, 10);
if (!isSupported(a, MatrixLibraries.SVD, MatrixLayout.SQUARE, Size.SMALL, EntryType.RANDN)) {
return;
}
a.randn(Ret.ORIG);
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例5: testSVDSquareRandLarge
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDSquareRandLarge() throws Exception {
if (!isTestLarge()) {
return;
}
Matrix a = createMatrixWithAnnotation(109, 109);
a.randn(Ret.ORIG);
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例6: testSVDFatSmall
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDFatSmall() throws Exception {
Matrix a = createMatrixWithAnnotation(4, 6);
if (!isSupported(a, MatrixLibraries.SVD, MatrixLayout.FAT, Size.SMALL, null)) {
return;
}
for (int r = 0, v = 1; r < a.getRowCount(); r++) {
for (int c = 0; c < a.getColumnCount(); c++) {
a.setAsDouble(v++, r, c);
}
}
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例7: testSVDTallSmall
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDTallSmall() throws Exception {
Matrix a = createMatrixWithAnnotation(6, 4);
if (!isSupported(a, MatrixLibraries.SVD, MatrixLayout.TALL, Size.SMALL, null)) {
return;
}
for (int r = 0, v = 1; r < a.getRowCount(); r++) {
for (int c = 0; c < a.getColumnCount(); c++) {
a.setAsDouble(v++, r, c);
}
}
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例8: testSVDTallLarge
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDTallLarge() throws Exception {
if (!isTestLarge()) {
return;
}
Matrix a = createMatrixWithAnnotation(140, 121);
for (int r = 0, v = 1; r < a.getRowCount(); r++) {
for (int c = 0; c < a.getColumnCount(); c++) {
a.setAsDouble(v++, r, c);
}
}
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例9: main
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
public static void main(String[] args) throws Exception {
// create a dense empty matrix with 4 rows and 4 columns
Matrix dense = DenseMatrix.Factory.zeros(4, 4);
// set entry at row 2 and column 3 to the value 5.0
dense.setAsDouble(5.0, 2, 3);
// set some other values
dense.setAsDouble(1.0, 0, 0);
dense.setAsDouble(3.0, 1, 1);
dense.setAsDouble(4.0, 2, 2);
dense.setAsDouble(-2.0, 3, 3);
dense.setAsDouble(-2.0, 1, 3);
// print the final matrix on the console
System.out.println(dense);
// create a sparse empty matrix with 4 rows and 4 columns
Matrix sparse = SparseMatrix.Factory.zeros(4, 4);
sparse.setAsDouble(2.0, 0, 0);
// basic calculations
Matrix transpose = dense.transpose();
Matrix sum = dense.plus(sparse);
Matrix difference = dense.minus(sparse);
Matrix matrixProduct = dense.mtimes(sparse);
Matrix scaled = dense.times(2.0);
Matrix inverse = dense.inv();
Matrix pseudoInverse = dense.pinv();
double determinant = dense.det();
Matrix[] singularValueDecomposition = dense.svd();
Matrix[] eigenValueDecomposition = dense.eig();
Matrix[] luDecomposition = dense.lu();
Matrix[] qrDecomposition = dense.qr();
Matrix choleskyDecomposition = dense.chol();
}
示例10: testSVDFatLarge
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
@Test
public final void testSVDFatLarge() throws Exception {
if (!isTestLarge()) {
return;
}
Matrix a = createMatrixWithAnnotation(123, 142);
if (!isSupported(a, MatrixLibraries.SVD, MatrixLayout.FAT, Size.LARGE, null)) {
return;
}
for (int r = 0, v = 1; r < a.getRowCount(); r++) {
for (int c = 0; c < a.getColumnCount(); c++) {
a.setAsDouble(v++, r, c);
}
}
Matrix[] svd = a.svd();
Matrix prod = svd[0].mtimes(svd[1]).mtimes(svd[2].transpose());
assertEquals(0.0, prod.minus(a).getRMS(), TOLERANCE);
if (a instanceof Erasable) {
((Erasable) a).erase();
}
}
示例11: train
import org.ujmp.core.Matrix; //导入方法依赖的package包/类
public void train(ListDataSet dataSet) {
System.out.println("training started");
Matrix x = Matrix.Factory.zeros(dataSet.size(), getFeatureCount(dataSet));
int i = 0;
for (Sample s : dataSet) {
Matrix input = s.getAsMatrix(getInputLabel()).toColumnVector(Ret.LINK);
for (int c = 0; c < input.getColumnCount(); c++) {
x.setAsDouble(input.getAsDouble(0, c), i, c);
}
i++;
}
System.out.println("data loaded");
mean = x.mean(Ret.NEW, ROW, true);
for (int r = 0; r < x.getRowCount(); r++) {
for (int c = 0; c < x.getColumnCount(); c++) {
x.setAsDouble(x.getAsDouble(r, c) - mean.getAsDouble(0, c), r, c);
}
}
std = x.std(Ret.NEW, ROW, true, true);
for (int r = 0; r < x.getRowCount(); r++) {
for (int c = 0; c < x.getColumnCount(); c++) {
x.setAsDouble(x.getAsDouble(r, c) / std.getAsDouble(0, c), r, c);
}
}
Matrix xtx = x.transpose().mtimes(x);
Matrix[] svd;
if (numberOfPrincipalComponents == 0) {
svd = xtx.svd();
} else {
svd = xtx.svd(numberOfPrincipalComponents);
}
u = svd[0];
System.out.println("training finished");
}