linear-algebra

    1热度

    1回答

    我有一个地形和一个物体应该在地形上移动,所以我做了一个函数,可以检测到哪个脸部是上面定位的物体的原点,所以我应该设置Y的对象从物体原点向量与脸部垂直的向量之间的交点的精确高度,我有脸部的三个顶点,因此我可以计算它的法线和它的原点以及最大和最小指针(有界框)。 enter image description here 编辑: 另一种形式的问题: 如果我知道一个点在3D空间中的三角形我怎么能知道它Ÿ如

    1热度

    1回答

    我使用spilu预处理矩阵,但是,要将此预处理器传入cg(内置共轭梯度法),需要使用LinearOperator函数,有人可以向我解释参数matvec,以及为什么我需要使用它。下面是我现在的代码 Ainv=scla.spilu(A,drop_tol= 1e-7) Ainv=scla.LinearOperator(Ainv.shape,matvec=Ainv) scla.cg(A,b,maxit

    1热度

    1回答

    我有这个函数,计算相关矩阵和预期工作,但我试图使它更有效率,摆脱循环,但我'这样做很麻烦。我下面的功能: def correlation(X): N = X.shape[0] # num of rows D = X.shape[1] # num of cols covarianceMatrix = np.cov(X) # start with covariance

    0热度

    1回答

    在Python numpy中,如何在主对角线上制作一个空心对称矩阵,其他位置的矩阵是零?例如 import numpy as np I = np.identity(8) # the identity matrix I array([[ 1., 0., 0., 0.], [ 0., 1., 0., 0.], [ 0., 0., 1., 0.], [ 0.,

    0热度

    1回答

    我无法弄清楚如何从我的对象的模型矩阵中提取围绕全局y轴的旋转。 我有我的对象作为glm :: mat4的当前模型矩阵,我需要取消围绕y轴的旋转。 glm中有没有什么功能我没注意到?

    -1热度

    1回答

    如何构造具有相同对角元素的两个正半定矩阵?假设M1和M2是两个psd矩阵。我想要对角线(M1)和对角线(M2)的矢量相同。非对角元素必须有所不同。 [a e f g] [a s t u] [e b h i] [s b v w] [f h c j] [t v c x] [g i j d] [u w x d] 认为第一个矩阵为M1,第二个矩阵为M2。请记住M1和M2都是PSD。 在此先感谢!

    2热度

    1回答

    我有线性系统来解决这个问题,它由大型稀疏矩阵组成。 我一直在使用scipy.sparse库和它的linalg子库来做到这一点,但我不能让一些线性求解器工作。 这里是一个工作例子再现问题对我来说: from numpy.random import random from scipy.sparse import csc_matrix from scipy.sparse.linalg import

    0热度

    1回答

    我有两个等级为3的张量,换句话说就是两个3D矩阵。我想要拿这两个矩阵的点积。我很困惑继续这个问题。帮助我用公式来做到这一点。

    1热度

    2回答

    在this question中找到了找到具有无限多解的非正方形线性系统的特定解法的解决方案。这导致了另一个问题: 如何找到所有的解决方案为一个非正方形线性系统与无限多的解决方案,与R?(见下面的无限集合的解决方案的一个可能的描述) 实施例:线性系统 x+y+z=1 x-y-2z=2 相当于A X = B与: A=matrix(c(1,1,1,1,-1,-2),2,3,T) B=matrix

    1热度

    1回答

    我有一个尺寸为3x1000的阵列。事实上每一列都是有趣的。我想用它来计算一个尺寸为3x3x1000的数组,其中平板i是原始数组的列i(换句话说,v %*% t(v))的外积。有没有一种干净的方式来做到这一点? 下面是一个示例输入矩阵和输出数组,在2x4矩阵的情况下。 mat_in <- cbind(c(1, 2), c(3, 4), c(5, 6), c(7, 8)) arr_out <- ar