通过距单元格的距离对矩阵中的向量求和(R)
通过距单元格的距离对矩阵中的向量求和(R)是一种矩阵运算方法，用于计算矩阵中向量之间的距离和。在这种方法中，矩阵中的每个向量都与其他向量进行比较，并计算它们之间的距离。距离可以使用不同的度量方法，如欧氏距离、曼哈顿距离或余弦相似度等。
这种方法的优势在于可以通过计算向量之间的距离来衡量它们之间的相似性或差异性。通过对矩阵中的向量进行距离求和，可以得到一个表示整个矩阵的总体相似性或差异性的值。
应用场景：
- 数据聚类：通过对矩阵中的向量进行距离求和，可以将相似的向量聚类在一起，从而实现数据的分类和分组。
- 特征提取：通过计算向量之间的距离，可以提取出矩阵中的重要特征，用于后续的数据分析和模型训练。
- 图像处理：在图像处理中，可以使用距离求和方法来比较不同图像之间的相似性，从而实现图像检索、图像匹配等功能。
腾讯云相关产品和产品介绍链接地址： 腾讯云提供了丰富的云计算产品和服务，以下是一些相关产品和介绍链接地址：
- 云服务器（CVM）：提供弹性计算能力，支持多种操作系统和应用场景。产品介绍链接
- 云数据库（CDB）：提供高可用、可扩展的数据库服务，支持多种数据库引擎。产品介绍链接
- 人工智能平台（AI Lab）：提供丰富的人工智能算法和工具，支持图像识别、语音识别、自然语言处理等应用。产品介绍链接
- 云存储（COS）：提供安全可靠的对象存储服务，支持海量数据的存储和访问。产品介绍链接
- 区块链服务（BCS）：提供一站式区块链解决方案，支持快速搭建和管理区块链网络。产品介绍链接
请注意，以上链接仅供参考，具体产品选择应根据实际需求进行评估和决策。
相关·内容
2回答
尽可能安排学生尽量减少作弊。
、、、、
给定一个m* n的矩阵，将k个学生置于这样的方式，这样考试中的作弊就能被最小化。
My方法：考虑m*n矩阵，将第一个学生放在X上(X=在给定矩阵中的任何单元格)。对于下一个学生，做一个BFS，对于每个空单元格，计算它的远程(R)。对于R，我们可以从值t=max(m，n)开始，在BFS中遇到的下一个单元格中，我们添加t级。存储这个单元格，如果R是最小的，到目前为止。最后将下一个学生放入R值最小的单元格中。对所有学生重复这一步骤(K)。
优化:如果我们到达一个R大于我们现在所拥有的细胞时，忽略它。
我们如何优化它?
有更好的算法来做这个吗?
有我的起始单元X和最终答案有关
例如，如果它的5x5垫子
浏览 0提问于2014-04-25得票数 0
1回答
在MATLAB中使用LIBSVM，需要预先计算直方图交集核矩阵.
假设x，y是两个向量。核函数为K(x，y) =和(min(x，y))。为了提高效率，大多数情况下的最佳做法是将操作向量化。
我想要做的是计算核矩阵，比如计算两个矩阵之间的欧几里德距离，比如pdist2(A，B，‘欧几里得’)。在定义函数'intKernel‘之后，我可以通过调用pdist2(A，B，intKernel)来计算交叉内核。
我知道函数'pdist2‘可能是一个选项。但是我不知道怎么写自定义的距离函数。然而，我不知道如何在一个浓缩表达式中编码向量(1×-M)和矩阵(M-by-N)之间的交集核。
“爬虫垫
浏览 4提问于2014-01-14得票数 7
回答已采纳
问题
我需要将大小为N×M的矩形放入n×m矩阵的自由区域，其中N=n*2-1，M=m*2-1。矩阵单元格如果是true，则被认为是自由的，如果是false，则被占用。中心单元格始终是true，并且由于矩形和矩阵的大小，总是位于矩形内。
另外的要求是矩形的左上角和中间单元格之间的距离必须尽可能小。
使用n=8和m=5的示例
其中灰色单元-占据，绿色中心单元，蓝色单元-解决矩形，红线距离之间的左上角矩形和中心单元格。
尝试
蛮力解具有O(N×M×n×m)的时间复杂度，这不是最优解。如果我对矩阵进行预处理，我可以消除某些单元格中的计算，但这仍然需要很长时间。
一开始，我认为我可以处理最大矩
浏览 4提问于2017-08-14得票数 7
回答已采纳
当我对一个矩阵进行求和时： import sympy
i = sy.symbols('i')
A=Matrix([1,i,3,3*i])
summation(A[i],(i,0,3)) 结果是： IndexError: Invalid index a[i] 我知道我可以用循环来解决这个问题。 难道没有支持渐近矩阵的求和函数吗？
浏览 13提问于2019-03-09得票数 1
我有一个关于旋转(R)和平移(T)矩阵的问题，我希望有人能帮助我，因为我们在机器人学中经常使用R和T来找到机器人的位置。
我有相对于相机帧的对象帧的R和T，以及关于同一相机的第二个对象的R和T。如果两个对象都在一个公共表面/平面上
问题#1：(在同一表面/平面上)这意味着两个物体相对于相机的R是相同的！这个假设是正确的吗？
问题#2:如何使用平移矩阵计算对象之间的距离(例如沿x或y轴)？
我有一个1x3的平移矩阵和一个3x3的旋转矩阵，我从R和T导出了一个4x4的变换矩阵=>
提前感谢
浏览 0提问于2012-01-28得票数 6
回答已采纳
1回答
素食者中的R大距离矩阵
、、、
我在一台拥有128 GB内存的机器上运行R3.2.3。我有一个123028行x 168列的大矩阵。我想在R中使用分层聚类算法，因此在这样做之前，我尝试用Bray方法在纯素包中使用vegdist()函数在R中创建一个距离矩阵。我得到了一个内存分配错误：
df <- as.data.frame(matrix(rnorm(20668704), nrow = 123028))
library(vegan)
mydist <- vegdist(df)
素食主义错误(Df)：.Fortran不支持长向量(参数4)
如果我使用pryr包来找出距离矩阵需要多少内存，我就会发现需要121 G
浏览 3提问于2015-12-14得票数 4
回答已采纳
我有一个简单的函数来计算两个向量之间的距离，这样距离=点积/两个向量中元素的和。
function d = simpleDistance(a,b)
d = dot(a,b)/ (sum(a) + sum(b));
end
例如：simpleDistance([1 2], [3 4]) = (3 + 8) / (3 + 7) = 11/10 = 1.1
给定这个小矩阵r，我想计算r中每两行之间的相似性(距离= simpleDistance)。
r =
1 2
5 0
3 4
而不是两个嵌套循环，我想使用pdist函数，因为它要快得多！
n = siz
浏览 0提问于2014-03-12得票数 0
回答已采纳
1回答
请通过我对LSI (机械工程背景)的理解来理解我：
在LSI中执行SVD之后，有3个矩阵：
U，S，V转座。
U将单词与主题进行比较，S是衡量每个特征强度的一种方法。将主题与文档进行比较。
U dot S dot Vt
返回SVD之前的原始矩阵。没有做太多(没有)深度代数，似乎：
U dot S dot **Ut**
返回一项一项的矩阵，该矩阵提供了两个术语之间的比较。也就是说，一个术语与其他术语的关系，是一种比较词而不是成分的DSM (设计结构矩阵)。我可能完全错了，但我尝试了一个样本数据集，结果似乎是有意义的。这可能只是偏见(我想让它起作用，所以我看到了我想要的)。我不能发布结果，因为
浏览 6提问于2013-02-04得票数 1
3回答
更多的是一个一般性的MATLAB问题，而不是寻找编程建议--如果我有：
y = cellfun(@(x)sum(x(:)), Z, 'un', 0);
在每个单元矩阵中有总和实数的组合时，当我对每个单元矩阵的所有元素求和时，我是否总是会得到NaN= NaN，因为其中有NaN，还是会忽略它们而只对实数求和。我问这个问题的原因是因为我得到了：
y = [NaN] [NaN] [NaN]
[NaN] [NaN] [NaN]
[NaN] [NaN] [NaN]
单元格矩阵(单元格元素)的示例如下：
x{1,1} = NaN 2 3
4 5
浏览 0提问于2012-07-02得票数 3
回答已采纳
3回答
比较布尔矩阵
、、、
我有个小问题。我需要比较一个由1和0填充的二维数组(让我们称之为矩阵A -零实际上表示空白点，1是足球运动员在球场上的位置)与许多其他填充不同的矩阵(但同样，只有1和0)，结果应该是一些指示哪个矩阵与矩阵A最相似。相似度我指的是球员在球场上的分布(或位置)的相似性-因此将选择与矩阵A最相似的球员位置的矩阵。
有人能帮我解决这个算法问题吗？
我是用c++编写的，但是伪代码就足够了。问题只是一个比较算法。最好是比较函数的输出类似于similarity coefficient，我可以将其存储在一个数组中，然后使用它选择最相似的矩阵。但我就是不能想出一些算法来进行这种相似性比较。
编辑:一些关于相似性
浏览 2提问于2012-01-21得票数 3
我有个大问题。
我得到了一个很大的栅格(rows=180，columns=480，number of cells=86400)，首先我对它进行了二值化(所以只有1和0)，然后我对聚类进行了标记。(相互连接的1个单元具有相同的标记。)
现在我需要计算所有单元格之间的距离，这些距离不是0。有很多是安静的，这是我最大的问题。我这样做是为了获得我感兴趣的单元格的坐标(获取单元格的位置(即单元格编号)，这些位置不是0)：
V=getValues(label)
Vu=c(1:max(V))
pos=which(V %in% Vu)
XY=xyFromCell(label,pos)
这个很好用。因此，XY是
浏览 0提问于2017-01-19得票数 0
1回答
我有一个矩阵distance = [d11,d12,d13,d14;d15;...;dn1,dn2,dn3,dn4,dn5];和一个向量index(n,1)。索引的值介于1和5之间。
如果一个向量R(1,5)，我想根据索引得到距离的和。
举个例子：
distance = [1,2,4,1,2 ; 4,5,6,1,6 ; 7,8,9,5,8] and index = [1;1;3]
所以，我想要R(1) = 1+4 = 5, R(2) = 0, R(3) = 9, and R(4) = R(5) = 0
条件是不使用带有if条件的1:5以上的循环，以便在有数十亿个点的情况下最小化执行时间。
也
浏览 11提问于2017-03-15得票数 1
回答已采纳
1回答
我正在寻找一个Java Matrix库来执行数据分析和实现聚类算法(如K-means或DBSCAN )。
我发现了Colt和Parallel Colt(在大数据集和小数据集上表现最好)，但显然它们不支持字符串矩阵。数据集条目应该只是双矩阵。
有什么建议吗？
提前感谢您的帮助。
浏览 5提问于2012-02-26得票数 1
回答已采纳
1回答
我有3*4相机矩阵，它是P = K[R | -RC]，其中K是内部相机矩阵，R是旋转矩阵，C是相机相对于世界原点的位置，-RC是平移。我可以将P分解为K、R和C。我希望如何使用此相机信息在openGL中渲染。如何在openGL中使用K、R、C；设置视图矩阵和投影矩阵，使其表示摄像机的实际位置？
浏览 7提问于2017-11-09得票数 0
我想在计算机程序中从3D体积中随机获得2D切片，以生成噪声。我决定使用旋转矩阵，我知道它必须是旋转组SO(3)的成员。如何创建此组的统一随机成员？
浏览 0提问于2012-05-25得票数 2
回答已采纳
1回答
我有一组不同像素的RGB矩阵。(N像素=> n行，RGB => 3列)。我必须计算出这个矩阵中任意两个像素之间的最小RGB距离。我尝试了循环方法，但是由于集合太大(比方说N=24000)，看起来程序需要很长时间才能完成。还有别的办法吗？我读过关于pdist的文章，但是RGB欧几里德距离不能和它一起使用。
k=1;
for i = 1:N
for j = 1:N
if (i~=j)
dist_vect(k)=RGB_dist(U(i,1),U(j,1),U(i,2),U(j,2),U(i,3),U(j,3))
k=k+1
浏览 4提问于2014-06-24得票数 1
回答已采纳
给定一个二进制矩阵，我希望生成一个列表，以逐行显示1s的位置。例如：
> M <- matrix(c(1,1,0,0,1,1,1,0,1), 3, 3)
> M
[,1] [,2] [,3]
[1,] 1 0 1
[2,] 1 1 0
[3,] 0 1 1
> apply(M==1, 1, which, simplify = FALSE)
[[1]]
[1] 1 3
[[2]]
[1] 1 2
[[3]]
[1] 2 3
但是，在4.1.0之前的R版本中，apply()函数会自动简化，这意味着结果有时不
浏览 9提问于2022-03-18得票数 1
回答已采纳
K-质心聚类函数以数值数据矩阵作为输入.然而，我这里只有距离矩阵，我认为k-质心算法确实适用于距离矩阵。
从正式文件中用尽
使用
library(flexclust)
kcca(x, k, family=kccaFamily("kmeans"), weights=NULL, group=NULL,
control=NULL, simple=FALSE)
参数
x A numeric matrix of data, or an object that can be coerced to such a matrix (such as a numeric vector
浏览 2提问于2013-07-29得票数 0
回答已采纳
1回答
目的是评估6个已知类群之间的相似性和差异性。
最初的数据从6个已知组和2,700+变量开始，所有变量的评分都在0到100之间。
我使用R中的2700+包中的dudi.pca函数将ade4变量简化为5个主成分，以下是这些组件的特征值：
eigenvalue variance.percent cumulative.variance.percent
Dim.1 998.3274 36.635867 36.63587
Dim.2 670.1278 24.591848 61.227
浏览 0提问于2017-11-14得票数 0
我有一个问题，我试图用mathematica来解决。我有一个列表，其中包含位置测量的x和y坐标(以及在每个点上测量的量的z值)。因此，我的列表从list={{-762.369, 109.998，0.915951}，{-772.412,109.993,0.923894}，{-777.39,109.998，0.918108}，...} (x，y，z)开始。出于某些原因，我必须将所有这些x、y和z值填充到一个矩阵中。如果我对每个y坐标都有相同数量的x坐标(假设是80)，那么这将很容易，然后我可以使用Partitionlist,80，它产生一个有80列的矩阵(以及一些行，其数量由具有相同值的y坐标的数
浏览 0提问于2011-07-26得票数 3
回答已采纳
1回答
计算移位材料的范数
、、、
披露：交叉贴在虚幻的答案上，认为对GameDev.SE来说是足够普遍的
我目前正在绘制一个位移纹理，另外，一个正常的纹理为我的海洋。是否有一种聪明的方法，我可以在最后的材料中计算出正常值。也就是说，非位移正规是(0，0，1)，那么用ddx/ddy也许可以求出真正的正规？看上去像这个问题是在尝试一个类似的问题。。我并不真正理解ddx/ddy节点/函数，但它们似乎至少是相关的。
使用ddx，ddy节点，我尝试的是计算每个给定点：
dx = ddx(WorldPosition.x)
dy = ddy(WorldPosition.y)
dz/dx = ddx(height) / dx
dz/dy =
浏览 0提问于2015-06-01得票数 2
回答已采纳
1回答
给定一个边长为L的正方形矩阵，如何在R中提取落入能够填充该矩阵的最大可能圆的所有值？ ? 我为C++找到了Filled circle in matrix(2D array)，但是如何测试矩阵中每个单元格的位置是否在等式中？如何在使用申请示例时知道每个单元格的X和Y？
浏览 26提问于2019-02-18得票数 3
回答已采纳
我正在尝试在MR上实现DBSCAN，并使用下面粘贴的链接中的直觉
http://codingwiththomas.blogspot.com/2012/05/distributed-dbscan-intuition.html
我的问题是如何计算初始距离矩阵。我不想仅仅为了计算距离矩阵而运行hadoop作业，并将其存储在内存中，因为我不会是一个好的设计。有什么建议吗。
浏览 0提问于2013-04-07得票数 1
回答已采纳
1回答
我想用一个相机来估计相对相机的姿态。
为实现这一目标，采取了以下步骤。
标定相机获得摄像机参数K
使用SIFT/SURF查找对应点
基本矩阵识别；使用findFundamentalMat()，并输入匹配点2D
用E=K‘’FK估计本质矩阵和修正奇异约束E
分解本质矩阵得到旋转，R= UWVt / UW'Vt (U和Vt是E的SVD )，得到平移t；
但这样做后，平移t太过small.assume，初始相机位置为(0，0，0).After摄像机向前移动20 is，下一姿态为(0，0，20)，实际平移应为20 is，两者的.The电平相差很大。
怎么了？如何从翻
浏览 2提问于2015-04-14得票数 1
3回答
我刚刚读了@David的简介，发现了一些有用的技巧来培养良好的R-编程技能/习惯，其中一条特别让我印象深刻。我一直认为R中的应用函数是使用dataframes的基石，但他写道：
如果您正在使用data.frames，请忘记有一个名为apply的函数--不管您做什么--不要使用它。特别是在边距为1的情况下(这个函数的唯一好用法是对矩阵列进行操作-边距为2)。
一些很好的选择：?do.call，?pmax/pmin，?max.col，？rowSums/rowMeans/等，令人敬畏的matrixStats包(用于矩阵)，?rowsum等等。
有人能给我解释一下吗？为什么应用函数不受欢迎？
浏览 0提问于2018-06-06得票数 7
回答已采纳
1回答
R应用返回子矩阵
、、、
鉴于这一守则：
test=matrix(c(1,2,3,4,5,6,7,8,9,10,11,12),4)
splitData=data.frame(first=c(1,3),second=c(2,4))
apply(splitData,1,function (x) {test[x[1]:x[2],]})
我得到这个矩阵：
[,1] [,2]
[1,] 1 3
[2,] 2 4
[3,] 5 7
[4,] 6 8
[5,] 9 11
[6,] 10 12
为什么不给我一个矩阵列表呢？
预期结果：
[[1]]
[
浏览 3提问于2015-01-26得票数 1
回答已采纳
3回答
迭代两个稀疏矩阵
、、、
我使用boost稀疏矩阵保存bool's，并尝试编写一个比较函数来将它们存储在map中。这是一个非常简单的比较函数。基本上，我们的想法是将矩阵视为二进制数(在被展平为向量之后)，并根据该数的值进行排序。这可以通过以下方式完成：
for(unsigned int j = 0; j < maxJ; j++)
{
for(unsigned int i = 0; i < maxI; i++)
{
if(matrix1(i,j) < matrix2(i,j) return true;
else if(matrix1(i,j) > matrix2(i
浏览 6提问于2009-08-25得票数 3
我有向量v和矩阵M。我希望通过将每一行与欧几里得距离( M v，np.linalg.norm)进行比较，并从最小距离到最大距离进行排序，从而按行排序。我研究过np.argsort，但不知道如何将它应用于这个实例，或者它在这里是否有效。任何帮助都将不胜感激！
示例输入：
k = 5
n = 3
M = np.random.randn(k, n)
v = [1, 2, 3, 4, 5]
浏览 7提问于2022-02-11得票数 -1
回答已采纳
2回答
我有一个名为A的cell(1, n)单元格数组，每个单元格条目都包含一个mxn矩阵。因此，实际上，我的单元格数组包含大小为mxn的n矩阵。
然后我有了另一个名为B的单元数组，其中存储了n pxm矩阵。
我需要做的是将两者相乘，例如：A[1] * B[1]，A[2] * B[2]，...，A[n] * B[n]。然后，我需要将结果存储为它们自己的单个矩阵，并对它们求和。
矩阵对于乘法是保角的，但是因为单元阵列B比单元阵列A包含的行数少，所以当我使用cellfun(@times A, B, 'UniformOutput', true)时，我得到了一个unequal matrices
浏览 0提问于2012-12-01得票数 6
回答已采纳
我有一个非常简单的问题。
假设我们有下面的代码来计算各向同性天线的方向性。
ThDeg = 0:5:180;
dtheta = 5*pi/180;
dphi = 5*pi/180;
Th = ThDeg*pi/180;
% the above are the angles at which the following matrix is acquired in practical case. In this case we take a matrix of all ones.
U_iso = ones(72, 37); % our matrix assumed
omega_iso =
浏览 3提问于2010-07-04得票数 0
我刚从R开始，遇到了一个问题。
我想hclust一个包含大量数据的表。
数据由一个矩阵组成，该矩阵具有：
> nrow(hell)
[1] 202397
> ncol(hell)
[1] 39840
矩阵由整数组成(大部分是0，当然也有更高的值)。
现在，我设法在一台48核、280 it的linux机器上读取了我的15 it文件，但如果我想要计算距离矩阵，它会失败，并显示以下错误：
d <- dist(as.matrix(hell))
Error in unlist(x,recursive, use.names):
resulting vector exceeds vecto
浏览 4提问于2013-05-29得票数 2
1回答
我是ganna，用均方位移(MSD)来表示3-dim矩阵.这个解计算(x，y)的MSD，但我有一个(x，y，z)的矩阵。说：
N维的推导：
MSD=2nDt。
这段代码用于n = 2
import numpy as np
r = np.sqrt(xdata**2 + ydata**2)
diff = np.diff(r) #this calculates r(t + dt) - r(t)
diff_sq = diff**2
MSD = np.mean(diff_sq)
因此，如果我想利用n = 3的这段代码，可以将其更改为以下代码：
im
浏览 2提问于2020-09-21得票数 0
我已经创建了一个由50个随机生成的(x，y)对组成的矩阵，我需要将它转换成一个矩阵，它显示从任意点(行)到所有其他点(列中)的距离。这应该导致对角线0，并且是一个50x50矩阵。
我将使用公式计算距离：
创建距离矩阵：
n = 50
x = round(runif(n)*1000)
y = round(runif(n)*1000)
coordinates = cbind(x,y)
head(coordinates)
x y
[1,] 266 478
[2,] 372 861
[3,] 573 438
[4,] 908 245
[5,] 202 71
[6,] 898
浏览 4提问于2019-09-24得票数 0
我有一个自定义函数来计算图像的两个像素之间的权重(表示图上的节点)。
function [weight] = getWeight(a,b,img, r, L)
ac = num2cell(a);
bc = num2cell(b);
imgint1 = img(sub2ind(size(img),ac{:}));
imgint2 = img(sub2ind(size(img),bc{:}));
weight = (sum((a - b) .^ 2) + (r^2/L) * abs(imgint2 - imgint1)) / (2*r^2);
其中a =
浏览 2提问于2015-04-10得票数 1
1回答
矩阵向另一个矩阵的旋转
、、、
我有两个4x4旋转矩阵M和N。M描述了我当前在空间中的物体姿态，N是期望的物体姿态。现在我想将M矩阵旋转到N，这样在接下来的迭代中，对象将慢慢旋转到所需的位置。你知道怎么解决这个问题吗？
浏览 0提问于2018-04-19得票数 0
1回答
我有一个带N个元素的向量V，每个元素都可以是介于1和N之间的整数。给定这个向量，我想构造一个N×n矩阵W，其中列I包含1到N之间整数的频率，如它们出现在子向量V1:i中。
例如，假设N=5和n=7以及V=c(3,1,4,1,2,1,4)。那么我的矩阵W就会有元素
0,1,1,2,2,3,3
0,0,0,0,1,1,1
1,1,1,1,1,1,1
0,0,1,1,1,1,2
0,0,0,0,0,0,0
因为整数1(第一行)在V1中出现0次，在V1:2中出现一次，在V1:3中出现一次，在V1:4中出现两次，在V1:5中出现两次，在V1:6中出现三次，在V1:7中出现三次，等等。
浏览 0提问于2018-05-29得票数 1
回答已采纳
我一直在努力想办法解决这个问题，但没有取得多大的成功。我正在使用一个带有二进制数据(0和1s)的表。我设法使用R函数dist(x,method="binary")从数据中估计出一个距离矩阵，但我不太清楚这个函数是如何准确地估计距离矩阵的。它是否使用Jaccard系数J=(M11)/(M10+M01+M11)？
浏览 3提问于2014-05-15得票数 1
回答已采纳
1回答
我想在MATLAB中对文本进行分层聚类。我有四句话，
I have a pen.
I have a paper.
I have a pencil.
I have a cat.
我想把以上四句话聚在一起，看看哪一句比较相似。我知道统计工具箱有像pdist这样的命令来测量成对的距离，linkage用来计算集群的相似性等等。一个简单的代码如下：
X=[1 2; 2 3; 1 4];
Y=pdist(X, 'euclidean');
Z=linkage(Y, 'single');
H=dendrogram(Z)
工作良好，并返回树状图。
我想知道我能不能像上面提到的
浏览 1提问于2010-09-05得票数 3
回答已采纳
1回答
假设我有一个2x3矩阵A
1 2 3
4 5 6
以及长度为4的向量y：
0 1 2 1
以及另一个4x2矩阵B
0 0
1 1
2 2
3 3
我希望通过添加来自A行的B来多次更新B的列。并给出了需要更新的A列的索引。
使用for循环，可以这样做：
for i in np.arange(4):
A[:,y[i]] += B[i,:]
我使用ufunc.at实现了这一点，如：
np.add.at(A.T,y,B)
然而，ufunc.at的性能几乎与使用for循环一样糟糕。
如何获得不同的矢量化实现？
使用A[:,y]+=B.T进行更新似乎只更新每一列一次。
浏览 5提问于2017-06-07得票数 1
回答已采纳
2回答
在JMP软件中，当行数大于2000时，可以使用“快速病房”方法。来自快速病房
应用一种算法，对大量行更快地计算Ward的方法。计算时间更短，因为该算法不需要计算距离矩阵。只要有超过2000行，就会自动使用。
也做同样的事..。“在使用病房链接方法创建的层次簇树中查找最多四个群集。指定'SaveMemory‘为'on’来构造集群，而不计算距离矩阵。否则，如果机器没有足够的内存容纳距离矩阵，则可能会收到内存不足的错误。”
我在Python中寻找类似的东西，但它们似乎都需要提前计算的距离矩阵(对于275 k行和10列的问题，这需要大量的内存)。但是，在JMP/Matlab中，它在内存
浏览 8提问于2020-02-11得票数 0
回答已采纳
1回答
创建相似矩阵
、、、
我有一个训练集和一个向量测试集。所有的向量都有标记。
对于测试集中的每个标记向量，训练集中有3个具有相同标记的向量。
我用余弦距离来计算向量中元素之间的相似性。在图中，我们可以看到在测试集的6个向量和训练集的18个向量的子集中应用余弦距离相似性的结果。
📷
现在，我想创建标签的相似性矩阵。因此，在这种情况下，我需要一个6x6维的矩阵，但我不知道如何将这个分数矩阵转换为一个相似矩阵。
浏览 0提问于2017-09-15得票数 1
回答已采纳
假设我有1,000万个条目，每个条目都有一个100维的实数向量(实际上它们是word2vec嵌入)。对于每个项目，我想得到(大约)前200个最相似的项目，使用余弦相似。我目前在Hadoop (hive)中作为UDF函数实现的余弦相似性标准需要花费大约1s来计算1项的余弦相似性，而其他项目则有1,000万项。这使得运行整个矩阵是不可行的。我的下一步是在Spark上运行它，使用更多的并行化，但是它仍然不能完全解决问题。
我知道有一些方法可以减少spars矩阵的计算。但我的矩阵是，而不是稀疏。
如何有效地为每个项目获得最相似的项目？是否有一个近似的余弦相似，将是更有效的计算？
浏览 2提问于2020-07-22得票数 0
我知道Matlab中的pdist(X，distance)用于获取点的(nxd)矩阵，并计算它们之间的成对距离。我还知道，如果矩阵包含向量而不是点，那么它有一个额外的选项来计算余弦距离。
我想做的是取一个(n X D X t)矩阵，它保存样本随时间t的变化位置，并高效/快速地计算所有对/所有帧之间的余弦方向距离，其中a v(t)定义为由p(t+1) - p(t)计算的方向，p(t)指的是行M(p，:，t)。
显然，如果可以帮助的话，我不想使用循环。有什么建议吗？
任何帮助都非常感谢。
浏览 0提问于2012-04-22得票数 0
因此，我有以下限制：
如何用MATLAB高效地编写这篇文章？输入是x_mn、M和N。集合B={1,...,N}和集U={1,...,M}
我这样做了(因为我把x写成了折叠向量)
x=[x_11, x_12, ..., x_1N, X_21, x_22, ..., x_M1, X_M2, ..., x_MN]
%# first constraint
function R1 = constraint_1(M, N)
ee = eye(N);
R1 = zeros(N, N*M);
for m = 1:M
R1(:, (m-1)*N+1:m*N) = ee;
e
浏览 2提问于2014-06-18得票数 1
回答已采纳
假设我有3x3的单元格，每个单元格包含一个9x9的矩阵，我该如何使用cellfun对整个单元格数组的每一行求和呢？
当我尝试在花括号中使用:时，我不断得到错误“错误的单元格引用”。
我不想把它转换成矩阵，然后再转换回单元格。
非常感谢你们的智慧！
浏览 0提问于2012-08-17得票数 0
回答已采纳
2回答
我在R中有一个包含1000列和4行的矩阵。矩阵中的每个单元格都包含1-4之间的整数。我想知道两件事：
1)以任意顺序包含"1“、"2”、"3“和"4”的列数是多少？理想情况下，我希望代码不要求我输入1,2,3,4的每个可能的组合来执行计数。
2)包含3个可能的整数，但不包含所有4个的列数是多少？
浏览 0提问于2016-06-08得票数 1
2回答
我很难弄清楚向量是如何格式化的。我需要在ISwR库的cystfibr包中找到参与者的平均身高。当打印整个高度数据集时，它看起来是一个带有高度值的21x2矩阵和一个表示性别的1或2。但是，ncol返回的值为NA，表明它是一个向量。尝试获取矩阵的特定索引(heightdata1，)也会返回错误的维数错误。
我希望只对向量中的高度值求和，但当我运行代码时，我得到了男性和女性整数的总和。(25)
install.packages("ISwR")
library(ISwR)
attach(cystfibr)
heightdata = table(height)
print(heightd
浏览 2提问于2020-05-24得票数 2
1回答
给定一个矩阵A，我想用手完成多维缩放，而不是使用任何给定的R函数。
因此，我使用以下代码计算了居中矩阵B：
n<-nrow(A)
id<-diag(n)
e<-diag(id)
H <- id - (1/n)*e %*% etranspose
B <- (-1/2)* H %*% A %*% H
我的问题是:如果没有cmdscale函数或类似的东西，我如何才能使用B矩阵来完成对A矩阵的多维缩放？
浏览 0提问于2015-11-16得票数 5
在下面的代码中，我有4个节点的x和y坐标。前两个节点，即1和2属于类别c(矩阵c)，而后两个节点，即3和4属于类别d(矩阵d)。"dist“表示所有节点之间的距离。在matlab中有没有什么命令可以用来判断c中的任何节点与d中的任何节点之间的距离是否小于R，例如
if distance of any of the nodes in C from any of the nodes in D > R
%do this
end
我可以使用for循环来做这件事，但是要寻找一种更短的方法。谢谢
x=[1 2 4 4];
y=[3 5 6 1];
range=R
dist=[0.7 1.6
浏览 2提问于2017-06-10得票数 0
2回答
我需要测量两个n维向量之间的距离。似乎Mahalanobis距离在这里是一个很好的选择，所以我想试一试。
我的代码看起来像这样：
import numpy as np
import scipy.spatial.distance.mahalanobis
x = [19, 8, 0, 0, 2, 1, 0, 0, 18, 0, 1673, 9, 218]
y = [17, 6, 0, 0, 1, 2, 0, 0, 8, 0, 984, 9, 30]
scipy.spatial.distance.mahalanobis(x,y,np.linalg.inv(np.cov(x,y)))
但是我得到了这
浏览 0提问于2013-06-19得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
