-
19 votesanswersviews
为什么多态性在haskell(GHC)中如此昂贵?
我问这个问题,请参考this问题 . Don Stewart接受的答案:第一行是"Your code is highly polymorphic change all float vars to Double ..",它提供了4倍的性能提升 . 我有兴趣在Haskell中进行矩阵计算,我应该养成编写高度单态代码的习惯吗? 但是有些语言很好地利用了ad-hoc多态来生成快速代码... -
6 votesanswersviews
如何在 SceneKit 中移动旋转的 SCNNode?
下图显示了一个旋转的框,应该在 X 和 Z 轴上水平移动。 Y 应该不受影响以简化方案。盒子也可以是相机的 SCNNode,所以我猜这个投影在这一点上没有意义。 所以我们要说我们想要沿着红色箭头的方向移动盒子。如何使用 SceneKit 实现这一目标? 红色箭头表示方框的-Z 方向。它还向我们展示了它与摄像机的投影或与网格显示为深灰色线条的全局轴不平行。 我的最后一种方法是平移矩阵和旋转矩阵的乘积... -
4 votesanswersviews
如何使用 4x4 矩阵变换或欧拉角度旋转 3d 矢量?
我有一个 3d 矢量我正在申请作为一个物理力量: let force = SCNVector3(x: 0, y: 0, z: -5) node.physicsBody?.applyForce(force, asImpulse: true) 我需要根据移动设备的位置旋转力,这可以作为 4x4 矩阵变换或欧拉角。 var transform :matrix_float4x4 - 摄像机在世界坐... -
5 votesanswersviews
sceneView.session.setWorldOrigin 转换后,ARKit 节点消失
我有一些代码,其中包括用于获取标题的委托方法和一个转换。我将标题转换为弧度,并使用角度围绕 y-axis 旋转: ┌ ┐ Y = | cos(ry) 0 sin(ry) 0 | | 0 1 0 0 | | -sin(ry) 0 cos(ry) 0 | ... -
1 votesanswersviews
在 ARKit 中将屏幕点取消投影到水平面
我正在尝试使用 ARKit 将屏幕上的某些点投影到世界空间中 ARPlaneAnchor 上的坐标上。本质上,我想根据相机指向的位置在地平面上绘制标记。 ARKit 提供了一些方法来执行此操作,但是它们是附加到 ARCamera 类的,但是不幸的是,在我的场景中,我只能访问摄像机的投影,视图和变换矩阵,而没有引用 ARFrame 或 ARCamera 对象本身。据我所知,我已经写了与 gluUnp... -
0 votesanswersviews
sympy矩阵元素?
我对宏做错了什么? from sympy import * from decimal import * def MyMatrixRound(A): m = A.shape[0] n = A.shape[1] for i in range(m): for j in range(n): A[i, j]=round(A[i, j]+0.2,... -
-7 votesanswersviews
如何从C中的矩阵中删除列和行? [等候接听]
存在二维矩阵a [i] [j],条件是元素的数量直接在代码中给出 . 例如:i = 1,2 ... n和j = 1,2 ... m . 该矩阵应由程序员自己填充,然后通过删除列和行进行更改 . 我需要删除此矩阵中的任何行和任何列,然后显示修改后的矩阵 . 如果您知道任务的解决方案,请参阅代码形式的示例 . -
0 votesanswersviews
如何在不制作矩阵的情况下从数据中删除子集或超集?
假设你有这个data.frame: d <- data.frame(ID=c(1,1,1,2,2,3,3,4,4,6,6,6,8,8),FP=-c(1,2,3,1,2,3,2,1,4,1,4,3,1,4)*100) 如果你做了应急矩阵: table(d$FP,d$ID) 1 2 3 4 6 8 -400 0 0 0 1 1 1 -300 1 0 1 0 1 0 ... -
0 votesanswersviews
cv2.estimateRigid转换最小点数?
cv2.estimateRigidTransform 所需的最小点数是多少? 据我所知 fullAffine=False 它有4个自由度,所以2点就足够了 . 然而: 使用2 numpy数组作为输入: src_pts_subset.shape (2, 2) tgt_pts_subset.shape (2, 2) type(src_pts_subset) <class 'numpy.ndarr... -
0 votesanswersviews
如何从Wishart Distrubtion生成随机协方差矩阵
我需要为项目生成一个n×n,正定的协方差矩阵 . 建议从Wishart分布中提取 . 如何在R中生成随机协方差矩阵,理想情况下也使用Wishart分布 . 我已经尝试过rwishart()获取值,但需要更多帮助 . 谢谢 -
31 votesanswersviews
线性索引,逻辑索引以及所有这些
我们习惯于在Matlab中使用不同形式的索引: standard(沿每个维度使用整数), 逻辑(使用逻辑值), linear(使用单个索引遍历具有多个维度的数组) . 乍一看,似乎这些形式是独占的:索引是标准的,逻辑的或线性的 . 然而,有时似乎在这些形式中的几种形式之间存在混合 . 例如, >> A = magic(3) A = 8 1 6 ... -
1 votesanswersviews
将第一行数据集分开并重新分配给新的matrix - python
我正在将文本文件读入python,其中包含x,y和z数组(描述坐标位置) . 我想从x,y,z数组中排除第一行,即源的坐标,并使用这些值创建一个新的数组xs,ys,zs . 这样我就可以单独绘制源位置,如代码的绘图分散部分所示 . 码: import numpy as np import matplotlib.pyplot as plt import matplotlib as mpl from ... -
0 votesanswersviews
翻转扁平矩阵
我被问到以下问题: 给定表示2D矩阵的1D阵列,以0(n)或更好的方式垂直翻转值 . 给出矩阵的高度和宽度 . 我这样做的方式如下: // Given width and height int w = 8; int h = 2; // Initialize the arrays float[] random = new float[w*h]; float[] output = new float... -
-1 votesanswersviews
用于解决“列表拼图”的递归矩阵函数 - 填写列表以与预定列表交换0
需要制作自动化的“拼图”解算器 . 例: matrix = [[0,0,0,0], [0,0,0,0], [0,0,0,0], [0,1,0,0]] piece1 = [[2,2], [2,0], [2,0], [2,0]] piece2 = [[3,3], ... -
1 votesanswersviews
如何在一个区域的扩展区域内找到一个点?
好的,所以我有一个表示为矩阵的 Map ,矩阵中包含值1的单元格表示我的区域,包含值2的单元格是从我的区域开始并且也在那里结束的路径 . 当你到达路径的尽头时,也就是当你回到你的区域时,路径内的区域和已经存在的区域应该与区域“合并”,以及路径本身 . 在下图中是具有路径的这种数学的示例 . 实现这一目标的最快方法是什么? 我考虑使用填充算法,但问题是我不知道如何在非常快速的方式中在包含值1和2(... -
2 votesanswersviews
双循环矩阵代替使用循环R.
我是使用apply / purrr函数的新手 . 我不知道何时以及如何使用此功能 . 我有矩阵的双循环,我想避免它们 . 有一种方法可以做到吗?我必须做这个操作(见下面的代码) NCols=4 NRows=4 set.seed(1234) myMat<-matrix(runif(NCols*NRows), ncol=NCols) myMat norm=matrix(0,NRows,NCo... -
0 votesanswersviews
如何解决涉及行列式的Maxima中的矩阵方程?
我需要在Maxima中从这个等式中得到alpha值,这涉及行列式 . 顺便说一句,M和K是8x8矩阵,它们是已知的,只有未知的是alpha . 有人可以提出解决方法吗? 我尝试了这段代码,但我需要添加一些命令来查找alpha矢量 . eqchar:ratsimp(determinant(ev(-lambda^2*(M_cond)+K_cond)))$ eq:eqchar=0; -
1 votesanswersviews
将多个.mat文件加载到矩阵中并对其执行操作
我有429个相同大小的数字矩阵(107行乘36列),存储在顺序命名为 .mat 的文件中(例如: subj_00000.mat ... subj_00428.mat ) . 这是我需要做的: 导入MATLAB工作区 . 导入后,平均所有矩阵以生成另一个矩阵,其尺寸也为107x36 . 最后,将平均矩阵的每列与每个原始429矩阵的每列线性相关,以生成429行和36列的新矩阵 . 到目... -
0 votesanswersviews
(IndexError)返回n x n矩阵的替换
我尝试使用Python实现Back-Substitution,从 Lx = b 解决 x 但是我得到了Error . 而且我不知道如何让迭代知道什么是错的 . 我误解了一个概念吗? import numpy as np def backSub(L: np.ndarray, b: np.ndarray)-> np.ndarray: length = len(L) x = np.... -
0 votesanswersviews
为什么在递归期间变量发生了变化(尝试计算矩阵行列式)
我正在尝试使用简单递归来计算矩阵(行主序)的行列式,但由于某种原因,变量 matrix 在最后一个循环期间发生了变化,我不会写任何改变变量值的代码所有 . import math def print_m(matrix): # Print the matrix for i in range(len(matrix)): print(matrix[i]) def determ... -
3 votesanswersviews
使用包含零的numpy创建nxn矩阵的最快方法[重复]
这个问题在这里已有答案: Create Multidimensional Zeros Python 4个答案 简短的问题 - 这是用python和numpy创建一个16x16(或nxn)矩阵的最快方法吗? a = np.matrix(np.zeros((16, 16), dtype = np.int)) -
4 votesanswersviews
如何使用矩阵在android中心点旋转位图
我使用矩阵来旋转我的图像,它正在工作但我对新图像位图有一些问题 . 即使我设置了枢轴,新的位图也有一些空间 . 我该怎么做才能旋转图像,但我不想要空间 . 我使用下面的代码来旋转位图: Bitmap bitmapOrg = BitmapFactory.decodeResource(getResources(),R.drawable.sticker01); int width = bitmapOrg... -
8 votesanswersviews
在SymPy中反转矩阵?
我想知道如何在Python中使用SymPy创建矩阵并计算其逆? 例如,对于此符号矩阵: -
4 votesanswersviews
在同情中得到矩阵乘法的元素方程
我有2个矩阵,其中第一个是稀疏的整数系数 . import sympy A = sympy.eye(2) A.row_op(1, lambda v, j: v + 2*A[0, j]) 第二个是象征性的,我在它们之间执行操作: M = MatrixSymbol('M', 2, 1) X = A * M + A.col(1) 现在,我想要的是获得元素方程式: X_{0,0} = A_{0,0} ... -
2 votesanswersviews
将Sympy矩阵替换为多项式
我有一个Sympy矩阵A和一个多项式表达式P,我想计算P(A) . 这是一个例子: x = Symbol('x') P = x**2 - 3*x + 5 A = Matrix([ [1,3], [-1,2] ]) P.subs(x, A) 我希望Sympy计算 A**2 - 3*A + 5*eye(2) (在这个例子中,结果是零矩阵) . 但是这失败并显示错误消息: AttributeError... -
3 votesanswersviews
如何用SymPy中的向量替换矩阵的对角元素?
我有一个矢量 X ,我这样创建: from sympy import * x1 = Symbol('x1') x2 = Symbol('x2') x3 = Symbol('x3') X = Matrix([x1, x2, x3]) 然后我还有一个矩阵 myMat ,它只包含一个: myMat = ones(3, 3) Matrix([ [1, 1, 1], [1, 1, 1], [1, 1... -
1 votesanswersviews
Matrix在SymPy中定义了未定义的函数
我正在寻找指定依赖于变量的矩阵量的可能性 . 对于如下工作的标量,使用未定义的函数: from sympy import * x = Function('f')(t) diff(x,t) 对于像Matrix这样的符号 x = MatrixSymbol('x',3,3) 我找不到相应的东西 . 有 i,j = Symbols('i j') x = FunctionMatrix(6,1,Lambd... -
7 votesanswersviews
如何在SymPy中扩展矩阵表达式?
在SymPy中,我试图执行矩阵乘法并在之后展开它 . 但是,SymPy似乎不支持矩阵表达式的扩展 . 例如,这是矩阵的4阶Runge-Kutta(RK4): from sympy import init_session init_session() from sympy import * A = MatrixSymbol('A', 3, 3) x = MatrixSymbol('x', 3, 1... -
3 votesanswersviews
如何使用Sympy导出Matrix元素
给定矩阵和向量的乘积 A.v 对于形状(m,n)和v的dim n,其中m和n是符号,我需要相对于矩阵元素计算导数 . 我还没有找到使用正确向量的方法,所以我从2 MatrixSymbol 开始: n, m = symbols('n m') j = tensor.Idx('j') i = tensor.Idx('i') l = tensor.Idx('l') h = tensor.Idx('h')... -
-3 votesanswersviews
如何在张量流中将k 2x2矩阵乘以k 2x2矩阵?
我在张量流中有2个3D张量,其中两个张量的形状为Kx2x2 . 张量代表一组2x2矩阵 . 有没有办法将第一张量中的2x2矩阵与第二张量中的相应矩阵相乘,以便最终得到Kx2x2张量?