在MATLAB中，可以使用hack来强制内置LU分解算法不使用置换矩阵(即，强制P=I)，如here所示。 python中有没有类似的黑客攻击？也就是说，有没有办法强制scipy.linalg.lu (或任何其他流行的LU算法)不使用置换矩阵？
浏览 16提问于2019-04-17得票数 2
1回答
鲤鱼C++ LU分解
、、、、
因为我要解不同的线性系统在A总是一样，b总是变化的情况下，我只想对A进行LU分解一次，并在不同的b中重用LU分解。不幸的是，我不知道如何在Armadillo中执行这种操作。我所做的只是A矩阵的LU分解：// ... fill the A matrix
浏览 1提问于2013-07-30得票数 2
我们知道消除大致需要1/3 n^3操作，如果我们使用存储在内存中的LU分解，它将简化为n^2操作。如果我们有一个带上对角线和下对角线的带状矩阵，我们可以跳过零，把它降到大约nw^2操作，如果我们使用LU分解，它可以在大约2nw操作中完成。在scipy.linalg中，我们有lu_factor和lu_solve，但它们似乎没有针对频带矩阵进行优化。我们也有solve_banded，但它直接解
浏览 42提问于2019-01-14得票数 0
1回答
我想计算A的LU分解。常用算法效率不高，因为它们必须循环矩阵的所有系数。是否有一种有效的算法可以利用CSR格式来计算LU分解？谢谢你的帮助!
浏览 3提问于2013-05-11得票数 0
在R核中没有LU分解函数。虽然这类分解是solve的一个步骤，但它并不是作为独立函数显式提供的。我们能为这个写一个R函数吗？它需要模拟LAPACK例程。这个函数对于教育和调试都是有用的。教育的好处是显而易见的，因为我们可以一列一列地说明因式分解/高斯消除。为了调试使用，以下是两个示例。
在中，有人问为什么R和Python中的LU
浏览 1提问于2018-08-04得票数 7
回答已采纳
3回答
根据我对LU分解的理解，这意味着对于下三角矩阵L和上三角矩阵U，矩阵A可以写为A= LU。然而，scipy中与LU分解(，，)相关的函数似乎涉及第三个矩阵P，使得A= PLU，P是置换矩阵( L，U与前面一样)。
这个排列矩阵的意义是什么？如果“真正的”LU分解总是可能的，为什么P不是单位矩阵呢？
浏览 0提问于2014-04-22得票数 3
2回答
我试图理解使用numpy和and库进行LU分解的必要性。所以，我在蟒蛇中厌倦了这个想法，用的是numpy和scipy。dev.7次运行中，每一次循环100000次)
lu, piv = linalg.lu_factor(A)编辑
浏览 0提问于2018-04-07得票数 0
传统智慧指出，如果您使用相同的A和不同的b多次求解Ax = b，则应该对LU使用LU分解。u, x) 这最终比仅仅使用 x = scipy.linalg.solve(A,b) 在循环中。scipy.linalg.solve()不是针对使用向前和向后替换的上对角线系统和下对角线系统进行了优化吗？或者，有没有可能存在一些编译技巧，其中python认识到它可以为scipy.linalg.solve部分进行LU分解<
浏览 45提问于2020-10-22得票数 2
回答已采纳
3回答
矩阵模板库矩阵求逆
、、、、
我好像找不到矩阵求逆的代码。我在谷歌上搜索了一些例子，它们似乎引用了上面的版本中似乎缺少的lu.h。有什么提示吗？那么，你知道如何在MTL4中进行矩阵求逆吗？
更新:我想我更了解马特了，我要去了。
浏览 6提问于2008-12-03得票数 3
回答已采纳
我一直试图开发一个程序来解决一个系统的Ax=b的平方矩阵A使用LU分解。但是，我意识到这种分解并不总是存在的(一种判断方法是，如果不需要行交换操作，那么就存在)。但是，我从许多来源中看到，这是计算Ax=b解决方案的一个很好的方法。
我的问题是:遇到没有LU分解的矩阵的频率是多少次？如果遇到这样的矩阵，他应该如何
浏览 1提问于2019-03-12得票数 0
回答已采纳
我希望我只是错过了一个简单的矩阵运算技巧，但问题是我所能访问的只是一个LU求解器(Matlab *或SuperLU)，我需要一个对称矩阵A的低密度脂蛋白分解，所以我想“没问题”，因为LU分解是唯一的，低密度脂蛋白分解是唯一的，那么D就是U的对角线。但问题是，它永远不会A= LU，而是PA = LU，PA不对称！所以我不知道如何从PA =
浏览 4提问于2017-10-12得票数 0
1回答
这个论坛的新成员。我试图运行八度的LU分解函数，其完全旋转如下：对于矩阵A，我一直得到这样的错误：
元素4是列排列的矩阵Q，发生了什么？
浏览 7提问于2020-01-14得票数 0
我想使用一些cula功能，如LU分解或矩阵逆，但我有一些关于指针输入的问题。例如，对于使用scikits.cuda.cula.culaDeviceSgetrf(m，n，a，lda，ipiv)进行LU分解，需要使用指针f "a“参数，但在python中没有指针(我知道python中的所有变量都是by ref)。那么在这种情况下我该怎么办呢？我应该使用ctype库来创建python吗？&quo
浏览 3提问于2012-08-31得票数 0
回答已采纳
我正在寻找最精确的C类实现来求解复杂的浮子型线性方程组(又名np.complex64)。都是为了LU分解。我想首先在numpy中实现它，确保一切顺利，然后将实现转换为C。这就是我现在得到的：import numpy as np
L, U = np.empty_like((la.lu_factor(A), b)np
浏览 1提问于2019-11-28得票数 1
回答已采纳
2回答
我正在将一些MATLAB代码转换为Python，并观察到\运算符和scipy.linalg.lstsq之间存在很大的数值差异，这显然是可互换的。在我的代码中，我计算了一些矩阵的LU分解，但是Python和Matlab对'L‘给出了稍微不同的答案。-1.0000 0; 0 -0.9350
浏览 5提问于2015-01-30得票数 0
回答已采纳
10回答
求4x4矩阵的逆矩阵
、、、、
我正在寻找一个关于如何反转4x4矩阵的示例代码实现。我知道有高斯元素，LU分解等，但我并没有详细研究它们，我只是在寻找代码来做这件事。语言理想情况下是C++，数据以列为主的16个浮点数组的形式提供。
浏览 8提问于2009-07-18得票数 90
回答已采纳
4回答
我有一个使用的程序，需要测试的是一个矩阵是否可以求逆。我知道我可以尝试一下并捕获一个异常，但这似乎不是一个好主意(将catch块作为“正常”代码路径的一部分似乎是不好的形式)。如果测试也返回逆操作(或者在比逆操作更好的O()中运行)，则更可取。
浏览 2提问于2009-01-27得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
