问题标签 [determinants]

当我在 Matlab 中为行列式运行这个拉普拉斯展开时，我确实得到“未定义的函数或变量 'A'”。如果有人可以帮助我纠正它并使其可行，我会很高兴。

谢谢

所以我采用矩阵的行列式，然后尝试使用 Sympy 中的 Solver 来求解特定变量的表达式。

我试图解决的决定因素非常复杂，但其中只有一个变量。作为一个长表达式，我不想将其全部粘贴，但我已将其缩短为给出相同结果的片段；

问题是，当我打印 Solutions 时，我得到了一个关于 Nperp 的表达式，而不是我想要的数值解。

我不确定问题是否在于 Sympy 无法处理多项式中的高幂，或者是否可能没有数值解，但我会感谢比我更有知识的人的一些想法。

谢谢！

编辑：代码未缩进

我有n行向量，我可以选择k向量来形成一个k*k矩阵。然后我取每个矩阵的行列式的绝对值。因此，我有(n choose k)可能的绝对值行列式。

我希望选择产生最高行列式的行向量。是否可以优化搜索以快速完成此操作？

我目前正在使用n ~= 10,000and k = 31，但我将使用各种设置。

我用 python 和 c++ 编写代码

我的程序由两个类组成。第一个创建一个二维数组并用用户输入填充它。第一个类工作正常，当我在 main 中调用它时，它能够创建和打印一个二维数组。但是，我试图将指向该二维数组的指针传递给第二个函数来计算矩阵的行列式。但是，在调用行列式函数后，我的程序不断崩溃。为什么我无法对这些数组元素进行乘法、加法或减法？

这是行列式类的实现文件：

以下是我调用该函数的 main 部分：

我观察到，当我通过 N < 10 阶矩阵的克莱默规则（两个行列式的商）求解线性方程组时，与 LAPACK 解相比，我得到相当大的残差。

这是一个例子：

现在linsolving这个，产生 LAPACK（来自 Intel MKL）：

x = [-0.000314947 -0.000589154 -0.00587876 0.0184799 0.01738 -0.0170484]

和克莱默规则（自己的实现）产生：

x = [-0.000314933 -0.000798058 -0.00587888 0.0184808 0.017381 -0.0170508]

注意x[1]区别。

我可以保证我的行列式计算是正确的。有没有人提出过类似的观察或可以说明这一点？

我需要帮助通过实验确定矩阵 nxn 行列式的计算复杂性

我的代码：

但是对于每个 n，我得到相同的时间，因此，计算复杂性：O(N) 这是不正确的，因为它应该是 O(N^3)。任何帮助将非常感激。

在 MATLAB 中使用行列式函数时，我观察到一个简单但特殊的行为，我想得到一些解释，因为我在函数帮助文档中没有找到任何关于它的内容。

我正在Q使用以下代码生成一个随机酉矩阵：

之后，我用det函数评估 Q 的行列式：

我希望结果是-1.000or 1.000。但是，格式似乎不是一成不变的。所以当我做这样的事情时：

detResults向量包含和1.000有时1。它只是打印格式问题还是由其他原因引起的？

我在让函数工作时遇到问题，该函数应该找到上三角矩阵的行列式。我的代码似乎返回了明显不正确的值，通常为零，我很确定这是由于我错误地定义了函数造成的。我怀疑这是我的一个基本错误，但在盯着它看了一段时间后，我还没有弄清楚。这是功能和打印代码：

和打印功能：

我将包含我迄今为止编写的完整代码以提供更多详细信息。目前它的基本应用是定义和 n 由 n 矩阵 (2

numpy.linalg.det正在发出一些我不明白的警告。这种情况并非每次都发生，但在下面的示例代码中会引发这种情况。

警告：“运行时警告：在 det r = _umath_linalg.det(a, signature=signature) 中遇到无效值

我四处寻找，但没有找到很多。这篇文章讨论了大型矩阵的溢出问题，但我怀疑这就是这里的问题。

我的代码仅适用于 3x3 矩阵如何进行一些修改以适用于 N x N 矩阵？`