问题标签 [ojalgo]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
158 浏览

java - Oj-Algo - 矩阵指数

我解决了以下等式:

在此处输入图像描述

为了解决它,我想使用矩阵指数:

在此处输入图像描述

我想了 3 种方法来做到这一点:

  1. 我可能会错过它,但 Oj-Algo 可以有一种简单的方法来计算 exp(A)(我在 MatrixStore javadoc 中没有找到它)
  2. 我从 EigenValue 方法 ([A] = [V][D][V] -1 ) 得到矩阵 D 和 V 然后我计算 然后首先出现的问题是我如何应用 x->exp(x*t) 函数到 D 的所有对角线元素?在此处输入图像描述
  3. 最后一个想法与 2 基本相同。但我之前将标量矩阵乘积存储在一个新矩阵 ([X] = [D]*(-t)) 中,然后我计算:

在此处输入图像描述

你能帮我找到我应该使用的最佳方式/方法/课程吗?谢谢

注意:这个问题是“后续问题”:初始问题

编辑:这是我现在尝试过的,这是最好的方法吗?

0 投票
1 回答
85 浏览

ojalgo - 使用 ojalgo 计算滚动平均值

使用 ojalgo 计算滚动平均值的最佳方法是什么?

First SampleSet 添加每个数字的次数与窗口大小一样多。

其次,代码意图并不明显。

我有:

哪个不是最好的...

0 投票
1 回答
77 浏览

ojalgo - 为什么 Double 不浮动

为什么推断类型是 Double 而不是 Float?

Primitive32Array.FACTORY或者换一种说法,和有什么区别Primitive64Array.FACTORY

0 投票
1 回答
327 浏览

java - ojAlgo:切片或提取子矩阵

在 java 库ojAlgo中,如何对矩阵进行切片或从现有矩阵中提取子矩阵?

例如, matrix A = [[1,2,3],[4,5,6],[7,8,9]]

我正在寻找一种方法,它看起来像slice(a,b,c,d)a,b 是行的开始和结束索引,c,d 是列的开始和结束索引。

例如,如果我调用A.slice(1,3,1,3),它应该返回[[5,6],[8,9]]

我尝试了中的slice方法SparseStore,但它没有按我预期的那样工作。

0 投票
1 回答
147 浏览

ojalgo - 从 OjAlgo 中的 MatrixStore 中删除行/列

是否有任何方便有效的方法可以从SparseMatrixojAlgo 中删除行(或列)?我试图RowView从我的矩阵中检索,然后将应该保持元素一个元素的行复制到新矩阵,但这不是非常有效的方法。

0 投票
1 回答
138 浏览

ojalgo - ojAlgo 中矩阵的秩

我目前正在使用 ojAlgo v45.1.0。我有一个关于如何获得矩阵的轨迹和总和的问题。由于我在 PrimitiveDenseStore 类中存储了一个矩阵,因此查看计算迹线和矩阵之和的方法是不可行的。那么,有谁知道如何在 PrimitiveDenseStore 类中获取矩阵的迹线和总和?先感谢您!

0 投票
1 回答
123 浏览

java - ojalgo 中的 QR/SVD 分解是否需要与列一样多的行?

在 ojalgo 中对 mxn 矩阵 A 进行 QR 或 SVD 分解时,我遇到了障碍。我的目的是找到列零空间的基础。如果 m >= n,一切正常。例如,秩为 2 的 5 x 4 矩阵 A 的转置 A' 的 QR 分解给出了一个 4 x 4 Q 矩阵,其最后两列跨越 A 的零空间。

另一方面,如果我从秩为 5 的 5 x 7 矩阵 A 开始(并对 A' 进行 QR 分解),我会得到正确的秩,但 Q 是 5 x 5 而不是 7 x 7,而且我没有'没有得到零空间基础。同样,具有相同矩阵 A 的 SVD 得到五个正奇异值(无零),Q2 矩阵为 5 x 7 而不是 7 x 7(无空向量)。

这是预期的行为吗?我找到了一个 n > m 矩阵的解决方法(将 nm 行的零添加到 A),但它很笨重。

0 投票
1 回答
72 浏览

kotlin - ojAlog - Kotlin 中的 ConvexSolver:二维数组

我正在尝试实现一些示例,因为我计划探索 ojAlgo 以进行优化。我的问题很简单。

在 Java 中,我可以轻松编写:

我尝试做与 kotlin 相同的事情:

但似乎 .rows 不能用我给出的论点来调用。

也许我做了一些愚蠢的事情,但我会很感激你的帮助。

谢谢你。

0 投票
1 回答
92 浏览

constraint-programming - 在 ojAlgo 中显示所有可能的解决方案

问题:是否可以调整上述内容以产生所有可能的结果而不是优化结果。

详细信息:给定一组装甲数据,我想生成一组满足我的约束的组合。

注意:在推荐库之前,请注意库必须与 Android 兼容。

0 投票
0 回答
49 浏览

java - 如何在 ojalgo 中获取求解器任务以返回正确的大小结果

我正在尝试 SolverTask 的各种实现者,目的是解决 A * x = b,在我的测试中,A 是 601 x 1100,b 是 601 x 1。到目前为止,我已经尝试过 GaussSeidelSolver 和 JacobiSolver,都通过直接调用 solve(A, b) 以及在 Gauss-Seidel 的情况下通过调用 preallocate(601, 1100, 1) 并将结果作为第三个传递给 solve()目的。无论我做什么,返回的矩阵都是 601 x 1 而不是 1100 x 1。

知道我做错了什么吗?