用稀疏矩阵和ojAlgo求解矩阵微分方程

时间:2018-03-28 06:26:54

标签: matrix sparse-matrix ojalgo

我正在使用ojAlgo开发一个java进化工具,我尝试解决以下等式:

this differential equation

其中A是稀疏矩阵(目前矩阵的维数为2000 x 2000,稍后将进行缩放),A不对称且仅使用实数值。

我进行了一些研究,我试图在github wiki / javadoc上找到解决这个等式(使用SparseStore)的方法,但我没有找到办法。你能帮我找到我应该使用的方法/课程吗? 谢谢

1 个答案:

答案 0 :(得分:1)

没有直接/特定的方法来解决ojAlgo中的微分方程。您必须知道如何操作(使用笔和纸)然后ojAlgo可以帮助您执行计算。

这里的主要问题是寻找特征对,对吧?

    Eigenvalue<Double> evd = Eigenvalue.PRIMITIVE.make(matrix);

    evd.decompose(matrix);

    Array1D<ComplexNumber> values = evd.getEigenvalues();
    MatrixStore<ComplexNumber> vectors = evd.getEigenvectors();
    Eigenpair pair = evd.getEigenpair(0); // One of the pairs