问题标签 [linear-algebra]

在 wpf 中做一些 3D 的东西 - 想要使用更简单的测试来查看一切是否正常（在移动到曲线之前）。

基本问题给出了两个点 x1,y1,z1 和 x2,y2,z2 我计算了点之间的距离。但是如何找到位于该线上某个距离的另一个点 (x3,y3,z3) 的坐标呢？

即，如果我的线在 -50,0,0 和 50,0,0 之间有 100 长，那么沿线 100 * 0.1 处的点的坐标是什么？

我认为这是一个简单的公式，但我还没有找到它......

如何在 lisp 中执行向量计算，例如向量的大小、向量的范数、距离（两点之间）、点积、叉积等。

谢谢。

问候，

我不确定这是否是一个愚蠢的问题。

假设我有 3 个 numpy 数组 A1、A2、A3 和 3 个浮点数 c1、c2、c3

我想评估 B = A1*c1+ A2*c2+ A3*c3

例如，numpy 会计算这个，

还是比这更聪明？在 c++ 中，我有一种巧妙的方法来抽象这种操作。

我定义了一系列通用的“LC”模板函数，用于线性组合的 LC，例如：

然后将其专门用于各种类型，

例如，对于一个数组，代码看起来像

从而避免必须创建新的中间数组。

这可能看起来很乱，但效果很好。

我可以在 python 中做类似的事情，但我不确定它是否需要。

提前感谢您的任何见解。-缺口

是的，以前已经做过一百万次了，但该死的，我还想再做一次。我正在为 C++ 编写一个简单的矩阵库，目的是做对。我遇到了一些在数学中相当明显的东西，但对于强类型系统来说并不那么明显——事实上，1x1 矩阵只是一个数字。为了避免这种情况，我开始沿着矩阵的毛茸茸的路径作为向量的组合，但也偶然发现了一个事实，即两个向量相乘可能是一个数字或一个对，这取决于两者的方向。

我的问题是，在 C++ 或 Java 等强类型语言中处理这种情况的正确方法是什么？

我正在尝试求解一组 Ax = 0 形式的方程。A 是已知的 6x6 矩阵，我使用 SVD 编写了以下代码来获得在一定程度上起作用的向量 x。答案大致正确但不足以对我有用，如何提高计算的精度？将 eps 降低到 1.e-4 以下会导致函数失败。

这个很棒的 SO 答案指向了一个很好的稀疏求解器Ax=b，但是我有这样的限制，即x每个元素x都是.>=0<=N

此外，A它很大（大约 2e6x2e6），但<=4每行元素非常稀疏。

有什么想法/建议吗？我正在寻找类似 MATLABlsqlin但具有巨大稀疏矩阵的东西。

我本质上是在尝试解决稀疏矩阵上的大规模有界变量最小二乘问题：

编辑： 在CVX中：

R 有一个qr()函数，它使用 LINPACK 或 LAPACK 执行 QR 分解（根据我的经验，后者快 5%）。返回的主要对象是一个矩阵“qr”，它包含在上三角矩阵 R（即R=qr[upper.tri(qr)]）中。到目前为止，一切都很好。qr 的下三角形部分包含“紧凑形式”的 Q。可以使用 qr 从 qr 分解中提取qr.Q()Q。我想找到 的倒数qr.Q()。换句话说，我确实有 Q 和 R，并且想把它们放在一个“qr”对象中。R 是微不足道的，但 Q 不是。目标是应用于它，这比在大型系统上qr.solve()要快得多。solve()

我一直在开发一个涉及等距空间的 Flash 游戏。我已经使用as3isolib了一段时间了，我对它的易用性印象不深。无论我是以错误的方式接近它还是它使用起来不是很好，这是另一个帖子的问题。

无论如何，我一直在想一种不同的方法来解决等距位置的问题，我想我有一个可行的想法。本质上，要渲染到等空间的每个对象都保持一个 3 坐标位置。然后，这些项目将使用相机注册，该相机根据此 Wikipedia 文章中的数学计算将该 3 坐标位置投影到屏幕上的 2 坐标点。然后，MovieClip 在该点被添加到舞台（或者可能是相机的 MovieClip），并在该点的 y 值的子索引处。这样，我认为离相机更近的对象将在更远的对象“上方”，并将在它们上方渲染。

所以我的问题有两个：

1. 你认为这个想法会像我想象的那样奏效吗？

2. 我应该查看任何现有的 3D 矩阵/矢量包吗？我知道 Flex 3 中有一个 Matrix3 类，但我们没有在这个游戏中使用 Flex。

谢谢！

我将 numpy 用于数值线性代数。例如，我怀疑如果我对执行某些计算的方式进行一些小修改，从而提高内存效率，我可以获得更好的性能。

我想知道 python 中是否有任何形式的工具来检测缓存和 TLB 未命中。有一个非常好的 api，PAPI，我在最近的课程中了解到，但它没有 Python 接口：

http://icl.cs.utk.edu/papi/overview/index.html

另外，一般有没有一种很好的方法来分析 numpy 或其他 python 数字代码？timeit 模块很难集成到代码中。mpi4py 有一种很好的方式来使用 MPE 库进行分析。演示代码片段 (demo/mpe-logging/cpilog.py)：

将创建一个可以以图形方式显示的日志文件。但这有点特定于 MPI。

谢谢。

R中是否有产生reduced row echelon form矩阵的函数？这个参考说没有。你同意？