线性代数是数学的核心体，研究这些空间之间的向量空间和线性映射。主要主题包括：矩阵，向量，线性方程的同时系统，对角线形式，特征值，特征向量，投影，点，线，平面，超平面，旋转和反射。
使用 scipy.optimize.linprog 最小化 Ax - b 的 L1 范数
众所周知，最小化 x 上 Ax - b 的 1-范数问题可以写成线性规划问题。我想使用 scipy.optimize.linprog 来解决这个最小化问题，...
为什么通过torch.linalg.eig计算出的特征对没有通过特征对测试？
运行以下代码： 进口火炬 火炬.manual_seed(47) m = torch.rand([8, 8]) eigvals, eigvecs = torch.linalg.eig(m) print(m.type(torch.complex64) @ eigvecs[0] - eigvals[0] * eigvecs[0]...
我正在寻找一个Python代码来创建一个5乘5的零一矩阵，行列式等于5。（这样的矩阵确实存在，但它是一些数学定理的结果。我喜欢尝试找到.. .
在 Julia 中，我想使用函数 qrfact()，我认为它是 LinearAlgebra 包的一部分。然而，当我尝试使用该函数时，Julia 声称函数 qrfact() 不是 de...
reverse_cuthill_mckee：如何避免对角线上的零以及如何取回原始矩阵？
我正在尝试使用reverse_cuthill_mckee来减少稀疏、对称、复杂矩阵的带宽，从而提高其条件数。但当我在更简单的矩阵上进行测试时，似乎
我想求解 3 个变量 Ax=b 的线性系统，其中 A 是 4x3 矩阵，b 是向量。 A 是矩阵 A=[1 1 1 ; 1 2 1 ; 1 2 3 ; 1 4 1] 而 b=[1 1 1 2] 完整的 4x4 矩阵 A|b 有 ra...
NumPy 中的右矩阵除法，还有比 np.linalg.inv() 更好的方法吗？
进行正确矩阵除法（即 pyhton numpy 中的 xA = B）的最快方法是什么？ A 和 B 是 (NxN) 方阵，并且 A 是可逆的，因此我想计算方程 x = BA^{-1}。 我问
在 Julia 中矩阵相乘时有什么方法可以减少分配次数吗？如果可能的话，我想减少这个数字。我已附上显示此内容的屏幕截图。
在 Julia 中，我想连接几个数组（也可以将它们相乘）。在我的程序中，我将其编写如下： [Uᵣ Qₐ]*Uₖ [VᵣQᵦ]*Vₖ 然而，这个数组串联非常expe...
使用“scipy.optimize.least_squares()”用可变数量的参数拟合两条曲线
我正在使用 scipy.optimize.least_squares() 函数编写代码来适应系列中相关函数对的参数，并且似乎测试参数没有传递给 fu...
我正在尝试使用 Whittaker-Eilers 来平滑噪声数据。使用本文创建一个方法来求解我所写页面上的方程 11 和 12 中的 A： 私人 LUDecomposition buildMatrices...
我想随机生成一个任意大小的数组，其中数组的所有元素都是绝对值为一的复数。在 Julia 中有什么办法可以做到这一点吗？
我有一个大网格，每个点都有一个关联的矩阵（3x3）。我想要执行对称和斜对称分解。我知道算法是 𝐴=0.5(𝐴+𝐴T)+0.5(𝐴−𝐴T) 我正在检查...
在 Julia 中，我想编写一个函数，多次提示用户输入矩阵，然后将其输入存储到数组中。到目前为止我已经尝试过以下操作： 函数接受列表（...
我想生成一个矩阵，每一行和每一列中只有一个非零条目，其中每个非零条目的模数为一。朱莉娅有什么办法吗？
生成随机正交矩阵的一种简单方法是首先对随机矩阵进行采样，然后应用奇异值分解 def random_orthonormal_matrix(n): 随机矩阵...
假设我有一个大小为 N=150,000 的稀疏对称实不定矩阵 A。该矩阵是在具有相同结构（非零元素的位置）的算法的每次迭代时生成的，但是
我有一个块对角矩阵，其中一些对角块为零。矩阵是一个运算符的表示，我想要这个运算符的逆（各个块有...
我想从 Qt C++ OpenGL ES 2.0 中的变换矩阵获取缩放、旋转和平移，以便为 Android 和 WebAssembly 制作具有线性插值的关键帧动画。我有...
我试图概括一个类来处理任何维度和类型（float/int/double）的向量。几乎所有常见操作都可以直接推广到更高维度，除了