大矩阵求逆方法
大矩阵求逆方法是指求解大矩阵的逆矩阵的算法。在数学中，矩阵的逆矩阵是指一个矩阵乘以它的逆矩阵等于单位矩阵。对于一个方阵，如果它是可逆的，那么它有一个唯一的逆矩阵。对于一个大矩阵，求解逆矩阵的方法有很多种，下面介绍几种常用的方法：
- LU分解法：LU分解是将一个矩阵分解为下三角矩阵L和上三角矩阵U的乘积。然后通过求解L的逆矩阵和U的逆矩阵，最后将它们相乘得到原矩阵的逆矩阵。
- 高斯消元法：高斯消元法是通过对矩阵进行行变换，将矩阵转换为阶梯形式，然后通过逆变换将矩阵的逆矩阵求解出来。
- SVD分解法：SVD分解是将一个矩阵分解为UΣV^T的形式，其中U和V是正交矩阵，Σ是对角矩阵。然后通过求解Σ的逆矩阵和V的逆矩阵，最后将它们相乘得到原矩阵的逆矩阵。
以上是常用的大矩阵求逆方法，但是对于特殊类型的矩阵，还有其他的求逆方法。需要根据具体情况选择合适的方法。
相关·内容
1回答
我必须对一个大的稀疏矩阵求逆。我无法摆脱矩阵求逆，唯一的捷径就是了解主对角线元素的概念，而忽略非对角线元素(我不希望这样做，但作为一种解决方案，它是可以接受的)。
我需要求逆的矩阵通常很大(40000 *40000)，并且只有一小部分非零对角线。我目前的方法是构建稀疏的所有内容，然后
posterior_covar = np.linalg.inv ( hessian.todense() )
这显然需要很长的时间和大量的内存。
有什么提示，或者这只是一个耐心的问题，或者让问题变得更小？
浏览 0提问于2013-02-28得票数 6
回答已采纳
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<- 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
回答已采纳
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
1回答
我有两个稀疏矩阵A(亲和力矩阵)和D(对角矩阵)，维数为100000*100000。我必须计算拉普拉斯矩阵L= D^(-1/2)*A*D^(-1/2)。我正在使用scipy CSR格式的稀疏矩阵。
我没有找到任何方法来求稀疏矩阵的逆。如何求稀疏矩阵的L和逆矩阵？还建议使用python来做这件事是否有效，或者我是否应该调用matlab函数来计算L？
浏览 1提问于2012-02-12得票数 3
1回答
用于加密的矩阵乘法求逆
、、、、
我正在尝试组合一个加密算法，但我陷入了以下问题，我甚至不知道它应该是这样的或不是这样的！
问题是：
我有16字节的矩阵要乘以16,16矩阵，结果是一个16字节的矩阵。
然后我应该将结果矩阵乘以逆，这里我想我应该得到原始的16字节矩阵(根据算法数据表)。
所以你能告诉我怎样才能找回原始矩阵吗？
提前感谢您的帮助。
致以敬意，
英语。Aws
浏览 2提问于2011-12-17得票数 0
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
回答已采纳
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
回答已采纳
3回答
矩阵模板库矩阵求逆
、、、、
我正在尝试用版本Boost boost_1_37_0和MTLmtl4-alpha-1-r6418来求逆矩阵。我好像找不到矩阵求逆的代码。我在谷歌上搜索了一些例子，它们似乎引用了上面的版本中似乎缺少的lu.h。有什么提示吗？
建议复制lu.h，但这似乎来自MTL2而不是MTL4。我在用VS05或更高版本的MTL2编译时遇到了问题。
那么，你知道如何在MTL4中进行矩阵求逆吗？
更新:我想我更了解马特了，我要去了。
浏览 6提问于2008-12-03得票数 3
回答已采纳
2回答
如何求距离矩阵的逆矩阵？
、、、、
我有一个巨大的距离矩阵。 示例：(10000 * 10000).. 有没有一种有效的方法来求逆矩阵？ 我试过numpy的Inv()，但它太慢了。 有没有更有效的方法？
浏览 30提问于2020-10-12得票数 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
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
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
回答已采纳
4回答
这是我试图解决的方程式：
H= (X'*X)^-1*X'*y
其中X是一个矩阵，y是一个向量((X' X) ^-1是X-转置时间X的逆)。我已经在Matlab中将其编码为：
h = (X'*X)\X'*y
我相信这是正确的。问题是X大约是10000x10000，而试图计算这个倒数会使Matlab崩溃，即使是我能找到的最强大的计算机(16核，24 16)。有没有办法把它拆分开，或者设计一个库来做这么大的倒排？
谢谢。
浏览 0提问于2013-03-13得票数 4
回答已采纳
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回答
该任务要求我生成一个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
2回答
我有一个nxn矩阵C，并使用numpy.linalg中的inv求逆以获得Cinverse。我的C矩阵具有10**4阶的元素，但是我的Cinverse矩阵具有10**12阶或更高阶的元素(不确定是否正确)。当我执行numpyp.dot(C,Cinverse)时，我得不到单位矩阵。为什么会这样呢？
我有一个向量x，我将它自身相乘，得到一个矩阵。
x=array([ 121.41191662, 74.22830468, 73.23156336, 75.48354975,
79.89580817])
c=np.outer(xvector,xvector)
这是一个5x5矩阵。
然
浏览 0提问于2017-02-15得票数 0
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
我试图用一些符号变量来求逆一个矩阵，但是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回答
实时求逆矩阵
、、、
我一直在开发一个实时工作的C语言控制软件。软件实现了被控系统的离散状态空间观测器等。对于观察者的实现，需要计算4x4维的矩阵的逆。逆矩阵计算必须每隔50微秒进行一次，值得一提的是，在这段时间内还将进行其他非常耗时的计算。因此，逆矩阵计算所消耗的时间必须远远少于50微秒。还有必要说明的是，所使用的DSP没有支持浮点运算的ALU。 我一直在寻找一些有效的方法来做到这一点。我的一个想法是准备计算矩阵4x4的行列式的通用公式和计算4x4矩阵的伴随矩阵的通用公式，然后根据下面给定的公式计算矩阵的逆矩阵。 ? 你对这种方法有什么看法？
浏览 46提问于2020-09-10得票数 1
我需要逆大矩阵，我想修改我当前的LAPACKE版本例程，以便利用GPU NVIDIA卡的强大功能。
事实上，我的LAPACKE例程对于相对较小的矩阵很好地工作，但对大型矩阵却不起作用。
下面是这个LAPACKE例程的实现：
#include <mkl.h>
// Passing Matrixes by Reference
void matrix_inverse_lapack(vector<vector<double>> const &F_matrix, vector<vector<double>> &F_outpu
浏览 8提问于2021-08-14得票数 0
回答已采纳
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回答
我有计算下三角矩阵求逆的代码。如何通过对下面的代码稍加修改来计算上三角矩阵的求逆？ 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回答
我正在通过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
回答已采纳
6回答
矩阵求逆
、、
我正在为class开发一个程序，它读取一个矩阵，然后用单位矩阵打印出来。然后，我必须一步一步地将其简化为其单位矩阵，同时将相同的计算应用于单位矩阵，从而得到该矩阵的逆。
我被代码的行减少部分卡住了。我正在尝试创建一个循环，它将一个nxn矩阵简化为一个单位矩阵，而不是为每一列创建一个。我需要一些帮助的循环，行减少矩阵一列一列。
任何人都有任何想法。(我是编程的初学者，我使用C语言，发布的解决方案也可以，但我正在寻找一种算法，可以一次减少所有行)。
浏览 0提问于2009-11-06得票数 0
2回答
我有很多需要在Matlab中求逆的大矩阵(大约5000 x 5000)。我实际上需要反转，所以我不能使用mldivide，这对于求解一个b的Ax=b来说要快得多。
我的矩阵来自一个问题，这意味着它们有一些很好的性质。首先，它们的行列式是1，所以它们绝对是可逆的。但是，它们是不可对角化的，或者我会尝试将它们对角化，颠倒它们，然后再将它们放回原处。它们的条目都是实数(实际上是有理的)。
我使用Matlab来获得这些矩阵，对于这些东西，我需要对它们的求逆进行处理，所以我更喜欢一种加速Matlab的方法。但是，如果我可以使用另一种更快的语言，请让我知道。我不知道很多其他的语言(只懂一点C和一点不懂Ja
浏览 2提问于2011-06-28得票数 13
回答已采纳
1回答
我有一个通用矩阵A，它是对称的，正定的，并且是稀疏填充的(它也很大，比如由几十到几十万行组成)。我想让一个神经网络学习如何找到一个矩阵B，它类似于A的逆(即，乘以A将使光谱半径的最小扩展)尽可能准确和快速。
我想知道是否有一些人工智能分支专注于这类问题，或者如果没有，我如何处理它。
我知道用来做这件事的经典算法，但我只想看看机器学习在这个问题上提供的可能性。
浏览 0提问于2016-10-19得票数 11
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回答
我正在用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
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
回答已采纳
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
回答已采纳
我正在寻找计算矩阵求逆的最有效的方法。我知道这与LU分解有关，或者是把它分成两个矩阵，上三角和下三角矩阵，以及这两个系统的求解，但我不能用MATLAB写出代码。我甚至不知道从哪里开始，也不知道如何在我的问题中使用系统的结果。我真的需要它，因为我必须求解一个10.000 x 10.000的矩阵，而这个简单的方法永远需要花费很长时间。
浏览 1提问于2014-03-23得票数 0
5回答
大矩阵求逆
、、、、
我正在考虑取一个大矩阵的逆矩阵，通常大小为1000x1000，但有时会超过100000 x 100000 (由于时间和内存的原因，它目前正在失败)。我知道正常的情绪是‘不要反其道而行之，找一些其他的方法来做’，但这在目前是不可能的。这样做的原因是因为使用了已经制作的软件，这些软件希望得到矩阵的逆矩阵。(注意:我正在寻找改变这一点的方法，但这将需要很长时间)
目前，我们正在使用来自数值副本的LU分解方法，我目前正在测试eigen库。eigen库似乎更稳定，速度更快，但我仍处于准确性测试阶段。我已经快速浏览了其他库，如ATLAS和LAPACK，但还没有对它们进行任何实质性的测试。
似乎eigen库
浏览 1提问于2012-06-27得票数 18
回答已采纳
我需要取一个矩阵的逆矩阵来拟合我在Rust中的数据的线性回归。我已经尝试过使用ndarray来做这件事，但是没有找到一个不使用BLAS后端的好方法。有什么建议吗?如何在不使用BLAS后端的情况下在Rust中找到矩阵的逆？
浏览 34提问于2020-12-18得票数 0
1回答
我正在试着求4x4矩阵的逆，我看过一些用几十次计算来计算行列式的例子，这对我来说似乎不起作用，我可能在某个地方搞砸了，但我一直在试图找到另一种方法来解决这个问题，我的答案是：
public Matrix inverse()
{
double[] array = new double[]{
1/m[0], m[4], m[8], 0,
m[1], 1/m[5], m[9], 0,
m[2], m[6], 1/m[10], 0,
-m[12], -m[13], -m[14],
浏览 2提问于2013-05-04得票数 0
回答已采纳
2回答
矩阵求逆R
、、
我想求一个平方对称正定矩阵的逆矩阵。我知道R中有两个函数solve()和chol2inv()，但它们的结果是不同的。我想知道为什么会这样？
谢谢。
浏览 1提问于2013-03-11得票数 2
回答已采纳
我有一个200000 x 200000大小的大矩阵。我需要计算它的逆。但是它在使用numpy.linalg.inv时出现内存错误。有什么方法可以计算大型矩阵的逆。
浏览 0提问于2019-03-23得票数 1
2回答
矩阵求逆GPU MATLAB
、、、、
我正在寻找一种方法来找到一个矩阵的伪逆，这样它就可以在GPU上完成。SVD/QR很难并行化，并且不被MATLAB的GPU支持，但似乎LU虽然可以并行运行，但也不被MATLAB的GPU支持。我比较了它的性能，它似乎比在单核CPU上运行要慢。
我正在寻找一个我可以使用的伪逆(甚至是方阵的正则逆)。根据Matlab，使用mldivide ()执行适用于GPU的高斯消除。
我尝试使用A\I，但不幸的是它不能在GPU上高效运行。
有没有人能帮我找到一个优化的并行逻辑单元或高斯消除的代码？
我听说过MAGMA包，但它似乎有很多工作要安装和编译，我真的需要这个简单的东西。
C++代码也是受欢迎的。
谢谢，吉
浏览 2提问于2013-06-22得票数 1
回答已采纳
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回答
我有一个非常大的吸收马尔可夫链(规模到问题大小--从10个州到数百万个州)，这是非常稀疏的(大多数州只能对其他4或5个状态做出反应)。
我需要计算这条链的基本矩阵的一行(每种状态的平均频率，给定一个启动状态)。
通常，我会通过计算(I - Q)^(-1)来做到这一点，但是我还没有找到一个好的库来实现稀疏矩阵逆算法！我看过几篇关于它的论文，大部分都是P.h.D级的。
我的大部分谷歌结果指向我的帖子，谈论如何不使用矩阵逆时，解决线性(或非线性)方程组.我觉得这没什么特别的帮助。基本矩阵的计算是否类似于求解方程组，而我只是不知道如何用另一种形式来表示一个？
因此，我提出两个具体问题：
：计算稀疏矩阵
浏览 6提问于2012-07-29得票数 14
回答已采纳
标题已经解释过了。如果我需要找一个矩阵的逆，有什么理由我应该使用A\eye(size(A))而不是inv(A)吗？在你问:是的，我真的需要逆，不只是计算。
PS：
isequal(inv(A), A\eye(size(A)))
ans =
0
那么哪一个更精确？
UPDATE：这个问题被关闭了，因为它看起来是"“问题的重复。这个问题在这里有很大的不同，因为它不是求解函数inv的速度，也不是函数的精度，而是inv和.\eye在计算矩阵的真逆方面的差异。
浏览 4提问于2016-01-27得票数 5
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
