matrix-inverse相关内容

如果十进制更长，则Fortran-逆矩阵结果不相同

我的真实数据是第一个输入的，但是结果的倒数却是如此之大.当您与第一和第二输入进行比较时，它们是相同的数据.只有十进制大小差异.为什么会有不同的结果?因为它们是相同的数据.他们怎么会有不同的结果?您可以看到结果和输入.太奇怪了. program test Implicit none double precision,allocatable,dimension(:,:) :: ..
发布时间：2020-05-07 19:32:38 其他开发

矩阵求逆R

我想对一个平方对称正定矩阵求逆.我知道R中有两个函数solve()和chol2inv()，但是它们的结果不同.我需要知道为什么会这样吗? 谢谢. 解决方案 对于solve，您需要给出原始矩阵，但对于chol2inv，则使用预先计算的cholesky分解: set.seed(1) a ..
发布时间：2020-05-07 19:06:45 其他开发

如何在不采用矩阵逆的情况下有效计算diag(X％*％solve(A)％*％t(X))?

我需要以下对角线: diag(X %*% solve(A) %*% t(X)) 其中，A是满秩方阵，X是矩形矩阵. A和X都是稀疏的. 我知道找到矩阵的逆是不好的，除非您真的需要它.但是，我看不到如何重写公式，以便用两个参数将solve(A)替换为solve，这样就可以在不显式求逆的情况下求解线性系统.有可能吗? 解决方案 也许在我真正开始之前，我要提一下，如果你这样做 ..
发布时间：2020-05-07 18:48:12 其他开发

用Python执行模块化矩阵求逆的最简单方法?

我想采用Python [[1,2]，[3,4]] mod 7矩阵的模逆.我已经看过numpy(它执行矩阵求逆，但不进行模块化矩阵求逆)，并且在网上看到了一些数字理论软件包，但是似乎没有什么东西可以做这种相对通用的过程(至少，对我来说似乎是相对通用的). /p> 顺便说一下，上述矩阵的逆是[[5,1]，[5,3]](mod 7).我希望Python可以帮我做. 解决方案 当舍入错误不成问题 ..
发布时间：2020-05-07 18:46:37 Python

Matlab:与滞后相关运算有关的混乱

我有一个时间序列模型y(t)= h^T y(t-1) + n(t)，其中n(t)是激发并驱动过程的高斯白噪声. y是t = 1,2,...的线性回归模型的输出，表示数据点的数量. 问题:如果Correlation矩阵为Ryy = E[y(t)*y(t)^T],，则可以计算滞后随机变量(例如 )的Correlation [E[y(t-1)*y(t-1)']] 通常，这些运算符和表达式 ..
发布时间：2020-05-06 15:36:17 其他开发

仅向矩阵的对角元素添加一个小值

我是matlab的新手，我正在尝试找出值非常小的矩阵逆.当我尝试找到逆时，我得到一个错误，说矩阵是奇异的.建议的解决方案之一是尝试向对角线元素添加一些元素.我知道我必须使用眼图和诊断方法，但是我无法提出正确的解决方案. 任何评论都会有所帮助. 解决方案 如果您只想将身份矩阵或其倍数添加到 square 矩阵中，则可以进行 A_new = A_old + k*eye(size(A ..
发布时间：2020-05-06 15:19:13 其他开发

R ginv和Matlab pinv产生不同的结果

与MATLAB pinv()函数相比，R中MASS包中的 ginv()函数产生的值完全不同.他们俩都声称会产生矩阵的Moore-Penrose广义逆. 我尝试为R实现设置相同的公差，但差异仍然存在. MATLAB默认tol:max(size(A)) * norm(A) * eps(class(A)) R默认值tol:sqrt(.Machine$double.eps) 复制: ..
发布时间：2020-05-06 14:25:45 其他开发

改善条件恶劣的矩阵

我有一个条件恶劣的矩阵，其 rcond() 是接近于零，因此，该矩阵的逆并不正确.我尝试使用 pinv() ，但这不能解决问题.这就是我采取的相反方法: X = (A)\(b); 我正在寻找针对此问题的解决方案，并找到了此链接(最后一种解决方案)以改善矩阵.那里的解决方案建议使用此方法: A_new = A_old + c*eye(size(A_old)); c > 0所在的位 ..
发布时间：2020-05-06 12:41:12 其他开发

矩阵求逆还是Cholesky?

我正在开发一种解决 Ax = b 的算法，其中已知 A 和 b . 有两种方法可以执行 x = A -1 b 或使用Cholesky.我知道矩阵将始终是正方形且为正定，尽管det( A )可能为零.在极少数情况下，我可以忽略它.但是从计算和效率的角度来看，创建逆矩阵是否效率太低? 解决方案 通常，您总是想使用求解器.实际求解器的运行速度应与乘以逆运算的速度一样快.与进行分解相比，计算 ..
发布时间：2020-05-06 11:09:32 其他开发

2x2矩阵的数值稳定逆

在我使用C进行的数值求解器中，我需要求一个2x2矩阵的倒数，然后在右边将其乘以另一个矩阵: C = B . inv(A) 我一直在使用以下2x2逆矩阵的定义: a = A[0][0]; b = A[0][1]; c = A[1][0]; d = A[1][1]; invA[0][0] = d/(a*d-b*c); invA[0][1] = -b/(a*d-b*c); invA[1] ..
发布时间：2020-05-06 11:08:22 Python

在3D theano张量上广播linalg.pinv

在下面的示例中，有一个大小为(4，3，3)的3d numpy矩阵+一种关于如何计算numpy中3 * 3矩阵中的4个矩阵的每一个的pinv的解决方案.我也尝试使用在numpy中使用的相同功能，theano希望实现相同，但是失败了.知道如何在theano中做到这一点吗? dt = np.dtype(np.float32) a=[[[12,3,1], [2,4,1], [2,4,2 ..
发布时间：2020-04-30 12:10:50 Python

逆矩阵的奇异性

作为数据，我得到一个矩阵A，但是在我的算法中，我需要对其求逆.我要做的是: C = inv(A) + B; 然后在另一行中更新A.在接下来的循环中，对于该算法，我同样需要(更新)逆.等等.在以后的周期中，我得到以下信息: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1 ..
发布时间：2020-04-30 12:10:38 其他开发

PHP矩阵逆

我看到了这个问题，并弹出了这个想法. 在PHP中是否有一种有效的方法? 编辑 最适合演示吗? 解决方案 您可以使用pear软件包 Math_Matrix . ..
发布时间：2020-04-30 12:08:22 PHP

Matlab逆运算和警告

不太清楚这意味着什么. “警告:矩阵对于工作精度至关重要." 我有一个3x4矩阵，称为矩阵bestM 矩阵Q是bestM的3x3，矩阵m是bestM的最后一列 我想做C =-Q *矩阵m的逆矩阵 我得到那个警告 和C = [Inf Inf Inf]这是不对的，因为我正在计算世界上的相机中心 bestM = [-0.0031 -0.0002 0.0005 0.9788; ..
发布时间：2020-04-30 12:06:44 其他开发

MATLAB中最有效的矩阵求逆

在MATLAB中计算某些方阵A的逆时，使用 Ai = inv(A) % should be the same as: Ai = A^-1 MATLAB通常会通知我这不是最有效的反相方法. 那么更有效率吗?如果我有一个方程式系统，则可能使用/，\运算符. 但是有时候我需要逆运算来进行其他计算. 最有效的反转方法是什么? 解决方案 我建议使用svd(除非您确实完全确定您的矩阵没 ..
发布时间：2020-04-30 12:05:32 其他开发

python中稀疏矩阵的伪逆

我正在处理来自神经影像的数据，由于数据量很大，我想在我的代码(scipy.sparse.lil_matrix或csr_matrix)中使用稀疏矩阵. 尤其是，我将需要计算矩阵的伪逆来解决最小二乘问题. 我已经找到了sparse.lsqr方法，但是它不是很有效.有没有一种方法可以计算出Moore-Penrose的伪逆(对于正常矩阵，它对应于pinv). 我的矩阵A的大小约为600'000 ..
发布时间：2020-04-30 12:04:22 Python