linear-algebra

0热度

1回答

使用dplyr语法，我正在寻找一个“优雅”的方式来mutate新变量到给定的数据帧，其中包含不同时间点的百分比变化。作为一个例子，考虑以下data.frame包含变量（var.a，var.b，var.c）在不同时间点（A，B）服用。 set.seed(123) df <- data.frame(replicate(6,sample(1:100,50,rep=TRUE))) names(df)

0热度

1回答

将点云投影到2D

半径为r的球体内有一个点云，这些点的共同系统位于球体的中心。这个想法是从这个云层的许多角度来“拍照”这个云。 “相机”位置取决于图像中显示的角度θ（方位角）和角度（高程）。我至少需要10000张图片或观点。我该如何处理？我做：继this link我投射点每架飞机，因为我还需要visulize它们3D。就像这样：所以我谁属于“照片”，但仍与原coordenate系统的平面投影点的coorde

2热度

1回答

当Mathematica轻松成功时，SymPy的零空间陷入困境

我一直在尝试使用命令A.nullspace在SymPy中找到符号矩阵A的零空间。现在，我输入的矩阵的计算没有完成（或比我等待的时间更长）。奇怪的是，Mathematica在几分之一秒内完成了这种计算。我给我使用数学作为输入代码的具体例子（应该运行） L = ({{0, -I a Sqrt[b] Conjugate[g] Conjugate[w], 0, I a Sqrt[b] g w, 2 g

2热度

1回答

sklearn截断svd不适用于复杂的矩阵

我试图使用sklearn的TruncatedSVD（尽管我实际上调用内部函数randomized_svd来获得实际的u，s，v矩阵）的随机版本。虽然它工作正常进行实矩阵，对于复杂的矩阵我不能回来，即使奇异值是完全正确的原始矩阵： >>> import numpy as np >>> from sklearn.utils.extmath import randomized_svd >>> N =

-1热度

1回答

计算给定3D旋转的向量向量

给定欧拉角的向量3，如何在数学上找到对象正面对的方向。换句话说，Unity如何计算'transform.forward'向量？

1热度

1回答

征：反转的方形矩阵（最低标准溶液）

我有很简单的问题，但不幸的是我无法找到本征文档中回答我有一个“胖”矩阵一个（行数少于数cols），我想找到矩阵的最小规范伪逆。理想情况下，我想通过最低标准找到它QR分解如this幻灯片中所述。根据幻灯片我可以用简单直接的方法，通过使用这个公式 A.transpose() * (A * A.transpose()).inverse() 要做到这一点，但我希望有在征更优雅的解决方案 PS对不

2热度

1回答

SVD-like矩阵变换（等级标准形式）

我正在MATLAB中实现一篇研究论文，并遇到了一个矩阵变换，我不知道如何在MATLAB中完成。这， P*L*Q = [I O] 其中P，Q是变换矩阵，L是给定矩阵，I，O是身份和零分别矩阵。任何人都可以帮助我通过一些函数或算法在MATLAB中做到这一点，所以我可以通过我的代码实现这个吗？

0热度

1回答

Python linprog最小化 - 单纯形法

我正在使用scipy.optimize.linprog库来计算使用simplex方法的最小化。我在我的教科书中正在处理这个问题，我希望有人能够指出我正确的方向，因为我没有得到我期望的输出。问题是： Minimize w = 10*y1 + 15*y2 + 25*y3 Subject to: y1 + y2 + y3 >= 1000 y1 - 2*y2 >= 0

0热度

1回答

确定旋转矩阵从轴和角度

我一直在使用4×4矩阵的3D旋转的角度。我遇到了很多很棒的信息，但我对一个主题缺乏了解。我还没有理解如何确定轴的角度。如果我们看一下here，我们会发现从轴和角度描述旋转矩阵维基页面。我知道轴是这两个向量的叉积。例如： Vector1: (1,0,0) Vector2: (0,0,1) axis = Cross(Vector1, Vector2) 但是，我不知道如何找到角度。如果有人有

0热度

1回答

求解具有特征的小齐次线性系统的最快方法

我需要求解许多形如Ax = 0的小（n = 4）齐次线性系统，其中A是一个奇异矩阵。我目前使用下面的代码： void solve(const matrix_t& A, vector_t& x){ auto svd = A.jacobiSvd(Eigen::ComputeFullU | Eigen::ComputeFullV); auto V = svd.matrixV();