linear-algebra相关内容

np.linalg.eig为同一矩阵提供不同的特征向量

以下内容基于Python 3.7.6. 我正在尝试使用名为PySCF的软件包来解决简单的计算化学问题.其中一项计算涉及对称为Fock矩阵的2D数组的评估.PySCF使用函数get_fock()[ 1 ，表示右特征值没有保证的阶数，特征向量具有单位长度.这给了我们两次健全性检查. D, W = np.linalg.eig(F)断言(np.all(np.abs(W @ W.T-np.eye( ..
发布时间：2021-05-29 20:54:44 Python

Jacobi迭代到高斯-赛德尔

我为Jacobi方法编写了以下函数，需要对其进行修改以执行Gauss-Seidel 功能[x，iter] = jacobi(A，b，tol，maxit)%jacobi 迭代％x = zeros(size(b));[n，m] = size(A);i = 1;iter = maxit;对于 i=1:maxit对于j = 1:ny(j)=(b(j)-A(j，1:j-1)* x(1:j-1)-A(j，j ..
发布时间：2021-05-29 20:54:39 其他开发

贤者:有限阿贝尔群的自同构群?

对于有限阿贝尔群 G，比如说， G = AbelianGroup((4,4,5))，我希望Sage返回G的自同构群.这实现了吗? 官方排列组文档.在这里你去，虽然发电机可能不是你是什么寻找: G = AbelianGroup((4,4,5))H =差距(G).AutomorphismGroup()PermutationGroup(gap_group = H.AsPermGroup())具 ..
发布时间：2021-05-29 20:54:37 其他开发

在Python中消除非正方形矩阵的线性相关列

我有一个矩阵A = np.array([[1,1,1]，[1,2,3]，[4,4,4]])，我只想在新数组中使用线性独立的行矩阵.答案可能是 A_new = np.array([1,1,1],[1,2,3]]) 或 A_new = np.array([1,2,3],[4,4,4]) 由于我的矩阵很大，因此需要将矩阵分解为较小的线性独立的满秩矩阵.有人可以帮忙吗? 解决方案 执行此操作 ..
发布时间：2021-05-29 20:54:32 Python

用R解欠定线性系统

R 可以解决欠定线性系统: A =矩阵((1:12)^ 2,3,4，T)B = 1:3qr(A)$ rank#3qr.solve(A，B)#解决方案将有一个零，不一定是一个零#0.1875 -0.5000 0.3125 0.0000resolve(qr(A，LAPACK = TRUE)，B)#0.08333333 -0.18750000 0.00000000 0.10416667 (它给出 ..
发布时间：2021-05-29 20:54:08 其他开发

如何在c中实现Laplace展开算法?

由于无法弄清楚问题的中间部分，我无法找出一种使该算法起作用的方法.到目前为止，这是我的代码: int det(int matrixSize，int matrix [] [matrixSize]){int行列式= 0，matrixValues [matrixSize * matrixSize]，matrixFirstRowValues [matrixSize * matrixSize]；for( ..
发布时间：2021-05-29 20:54:05 其他开发

在Eigen3中实现Bartels–Stewart算法?

过去当我需要解 Sylvester 方程时，AX + XB = C，我使用了 scipy 的函数，solve_sylvester [1]，显然是通过使用Bartels-Stewart算法将事物分解为上三角形式，然后使用 lapack 求解方程来实现的. 我现在需要使用 eigen 求解方程. eigen 提供了一个函数 matrix_function_solve_triangular_syl ..
发布时间：2021-05-29 20:54:01 其他开发

如何使用相同维度的两个矩阵执行逐元素的自定义函数

未能找到任何有关此的信息.如果我有两个相同维度的m x n矩阵，有没有办法对它们应用numpty的逐个元素函数?为了说明我的意思: 自定义函数为F(x，y) 第一矩阵: array([[a，b]，[c，d]，[e，f]]) 第二矩阵: array([[g，h]，[i，j]，[k，l]]) 是否可以使用numpy中的上述两个矩阵在下方获得所需的输出 array([[F ..
发布时间：2021-05-29 20:53:58 Python

使用scipy.sparse.linalg线性系统求解器的问题

我要解决的线性系统是由大型稀疏矩阵组成的. 我一直在使用 scipy.sparse 库及其 linalg 子库来执行此操作，但是我无法获得一些线性求解器工作. 这是一个有效的示例，可以为我重现此问题: 来自numpy.random的 随机导入从scipy.sparse导入csc_matrix从scipy.sparse.linalg导入溶解度，分钟N = 10A = csc_matr ..
发布时间：2021-05-29 20:53:55 Python

如何找到具有不同特征值的两个矩阵的公共特征向量

我正在寻找或更确切地说在两个矩阵 A 和 B 之间建立公共特征向量矩阵X，例如: AX = aX带有"a"特征值对应的对角矩阵BX ＝ bX，其中"b"为0.特征值对应的对角矩阵 其中 A 和 B 是正方形和可对角线化的矩阵. 我看过类似的帖子，但没有设法得出结论，即当我构建由以下内容定义的最终想要的同构异形 F 时具有有效的结果: F = PDP ^ -1 我还阅读了维基百科 ..
发布时间：2021-05-29 20:53:54 其他开发

寻找具有无限多个解的非平方线性系统的所有解

在此问题中找到了找到特定问题的解决方案具有无限多个解的非平方线性系统的解.这导致了另一个问题: 如何找到带有R的具有无限多个解的非平方线性系统的所有解?解决方案) 示例:线性系统 x + y + z = 1x-y-2z = 2 等效于 A X = B ，其中: A = matrix(c(1,1,1,1，-1，-2)，2,3，T)B =矩阵(c(1,2)，2,1，T)一种[， ..
发布时间：2021-05-29 20:53:49 其他开发

使用Power方法从Python 3x3矩阵中获取特征值

我正在尝试使用Python中的Power方法从3x3矩阵中获取所有特征值.但是由于某种原因，我的方法从正确的特征值中返回了不同的特征值. 我的矩阵:A = [[1、2、3]，[2、4、5]，[3、5，-1]] 正确的特征值:[8.54851285，-4.57408723，0.02557437] 我的方法返回的特征值:[8.5485128481521926，4.5740872291 ..
发布时间：2021-05-29 20:53:45 Python

numpy的基数变化

给出numpy中的两个基本矩阵 basis_old 和 basis_new ，库中某处有一个函数来获取转换矩阵以转换矢量 vec basis_old 中的code>到在 basis_new 中的表示形式? 例如，如果我在标准基础上有一个向量 vec = [1,2,3] [1,0,0]，[0,1,0]，[0,0,1] ，如何将其转换为另一个基础，例如 e1 = [1 0 0]e2 = [0 ..
发布时间：2021-05-29 20:53:44 Python

用于进行高斯消去的BLAS/LAPACK例程

我是BLAS/Lapack的新用户，我只是想知道是否有一个例程可以消除高斯，甚至消除高斯-乔丹?我在Google上搜索并查看了他们的文档，但仍然找不到它们. 非常感谢您的帮助！ 解决方案 高斯消除与LU分解基本相同.例程xGETRF计算LU分解(例如，对于真正的双精度矩阵， DGETRF ).U因子对应于高斯消除后的矩阵.U因子在出口处存储在矩阵A的上三角部分(包括对角线)中. ..
发布时间：2021-05-29 20:53:41 其他开发

Scipy - 找到矩阵列空间的基数

我正在尝试编写一个简单的Simplex算法，其第一步是找到一个基本可行的解决方案: 选择A的线性独立列的集合B 将与不在B中的列相对应的x的所有分量设置为零. 求解m个所得方程，以确定x的分量.这些是基本变量. 我知道解决方案将涉及使用 scipy.linalg.svd (或 scipy.linalg.lu )和一些 numpy.argwhere / numpy.where 魔术， ..
发布时间：2021-05-29 20:53:40 Python

c(...％*％...)和sum(... * ...)之间的差异

此问题是 解决方案 我将这些注释编译为答案，如果我错过任何内容，其他人都可以加入. ％*％和 * 是完全不同的运算符: * 执行逐元素乘法，而％*％执行线性代数矩阵乘法.这些是非常不同的操作，演示如下: 1:4 * 2:5#[1] 2 6 12 201:4％*％2:5#[，1]#[1，] 40总和(1:4 * 2:5)#[1] 40 如果您要从两个向量相乘中查找单个汇总统计信 ..
发布时间：2021-05-29 20:53:23 其他开发

我如何找出A * B是Numpy中的Hadamard或Dot产品?

如果我在导入numpy的python代码中看到以下行: c = a * b 确定此操作是以Hadamard(元素方式)还是点积(点方式)操作执行的最简单，最实用的方法是什么? 对于Hadamard产品而言，A和B的列和行大小必须相同是正确的.对于点积，仅A的列大小必须与B的行大小相同，对吗?所以我可以查找两者的形状并找出使用了哪个操作? 解决方案 这在 ..
发布时间：2021-05-29 20:53:22 Python

Three.js-向上向量如何与lookAt()一起使用?

我试图了解up.vector和lookAt()在three.js中如何协同工作.我正在设置该axisHelper的向上矢量，以便Y轴始终指向目标地理位置，这标志着向上矢量的位置.它绕X轴旋转，如X和Y所预期的那样工作.当我尝试调整向上矢量的Z值时，我希望轴绕X轴旋转，但是什么也没发生. http://jsfiddle.net/68p5r/4/[编辑:我添加了地理信息以显示向上目标位置.] ..
发布时间：2021-05-29 20:53:18 其他开发

左矩阵除法和Numpy解决

我正在尝试将包含\运算符的代码从Matlab(八度)转换为Python.示例代码 B = [2; 4]b = [4; 4]\ 这有效并产生1.2作为答案.使用此网页 http://mathesaurus.sourceforge.net/matlab-numpy.html 我翻译成: 将 numpy 导入为 np将 numpy.linalg 导入为 linB = np.arra ..
发布时间：2021-05-29 20:53:16 Python

O(mn)是否在O(n ^ 2)中?

一个简单的问题.使用m x n矩阵，我正在执行一些O(mn)运算.我的问题是O(mn)是否在O(n ^ 2)中.在大O上查看Wikipedia，我会这样想，但是我在复杂性方面一直很糟糕，所以我希望有人能澄清一下. 对于 m x n 矩阵， 解决方案 O(mn)意味着您要对矩阵的每个值进行恒定的工作. O(n ^ 2)表示，对于每一列，您正在做的工作是O(#列).请注意，此运行时间随行数 ..
发布时间：2021-05-29 20:53:13 其他开发