linear-algebra相关内容

这些函数是列优先还是行优先?

我正在使用矩阵比较用于 3D 图形的两个不同的线性数学库.以下是来自两个库的两个相似的 Translate 函数: 静态Matrix4翻译(T x, T y, T z){矩阵4米;m.x.x = 1;m.x.y = 0;m.x.z = 0;m.x.w = 0;m.y.x = 0;yy = 1;m.y.z = 0;m.y.w = 0;m.z.x = 0;m.z.y = 0;m.z.z = 1;m. ..
发布时间:2021-11-25 01:46:35 C#

返回 CATransform3D 以将四边形映射到四边形

我正在尝试导出一个 CATransform3D,它将一个带有 4 个角点的四边形映射到另一个带有 4 个新角点的四边形.我花了一点时间研究这个,似乎这些步骤涉及将原始 Quad 转换为 Square,然后将该 Square 转换为新的 Quad.我的方法看起来像这样(从这里借用的代码): - (CATransform3D)quadFromSquare_x0:(float)x0 y0:(float ..
发布时间:2021-11-25 01:06:16 移动开发

如何计算相机的 z 距离以在 3D 空间中以 100% 的原始比例查看图像

如何计算相机与 3D 空间中的对象(在本例中为图像)的距离,以使图像处于其原始像素宽度. 考虑到相机的纵横比、fov 和图像的原始宽度/高度(以像素为单位),我是否认为这是可能的? (如果相关,我在这个特定实例中使用 THREE.js). 感谢任何可以帮助或引导我走向正确方向的人! 解决方案 感谢大家的投入! 在进行了一些挖掘,然后弄清楚这一切如何适合我试图用 TH ..
发布时间:2021-11-25 00:52:50 C#

在 3D 空间中旋转矢量

我正在 opengl es 中制作一个 android 项目,该项目使用加速度计来计算特定轴的变化,我的目标是旋转我的类似航天器的物体的运动矢量.问题是我无法理解旋转矩阵背后的数学原理.默认的运动向量是 0,1,0 ,表示 +y,所以物体一开始是向上看的.我正在尝试旋转它的运动向量,这样我就可以将对象移动到它指向的地方.我可以收集手机的旋转变化.x轴:旋转[0],y轴:旋转[1],z轴:旋转[2] ..
发布时间:2021-11-25 00:34:11 C#

如何计算垂直线段的端点?

我知道线段的端点以及我想创建的垂直端盖的距离/大小,但我需要计算垂直线的端点.我一直在使用 45-45-90 三角形和点积将头撞在墙上,但我似乎无法将它们组合在一起. 我知道蓝色的点和红色的点的距离,我需要找到红色的点. 在标记为重复之前,我尝试了 这个问题 但它导致端盖总是垂直的. http://rauros.net/files/caps.png http://rauros.n ..
发布时间:2021-11-24 22:00:16 C#

二维交叉积定义

在确定线段是否在多边形内中,我注意到接受的答案有一个不寻常的二维交叉定义: (u1, u2) x (v1, v2) := (u1 - v2)*(u2 - v1) 我从来没有遇到过这样的二维叉积的定义.任何人都可以启发我了解这种定义的起源吗? 解决方案 我建议你看看 外代数.它概括了叉积和行列式的概念.描述平面区域的“动机示例"部分完全回答了您的问题. 它适用于任何维度.3D 是一 ..
发布时间:2021-11-24 21:41:58 C#

numpy 中的快速外张量积

我有两个 numpy 数组: x 形状 ((d1,...,d_m))y 形状 ((e_1,...e_n)) 我想形成外张量积,即numpy数组 z 的形状 ((d1,...,d_m,e_1,...,e_n)) 这样 z[i_1,...,i_n,i_{n+1}...,i_{m+n}] == x[i_1,...i_m]*y[i_{m+1},...,i_{m+n}] 我必须多次执行上述外部乘 ..
发布时间:2021-11-18 03:44:31 Python

线性索引上三角矩阵

如果我有矩阵的上三角部分,在对角线上偏移,存储为线性数组,如何从线性数组中提取矩阵元素的 (i,j) 索引数组的索引? 例如线性数组[a0, a1, a2, a3, a4, a5, a6, a7, a8, a9 就是矩阵的存储 0 a0 a1 a2 a30 0 a4 a5 a60 0 0 a7 a80 0 0 0 a90 0 0 0 0 而且我们想知道与线性矩阵中的偏移量相对应的数组中的 ..
发布时间:2021-11-18 03:16:56 C/C++开发

大内存映射数组的高效点积

我正在处理一些相当大的、密集的 numpy 浮点数组,这些数组当前驻留在 PyTables CArray 中的磁盘上.我需要能够使用这些数组执行高效的点积,例如 C = A.dot(B),其中 A 是一个巨大的 (~1E4 x 3E5 float32) 内存映射数组,B 和 C 是驻留在核心内存中的较小的 numpy 数组. 我目前正在做的是使用 np.memmap 将数据复制到内存映射的 ..
发布时间:2021-11-18 02:58:07 Python

为什么 numpy.linalg.solve() 提供比 numpy.linalg.inv() 更精确的矩阵求逆?

我不太明白为什么 numpy.linalg.solve() 给出了更精确的答案,而 numpy.linalg.inv() 有点崩溃,给出 (我相信是)估计. 举一个具体的例子,我正在求解方程 C^{-1} * d 其中 C 表示一个矩阵,d是一个向量数组.为了便于讨论,C 的维度是形状 (1000,1000) 和 d 是形状 (1,1000). numpy.linalg.solve(A ..
发布时间:2021-11-18 02:38:55 Python

我如何翻译抛物线?

我正在开发带有一些简单动画的 iPhone 应用. 我有一个要翻译的视图,但不是沿着一条线.我想用抛物线翻译它.想象一下,我正在为一辆沿着弯曲道路行驶的汽车制作动画. 我知道我可以将转换正确设置为 CGAffineTransform 的实例 问题是,我不知道如何创建转换.我知道如何缩放、翻译等,但我如何进行抛物线翻译?甚至有可能吗? 解决方案 要沿平滑曲线设置动画,您需要 ..
发布时间:2021-11-10 05:44:45 移动开发

生成与一组向量线性无关的随机向量

我正在尝试提出一种算法,该算法将允许我生成一个随机的 N 维实值向量,该向量与一组已生成的向量线性无关.我不想强迫它们正交,只是线性无关.我知道 Graham-Schmidt 存在于正交化问题中,但是否有一种更弱的形式只给你线性无关的向量? 解决方案 Step 1. Generate random vector vr. 步骤 2. 将 vr 复制到 vo 并更新如下:对于 v1、v2 ..
发布时间:2021-09-17 19:16:15 其他开发

JS 的线性变换和矩阵乘法失败

我正在学习线性代数并尝试使用基本的线性变换(旋转、缩放、平移)制作一个小程序. 以下是完整的示例: https://codesandbox.io/embed/determined-diffie-t2iy5?fontsize=14&hidenavigation=1&theme=dark 我编写了用于为每个变换生成每个矩阵的函数以及计算它们的函数(将矩阵与一个点相乘,将矩阵相乘). ..

sympy 符号矩阵平方根

是否可以使用 sympy sqrtm 函数来计算符号矩阵的平方根?例如,如果我有矩阵 [[1, a], [a, 1]] 且 a > 0,我可以使用 sympy 来计算平方根吗? 解决方案 from sympy import *a=symbols('a')m=矩阵([[1, a], [a, 1]])米**(1/2) 给予 矩阵([[ (-a + 1)**0.5/2 + (a + 1)** ..
发布时间:2021-09-03 20:12:55 Python

如何在python中存储大于1 * 10 ^ 310的浮点数?

我正在开发一个输出大矩阵的条件数的程序,因此我使用了幂方法来获得最大特征值,但是这些值是大于1 * 10 ^ 310的大数(浮点数),并且最后值变成“无穷大",我尝试了十进制模块,但它是相同的.如何存储那些较大的浮点值?还是另一种使用较短值的方法?(我不允许使用任何将过程明确地帮助为Numpy的模块) 解决方案 您要使用十进制模块: 从十进制导入十进制x =小数('1.345e1310' ..
发布时间:2021-05-29 20:55:01 Python

如何对2D矩阵的每一行对外部乘积矩阵的求值向量化?

我正在尝试加快评估外部产品矩阵的过程.我有一个名为 a 的4 * n矩阵.我想通过公式 计算 a 的每一行的外部乘积矩阵: K = a * a'; 如果我使用 for 循环对该过程进行编码,则如下所示: K = zeros(4,4,size(a,2));对于i = 1:size(a,2)K(:,:,i)= a(:,i)* a(:,i)';结尾 我发现了另一种使用 cellfun ..
发布时间:2021-05-29 20:54:54 其他开发