3回答
我在想一个算法来解决与p素数的同余ax = 1 mod p。我在考虑用费马定理。因为我知道
a ^ (p-1) = 1 mod p
那就是
a ^ (p-1) = a * (a ^ (p-2))
这意味着a ^ (p-2) mod p是解决方案。不幸的是，这个解决方案，虽然在数学上是正确的，但对计算机来说并不好，因为对于大素数，我必须做a ^ (p-2)，这通常是不可计算的。
哪种算法对计算机科学有好处？
浏览 0提问于2012-12-31得票数 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
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回答
我正试图编写一个代码，使用翻版库在C中反演一个复杂的矩阵，但是，我遇到了一个分割错误，这似乎取决于矩阵的大小N。更重要的是，每次我编译程序或触摸任何东西时，发生分段错误的大小都是不同的。这使我认为代码试图访问分配不良或被禁止的内存的某个地方。不幸的是，我不明白这是如何发生的，因为它似乎与LAPACKE函数主题相关。实际上，当函数/*MatrixComplexInv(invA,A,N);*/ (其中调用LAPACKE函数进行反演)被注释时，分割错误就不会发生。下面是可以自己编译和运行的工作代码。
#include <stdio.h>
#include <lapacke.h&g
浏览 7提问于2019-02-23得票数 1
回答已采纳
3回答
我想取整数的模逆(k≥1)，然后将结果乘以另一个整数，如下表达式所解释：
result=((x^(-k)))*y mod z
我如何实现这个表达式，其中k≥1？
浏览 0提问于2019-05-28得票数 3
回答已采纳
3回答
我想为大整数找到(n选择r)，并且我还必须找出该数字的mod。
long long int choose(int a,int b)
{
if (b > a)
return (-1);
if(b==0 || a==1 || b==a)
return(1);
else
{
long long int r = ((choose(a-1,b))%10000007+(choose(a-1,b- 1))%10000007)%10000007;
return r;
}
}
我正在使用这段代码，但我
浏览 0提问于2014-06-27得票数 2
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
回答已采纳
2回答
是否有一个线性О(n+m)时间算法仅求给定图G(V,E)的最小生成树的值r？我们不想找到MST，只是它的边之和。
我一直在寻找问题的解决方案，但是Kruskal和Prim的算法由于其使用的比较结构(UnionFind(Kruskal)和PQ(Prim))而具有更高的复杂性。此外，他们还找到了MST，这是不想要的，也许还有更快的方法只找到r。
浏览 6提问于2016-02-04得票数 1
回答已采纳
2回答
我有很多需要在Matlab中求逆的大矩阵(大约5000 x 5000)。我实际上需要反转，所以我不能使用mldivide，这对于求解一个b的Ax=b来说要快得多。
我的矩阵来自一个问题，这意味着它们有一些很好的性质。首先，它们的行列式是1，所以它们绝对是可逆的。但是，它们是不可对角化的，或者我会尝试将它们对角化，颠倒它们，然后再将它们放回原处。它们的条目都是实数(实际上是有理的)。
我使用Matlab来获得这些矩阵，对于这些东西，我需要对它们的求逆进行处理，所以我更喜欢一种加速Matlab的方法。但是，如果我可以使用另一种更快的语言，请让我知道。我不知道很多其他的语言(只懂一点C和一点不懂Ja
浏览 2提问于2011-06-28得票数 13
回答已采纳
1回答
我有一个矩阵M，它有576行，15 columns.First列都是1。第二列从1开始，逐个增加到576。第三列是第二列的正方形。第四列是第二列的三次方。第五列是第二列的四次方。它是这样的。最后，第15列是第二列的14次方。该操作是将矩阵M的转置与其自身相乘，并取其结果的逆。然而，Matlab给出了这样的警告:警告:矩阵接近于单数或缩放不良。结果可能不准确。RCOND = 1.889839e-79。这是一种与之相反的方式吗？我也尝试了这个操作符，\但同样的警告也出现了。
浏览 1提问于2016-10-14得票数 0
1回答
我试图用GMP库找到x的乘积逆；我知道有一个内置函数，但我想写我自己的函数。
这是我扩展的欧几里得算法实现：
void extended_euclides(mpz_t r,mpz_t x,mpz_t y,mpz_t a, mpz_t b){
mpz_t r_2,r_1,x_2,x_1,y_2,y_1,q,tmp;
mpz_inits(r,x,y,NULL);
mpz_init_set(r_2,a);
mpz_init_set(r_1,b);
mpz_init_set_ui(x_2,1);
mpz_init_set_ui(y_2,0);
浏览 4提问于2015-11-10得票数 1
回答已采纳
2回答
下面的级数的总和是无穷大的：
1 1/2 1/3 1/4 1/5 ...
根据某位科学家的解释，无穷大是任何一点都不存在的点，即inf +x= inf或inf ~ (inf + x) = 0。因此，在这一理论的基础上，使用了以下算法：
float sum=0.0;
for(int i=0;;i++){
if((sum+(1.0/i))==sum)
break;
sum+=(1.0/i);
}
/* print the value of sum */
该算法在C和JAVA中运行，并以inf的形式给出输出。用C和Java编写的打印语句分别为
printf("
浏览 3提问于2016-03-08得票数 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
回答已采纳
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回答
对于一个主题，我使用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
我试图重新创建一个RSA加密程序，我已经在Java中制作成Erlang。不过，我不知道如何生成私钥。我最初的java代码：
privateKey = publicKey.modInverse(phi);
我找不到任何通用的算法来在网上找到'd‘或私钥。大多数功能小的简单方程不能在更大的问题上实现，或者本教程只给出私钥而不解释过程。有人能指点我的方向吗?这样我就能学会生成私钥了？或者，如果Erlang中有一个模反函数，请说明它是什么。
提前谢谢你
编辑:实际要解的方程是e*d mod ( phi ) = 1，其中e是公钥，d是私钥，phi=p1。当所有其他变量都已知时，我很想为d
浏览 3提问于2013-10-31得票数 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
回答已采纳
我需要逆大矩阵，我想修改我当前的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
回答已采纳
我正在尝试对大于10,000 x 10,000的矩阵执行求逆。
InverterTask<Double> matrixInverter = InverterTask.PRIMITIVE.make(storeM);
try{
storeI = matrixInverter.invert(storeM);
}catch (RecoverableCondition e){
throw new RuntimeException(e);
}
storeM是一个10,000 x 10,000大小的矩阵。
但是，我遇到了以下错误：
Exception in thr
浏览 4提问于2019-05-03得票数 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
回答已采纳
1回答
我正在尝试学习凯克卡克哈希函数的密码分析。关于零和判别器的一篇论文讨论了逆度为3的Keccak。
我无法弄清楚:你怎么得到3的逆度？例如，Keccak中的一个排列Chi在GF(2)中进行乘法。所以我猜这意味着，Keccak的图像变换有2级，但我不知道它们是如何推导出3级的，用于前像、逆或者向后。
如果有人能向我指出一个资源，这将是很大的帮助，因为我找不到一个。
浏览 0提问于2013-06-25得票数 6
回答已采纳
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回答
如何求几乎奇异矩阵的逆？
、、、
我正在用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回答
求满足模数的最小值
、、、、
我遇到的问题是x= (16807 X)% 65536
ie 16807k≡x (mod 65536)
我需要计算，知道x。到目前为止，我最大的努力是一种野蛮的力量。有计算k的数学方法吗？如果不是，我希望对我当前的代码进行任何优化。
t = x;
while ( t += 15115 ) // 16807k = 65536n + x - this is the n
{
if (t%16807 == 0)
return t/16807;
}
return x;
编辑:将+=更改为15115
浏览 0提问于2014-02-20得票数 1
回答已采纳
找到一个更简单的函数的反演是很简单的。我通过翻转方程中的x和y来找到这样做的方法，并对y进行求解，但我被困在了某个部分。
Y= (6*x) mod 13
X= (6*y) mod 13
浏览 1提问于2018-05-06得票数 0
回答已采纳
标题已经解释过了。如果我需要找一个矩阵的逆，有什么理由我应该使用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
回答已采纳
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回答
我指的是InterviewBit上的问题“带重复的排序排列”，我的解决方案可以产生正确的输出，除了长字符串值。这通常是由大量的阶乘溢出造成的。
我已经通过使用Java的BigInteger数学类来解决问题，但是解决方案提示建议使用“模乘逆”替代计算(N-1)！/ (p1！* p2！* p3！！)其中p1、p2和p3是字符串中重复字符的频率。
那么，我的问题是，“模乘逆”如何帮助解决不符合整数基元类型的大阶乘值，以及它背后的数学直觉是什么？我知道如何解决这个编程问题，但唯一阻止成功提交的部分是长字符串值。
非常感谢您对此的任何解释！我的解决方案如下所示，无需使用BigInteger类。
publ
浏览 1提问于2016-10-01得票数 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
回答已采纳
3回答
矩阵模板库矩阵求逆
、、、、
我正在尝试用版本Boost boost_1_37_0和MTLmtl4-alpha-1-r6418来求逆矩阵。我好像找不到矩阵求逆的代码。我在谷歌上搜索了一些例子，它们似乎引用了上面的版本中似乎缺少的lu.h。有什么提示吗？
建议复制lu.h，但这似乎来自MTL2而不是MTL4。我在用VS05或更高版本的MTL2编译时遇到了问题。
那么，你知道如何在MTL4中进行矩阵求逆吗？
更新:我想我更了解马特了，我要去了。
浏览 6提问于2008-12-03得票数 3
回答已采纳
1回答
所以我在循环中嵌入了一个循环：
int a,b,n;
for (a = 1; a <=n; a++) {
for (b = 0; b < n; b+=a)
cout << "hey" << endl;
}
N是2的幂
我试着去理解如何计算时间复杂度，但是我很难弄清楚这其中的大θ符号。
我知道外部循环在O(n)时间内运行，但是由于b+=a的原因，我不确定内环，我知道如果我有两个循环的时间，我可以把它们乘以得到函数的大θ时间，但是我不知道内环运行在什么位置。
当我插入样本n's时。2，4，8，16)，然后内环分别环
浏览 2提问于2014-09-16得票数 2
回答已采纳
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回答
我有两个稀疏矩阵A(亲和力矩阵)和D(对角矩阵)，维数为100000*100000。我必须计算拉普拉斯矩阵L= D^(-1/2)*A*D^(-1/2)。我正在使用scipy CSR格式的稀疏矩阵。
我没有找到任何方法来求稀疏矩阵的逆。如何求稀疏矩阵的L和逆矩阵？还建议使用python来做这件事是否有效，或者我是否应该调用matlab函数来计算L？
浏览 1提问于2012-02-12得票数 3
我需要反转一个BigInteger。
假设我有BigInteger x;，并且我需要计算x.modPow(new BigInteger("-1"), p)。
我收到以下错误：java.lang.ArithmeticException: BigInteger not invertible。
浏览 0提问于2011-11-30得票数 5
回答已采纳
1回答
我在Google上搜索，找不到用GSL计算矩形矩阵逆的函数。由于很难找到答案，当其他人需要寻找矩形矩阵的逆时，这里的答案将对他们有所帮助。如果不可能使用GSL，那么请提供一些易于使用和提供矩形矩阵逆的可选库。
浏览 8提问于2013-08-14得票数 5
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
回答已采纳
我做了重量有向图，就像这样
6
0 3 4 INFINITY INFINITY INFINITY
INFINITY 0 INFINITY 7 INFINITY INFINITY
INFINITY 3 0 5 11 INFINITY
INFINITY INFINITY INFINITY 0 6 3
INFINITY INFINITY INFINITY INFINITY 0 4
INFINITY INFINITY INFINITY INFINITY INFINITY 0
首先，我使用一些整数值来表示无穷大，比如99或20000。但是当我发现它是错误的时，v5 -> v4必须表示无穷大，但是表
浏览 0提问于2011-11-16得票数 4
回答已采纳
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回答
用于加密的矩阵乘法求逆
、、、、
我正在尝试组合一个加密算法，但我陷入了以下问题，我甚至不知道它应该是这样的或不是这样的！
问题是：
我有16字节的矩阵要乘以16,16矩阵，结果是一个16字节的矩阵。
然后我应该将结果矩阵乘以逆，这里我想我应该得到原始的16字节矩阵(根据算法数据表)。
所以你能告诉我怎样才能找回原始矩阵吗？
提前感谢您的帮助。
致以敬意，
英语。Aws
浏览 2提问于2011-12-17得票数 0
5回答
我编写了一个函数，它接受一个长值的n，并使用它作为要遍历的迭代次数。函数应该给出π的一个很好的估计，然而，大n的所有值都趋向于3.000，而不是3.1415，所以我不知道发生了什么？
我做错什么了吗？
这是我的密码：
double estimate_pi(long long n){
double randomx, randomy, equation, pi;
long long i, incircle = 0;
for(i = 0; i < n; i++){
randomx = (double)(rand() % (1+1-0) + 0);
浏览 1提问于2019-03-13得票数 7
回答已采纳
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
回答已采纳
5回答
大矩阵求逆
、、、、
我正在考虑取一个大矩阵的逆矩阵，通常大小为1000x1000，但有时会超过100000 x 100000 (由于时间和内存的原因，它目前正在失败)。我知道正常的情绪是‘不要反其道而行之，找一些其他的方法来做’，但这在目前是不可能的。这样做的原因是因为使用了已经制作的软件，这些软件希望得到矩阵的逆矩阵。(注意:我正在寻找改变这一点的方法，但这将需要很长时间)
目前，我们正在使用来自数值副本的LU分解方法，我目前正在测试eigen库。eigen库似乎更稳定，速度更快，但我仍处于准确性测试阶段。我已经快速浏览了其他库，如ATLAS和LAPACK，但还没有对它们进行任何实质性的测试。
似乎eigen库
浏览 1提问于2012-06-27得票数 18
回答已采纳
13回答
求矩阵的逆
、、
不知怎么的，我们还没有找到任意大小的方阵的逆的挑战，尽管有3x3和4x4的逆矩阵和更复杂的版本矩阵。
您的任务是，给定一个平方n\times n非奇异矩阵M，输出满足以下条件的矩阵M^{-1}
MM^{-1} = I_n
计算M^{-1}有许多方法和公式，但最著名的方法之一是
M^{-1} = \frac1{\det(M)}\text{ adj}(M)
其中\det表示行列式，\newcommand{\adj}{\text{adj}}\adj表示副词
一些定义：
I_n：n\times n 单位矩阵，即一个n\times n矩阵，其中主导对角线完全由1s和其余的0s组成。
非奇异性:保
浏览 0提问于2020-10-20得票数 29
尝试使用AugPake协议在java中实现BigInteger。
我在计算K=Y^z \mod p时遇到了一些困难，因为z总是用于z={1\over x+(w*r)} \mod q的0。作为1和一个非常大的数字之间的比率，使用BigInteger，z总是为零。
能帮我理解如何计算z吗？也许我误解了一些文件。
谢谢你的回应。
要计算Z，我使用扩展的欧几里德算法: ax + qy = gdc(a，b)
非常感谢
浏览 0提问于2018-01-14得票数 2
2回答
我正在尝试将我的程序从MATLAB转换成Python。在我尝试求逆我的一个矩阵(22x22)之前，一切似乎都运行得很好。对于逆矩阵，我得到不同的值。在MATLAB中，我使用了函数inv(J)和A\b来求逆矩阵，而在python中，我使用了J.I和np.linalg.inv(J)。虽然python中的两个函数返回相同的结果，但它们与我在MATLAB中得到的结果不同。即使行列式的值对于这两种情况也是不同的。MATLAB返回值为10^23，python返回值为10^20。一个矩阵有可能有多个逆数吗？或者这与计算的精度有关吗？
在这个帖子中有一个类似的讨论：。然而，在那次讨论中，我没有找到任何确凿的答
浏览 27提问于2018-07-28得票数 2
3回答
我正在学习“大O”( Big )，并且被一个问题困住了。
问题-
int sample_function (int x) {
if (x<1){
return x;}
int y = x/2;
return sample_function (y) + sample_function (x-y);
}
如果我们调用(n/2)的2个递归调用，那么大O会是什么？我知道除法递归函数(n/2)的大O是O(log(n))，但不确定这个问题。
浏览 3提问于2020-12-11得票数 2
回答已采纳
2回答
据我所知，分解/分解(LU，QR，Cholesky等)是必需的，然后是基于因式分解的矩阵逆计算。有没有其他方法可以绕过它(我想知道在的试用版中，我是否可以继续使用免费提供的6个函数)？提前感谢你的评论。
浏览 3提问于2009-11-05得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
