在R中创建矩阵，有没有更快的方法？
在R中创建矩阵的常见方法是使用函数matrix()，但如果需要更快的方法，可以考虑使用以下两种方式：
- 使用
array()函数：array()函数可以创建多维数组，包括矩阵。相比于matrix()函数，array()函数在创建矩阵时更快。以下是使用array()函数创建矩阵的示例代码：
n <- 1000 # 矩阵的维度
m <- 1000
mat <- array(0, dim = c(n, m)) # 创建一个n x m的矩阵，初始值为0- 使用矩阵运算函数：R中的矩阵运算函数通常比直接创建矩阵更快。例如，可以使用
outer()函数创建一个矩阵，该矩阵的元素是两个向量的运算结果。以下是使用outer()函数创建矩阵的示例代码：
n <- 1000 # 矩阵的维度
m <- 1000
vec1 <- 1:n
vec2 <- 1:m
mat <- outer(vec1, vec2, FUN = "*") # 创建一个n x m的矩阵，元素为vec1和vec2对应位置的乘积这些方法都可以更快地创建矩阵，具体选择哪种方法取决于具体的使用场景和需求。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云云服务器（CVM）：https://cloud.tencent.com/product/cvm
- 腾讯云云数据库 MySQL 版：https://cloud.tencent.com/product/cdb_mysql
- 腾讯云对象存储（COS）：https://cloud.tencent.com/product/cos
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云物联网通信（IoT Hub）：https://cloud.tencent.com/product/iothub
- 腾讯云移动推送（TPNS）：https://cloud.tencent.com/product/tpns
- 腾讯云区块链服务（BCS）：https://cloud.tencent.com/product/bcs
- 腾讯云元宇宙（Tencent Real-Time Rendering）：https://cloud.tencent.com/product/trr
相关·内容
1回答
有没有人知道一种快速创建矩阵的方法，就像在R中创建如下矩阵一样。
[,1] [,2] [,3] [,4]
[1,] 1 1 1 1
[2,] 1 2 2 2
[3,] 1 2 3 3
[4,] 1 2 3 4
上面的矩阵是4x4，我想创建类似10000x10000的东西。
浏览 0提问于2013-07-10得票数 3
回答已采纳
2回答
我正在自己学习R，在尝试使用markovchain包在Rstudio中构建转移概率矩阵时遇到了一些问题。首先，我尝试计算DNA序列的转移概率。
ATTCAACACATCCAGCCACATGCTCCGAGAGGAGGCAGAGGGCCCCCGGAATGATGCTTACCGAGATTCTTGTTTTTATCCTCGTGGTTGTTTAAAAACGAGTTGAAACTGACGGCATGTCGGACTATAAGCTACTTACTCACCATAGACGTGACCATAGGCCCTAAAACGTTACCGAGATATTCACTTCTAATAACAGTTGTCGGCAGAGCCAAAAGGCCGGGTGAT
浏览 3提问于2017-11-16得票数 1
2回答
在处理大型SciPy CSR稀疏矩阵时，我注意到将矩阵切片以从矩阵中获得单行的速度非常慢，因为它似乎是在复制。
有没有办法创建一个稀疏矩阵，它引用现有的行而不是复制它，也许有比CSR矩阵更合适的实现？
我的实现需要的是元素和行的快速查找，以及向量的所有非零索引的快速查找。我永远不需要以任何方式更改矩阵或对矩阵执行其他操作。
浏览 0提问于2017-07-03得票数 1
我有一个包含笛卡尔坐标点的3D矩阵。在下面显示的示例中，给出了尺寸为10x10x10x3的这种矩阵的缩小版本(a)。然后我有一个大小为3x3的旋转矩阵R，我想对矩阵a中给出的每个3x1位置向量做一个点积，例如a[0,0,0]是最左上角坐标点的位置向量。 一种简单的方法是使用for循环- import numpy as np
a = np.ones([10, 10, 10, 3])
R=np.ones([3,3])
a=np.reshape(a,(1000,3))
b=np.array([np.dot(R,xyz) for xyz in a])
b=np.reshape(b, (10,10,
浏览 16提问于2021-07-15得票数 0
1回答
我正在尝试计算数据集中unigram的互信息。在尝试这样做的时候，我试图提高在numpy ndarray中循环的速度。我有以下代码，其中我使用了一个已经创建的矩阵'C‘，它具有6018行和27721列，以便计算PMI矩阵。有什么想法可以提高for循环的速度(目前它几乎需要4小时才能运行)？我在其他一些关于使用Cython的文章中读到过，但是有没有其他选择呢？提前说一句，谢谢你的帮助。
# MAKE MUTUAL INFO MATRIX, PMI
print "Creating mutual information matrix"
N = C.sum()
invN = 1
浏览 1提问于2015-01-26得票数 0
1回答
我想要获取data.framework/矩阵的列，并在dataframe的每个单元格([i, j])之间应用一个函数，其中i和j是沿着data.frame列的序列。基本上，我希望以与cor函数与data.frame相同的方式填充单个单元格的矩阵。
这是一个相关的问题:然而，，我在随机化测试中使用这个方法，并多次重复操作(创建许多矩阵)。我在找最快的方法来做这个手术。我已经用并行处理加快了速度，但是我仍然对这个速度不满意。也不能假设矩阵输出是对称的，即cor生成对称矩阵的方式(我的例子将反映这一点)。
我今天在data.table网页()上看到了以下内容：
500+比DF[i,j]<-v
浏览 4提问于2014-05-22得票数 3
回答已采纳
2回答
给定两个矩阵，我想创建每行的平方差和的新数组，但我似乎找不到一种方法。 为了更清楚地理解我的意思，让我们举个例子。我想在numpy矩阵计算中执行以下for循环： a = np.array([[0.0, 0.0, 0.0], [1.0, 1.0, 1.0]])
b = np.array([[1.2, 2.3, 3.4], [4.5, 5.6, 7.8], [9.10, 10.11, 11.12]])
summed = np.ones((2,3))
for i, aSample in enumerate(a):
for j, bSample in enumerate(b):
浏览 38提问于2021-09-13得票数 1
回答已采纳
我会先说我是一个R菜鸟，我认为这可能有一个简单的解决方案，但我正在努力寻找它。
我有一个2列1,000行的矩阵。保持行固定，我想创建一个新的变量，从2列中随机选择一个元素。例如，制作一个简单的矩阵：
matrix(c(1,1,4,6,1,3,2,1,1,7), ncol=2)
[,1] [,2] [,3]
[1,] 1 3 3
[2,] 1 2 1
[3,] 4 1 4
[4,] 6 1 1
[5,] 1 7 7
在上面的简化矩阵中，第三列(我刚刚
浏览 0提问于2011-11-08得票数 3
1回答
我缺乏对基本面的深入了解，这对解决这些问题的挑战造成了影响。
HackerRank矩阵旋转问题是一个非常有趣的问题。我推荐那些试图丰富他们的编码技能的人使用hackerrank ()
问题的总结是，给出了整数的R矩阵，其中R和C的最小值必须是偶数。你必须逆时针旋转矩阵x次。旋转适用于矩阵的元素，而不是矩阵的维数，以防不清楚。
所以我用两种算法解决了这个问题。它们都是非常相似的，因为你可以想象像洋葱层一样的矩阵，在每一层中循环，并旋转该层中的元素。旋转的次数仅仅是x%(该层中元素的计数)，所以如果给您x=1,000,000，那么重复完全旋转是没有意义的。
第一个最快的是：
第二种方法不是循环循环
浏览 9提问于2016-03-07得票数 1
回答已采纳
在MATLAB中，我尝试在多个X和Y上使用内置函数[r, p] = corr(X1,Y1);的输出创建一个矩阵，然后将所有r和P合并到它们各自的矩阵R和P中。
R = [corr(X1,Y1), corr(X2,Y2); (...)
corr(X3,Y3), corr(X4,Y4)];
默认情况下，corr返回r值。有没有办法也能在p上做到这一点？下面是我做的很长的方法，我只是想知道有没有像上面这样更短更容易的方法。
首先找到每个r和p
[r1, p1] = corr(X1,Y1);
[r2, p2] = corr(X2,Y2);
[r3, p3] = corr(X3,Y3);
..
浏览 0提问于2013-02-22得票数 1
3回答
我需要从权重矩阵的值开始创建一个矩阵。在创建和迭代矩阵时，哪种结构在速度上是保持矩阵的最佳结构？我在考虑一个列表或一个numpy 2D数组，但它们在我看来都很慢。我需要的是：
numpy array
A = np.zeros((dim, dim))
for r in range(A.shape[0]):
for c in range(A.shape[0]):
if(r==c):
A.itemset(node_degree[r])
else:
A.itemset(arc_weight[r,c])
或
list
浏览 1提问于2014-05-29得票数 1
回答已采纳
2回答
我正在写R代码，它使用编译后的C代码。
从“编写R扩展”文档中，我了解到有许多R可执行文件/DLL可以从C代码中调用。头文件‘Rmath.h’列出了许多可用的函数，其源代码列在此网站上：
我需要计算许多矩阵的奇异值分解，但是我在上面的网站上找不到这样做的子例程。(所以我假设Rmath.h不包含SVD子例程)有没有在C- in -R代码中进行特征值计算的简单方法？
非常感谢。
浏览 0提问于2013-01-05得票数 5
回答已采纳
1回答
我有这行R代码：
croppedDNA <- completeDNA[,apply(completeDNA,2,function(x) any(c(FALSE,x[-length(x)]!=x[-1])))]
它所做的是识别DNA序列矩阵(1行=1序列)中不通用(信息)的位点(cols)，并从矩阵中对它们进行子集，以形成新的“裁剪矩阵”，即去掉其值相同的所有列。对于大数据集，这需要大约6秒。我不知道我是否可以在C++ (C++的初学者)中做得更快，但这对我来说是很好的尝试。我的想法是使用Rcpp，循环遍历CharacterMatrix的列，拉出列(站点)作为CharacterVector
浏览 2提问于2013-05-15得票数 6
回答已采纳
给定一个2x2矩阵r = np.array([[1, 0], [0, 0]], dtype=complex)，我如何创建一个块矩阵，使我获得N个r作为它的条目，以便我得到
np.array([r, r,..., r, r])
有没有像np.something(r, 5)一样的东西，所以它提供了np.array([r, r, r, r, r])
浏览 0提问于2019-07-21得票数 0
1回答
我需要为我的工作生成大量的随机均值不变正交矩阵。均值不变矩阵具有A*1_n=1_n性质，其中1_n是标量1基本np.ones(n)的大小为n的向量.我使用Python，特别是Numpy来创建我的矩阵，但是我想确保我的方法是正确的，也是最有效的。另外，我还想介绍我尝试过的三种不同的正交化方法的结果，并希望能解释为什么一种方法比另一种方法更快。在这篇文章的结尾，我问了四个关于我的发现的问题。
一般情况下，为了创建一个均不变的随机正交矩阵A，需要创建一个随机平方矩阵M1，将其第一列替换为1列，并将该矩阵正交化。然后，再用另一个矩阵M2这样做，最终的均值不变随机正交矩阵A= M1*(M2.T)。这一过
浏览 20提问于2021-02-26得票数 0
2回答
我已经为这个结果做了一些挖掘，但这里的大多数问题都有关于cbind函数和基本矩阵连接的信息。我想要做的事情有点复杂。
例如，假设我有一个NxM矩阵，它的第一列是每一行的惟一标识符(幸运的是，在本例中是按该标识符排序的)。出于与这个查询无关的原因，我将这个矩阵的行分成( n_i )个xM矩阵，使得n_i的和= N。
我打算对每个子矩阵分别运行分析，然后使用唯一标识符将数据再次组合在一起。
举个例子:假设我有一个10xM的矩阵data。拆分后，我将收到矩阵subdata1和subdata2。如果您要查看矩阵的内容：
data[,1] = 1:10
subdata1[,1] = c(1,3,4,6,
浏览 0提问于2014-06-25得票数 0
2回答
大尺寸多矩阵
、、
我想尝试在matlab中计算三个矩阵的乘积。矩阵的形成说明如下：
L = D^(-1/2) * A * D^(-1/2);
D、A和L是n*n矩阵。A和L是，不是对角线或稀疏，而是D是对角线。在这种情况下，n = 16900。当我在matlab中计算L时，需要很长时间，大约4个小时！
我的问题是:有没有更有效的方法来计算L？
浏览 1提问于2014-02-13得票数 0
回答已采纳
1回答
无环的三维变维对象的矩阵和
、、、、
我在R中有一个三维对象，它包含n个方阵.例如：
myObject[,,1] # returns a square matrix
myObject[,,2] # returns a square matrix of the same size
...
对象中的所有矩阵都是相同大小的。我想把所有矩阵加在一起，没有循环。如果我知道对象中有多少个矩阵，这就足够简单了。例如：
matrixSum <- myObject[,,1] + myObject[,,2] + myObject[,,3]
问题是，我需要对几千个这样的对象这样做，并且每个对象中有一个可变的矩阵数。我有没有办法不用循环就能做到这一
浏览 1提问于2014-10-31得票数 1
回答已采纳
2回答
我有一个实数的M乘N矩阵R。输入Rm，n给出了m个月的股价回报和n个月的模拟，我想从所有的模拟中找出最早的一个月出现负回报。
我可以为每一列找到第一次出现负返回，然后找到其中的最小值。但有没有更有效的方法？
问题：在R中最有效的方法是什么
浏览 1提问于2015-05-03得票数 0
我正在研究概率矩阵分解。我想将评分矩阵分为训练和验证。我看到的代码不是随机样本。
我一直在努力想办法做到这一点。
举个简单的例子，原始矩阵R是(10,10)。我正在寻找一个训练矩阵，它是(10,10)，但只有矩阵R的值的70%随机采样，其他值都是0。并且验证矩阵也是(10,10)，但是具有剩余的30%以上的值，其余的值为0。
这里的关键是保持矩阵的大小不变。
有人知道怎么做吗？
编辑:考虑更多: TrainingMatrix + ValidationMatrix =矩阵R，如果你做元素加法
浏览 10提问于2020-05-15得票数 0
我想要创建三维数组，它是2d数组的函数，并在每个2D数组上应用矩阵操作。现在，我正在使用for循环来创建一系列2d数组，如下所示：
for i=1:50
F = [1 0 0; 0 i/10 0; 0 0 1];
B=F*F';
end
有没有一种方法可以不使用for循环来完成此操作？我试过这样的事情：
F(2,2) = 0:0.1:5;
以及：
f=1:0.1:5;
F=[1 0 0; 0 f 0; 0 0 1];
如果不使用循环来创建它们，但是两者都会产生维度不一致的错误。
我还想在代码中对F执行矩阵运算，例如
B=F*F';
并想把F的某些部分画成其他函
浏览 2提问于2015-08-27得票数 1
回答已采纳
我目前正在将R代码转换为Julia (0.6)，我想知道是否有Julia函数可以让我创建一个可与ngCMatrix类的R稀疏矩阵相媲美的对象。它们是非常有效的，因为它们只承载非零数的位置，而不必实际存储非零数。
sparseMatrix(1:10,1:10)
10 x 10 sparse Matrix of class "ngCMatrix"
[1,] | . . . . . . . . .
[2,] . | . . . . . . . .
[3,] . . | . . . . . . .
[4,] . . . | . . . . . .
[5,] . . . .
浏览 9提问于2018-07-24得票数 3
4回答
我有一个矩阵，我想创建一个新的矩阵，它将是旧的矩阵，但没有第一行和第一列。有没有办法在不使用循环的情况下做到这一点？
浏览 0提问于2010-08-18得票数 1
回答已采纳
我有一个值的向量v和一个索引的向量r。我想将这些值存储在矩阵m中，如下所示：
for i = 1:length(v)
m(i, r(i)) = v(i);
end
以向量化的方式完成这项工作的最快方法是什么？
浏览 0提问于2012-01-09得票数 0
回答已采纳
2回答
关于R中的矩阵，我有一个非常愚蠢的问题:如何检查一个矩阵是否满R？
例如：
m1 = matrix(NA, 2, 2)
m1[,1] = c("a", "b")
因此此时，这个矩阵应该如下所示：
a NA
b NA
我试着使用is.null()检查它是否为空，并返回False，这正是这里所期望的。然而，我不知道R中是否有一个函数，它可以判断矩阵中的每个条目是否都有赋值(不管它是什么，整数，字符，只要它有赋值)？
提前感谢！
浏览 6提问于2017-06-28得票数 0
回答已采纳
R版本的igraph包包含distances函数，根据文档，该函数返回传递给它的图的距离矩阵。 但是，我在这个包的python版本中找不到这个函数或类似的函数。我能找到的唯一提供的函数是Point.distance，在我对文档的理解中，它只计算2D平面中两点的距离，而不是网络图。 在这个包的python版本中，有没有一个函数可以返回图的距离矩阵？ 从本质上讲，我正在寻找的是以下内容： # create the graph
graph = Graph(n=5, edges=[(0,1), (1,4), (2,3), (2,4), (2,5) , (3,5), (4,0)], directed=
浏览 56提问于2021-11-23得票数 0
回答已采纳
假设我想从头开始创建一个3D环境，没有任何类库。我的直觉是创建一个3D矩阵(一个向量，因为我使用的是AS3)，然后操纵对象和相机在矩阵中的位置。然后我会计算活动摄像头的位置以及它应该能看到的矩阵的多少，并将矩阵的这些部分(对象)渲染到视图中。
当然，我会有一个单独的矩阵来保存对象的位置，以便有效地访问。
这是有效的吗？处理器是否能够在单线程AS3应用程序中处理这些计算？如果不是，什么是更好的方法？
浏览 1提问于2014-02-07得票数 0
1回答
嗨，我需要在R中绘制一个箱线图。我有两个矩阵a和b。我为a创建了一个箱图，并想在a的同一个图上为b创建箱图。b矩阵的框图应该位于a的框图的胡须上。
有没有办法在R里做?？
浏览 0提问于2012-04-02得票数 6
回答已采纳
有没有一种方法可以在R中创建一个矩阵，而该矩阵中的每个元素都是另一个矩阵？我以前用Python做过，但是当我做的时候
X <- matrix(rep(0,200),nrow=200,ncol=1)
for (i in 1:200){ X[i,] <-matrix(rep(0,32),nrow=8,ncol=4)}
它不适用于R。
谢谢!
浏览 2提问于2018-11-29得票数 5
回答已采纳
1回答
有人知道如何在Matlab中实现单元格(2，2)的R当量吗？
在Matlab中，这创建了一个2x2“矩阵”，其中每个单元可以是任何类型的数据，比如另一个矩阵或其他什么。
因此，基本上，它可以是一个矩阵，如果这是用户想要的。
有没有办法在R中做到这一点？
浏览 2提问于2014-08-12得票数 3
回答已采纳
4回答
我正在尝试创建R中的块循环矩阵。块循环矩阵的结构如下所示。
C0 C1 ... Cn-1
Cn-1 C0 C1 ... Cn-2
Cn-2 Cn-1 .... Cn-3
and so on
我拿到积木了
C0 .... Cn-1
创建矩阵的最简单方法是什么？有没有已经可用的功能？
浏览 1提问于2013-07-04得票数 5
回答已采纳
1回答
有人知道有什么好的稀疏矩阵库吗？我需要它在大型稀疏矩阵(10,000 x 10,000)上做kronecker积和乘法。
现在我们正在使用R，它很好地处理它们。不幸的是，我们不能把R投入生产，所以寻找替代品。我一直在玩Math.net数字，但看起来有点慢。
浏览 4提问于2012-08-01得票数 1
我想找出一个无向图的两个顶点之间固定长度"K“的不同可能路径的数量，我看到了许多使用邻接矩阵的实现，但有没有办法在没有邻接矩阵的情况下做到呢？如果有办法，请给点提示。
浏览 13提问于2021-05-02得票数 0
回答已采纳
在Matlab中，不需要任何编码，我就可以创建一个矩阵，打开它的电子表格，从Excel中复制多列数值并将它们粘贴到电子表格中。然后，我可以右键单击这个矩阵，并立即绘制它。
我试着在googling上搜索如何在R中做同样的事情，似乎所有的事情都涉及到创建一个函数，用for循环迭代每个值。这看起来有点麻烦，在RStudio中有没有同样简单的方法来做到这一点呢？
谢谢。
浏览 1提问于2015-02-15得票数 0
1回答
无环矩阵的优化重组
、、、、
TL；博士:我正试图在Matlab中优化下面的简短代码。因为它涉及大矩阵上的循环，所以太慢了。
for i = 1:sz,
for j = 1:sz,
if X(j) == Q(i) && Y(j) == R(i),
S(i) = Z(j);
break
end
end
end
Specifics：基本上，我从x、y和z数据的三个向量开始，我想把它们绘制成一个曲面。生成x和y数据的网格，然后为相应的z值创建一个矩阵。
[X, Y] = meshgrid(x, y);
Z = grid
浏览 0提问于2019-03-15得票数 0
回答已采纳
1回答
出于速度的原因，我正在考虑将R脚本转换为C代码，并将其打包为.exe。我是C语言的新手。
我的问题是，在C中它会明显更快吗？速率限制步骤是一种排序算法，它必须多次应用于大向量。我不确定R中的矢量化功能是否会对此有所帮助，还是会减慢速度。我还读到过for循环在R中的效率很低。
如果我应该用C语言来做这件事，什么库可以帮助我模仿R的一些数据处理函数，比如基本的矩阵操作？我应该从哪里开始呢？现在我甚至不知道如何将我的数据读入C(逗号分隔的文本文件)。
浏览 0提问于2012-04-13得票数 10
回答已采纳
2回答
寻找图像部分图像的坐标
、、、
我正在寻找一种方法，以找到我的屏幕上的一个按钮的位置使用matlab。现在，我已经为此创建了以下代码:用这段代码，我对我的计算机做了一个屏幕截图，在matlab中加载它，并使它加倍。然后对目标图像做同样的操作，需要在我的屏幕上找到的图像，然后尝试匹配这个矩阵，以找到矩阵和总矩阵，并显示该图像左上角的坐标。问题是，用这种方法和1920*1080的屏幕分辨率，大约需要15分钟才能通过这个完整的过程。
因此，我的问题是，是否有一种更容易/更快的方法在我的屏幕上或在另一张图像中找到一个图像的中心？到目前为止，我还没有发现，即使在使用图像处理工具箱。
clear all
close all
robo =
浏览 3提问于2014-03-03得票数 0
回答已采纳
1回答
我有一个矩阵A = [1 2 4 4 8 8 8 4 1 7 7 8 8 9]，我想要创建一个新的矩阵，它在A中有相同的元素数。
我有2 1，1 2，3 4，2 7，5 8和1 9。
我的新矩阵应该是数字；每个数字的数量。
newmatrix = [1 2 4 7 8 9; 2 1 3 2 5 1]
如何从A创建新的矩阵？
浏览 4提问于2014-08-18得票数 2
回答已采纳
我有一个存储在350x130矩阵中的值集合，我正在用它在Matlab中创建一个简单的图：
figure;
plot(r)
有没有办法用另一个350x1矢量，v来给r生成的线条着色？V包含一个整数集合，在本例中，这些整数从1到4不等。
我想象这个调用应该是这样的，其中r是350x130，v是350x1：
figure;
plot(r,v)
浏览 0提问于2018-02-07得票数 1
回答已采纳
1回答
我问题的基本目的是如何使用R包在Matrix中实现矩阵操作的最佳性能。特别是，我想并行化运算(乘法)，并使用CUDA GPU上的计算处理稀疏矩阵。
详细信息
根据Matrix包在R cran中的文档
矩阵类的一种丰富的层次结构，包括三角矩阵、对称矩阵和对角线矩阵，既密集又稀疏，具有模式、逻辑和数字条目。使用'LAPACK‘和'SuiteSparse’库对这些矩阵进行大量的方法和操作。
看来，多亏了SuiteSparse，我应该能够使用GPU (CUDA)对稀疏矩阵执行基本操作。特别是，SuiteSparse的文档列出了以下内容：
SSMULT和SFMULT:稀疏矩阵乘
浏览 10提问于2017-11-27得票数 2
回答已采纳
2回答
我希望从我的矩阵中删除一个星号(*)，并将该矩阵写入文本文件，其余的元素将在没有空格或任何类型的分隔符的情况下相互连接。我写了这段代码
for b = 1 : length(out7num_r7_nt_back)
if ~(out7num_r7_nt_back(b) == '*')
out7num_r7_back(b) = '';
end
end
disp(out7num_r7_nt_back);
dlmwrite('my_data.txt',out7num_r7_nt_back, '');
我收到
浏览 2提问于2012-11-12得票数 3
1回答
解复矩阵方程
、、、
我有一个涉及矩阵的复杂方程式：
R = expm(X)*A + (expm(X)-I)*inv(X)*B*U;
其中R、B和U是已知矩阵。
I是一个单位矩阵。
我需要解决X的问题。有没有办法在MATLAB中解决这个问题？
浏览 0提问于2019-04-28得票数 1
1回答
我正在尝试计算P^100，其中P是我的转移矩阵。我想通过对角化P来实现这一点，这样我们就有了P = Q*D*Q^-1。 当然，如果我可以得到这种形式的P，那么我就可以很容易地计算出P^100 = Q*D^100*Q^-1 (其中*表示矩阵乘法)。 我发现，如果你只是做P^5，你得到的所有回报都是一个矩阵，其中你的每个P项都被提升到5次方，而不是矩阵的5次方(P*P*P*P*P)。 我在这里发现了一个问题，问的是如何检查矩阵是否可对角化，而不是如何显式构造矩阵的对角化。在MATLAB中，这非常简单，但是，我使用的是R，而不是MATLAB。
浏览 20提问于2020-02-03得票数 0
回答已采纳
我正在制作一个包含两个Rcpp函数的包。第一个函数用于创建一个矩阵，该矩阵将被第二个函数多次使用。矩阵存储在R的全局环境中，在调用这两个函数之间。
M <- myFirstRcpp(X)
P <- mySecondRcpp(M)
根据输入参数，第二个函数将对输入矩阵(由第一个函数创建)进行更改，然后再从输入矩阵中计算向量(aFunction是mySecondRcpp()中的C++ )：
IntegerVector aFunction( SEXP Qin, SEXP param ) {
NumericMatrix Q(Qin);
// Some changes mad
浏览 5提问于2014-01-22得票数 2
回答已采纳
我正在尝试实现SIFT算法。我有一个矩阵，比方说D，维度是(3，m，n)。对于矩阵中间层中的每个点，即A(2，:，:) (Matlab是1索引的)，我需要在该点的3x3邻域中找到局部最小值。我这样做的方式非常幼稚，如下所示：
img = D(2, :, :);
for r = 2:size(img, 1)-1
for c = 2:size(img, 2)-1
nbr = D(:, (r-1):(r+1), (c-1):(c+1));
if abs(D(2,r,c) - min(nbr(:))) <= eps(0.5) || abs(D(2,r,c)
浏览 2提问于2016-09-24得票数 0
1回答
我想将读入Stata。例如，我有一个由Python命令sio.mmwrite()创建的稀疏矩阵auth.mfx，我希望将其作为空间权重矩阵。相应的命令spatwmat需要一个.dta文件：
insheet using "auth_matrix.mtx", clear
只是将文件作为文本文件读取，而不是矩阵。有没有办法将矩阵市场文件读入Stata？显然，它们是R、Python、Matlab、Fortran等语言中公认的格式。
浏览 0提问于2015-05-21得票数 0
有时，当我想得到一个关于某些变量的分析表达式时，我会感到困惑。在以下情况下，一个是Numpy数组(T)，另一个是Sympy矩阵(X)。我知道直接乘它们不是个好主意，所以我决定把T转换成Sympy矩阵。然而，要想得到这个大型矩阵的结果需要很长时间。有没有更有效的方法？谢谢。
import numpy as np
import sympy as sp
T = np.random.rand(100, 5000)
x = sp.symbols('x:'+str(5000))
X = sp.Matrix(x)
W = sp.Matrix(T)
V = W * X
浏览 13提问于2022-05-23得票数 0
1回答
R基函数glm()采用Fishers评分法求解MLE，而glmnet则采用坐标下降法求解同样的方程。坐标下降比Fisher评分更省时，因为Fisher评分除了计算其他一些矩阵运算外，还计算二阶导数矩阵。这使得执行起来很昂贵，而坐标下降可以在O(np)时间内完成同样的任务。
为什么R基函数使用Fisher评分？这种方法比其他优化方法有优势吗？坐标下降和Fisher得分比较如何？我是比较新的做这个领域，所以任何帮助或资源将是有帮助的。
浏览 0提问于2014-07-03得票数 11
1回答
假设我有一个矩阵，每行都是一个标准的基向量，即每行恰好包含一个1，其他列都是0。
有没有一种方便的方法来创建这样的矩阵(即给出每一行中位置的向量)？
另外，有没有一种我应该表示这样的矩阵的方法，以便在倍频程中更有效地完成与它的乘法？
浏览 0提问于2011-11-18得票数 1
回答已采纳
3回答
在R中使用多核计算SVD
、、、、
我想在一个大的稀疏矩阵(17k X 2m)上运行R中的svd()，并且我可以访问一个集群。有没有一种简单的方法来计算使用多核的R的SVD？
RScaLAPACK包()似乎使这成为可能，但它似乎不再受到积极支持()，我认为这是有原因的。
浏览 5提问于2013-05-09得票数 5
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
