我正在寻找计算矩阵求逆的最有效的方法。我知道这与LU分解有关，或者是把它分成两个矩阵，上三角和下三角矩阵，以及这两个系统的求解，但我不能用MATLAB写出代码。我甚至不知道从哪里开始，也不知道如何在我的问题中使用系统的结果。我真的需要它，因为我必须求解一个10.000 x 10.000的矩阵，而这个简单的方法永远需要花费很长时间。
浏览 1提问于2014-03-23得票数 0
1回答
如何求大稀疏矩阵的伪逆？
、、、、
我必须反演一个大型稀疏矩阵(50000 X 12000)。它最初以numpy.ndarray的形式存储，矩阵的大小约为3.5GB。我已经尝试过使用numpy.linalg.pinv来转换这个矩阵，但是它崩溃了jupyter笔记本内核。将此numpy.ndarray转换为scipy.sparse.csr_matrix (稀疏矩阵格式)是可行的，但我不知道有任何函数可以计算csr_matrix的伪逆。
如何求大型稀疏矩阵的伪逆？
浏览 4提问于2022-09-06得票数 0
1回答
我有两个稀疏矩阵A(亲和力矩阵)和D(对角矩阵)，维数为100000*100000。我必须计算拉普拉斯矩阵L= D^(-1/2)*A*D^(-1/2)。我正在使用scipy CSR格式的稀疏矩阵。
我没有找到任何方法来求稀疏矩阵的逆。如何求稀疏矩阵的L和逆矩阵？还建议使用python来做这件事是否有效，或者我是否应该调用matlab函数来计算L？
浏览 1提问于2012-02-12得票数 3
3回答
矩阵模板库矩阵求逆
、、、、
我正在尝试用版本Boost boost_1_37_0和MTLmtl4-alpha-1-r6418来求逆矩阵。我好像找不到矩阵求逆的代码。我在谷歌上搜索了一些例子，它们似乎引用了上面的版本中似乎缺少的lu.h。有什么提示吗？
建议复制lu.h，但这似乎来自MTL2而不是MTL4。我在用VS05或更高版本的MTL2编译时遇到了问题。
那么，你知道如何在MTL4中进行矩阵求逆吗？
更新:我想我更了解马特了，我要去了。
浏览 6提问于2008-12-03得票数 3
回答已采纳
3回答
我想要得到一个倍频程的可逆矩阵，但作为整数矩阵，所以：
x = [9,15;19,2];
inv(x)
这就是我得到的：
[-0.0074906, 0.0561798; 0.0711610, -0.0337079]
但是我想让[22,17;25,21]知道如何求矩阵的逆？
浏览 0提问于2010-01-17得票数 7
回答已采纳
1回答
求CSS3变换矩阵的逆
、、、
如果我有一个CSS3变换矩阵，我可以如何求逆这个矩阵？
例如，如果我有这个矩阵，旋转45度，平移10px，10px
矩阵( 0.7071067811865475，0.7071067811865476，-0.7071067811865476，0.7071067811865475，-0.0000000000000008881784197001252，14.142135623730951)
如何计算逆数？
矩阵( 0.7071067811865475，-0.7071067811865476，0.7071067811865476，0.7071067811865475，-14.1421356237309
浏览 2提问于2012-03-20得票数 1
回答已采纳
1回答
我试着求出协方差矩阵的逆。到目前为止，我有：
A<- c(1, 3, 5, 7)
B<- c(8, 4, 1, 4)
C<- c(9, 2, 4, 6)
D<- c(1, 4, 2, 0)
M<- matrix(c(A, B, C, D), 4)
> M
[,1] [,2] [,3] [,4]
[1,] 1 8 9 1
[2,] 3 4 2 4
[3,] 5 1 4 2
[4,] 7 4 6 0
>
> M<- cov(M)
浏览 34提问于2018-02-27得票数 1
回答已采纳
1回答
我有一个矩阵M，它有576行，15 columns.First列都是1。第二列从1开始，逐个增加到576。第三列是第二列的正方形。第四列是第二列的三次方。第五列是第二列的四次方。它是这样的。最后，第15列是第二列的14次方。该操作是将矩阵M的转置与其自身相乘，并取其结果的逆。然而，Matlab给出了这样的警告:警告:矩阵接近于单数或缩放不良。结果可能不准确。RCOND = 1.889839e-79。这是一种与之相反的方式吗？我也尝试了这个操作符，\但同样的警告也出现了。
浏览 1提问于2016-10-14得票数 0
3回答
求高条件数矩阵的逆
、、、、
我试图找到的反义词，这样我就可以和一起使用它。然而，我从得到的结果是一个充满1.02939420e+16的矩阵。我一直试图找出原因，考虑到Wolfram会给我正确的答案，这似乎与矩阵的条件数有关，在这种情况下是3.98290435292e+16。
虽然我想了解这背后的数学，但我现在真正需要的只是一个解决这个问题的方法，这样我就可以继续实现了。有没有办法找出这样的矩阵的逆？或者是否有可能直接从数据中找到逆协方差矩阵？
编辑：矩阵数据(与pastebin链接相同)
[[ 0.46811097 0.15024959 0.01806486 -0.03029948 -0.12472314 -0.11
浏览 5提问于2016-04-22得票数 0
回答已采纳
1回答
我有计算下三角矩阵求逆的代码。如何通过对下面的代码稍加修改来计算上三角矩阵的求逆？ function L = L_inv(A)
[n,n] = size(A);
L = zeros(n);
for i=1:n
L(i,i) = 1/A(i,i);
for j=i+1:n
L(j,i)=-A(j, i:j-1)*L(i:j-1,i)/A(j,j);
end
end
浏览 220提问于2021-01-05得票数 0
回答已采纳
2回答
如何求几乎奇异矩阵的逆？
、、、
我正在用C++和CUDA实现一个算法。但是当我试图找到一个特殊矩阵的逆时，我遇到了麻烦。该矩阵具有以下特点：
它是一个方阵(假设：(m+3)x(m+3)，m>0)；
它的转置矩阵是它的自我；
它的主要对角线必须是零；
它必须在右下角有一个3x3.0矩阵；
你可以用这种形式考虑这个矩阵:h=A，B ;B‘,0；
我尝试过一些方法，但都失败了：
伪逆矩阵：
我首先使用matlab，当我尝试使用inv(H'*H)：警告时，得到了错误或警告:矩阵与工作精度是奇异的，或矩阵接近奇异或严重缩放的。
几种逼近方法：
参考资料如下：i发现了两种方法:G
浏览 3提问于2014-01-05得票数 3
回答已采纳
1回答
对于一个主题，我使用numpy实现了来自的Hill Cipher示例。到目前为止我的代码如下：
import numpy as np
msg = np.array([0, 2, 19]) # CAT
key = np.matrix([
[6, 24, 1],
[13, 16, 10],
[20, 17, 15]
])
ciphertext = np.mod(np.dot(key, msg), 26) # [5, 8, 13] # POH
# And here things go wrong
invKey = key.I
invKey2 = np.linalg.in
浏览 1提问于2018-08-19得票数 0
2回答
如何求距离矩阵的逆矩阵？
、、、、
我有一个巨大的距离矩阵。 示例：(10000 * 10000).. 有没有一种有效的方法来求逆矩阵？ 我试过numpy的Inv()，但它太慢了。 有没有更有效的方法？
浏览 30提问于2020-10-12得票数 0
2回答
我通过将向量与模型视图矩阵相乘来转换向量。
glGetDoublev( GL_MODELVIEW_MATRIX, cam );
newX=cam[0]*tx+cam[1]*ty+cam[2]*tz+cam[3];
newY=cam[4]*tx+cam[5]*ty+cam[6]*tz+cam[7];
newZ=cam[8]*tx+cam[9]*ty+cam[10]*tz+cam[11];
如何获得反向变换的模型视图矩阵的逆矩阵？
浏览 0提问于2014-04-21得票数 0
3回答
我正在使用一个基于ARM Cortex-A8的处理器，并且我有几个地方可以计算3x3矩阵的逆操作。
由于Cortex-a8处理器有一个霓虹灯SIMD处理器，我有兴趣将这个协处理器用于3x3矩阵求逆，我看到了几种4x4实现(Intel SSE和)，但没有看到使用SIMD指令执行3x3矩阵求逆操作的地方。freevec使用的方法是将4x4矩阵拆分为4，2x2矩阵，然后对每个矩阵进行微小操作，最后将结果重新组合，得到4x4矩阵的逆。我在3x3矩阵中看不到这样的机会。
有谁知道如何使用SIMD指令进行3x3矩阵求逆，这对我会很有帮助吗？
关于Vikram
浏览 2提问于2010-07-26得票数 5
回答已采纳
我需要取一个矩阵的逆矩阵来拟合我在Rust中的数据的线性回归。我已经尝试过使用ndarray来做这件事，但是没有找到一个不使用BLAS后端的好方法。有什么建议吗?如何在不使用BLAS后端的情况下在Rust中找到矩阵的逆？
浏览 34提问于2020-12-18得票数 0
1回答
用于加密的矩阵乘法求逆
、、、、
我正在尝试组合一个加密算法，但我陷入了以下问题，我甚至不知道它应该是这样的或不是这样的！
问题是：
我有16字节的矩阵要乘以16,16矩阵，结果是一个16字节的矩阵。
然后我应该将结果矩阵乘以逆，这里我想我应该得到原始的16字节矩阵(根据算法数据表)。
所以你能告诉我怎样才能找回原始矩阵吗？
提前感谢您的帮助。
致以敬意，
英语。Aws
浏览 2提问于2011-12-17得票数 0
2回答
对于循环求逆矩阵
、、
我正在尝试在matlab中求一个矩阵的逆，但是我正在努力。
我假设我需要以某种方式使用for循环来得到一个反3x3矩阵，每个矩阵包含801个点。
inv(A11(1) A12(1) A13(1);A21(1) A22(1) A23(1);A31(1) A32(1) A33(1))
例如，这个倒数将给出矩阵的801个点中的第一个
浏览 2提问于2016-04-27得票数 0
1回答
我试图乘两个矩阵，然后计算J中结果矩阵的逆，我的矩阵如下所示：
x =: i. 2 3
0 1 2
3 4 5
y =: i. 3 2
0 1
2 3
4 5
我会用什么命令将这两者相乘，然后计算逆？
如果矩阵没有逆，那么处理这个问题的命令是什么？
浏览 4提问于2017-09-29得票数 1
回答已采纳
1回答
该任务要求我生成一个50列50行的矩阵，其中种子1007092020的随机库的范围为0,1。 import numpy as np
np.random.seed(1007092020)
A = np.random.randint(2, size=(3,3)) 然后我必须找到A的一个具有LU分解的逆矩阵。我不知道该怎么做。
浏览 33提问于2021-09-19得票数 0
1回答
我在Google上搜索，找不到用GSL计算矩形矩阵逆的函数。由于很难找到答案，当其他人需要寻找矩形矩阵的逆时，这里的答案将对他们有所帮助。如果不可能使用GSL，那么请提供一些易于使用和提供矩形矩阵逆的可选库。
浏览 8提问于2013-08-14得票数 5
我试图用一些符号变量来求逆一个矩阵，但是Matlab只返回了'FAIL‘。我使用inv(K)。这是Matlab经过一些计算后输出的K：
K =
[ 11/80, 7/80, -11/80, -7/80, 0, 0]
[ 7/80, 11/80, -7/80, -11/80, 0, 0]
[ -11/80, -7/80, 11/80, 7/80, 0, 0]
[ -7/80, -11/80, 7/80, 11/80, 0, 0]
[ 0, 0, 0, 0, 0, 0]
[ 0, 0,
浏览 8提问于2012-03-29得票数 1
回答已采纳
1回答
加速数值矩阵求逆
、、、、
我正在使用Numpy/Scipy来求一个20k矩阵的逆矩阵，速度很慢。我试过了：
(1) M_inv = M.I
(2) Ident = np.Identity(len(M))
M_inv = scipy.linalg.solve(M, Ident)
(3) M_inv = scipy.linglg.inv(M)
但没有看到任何加速。
有没有其他方法可以加速这一过程？
浏览 0提问于2012-05-03得票数 3
回答已采纳
1回答
我有一个通用矩阵A，它是对称的，正定的，并且是稀疏填充的(它也很大，比如由几十到几十万行组成)。我想让一个神经网络学习如何找到一个矩阵B，它类似于A的逆(即，乘以A将使光谱半径的最小扩展)尽可能准确和快速。
我想知道是否有一些人工智能分支专注于这类问题，或者如果没有，我如何处理它。
我知道用来做这件事的经典算法，但我只想看看机器学习在这个问题上提供的可能性。
浏览 0提问于2016-10-19得票数 11
3回答
我正在用python编写一些代码，这些代码需要频繁地对大方阵(100-200行/列)求逆。
我达到了机器精度的极限，所以我开始尝试使用mpmath来进行任意精度的矩阵求逆，但即使使用gmpy，它也非常慢。
以精度30 (十进制)对大小为20、30、60的随机矩阵求逆需要大约0.19秒、0.60秒和4.61秒，而在mathematica中相同的操作需要0.0084、0.015和0.055秒。
这是在arch linux机器上使用python3和mpmath 0.17 (不确定gmpy版本)。我不知道为什么mpmath会慢这么多，但是有没有开源库可以达到mathematica的速度(即使是1/2的
浏览 5提问于2013-03-10得票数 3
1回答
我必须对一个大的稀疏矩阵求逆。我无法摆脱矩阵求逆，唯一的捷径就是了解主对角线元素的概念，而忽略非对角线元素(我不希望这样做，但作为一种解决方案，它是可以接受的)。
我需要求逆的矩阵通常很大(40000 *40000)，并且只有一小部分非零对角线。我目前的方法是构建稀疏的所有内容，然后
posterior_covar = np.linalg.inv ( hessian.todense() )
这显然需要很长的时间和大量的内存。
有什么提示，或者这只是一个耐心的问题，或者让问题变得更小？
浏览 0提问于2013-02-28得票数 6
回答已采纳
2回答
我在python中做矩阵求逆，我发现结果因数据规模的不同而不同，这非常奇怪。
在下面的代码中，预期A_inv/B_inv = B/A。但是，它显示A_inv/B_inv和B/A之间的差异随着数据规模的不同而变得越来越大……这是不是因为Python不能精确地计算大值矩阵的逆矩阵？
此外，我检查了B的条件数，无论规模是多少，它都是一个常数~3.016。
谢谢！
import numpy as np
from matplotlib import pyplot as plt
D = 30
N = 300
np.random.seed(10)
original_data = np.random.sa
浏览 1提问于2017-08-31得票数 1
3回答
我正在尝试对iPhone加速度计输出的数据应用卡尔曼滤波器。我需要尽可能快地执行矩阵乘法和求逆，所以我对使用GPU执行这两项任务的可能性感到好奇。到目前为止，我只找到一个矩阵乘法的引用：
float mBone01[16] = { ... }
float mBone02[16] = { ... }
float mResult[16];
glMatrixMode ( GL_MODELVIEW );
glLoadIdentity( );
glLoadMatrix ( mBone01 );
glMultMatrix ( mBone02 );
glGetMatrix ( GL_MODELV
浏览 0提问于2010-01-10得票数 3
3回答
我正在尝试用boost做一个简单的矩阵求逆运算。但是我得到了一个错误。基本上，我试图找到的是inversted_matrix =逆(trans(矩阵)*矩阵)，但我得到了一个错误
Check failed in file boost_1_53_0/boost/numeric/ublas/lu.hpp at line 299:
detail::expression_type_check (prod (triangular_adaptor<const_matrix_type,
upper> (m), e), cm2)
terminate called after throwing
浏览 0提问于2013-06-22得票数 5
回答已采纳
6回答
矩阵求逆
、、
我正在为class开发一个程序，它读取一个矩阵，然后用单位矩阵打印出来。然后，我必须一步一步地将其简化为其单位矩阵，同时将相同的计算应用于单位矩阵，从而得到该矩阵的逆。
我被代码的行减少部分卡住了。我正在尝试创建一个循环，它将一个nxn矩阵简化为一个单位矩阵，而不是为每一列创建一个。我需要一些帮助的循环，行减少矩阵一列一列。
任何人都有任何想法。(我是编程的初学者，我使用C语言，发布的解决方案也可以，但我正在寻找一种算法，可以一次减少所有行)。
浏览 0提问于2009-11-06得票数 0
1回答
实时求逆矩阵
、、、
我一直在开发一个实时工作的C语言控制软件。软件实现了被控系统的离散状态空间观测器等。对于观察者的实现，需要计算4x4维的矩阵的逆。逆矩阵计算必须每隔50微秒进行一次，值得一提的是，在这段时间内还将进行其他非常耗时的计算。因此，逆矩阵计算所消耗的时间必须远远少于50微秒。还有必要说明的是，所使用的DSP没有支持浮点运算的ALU。 我一直在寻找一些有效的方法来做到这一点。我的一个想法是准备计算矩阵4x4的行列式的通用公式和计算4x4矩阵的伴随矩阵的通用公式，然后根据下面给定的公式计算矩阵的逆矩阵。 ? 你对这种方法有什么看法？
浏览 46提问于2020-09-10得票数 1
在我的损失函数中的某处，我对一个大小为64*64的复数矩阵求逆。尽管torch.tensor支持复数矩阵求逆，但无法在训练循环中计算梯度，因为我得到了这个错误：
RuntimeError: inverse不支持对具有复杂类型的输出进行自动区分。
有没有人有解决这个问题的办法？也许是一个自定义函数而不是torch.inverse？
浏览 5提问于2021-01-14得票数 3
如何计算一个巨大( 14000*14000)矩阵的逆矩阵？我正在使用R进行尝试，错误消息显示为：
Error: cannot allocate vector of size 762.9 Mb
我知道这表示无法获取内存。任何帮助都将不胜感激。
即使我试图求一个标量单位矩阵的逆矩阵，R也不能解决它。标量单位矩阵的逆矩阵将是相同的矩阵，但R也不能求解它。
> solve(diag(1,10000,10000))
Error: cannot allocate vector of size 762.9 Mb
浏览 1提问于2014-05-22得票数 0
简单的问题。它有一个运算吗，还有，你如何创建一个单位矩阵？我已经检查了文档，但找不到任何操作。
浏览 66提问于2018-08-12得票数 2
回答已采纳
1回答
我正在试着写一个代码来解密希尔密文。我们被教导如何在使用线性代数的情况下手动完成这项工作。我的代码需要我找到一个方阵的逆矩阵。 假设我正在为一个2x2矩阵手动执行此操作。如果我的矩阵是： |1 4 |
|11 11| 如果我找到行列式(mod 26，因为它毕竟是用于密码的)和矩阵的伴随项，我将得到： |121 44|
|-121 11| 我可以在以后找到mod 26的元素，但我想说的是我将得到整数。 现在为了编写代码，我遇到了时间紧张，所以我想我可以使用numpy来求这个矩阵的逆，而不是写单独的代码来求逆和伴随。这就是我所做的： a=[[1,11],[4,11]]
print(linal
浏览 25提问于2019-06-27得票数 0
1回答
我有一个矩阵，它没有通过我正在计算的朴素贝叶斯分类器的奇异测试。我正在处理等式中的ln(det(sigma))部分。
if np.linalg.cond(covarianceMatrix) < 1/sys.float_info.epsilon:
return np.log(np.linalg.det(covarianceMatrix))
else:
return a pseudo determinant
当协方差矩阵是奇异的时，我必须找到伪行列式。我该怎么做呢？
浏览 2提问于2016-03-21得票数 5
我刚刚打印了大小为n*n的矩阵
代码是(假设已经定义了所有矩阵)：
public class ReadContents {
public static void main(String args[]) {
System.out.println();
System.out.println("The Matrix Is:");
double mat[][]=new double[col][col];
for(k=0;k<col;k++) {
for(p=0;p<col;p
浏览 2提问于2011-05-06得票数 1
2回答
我实现了一个类似AES的分组密码。但是我不能解密的原因是我不能得到逆MDS矩阵。我使用的MDS矩阵是GF(2^8) \implies GF(2^8)类AES上的3x3mds矩阵。
\begin{bmatrix} 1 & 2 &2 \\ 2 & 2 & 1\\ 2 & 1 & 2\\ \end{bmatrix}
加密过程与AES混合列相同，所使用的MDS矩阵如下所示。
我需要解密所需的MDS逆矩阵。在GF(2^8)上这个MDS矩阵的逆是什么？
我寻找欧几里得-沃利斯算法或扩展的欧几里德算法，但我不太理解它。
如何求3x3MDS矩阵的逆？
浏览 0提问于2020-05-16得票数 4
1回答
我正在上计算机图形学的速成课程，我们刚刚介绍了2D变换，我的导师指出，计算矩阵的逆比转置计算更密集，这就是为什么正交矩阵在计算中如此有用(因为正交矩阵的转置也是它的逆矩阵)。由于时间的限制，以及课程的性质，他没有详细说明为什么会这样，我想知道这里是否有人能够做到这一点。
我特别感兴趣的是每个进程中涉及的CPU/GPU指令的差异，或者如果我弄错了，那么是堆栈中出现瓶颈的最低级别。我也会对学习任何资源、书籍、网站等感兴趣，在那里我可以了解更多关于这些效率/瓶颈的知识。
浏览 2提问于2011-09-16得票数 1
回答已采纳
2回答
我正在通过Cholesky分解使用不同的库来计算方阵的逆矩阵。然而，我的结果并不像我预期的那样。我不是数学专家，但我希望能得到一个更接近的结果。
我正在使用MLK，magma和CULA库在CPU和GPU中计算矩阵的逆。在做完计算之后，我注意到这些库的结果总是在一个元素上不同。假设我想要计算A= [0.237306,0.000458;0.000458,0.238497]的倒数
A[0] = 0.237306
A[1] = 0.000458
A[2] = 0.000458
A[3] = 0.238497
我得到的结果是：
inv(A)[0] = 4.213983
inv(A)[1] = -
浏览 3提问于2012-07-20得票数 0
回答已采纳
1回答
我想求一个矩阵q+1e-5*np.ye( d) (size d×d)的逆，并使用下面的代码来获得近似结果。
Q = X.dot(X.T) # X is a large sparse matrix, Q is singular
P = np.linalg.inv(Q+1e-5*np.eye(d))
但我得到了这个：
P=[[ nan nan nan ..., nan nan nan]
[ nan nan nan ..., nan nan nan]
[ nan nan nan ..., nan nan nan]
...,
[ nan nan na
浏览 0提问于2016-11-10得票数 0
2回答
给定一个方阵A作为数组，如下所示
import numpy as np
A = np.array(
[
[1, 2, 3],
[3, 4, 6],
[7, 8, 9],
]
)
当出现以下情况时，NumPy的在内部使用哪种算法
np.linalg.inv(A)
来计算A的矩阵求逆。
特别是，由于矩阵求逆在数值上可能是不稳定的(取决于矩阵的 )，是否存在取决于某些矩阵属性的特殊情况？
浏览 30提问于2021-05-12得票数 1
回答已采纳
我想求一个矩阵的逆矩阵。
我知道这首先涉及LU分解，然后是求逆步骤，但通过搜索apple的10.7文档，我找不到所需的函数！
这似乎是一个有用的post ，指出应该使用sgetrf_和sgetri_函数。但是，在搜索这些术语时，我在Xcode文档中什么也没有找到。
有没有人有这个矩阵运算的样板代码？
浏览 1提问于2012-07-01得票数 5
回答已采纳
1回答
避免矩阵求逆
、、
我想避免计算矩阵的逆。我是用C语言写代码的，所以我不能使用MATLAB命令。我想执行以下计算：
v = bk/G;
bk是2×1矩阵，G是2×2矩阵，v是浮点变量。
我目前正在这样做：
V = G^-1 * bk;
这样，我必须使用除法运算符来查找G^-1。我想避免使用除法。
浏览 3提问于2015-09-25得票数 0
2回答
据我所知，分解/分解(LU，QR，Cholesky等)是必需的，然后是基于因式分解的矩阵逆计算。有没有其他方法可以绕过它(我想知道在的试用版中，我是否可以继续使用免费提供的6个函数)？提前感谢你的评论。
浏览 3提问于2009-11-05得票数 2
回答已采纳
1回答
作为用于数值计算的python代码的一部分，我必须多次反演一些较大的(稀疏)矩阵(~100x100)。我真的很想加快程序的速度，我建议的方法之一就是调用C中的一个子程序来进行矩阵反演。
对于此任务，是否有任何推荐的、高效的、经过良好测试的C例程？
谢谢。
浏览 2提问于2012-06-25得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
腾讯会议活动推荐
腾讯云开发者
