首页 文章
  • 9 votes
     answers
     views

    矩阵变换:将SVG路径坐标转换为Leaflet坐标系

    简短版本:如何将SVG路径添加到Leaflet Map ,以便在 Map 坐标更改时路径将更新(例如,在缩放更改或幻灯片上)? 长版:你好,我有一个包含建筑轮廓的地形image . 在对图像进行地理校正后,我使用Photoshop将栅格数据转换为SVG . 我知道描述SVG边界的边界框的地理坐标,并知道SVG路径元素的内部坐标 . 我想知道现在将SVG的路径元素中描述的建筑物添加到Leaflet ...
  • 1 votes
     answers
     views

    绘制矢量到一条线

    我正在进行模拟,我有一架飞机,我需要能够飞到一条直线的起点 . 到达该点时,需要与线的角度对齐 . 起点可以是线上的任一点 . 它类似于模拟飞机降落在跑道上,但我不需要考虑高度 . example 我有以下信息: 飞机矢量 纬度/经度 Headers 速度 目的地线(两分) 点1纬度/经度 点2纬度/经度 飞机位置每0.5秒更新一次,限制为每秒3度转弯率 . 我目前...
  • 0 votes
     answers
     views

    从线性系统中消除冗余方程

    我有一个或多个方程和3个未知数 . 我想将它减少到可能的最小方程数 . 可以有0,1或无限多个解决方案 . 在0解决方案的情况下,我不需要缩减系统,只需要知道它正在发生 . 高斯消除可以,但每个人都说它在数值上不稳定 . 也许这对于这么小的系统来说无关紧要,只要你使用旋转?我也不需要行梯形式,所以它有点矫枉过正 . 他们说SVD更稳定,但我不知道如何从它产生的U,Sigma和V矩阵中获得减少的方程...
  • 4 votes
     answers
     views

    用于java的线性代数库[关闭]

    我有一组具有无限解决方案的线性方程(解决方案集)但我似乎无法找到任何可以生成解决方案集的Java库(他们都抱怨矩阵是“单数”,这会产生无限的解决方案?) . 我可能无法理解某些内容,但我已经将一些库调查为JAMA和la4j但是我无法为我生成解决方案集 . BR Tomas A. 编辑:我尝试使用在线“解算器”解决方程组,它为我生成解决方案集,所以我知道它应该是可能的 .
  • 0 votes
     answers
     views

    线性代数方程Ax = By的最小二乘解在特征C中

    我有一组矩阵形式的线性代数方程, Ax=By . 其中 A 是 36x20 的矩阵, x 是大小为 20 的向量, B 是 36x13 , y 是 13x1 . Rank(A)=20 . 因为系统是超定的(方程的数量多于变量),所以最小二乘解是可能的,i,e; x = (A^TA)^-1A^TBy . 我想要解决方案,以便最小化剩余误差e = Ax-By . 使用 Eigen/Den...
  • 4 votes
     answers
     views

    matlab中相同算法的结果不同

    我正在进行线性代数的分配,以比较QR分解算法Gram-Schmidt和Householder的性能和稳定性 . 在计算下表时我怀疑: 通过将Gram-Schmidt和住户应用于Hilbert矩阵A,矩阵Q和R是QR分解的结果矩阵,I是维数N的单位矩阵;和|| * ||是弗罗贝尼乌斯的常态 . 当我在不同的计算机上进行计算时,在某些情况下我会得到不同的结果,可能是由于这个原因?上表对应于在32位计...
  • 2 votes
     answers
     views

    当Mathematica轻松成功时,SymPy的零空间挣扎

    我一直在尝试使用命令 A.nullspace 在SymPy中找到符号矩阵 A 的零空间 . 现在,对于我输入的矩阵,计算没有完成(或者比我等待的时间更长) . 奇怪的是,Mathematica在几分之一秒内完成了这项计算 . 我将给出我作为Mathematica输入代码使用的具体示例(应该运行) L = ({{0, -I a Sqrt[b] Conjugate[g] Conjugate[w], 0...
  • 1 votes
     answers
     views

    具有多个解的线性方程组(R,rSympy)

    我需要求解一个线性方程组,它具有比方程更多的变量,因此不一定只有一个解 . 我可以将问题表示为具有已知行和列总和的未知变量矩阵,如下所示: X1 X2 X3 Y1 V11 V12 0 Y2 V21 V22 V23 Y3 V31 V32 V33 Xs和Ys分别是已知的列和行和,Vs是未知的变量值,它们...
  • -1 votes
     answers
     views

    你好 . 我想要一些帮助 . 你能帮我解决这个问题吗? [等候接听]

    叠加原理指出,如果两个或多个波在同一介质(空气,水等)中行进,则通过将各个波的位移相加来找到所得到的波 . 例如,对于波f和g,得到的波是y = f(x)g(x) . 一个 . 在0≤x≤2π的区间内,图f(x)= cos2x和g(x)= cos(2x-1) . f和g的图表有何不同?湾图f(x)= f(x)g(x)在与f和g相同的屏幕中 . 描述h图的行为 . C . 比较f,g和h的周期 ...
  • 1 votes
     answers
     views

    将自定义函数应用于sklearn中的稀疏矩阵

    我正在研究一个情绪分析项目 . 我在sklearn上用了一袋文字 . 我想将行正常化,如下所示: a = [1, 0, 1, 0, 0, 0, 0, 1, 0, 1] #a ->  4 non zero values 规范化的向量应该是: [0.25, 0.0, 0.25, 0.0, 0.0, 0.0, 0.0, 0.25, 0.0, 0.25] 对于b - > 8个非零值: b...
  • 1 votes
     answers
     views

    三维空间中三个三维点的线性插值

    我有三个3D点,如 p1(x1,y1,z1) , p2(x2,y2,z2) , p3(x3,y3,z3) . 我还有另外一点,但我只知道 x , y 这一点的值,如 p4(x4,y4,Z) ,其中 Z 是我想要计算的值 . 我确信 p4(x4,y4) 点是在 p1(x1,y1) , p2(x2,y2) , p3(x3,y3) 形成的三角形内部,通过delaunay三角测量方法检查 . 如何计算点...
  • 0 votes
     answers
     views

    用于大型系统的快速线性系统求解器

    如何解决系统矩阵带状对角线的大型线性系统?只使用 \ 需要的时间比我想要的要长 . 我愿意使用迭代方法,如果这是唯一的选择,但是是否已经实现了任何已知的高效m文件?我应该寻找什么?非常感谢!
  • 5 votes
     answers
     views

    求解非方阵中的线性方程组

    我有一个线性方程组,它构成了我需要解决的 NxM 矩阵(即非方形) - 或者至少试图解决以表明系统没有解决方案 . (更可能的是,没有解决方案) 据我所知,如果我的矩阵不是正方形(过度或不足),那么就找不到 exact solution - 我是否正确地想到这一点?有没有办法将我的矩阵转换为方形矩阵以计算确定性,应用高斯消元法,克莱默规则等? 值得一提的是,我的未知数的系数可能为零,因此在某...
  • 11 votes
     answers
     views

    相当于Python中2D多项式的“polyfit”

    我想为 a 系数找到最小二乘解 z = (a0 + a1*x + a2*y + a3*x**2 + a4*x**2*y + a5*x**2*y**2 + a6*y**2 + a7*x*y**2 + a8*x*y) 给定长度为20的数组 x , y 和 z . 基本上我正在寻找相当于numpy.polyfit但是对于2D多项式 . This question类似,但解决方案是通过MAT...
  • 11 votes
     answers
     views

    3D基本矩阵的对应关系

    在MATLAB中,我使用标准化的Eight point algorithm计算了Fundamental matrix(两个图像) . 从那我需要triangulate 3D空间中的相应图像点 . 根据我的理解,要做到这一点,我需要旋转和翻译图像的相机 . 当然最简单的方法是calibrate the cameras然后拍摄图像,但这对我的应用来说太紧张了,因为它需要这个额外的步骤 . 所以这让我留...
  • 0 votes
     answers
     views

    在简化的相机模型中解决3 DOF旋转问题

    我试图在受约束的投影相机模型中获得3x3旋转矩阵和焦距参数: 将像素坐标与真实世界坐标相关联的公式: 描述内在和外在参数的3x4矩阵 所以在我的约束相机模型中,u0和v0设置为1,Tx,Tv,Tz设置为0,因此DOF减小了4 DOF(旋转矩阵3 DOF为焦距(au,av) )) . 我的问题是,虽然我知道旋转矩阵是正交的,并且每个列/行彼此相关,因此3x3矩阵虽然有9个元素但只有3个DOF,但我想...
  • 6 votes
     answers
     views

    python中稀疏矩阵的伪逆

    我正在使用来自神经成像的数据,并且由于数据量很大,我想为我的代码使用稀疏矩阵(scipy.sparse.lil_matrix或csr_matrix) . 特别是,我需要计算矩阵的伪逆以解决最小二乘问题 . 我找到了sparse.lsqr方法,但效率不高 . 有没有一种方法来计算Moore-Penrose的伪逆(对应于正常矩阵的pinv) . 我的矩阵A的大小约为600'000x2000,并且在矩阵...
  • 6 votes
     answers
     views

    代表旅行推销员作为线性表达

    我在网上看到,人们可以将旅行商问题写成线性表达式,并使用CPLEX for java等软件进行计算 . 我有1000个城镇,需要找一小段距离 . 我计划将这1000个城镇划分为约100个城镇的集群,并在这些单独的集群上执行一些线性规划算法 . 我的问题是,我究竟如何将其表示为线性表达式 . 所以我有100个城镇,我相信每个人都知道TSP是如何工作的 . 我完全不知道如何编写满足TSP的线性约束,目...
  • 4 votes
     answers
     views

    在校准的立体视觉装备中,如何获得实现3D三角测量算法所需的“相机矩阵”?

    我试图实现Hartley&Zisserman的"Multiple View Geometry"(第12.2节)中的(相对简单的)线性齐次(DLT)3D三角测量方法,目的是在未来实现它们的全部,"optimal algorithm" . 现在,基于this question,我试图让它在Matlab中运行,稍后将其移植到C和OpenCV中,测试整合过程中的符合...
  • 1 votes
     answers
     views

    cv2.projectPoints如何实际工作?

    我正在尝试编写自己的cv2.projectPoints实现,仅用于练习,而我的实现并没有按预期运行 . 我假设没有失真的针孔相机,所以我所做的就是将以下等式应用于我的所有平面物点 . 让我们说在一些世界空间(飞机和摄像机都存在),飞机距离2个单位 . 因此,在世界坐标空间中,平面上的所有点都有 z=2 . 但是,在平面自己的坐标空间中,平面上的所有点都有 z=0 . 所以,矢量 [X Y Z...
  • 1 votes
     answers
     views

    张量流中补丁的点积

    我有两个相同大小的方形矩阵和一个方形补丁的尺寸 . 我想计算每对补丁之间的点积 . 基本上我想实现以下操作: def patch_dot(A, B, patch_dim): res_dim = A.shape[0] - patch_dim + 1 res = np.zeros([res_dim, res_dim, res_dim, res_dim]) for i in xr...
  • 31 votes
     answers
     views

    TensorFlow中矩阵和向量的有效元素乘法

    什么是倍增(逐元素)2D张量(矩阵)的最有效方法: x11 x12 .. x1N ... xM1 xM2 .. xMN 通过垂直向量: w1 ... wN 获得一个新的矩阵: x11*w1 x12*w2 ... x1N*wN ... xM1*w1 xM2*w2 ... xMN*wN 为了给出一些上下文,我们可以批量处理 M 数据样本,并且每个 N -element样本必须乘以存储在变量中的...
  • 22 votes
     answers
     views

    如何计算scipy稀疏矩阵行列式而不将其变为密集?

    我试图找出在python中找到稀疏对称和实矩阵行列式的最快方法 . 使用scipy sparse 模块,但真的很惊讶没有行列式功能 . 我知道我可以使用LU分解来计算行列式,但是没有看到一个简单的方法来执行它,因为 scipy.sparse.linalg.splu 的返回是一个对象并且实例化一个密集的L和U矩阵是不值得的 - 我不妨做 sp.linalg.det(A.todense()) wh...
  • 5 votes
     answers
     views

    解决稀疏的上三角系统

    我试图弄清楚如何有效地解决稀疏的三角形系统,在scipy稀疏中Au * x = b . 例如,我们可以构造一个稀疏的上三角矩阵,Au和一个右手边b,其中: import scipy.sparse as sp import scipy.sparse.linalg as sla import numpy as np n = 2000 A = sp.rand(n, n, density=0.4) + ...
  • 1 votes
     answers
     views

    SVM中的决策边界计算

    一旦得到模型的系数,我无法理解如何计算决策边界 . 这是我所指的链接:http://scikit-learn.org/stable/auto_examples/svm/plot_svm_margin.html 这是代码 # get the separating hyperplane w = clf.coef_[0] a = -w[0] / w[1] xx = np.linspace(-5, 5) ...
  • 4 votes
     answers
     views

    BLAS矩阵通过矩阵转置乘法

    我必须以 A'A 或更一般 A'DA 的形式计算一些产品,其中 A 是一般 mxn 矩阵, D 是对角线 mxm 矩阵 . 两者都是满级;即 rank(A)=min(m,n) . 我知道你可以节省大量时间就是这样的对称产品:鉴于 A'A 是对称的,你只需要计算产品矩阵的下部 - 或上部 - 对角线部分 . 这增加了要计算的 n(n+1)/2 个条目,这大约是大型矩阵的典型 n^2 的一半 . 这...
  • 2 votes
     answers
     views

    没有部分旋转的高斯消除

    晚上好, a=[[1, 2, -1], [-2, 4, 1], [2, 2, 3]] b=[-4,6, 1] n=3 def SPP(): for i in range(1,n): s=[0]*3 p=[0]*3 for j in range(1,n): s[i]+=max(s[i], abs(a[i][j...
  • 2 votes
     answers
     views

    当表示为对象的一维向量时,有效地旋转NxM矩阵(C)

    到目前为止,我已经想出顺时针旋转NxM(N不一定等于M)矩阵的唯一方法(当它表示为一维向量,高度和宽度变量分别存储)如下: struct matrix { vector<int> data; int height; int width; void rotate_90() { vector<int> newdata(height*width); ...
  • 7 votes
     answers
     views

    使用boost c ublas矩阵和向量类型的元素操作

    我想在增强矩阵和矢量类型上执行逐元素函数,例如:取每个元素的对数,取幂每个元素,应用特殊函数,如gamma和digamma等(类似于matlab对这些函数的处理应用于矩阵和向量 . ) 我想写一个帮助函数,对每个所需的函数强制执行这个就足够了,但这看起来很浪费 . 同样,boost wiki为vectorize standard functions提供了一些代码,但这看起来相当复杂 . valar...
  • 2 votes
     answers
     views

    在同义词中对矩阵进行逐列(或逐行)操作

    是否有办法在同义词中对矩阵进行逐列(或逐行)操作?例如,将矩阵的每一列除以其范数,或将矩阵的每一行乘以其范数?

热门问题