31得票5回答
谷歌的Android OpenGL教程是否教授了错误的线性代数知识?

在帮助另一个用户解决关于响应触摸事件安卓教程的问题后,我下载了源代码,但看到的内容让我感到困惑。该教程似乎不能确定它想使用行向量还是列向量,对我来说看起来混乱不堪。 在安卓矩阵页面上,他们声称他们的约定是列向量/列主元,这是OpenGL的典型做法。 我是正确的,还是有什么我没注意到的地方?...

30得票2回答
如何在NumPy中尝试捕获由奇异性引起的非法矩阵操作

在NumPy中,我正在尝试使用linalg在牛顿-拉弗森(NR)方案的每个步骤中计算矩阵的逆(问题大小是有意为之的,以便我们可以解析地计算Hessian矩阵的逆)。然而,在接近收敛时,Hessian矩阵变得接近奇异。 在NumPy中是否有任何方法可以测试矩阵是否被视为奇异(计算行列式不够健壮...

30得票2回答
在三维空间中拟合一条直线

有没有算法可以从一组3D数据点返回一条直线的方程?我可以找到很多关于从2D数据集得到直线方程的资料,但是没有关于3D数据集的。

28得票5回答
沿一个轴进行卷积

我有两个二维数组,第一个轴的大小相同。在Python中,我想仅沿第二个轴对这两个矩阵进行卷积运算。我希望获得下面的C,但不想同时计算第一个轴上的卷积。 我有两个二维数组,第一个轴的大小相同。在Python中,我想仅沿第二个轴对这两个矩阵进行卷积运算。我希望获得下面的C,但不想同时计算第一个...

28得票4回答
如何在不将scipy稀疏矩阵转换为密集矩阵的情况下计算其行列式?

我正在尝试找到在 Python 中计算稀疏对称实矩阵行列式的最快方法,使用 Scipy 的 sparse 模块,但惊讶地发现没有行列式函数。我知道可以使用 LU 分解来计算行列式,但是不知道有什么简单的方法,因为 scipy.sparse.linalg.splu 的返回值是一个对象,并且实例化...

28得票2回答
scipy.linalg.eig返回协方差矩阵的复特征值?

协方差矩阵的特征值应该是实数且非负,因为协方差矩阵是对称的并且是半正定的。 然而,看一下使用scipy进行的以下实验: >>> a=np.random.random(5) >>> b=np.random.random(5) >>> ab...

28得票3回答
BLAS: gemm与gemv的区别

为什么BLAS有一个用于矩阵-矩阵乘法的gemm函数,以及一个用于矩阵-向量乘法的单独的gemv函数呢?难道矩阵-向量乘法不是矩阵-矩阵乘法的特殊情况,其中一个矩阵只有一行/列吗?

28得票2回答
Haskell线性代数?

我开始尝试使用Haskell进行线性代数测试。 有没有关于这个目的最好的软件包的推荐?还有其他用Haskell进行基本矩阵操作的好资源吗? Haskell维基网站列出了几个相关资源。 我目前的关注点是 hmatrix 和bindings-gsl,两者看起来都很有前途。

28得票4回答
有没有快速的矩阵指数计算方法?

除了简单的分治算法,有没有更快的矩阵指数计算方法来计算Mn(其中M是一个矩阵,n是一个整数)?

27得票2回答
为什么我的翻译矩阵需要转置?

我正在使用OpenGL开发一个小型图形引擎,我的翻译矩阵存在一些问题。我正在使用OpenGL 3.3、GLSL和C++。 情况是这样的:我已经定义了一个小立方体,我想在屏幕上渲染它。该立方体使用自己的坐标系,因此我创建了一个模型矩阵来对其进行转换。为了使自己更容易入手,我首先使用平移矩阵作为...