1回答
我有一个字符矩阵：
mat1
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] "0" "B" "A" "C" "D" "D"
[2,] "0" "0" "B" "C" "C" "C"
[3,] "0" "0" "0" "D" "D" "C"
[4
浏览 7提问于2020-07-01得票数 0
回答已采纳
2回答
假设我有一个矩阵。假设我有一个下三角矩阵的列表。如何计算矩阵中元素的数量？
例如,
Matrix <- c(0, 4, 3, 1, 2,
0, 0, 3, 2, 1,
0, 0, 0, 2, 1,
0, 0, 0, 0, 1,
0, 0, 0, 0, 0)
Matrix <- matrix(Matrix, 5, 5)
> Matrix
[,1] [,2] [,3] [,4] [,5]
[1,] 0 0 0 0 0
[2,] 4
浏览 0提问于2018-05-16得票数 4
回答已采纳
1回答
BLAS (基本线性代数子程序)提供了许多其他编程语言，比如我使用的Matlab，它有快速的例程来做矩阵乘法之类的事情。
然而，当将多个矩阵相乘时，存在一个“括号”矩阵的最佳顺序。摘自
例如，假设A是10×30矩阵，B是30×5矩阵，C是5×60矩阵。然后,
(AB)C = (10×30×5) + (10×5×60) = 1500 + 3000 = 4500次运算
A(BC) = (30×5×60) + (10×30×60) = 9000 + 18000 = 27000次运算。
本文接着讨论了这种乘法的最优顺序的求解方法。我的问题是，这些优化程序中是否有在BLAS中使用过？如果没有，如
浏览 7提问于2016-05-12得票数 4
访问Armadillo矩阵的上三角元素或下三角元素的最有效(即平衡内存和速度)方法是什么？我知道我可以为元素提供一个整数向量，但当矩阵变得非常大时，我希望避免携带另一个较大的向量。或者有没有一种有效的方法来快速创建下/上三角索引？
例如，使用5x5矩阵
// C++11 Initialization
arma::mat B = { 1, 2, 3, 4, 5,
6, 7, 8, 9, 10,
11, 12, 13, 14, 15,
16, 17, 18, 19, 20,
浏览 0提问于2015-03-11得票数 2
1回答
我有3个矩阵，它们都是稀疏的，A、B和C。
我需要取AB的矩阵乘积，这会产生一个稠密的矩阵。在此之后，我需要AB (element wise *) C. C的元素明智的乘积是稀疏的，因此元素向乘法将使大部分稠密积AB为零，从而再次形成稀疏矩阵。
知道了这一点，我正试图找出一种策略，以避免实现所有密集的AB组件。
如果C_{i，j}为0，则不应实现AB_{i，J}。这意味着我可以跳过A_{row i}和B_{col }的点积。但是，在A行上编写一个for循环来选择我想要实现的行似乎效率很低。
还有另外一种方法来聪明地进行这种乘法吗？
这里有一个R中的数据生成器示例，尽管我拥有的实际产品AB比这个
浏览 4提问于2020-02-26得票数 1
1回答
从3D投影到屏幕PHP
、、、、
我正试图在屏幕上获得三维物体的位置。
我已经研究过这个wiki atricle (和许多其他的)，但似乎我得到了错误的答案。
$center =array(0,0,0);
$point = array(0, 30, 30);
$rot = array(90,0,0);
$A=array(
array( 1, 0, 0),
array( 0, cos($rot[0]), sin($rot[0])),
array( 0, -sin($rot[0]),
浏览 0提问于2015-10-27得票数 0
回答已采纳
1回答
我有一个非常大的矩阵，需要一些计算。由于for循环在R中的速度是出了名的慢，我想用一些更智能的函数来代替它。
这是我用一个小的示例矩阵编写的for循环。
d <- matrix(c(1,1,0,0,1,1,1,0,0), 3,3)
for (i in 1:nrow(d)) {
for (j in 1:ncol(d)) {
if (d[i,j] == 1) {
d[j, i] =1
} else {d[j,i] = 0}
}
}
这段代码很好地替换了想要的值，产生了一个对称矩阵，其中di，j= dj，i。但是，当矩阵很大时，它会占用大量的时间和内存。做这
浏览 0提问于2016-10-14得票数 0
2回答
我为每0.1秒的时间步骤创建了距离矩阵，每隔60秒一次。对于每个时间步骤，矩阵看起来都是这样的，它们的距离值是填充它们的：
时间= 0.1
a1 b2 c3 d4
a1 0 5.4 9.1 10.1
b2 5.4 0 5.0 3.2
c3 9.1 5.0 0 6.6
d4 10.1 3.2 6.6 0
时间= 0.2
a1 b2 c3 d4
a1 0 2.4 9.1 12.1
b2 2.4 0 6.7 3.6
浏览 6提问于2022-05-09得票数 0
回答已采纳
1回答
我想在C中转置一个矩阵，使用尽可能少的内存。因此，不能简单地创建矩阵的副本，填充它，然后删除另一个矩阵。因此，我希望在矩阵本身上进行每一次交换，并使用一个temp变量进行重用。矩阵元素使用指针算法进行索引，但据我所知，这些都是正确的。现在的问题是，该函数不会导致任何更改。矩阵结构是动态分配的，具有如下结构；
typedef struct {
int rows;
int cols;
int** data;
} matrix;
我现在编写的转置函数如下所示；
void transpose(matrix* m){
int i,j,temp;
for(i=0;i<m->row;i+
浏览 2提问于2012-03-30得票数 1
回答已采纳
我在一个单元格数组中有三个序列：
Input_cell= {'ABCD','ACD', 'ABD'}
S1= 'ABCD' % which means A<B<C<D
S2= 'ACD' % which means A<C<D % missing B in the full string of 'ABCD'
S3= 'ABD' % which means A<B<D % missing C in the full string
浏览 1提问于2015-10-08得票数 0
回答已采纳
2回答
我有一个有20行的数据，其中每一行是一个带对角的10×10上三角协方差矩阵，我想把它转换成一个对称的10X10矩阵。我一直在跟踪这个帖子：，但是当我这样做的时候，它给了我一个列表，而不是矩阵。
另外，我希望对每一行进行递归处理，并将得到的20个矩阵存储在一个列表中。对于可重复的例子，这里是一个20行的数据格式。每一行为3x3上三角矩阵，具有对角线。
g <- data.frame(a=runif(20),b=runif(20),c=runif(20),d=runif(20),e=runif(20),f=runif(20))
g1<-c(g[1,])
g2= matrix(0, 3,
浏览 0提问于2017-11-02得票数 2
回答已采纳
我试图用chol()函数得到R中以下矩阵的下三角Cholesky分解。但是，它继续返回上三角分解，我似乎无法找到获得下三角分解的方法，即使在查看了文档之后也是如此。以下是我使用的代码-
x <- matrix(c(4,2,-2, 2,10,2, -2,2,5), ncol = 3, nrow = 3)
Q <- chol(x)
Q
# [,1] [,2] [,3]
# [1,] 2 1 -1.000000
# [2,] 0 3 1.000000
# [3,] 0 0 1.732051
我基本上需要找到矩阵Q，这样的
浏览 3提问于2016-11-09得票数 4
回答已采纳
2回答
我有个向量说
a = c(1,2,3,4,5,6)
我想把它们组织成一个上三角形矩阵(不考虑对角线元素，它们都是零)逐行。我的目标是得到以下矩阵：
[,1] [,2] [,3] [,4]
[1,] 0 1 2 3
[2,] 0 0 4 5
[3,] 0 0 0 6
[4,] 0 0 0 0
但是下面的方法是用这个向量替换对角线元素，但是按列赋值。例如,
b= matrix(0, 4, 4)
b[upper.tri(b, diag=FALSE)]=a
它将给我以下矩阵
[,1] [
浏览 4提问于2015-06-11得票数 12
回答已采纳
1回答
我有一个很大的块状数组..。
power = ...
print power.shape
>>> (3, 10, 10, 19, 75, 10, 10)
这是对称的w.r.t。10x10部分，即以下二维矩阵是对称的。
power[i, :, :, j, k, l, m]
power[i, j, k, l, m, :, :]
对于i，j，k，l，m的所有值
我能利用这个因子-4增益吗？例如，当将矩阵保存到文件中时(用savez_compressed保存50 mb )
我的尝试：
size = 10
row_idx, col_idx = np.tril_indices(size)
浏览 0提问于2014-05-15得票数 4
回答已采纳
我有一个对称矩阵存储在一个以制表符分隔的文件中，该文件有3列，其中前两列是行和列的位置，第三列是数据。因为它是对称矩阵，所以只有一个三角形被表示，所以它看起来像这样：
A A 0.2
A B 0.1
B B 1.2
A C 0.9
B C 2.3
C C 3.4
A D 2.1
B D 4.3
C D 0.8
D D 1.0
...
有没有什么方便的函数可以方便地将其加载到pandas DataFrame中？
浏览 30提问于2020-02-23得票数 1
回答已采纳
2回答
我有一个数据框架，它有三个变量: a，b，c。数据集中的前两列是所有可能组合的两个变量的配对，第三列是它们之间的相关性。如下所示。
> var1 <- c("a","a","b")
> var2 <- c("b","c","c")
> r <- c(.55,.25,.75)
> as.data.frame(cbind(var1,var2,r))
var1 var2 r
1 a b 0.55
2 a c 0.25
3
浏览 9提问于2017-01-09得票数 2
3回答
我正在寻找一个项，可以用于矩阵的对角线等于零和上下三角形是对称的。
例如，如果我的矩阵大小为5乘5，则矩阵将如下所示：
0 1 2 3 4
1 0 5 6 7
2 5 0 8 9
3 6 8 0 10
4 7 9 10 0
这个矩阵有术语吗？
浏览 5提问于2017-07-07得票数 1
回答已采纳
1回答
是否有一种方法仅用一个变量(第一个加载矩阵)在C++中确定给定矩阵的行列式，而在下一个递归函数中仅使用该矩阵的一个引用？
如何利用矩阵中元素的坐标来确定给定矩阵的子矩阵的行列式，而不把它们作为矩阵，只使用第一矩阵中的元素及其坐标？可以使用递归完成吗?还是不应该使用递归？
浏览 1提问于2018-05-30得票数 1
回答已采纳
在Opencv 3.0中，我使用以下函数获得了基本矩阵：
E = findEssentialMat(points2, points1, focal, pp, RANSAC, 0.999, 1.0, mask);
从掩码值来看，内部节点的数量确实很高。
在我使用recoverPose函数之后：
recoverPose(E, points2, points1, R, t, focal, pp, mask);
但我注意到，传递的inliers数量(查看掩码)非常少(通常为零)。
我的问题是:有没有一种方法可以从recoverPose函数中获得更多的内部函数？如果不是，是否有其他方法可以获得R和t(不
浏览 3提问于2021-01-18得票数 0
2回答
生成所有可能的无向图
、、、
使用关联矩阵生成所有可能的图的有效解决方案是什么？这些问题等价于生成所有可能的二元三角矩阵。我的第一个想法是在迭代工具中使用python。例如，生成所有可能的4x4矩阵。
for b in itertools.combinations_with_replacement((0,1), n-3):
b_1=[i for i in b]
for c in itertools.combinations_with_replacement((0,1), n-2):
c_1=[i for i in c]
for d in itertools.combinati
浏览 0提问于2015-01-13得票数 1
我试图得到从一个三角形到另一个三角形的变换矩阵。我主要关心的是2D。稍后我们将看到3D (但对解决方案开放)。
我在读。
现在请纠正我:如果我的第一个等边三角形的边长为1，其最左边的点位于原点，则它的变换应该是恒等矩阵。
所以，读上面的解，我正在寻找的变换矩阵，应该是三角形B矩阵*inv(恒等式) == B矩阵。
我无法相信这一点，因为这似乎是错误的。在给定的图像中，我想把蓝色转换成红色，蓝色在原点有一个点。选择的红色三角形的变换点总是最接近原点的。最后一个丢失的图像是最后一个将蓝色合并为红色的变换(相同的三角形)。
我已经在使用特征(c++)进行转换和计算了。
问题：
，，在这个过程中，我遗
浏览 7提问于2022-06-20得票数 0
我有一个矩阵$A$和矩阵$B_i=L_i^\mathsf{T} M_i L_i$。我想估计一下$i\in 1，n$，
\begin{eqnarray}
m_i^\mathsf{T}(A-B_i)^{-1}m_i
\end{eqnarray}
我使用lapack/blas来获得$B_i$s，但是当我想要计算线性系统$(A-B_i)x = m_i$时，使用lapack，例程dpotrf和dpotrs，它们要求矩阵是下三角形形式或上三角形形式。如何使用lapack/blas计算上述计算？有没有用C或C++写的例子？从一个完整的矩阵复制到下三角形形式的成本很高，有没有什么方法可以直接在下三角形形式中对
浏览 1提问于2017-10-28得票数 1
1回答
我有一个向量，其中元素的顺序很重要，例如：
v<-c(1,2,3,4,5,6,7,8,9,10)
我想将我的向量排列成一个具有特定顺序的下/上三角矩阵：
1 2 4 7
0 3 5 8
0 0 6 9
0 0 0 10
或
1 0 0 0
2 3 0 0
4 5 6 0
7 8 9 10
我想我能做得到，但我不知道该怎么做，空格可以填上NAs或0。
谢谢
浏览 0提问于2019-05-31得票数 1
回答已采纳
我想用cusolverDnDpotrfBatched找到3x3矩阵的Cholesky分解，但是我没有得到下三角矩阵中应该存在的零点。这是我想要计算cholesky分解的矩阵，它是这样的吗？我遗漏了什么？我知道这个帖子 .Here是我的代码：
#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <cuda_runtime.h>
#include <cusolverDn.h>
#include <iostream>
void printMat
浏览 19提问于2021-08-06得票数 0
回答已采纳
1回答
我将尽量不用我的解释使事情变得太复杂，但我对如何最好地用从另一个包中获得的现有相关值填充没有重复值的三角相关矩阵感到困惑。这涉及到从文本文件列表中提取特定值。这就是我到目前为止所做的：
# read in list of file names (they are named '1_1', '1_2' .. so on until '47_48' with no repeat values generated)
filenames <- read_table('/home/filenames.txt', col_names =
浏览 0提问于2020-02-25得票数 0
1回答
好的，我编写了一些代码来矢量化一个对称矩阵，它只取唯一的元素并将它们转化为一维向量，同时还用root2将非对角元素相乘：
def vectorize_mat(mat):
assert mat.shape[0] == mat.shape[1], 'Matrix is not square'
n = int(mat.shape[0])
vec_len = 0.5*n*(n+1)
weight_mat = (np.tri(n,k=-1)*np.sqrt(2))+np.identity(n)
mask_mat = np.tri(n).astyp
浏览 0提问于2018-11-15得票数 1
回答已采纳
1回答
我在下三角矩阵的高斯消元法上遇到了麻烦，我不能想象这里的循环是如何工作的。我试着向后运行循环，但没有帮助。现在我得到的只是上三角矩阵的高斯消元法。
For[k = 1, k <= size - 1, k++,
For[i = k + 1, i <= size, i++,
If[tab[[k]][[k]] != 0,
help = tab[[i]][[k]]/tab[[k]][[k]];
For[j = k, j <= size, j++,
tab[[i]][[j]] = tab[[i]][[j]] - help
浏览 0提问于2016-12-19得票数 0
我有一个总共包含55个数值的列表。我想创建一个10x10矩阵，其中只填充下(或上)三角形矩阵(带有对角线本身)。我知道我可以使用lower.tri()创建一个下三角矩阵，但是，当我使用这个函数时，似乎数据不是按行填充的。如果我使用matrix(v, nrow= 10, ncol= 10, byrow= TRUE)，那么将填充整个矩阵，而不仅仅是下对角线。我见过类似问题的解决方案()，但在这个例子中，他们只使用了6个变量，而我有10个变量，而且这个解决方案对我来说是扭曲的。
v <- 1:55
m <- diag(10)
浏览 2提问于2019-02-07得票数 0
2回答
我有很多需要在Matlab中求逆的大矩阵(大约5000 x 5000)。我实际上需要反转，所以我不能使用mldivide，这对于求解一个b的Ax=b来说要快得多。
我的矩阵来自一个问题，这意味着它们有一些很好的性质。首先，它们的行列式是1，所以它们绝对是可逆的。但是，它们是不可对角化的，或者我会尝试将它们对角化，颠倒它们，然后再将它们放回原处。它们的条目都是实数(实际上是有理的)。
我使用Matlab来获得这些矩阵，对于这些东西，我需要对它们的求逆进行处理，所以我更喜欢一种加速Matlab的方法。但是，如果我可以使用另一种更快的语言，请让我知道。我不知道很多其他的语言(只懂一点C和一点不懂Ja
浏览 2提问于2011-06-28得票数 13
回答已采纳
我有一个对称矩阵，我想把它转换成R中的上三角/下三角矩阵，有没有这样做的方法？
我不能使用upper.tri和lower.tri做到这一点。使用它们可以得到一个矩阵，其中的项可以是TRUE或FALSE。
浏览 2提问于2014-10-15得票数 28
回答已采纳
设置以重现我的最小工作示例 我有以下矩阵 K <- matrix(c(1.250000e+00, 3.366892e-07, 4.641930e-10, 1.455736e-08, 1.049863e-06,
3.366892e-07, 1.250000e+00, 5.482775e-01, 8.606555e-01, 9.776887e-01,
4.641930e-10, 5.482775e-01, 1.250000e+00, 8.603413e-01, 4.246732e-01,
1.455736
浏览 48提问于2019-12-24得票数 4
回答已采纳
1回答
我目前正在实现一组协方差矩阵的在线估计算法，由于数字错误，我的协方差矩阵碰巧不是始终对称的。
我目前通过将上三角矩阵镜像到下三角矩阵来解决这个问题：
for i=1:T
tempSigma(:, :, i) = eye(D)/sigmaT(:,:,i);
temp = triu(tempSigma(:,:,i), 1);
tempSigma(:, :, i) = triu(tempSigma(:, :, i));
tempSigma(:, :, i) = tempSigma(:, :, i) + temp';
en
浏览 1提问于2013-03-20得票数 2
回答已采纳
我有一个上对角线三角形矩阵和一个下对角线三角形矩阵，我想合并这两个矩阵，但也要添加一个1的对角线。 这是我的资料..。 upper = np.array([[0.90858992, 0.96590361, 0.95616282],
[ np.nan, 0.90055735, 0.8956839 ],
[ np.nan, np.nan, 0.94590614]]) 下面的..。 lower = np.array([[0.90858992, np.nan, np.nan ],
浏览 12提问于2021-07-08得票数 2
回答已采纳
3回答
我正在寻找一种方法来执行对称稀疏矩阵-矩阵乘法:X=A，其中稀疏矩阵A以前存储在CSR3格式(上三角)，而矩阵B是密集非对称矩阵。在MKL库中是否有一个例程可以这样做？还是它们都需要CSR格式的完全稀疏矩阵(从CSR格式获取句柄)而不是三角矩阵(我构建了三角矩阵，因为我需要在MKL Pardiso中使用它)？我知道mkl_sparse_d_mv(.)例程，但是我找不到从对称稀疏矩阵获得稀疏矩阵句柄的方法，这个对称稀疏矩阵以前是像CSR格式的上三角矩阵一样存储的。谢谢你，丹妮尔
浏览 7提问于2020-06-24得票数 1
1回答
我有(n×n)个下三角矩阵，那么它的复杂度是多少？它说这将是O(n^2)。那么如果我使用一维数组来表示这个下三角矩阵，那么这次的复杂度是O(N)还是其他的呢？比较每一轮的两种复杂性。对于在c++上运行的操作。我需要每一轮的数组大小。
浏览 1提问于2021-11-14得票数 0
1回答
我有一个不对称的矩阵，基本上我想比较对角反对的元素如下：
如果对角反对的元素在符号上相等但相反，则保留元素的绝对值，而对角反对的值为零。
如果不是这样的话，那么这两个元素中的一个就是0(但我们不知道是哪一个)，所以取两者的绝对值。
完成后，转置矩阵的下三角形并将其添加到上三角形。
我提出了以下python循环：
for i in range(0, number_files):
for j in range(0, len(Identifier)):
for k in range(0,len(Identifier)):
if Matr
浏览 3提问于2014-11-28得票数 2
回答已采纳
1回答
在R中，我使用cov2cor()计算一个相关矩阵，如：
A,B,C,...
A 1,0.5,0.2,...
B 0.5,1,0.4,...
C 0.2,0.4,1,...
...
我如何重塑矩阵，以便将列按行堆叠，如：
X,Y,Correlation
A,B,0.5,
A,C,0.2,
...
B,C,0.4,
...
提醒A,A被排除在外，A,B B,A被当作重复处理，这样就排除了一个副本。
有什么简单的方法来实现这一点吗？
浏览 4提问于2014-02-05得票数 3
回答已采纳
我有一个16000*16000矩阵，想要找到最小的条目。这个矩阵是一个距离矩阵，所以它对角是对称的。为了每次精确地得到一个最小值，我将下三角形和对角设置为np.inf。下面是一个5*5矩阵的例子：
inf a0 a1 a2 a3
inf inf a4 a5 a6
inf inf inf a7 a8
inf inf inf inf a9
inf inf inf inf inf
我只想在上三角形中找到最小项的索引。然而，当我使用np.argmin()时，它仍然会遍历整个矩阵。有没有办法“忽略”下三角，提高速度？
我尝试过许多方法，例如：
使用蒙面数组使用triu_indices()提取
浏览 0提问于2021-03-19得票数 2
6回答
在内存中存储对称矩阵的最佳方式是什么？
这将是很好的节省一半的空间，而不会牺牲太多的速度和结构的复杂性。这是一个与语言无关的问题，但是如果您需要做一些假设，那么就假设它是一种很好的旧的普通编程语言，比如C或C++。
这似乎是一个有意义的东西，只要有一种方法可以让事情变得简单，或者只是在矩阵本身非常大的时候，我说的对吗？
为了正式起见，我的意思是这个断言对于我想要存储的数据总是正确的
matrix[x][y] == matrix[y][x]
浏览 1提问于2010-07-06得票数 22
回答已采纳
1回答
每当我尝试通过特征库使用任何矩阵或稀疏矩阵的selfadjointView属性时，我都会不断收到错误消息。下面是检查这一点的简单代码。在我的程序中，我尝试使用自伴矩阵：
#define EIGEN_YES_I_KNOW_SPARSE_MODULE_IS_NOT_STABLE_YET
#include <Eigen/Sparse>
#include <Eigen/Dense>
#include <Eigen/Core>
#include <iostream>
using namespace Eigen;
int main ()
{
SparseM
浏览 4提问于2011-08-25得票数 3
回答已采纳
1回答
我已经成功地计算了旋转，平移与本征相机矩阵的两个相机。我还得到了纠正图像从左和右相机。现在，我想知道如何计算一个点的三维坐标，一个图像中的一个点。在这里，请看绿点。我看了一下方程，但它需要基线，我不知道如何计算。你能告诉我用给定的信息(R，T和内禀矩阵)计算绿点的三维坐标的过程吗？
我也有一个基本矩阵和基本矩阵，以防我们需要它们。2.原始图像大小为960 x 720。校正图为925 x 669 3，左边的绿点为(562,185)，右边的为(542,185)。
浏览 2提问于2018-09-15得票数 3
回答已采纳
我试图使用相同名称的R包中的“plspm”函数来评估SEM。我有两种相互竞争的模式：
模式1：
A<-c(0,0,0)
B<-c(1,0,0)
C<-c(0,1,0)
sat_path=rbind(A,B,C)
innerplot(sat_path)
对比
模式2：
A<-c(0,1,0)
B<-c(0,0,0)
C<-c(0,1,0)
sat_path=rbind(A,B,C)
innerplot(sat_path)
因此，区别仅在于A和B之间箭头的方向。
当我运行代码(参见下面的代码)时，第一个模型(模型1)返回一个答案，而第二个模型(模型2)返回以下错
浏览 2提问于2017-05-04得票数 1
4回答
我有一个由相等的行和列组成的数据基/矩阵。我只想提取上或下三角形。
x<-data.frame(matrix(1:25,nrow=5))
colnames(x)<-LETTERS[1:5]
rownames(x)<-LETTERS[1:5]
x[upper.tri(x,diag=F)]
从这个结果中，不可能说明值来自哪个列和行的组合。因此，我希望在结果中包含行和列属性。就像这样：
Col Row Val
B A 6
C A 11
C B 12
...
我需要对一个大的相关矩阵这样做。谢谢。
浏览 3提问于2013-08-08得票数 6
回答已采纳
我有一个类似于这个元素的向量：
Sentences -> c("Sentence_1", "Sentence_2", "Sentence_3", "Sentence_4")
我还有一个矩阵，其中行和列由Sentences中的每个元素组成。
Sentence_1 Sentence_2 Sentence_3 Sentence_4
Sentence_1 1 0.25 0.74 0.35
Sentence_2 0.25
浏览 2提问于2020-09-05得票数 0
回答已采纳
2回答
我有一个对称矩阵mat
A B C
A 1 . .
B . 1 .
C . . 1
我想计算出它的两个最高元素。现在，由于它是对称矩阵，所以我想使用upper.tri，如下所示：
mat.tri<-upper.tri(mat) # convert to upper tri
mat.ord<-order(mat.tri,na.last=TRUE,decreasing=TRUE)[1:2] # order by largest
a.ind<-which(mat%in%mat.tri[mat.ord]) # get absolute indices
r.
浏览 5提问于2014-01-23得票数 1
回答已采纳
所以基本上，这是我需要创建的矩阵 sig = matrix(c(0.56^2, 0.56*7.77*-0.07, 0.56*13.48*-0.095, 0.56*16.64*-0.095,
0.56*7.77*-0.07, 7.77^2, 7.77*13.48*0.959, 7.77*16.64*0.936,
0.56*13.48*-0.095, 7.77*13.48*0.959, 13.48^2, 13.48*16.64*0.997,
0.56*16.64*-0.095, 7.77*16.64*0.9
浏览 25提问于2021-09-19得票数 0
回答已采纳
1回答
我将一个距离矩阵转换成一个两列的data.frame，我想过滤掉“对角线”和">0“值，保持相同的对象。
我有一个距离矩阵
mat<-c(0,0.5,0,0.5,0,0.3,0,0.3,0)
dim(mat)<-c(3,3)
rownames(mat)<-c("A","B","C")
colnames(mat)<-c("A","B","C")
我试过了
df = as.data.table(mat)
df1 <- stack(df)
setDT(df1)
浏览 0提问于2019-01-29得票数 0
回答已采纳
我仍然在学习使用data.table (从包)，即使在网上寻找帮助和帮助文件之后，我仍然在努力做我想做的事情。我有一个包含60多个列的大数据表(前三列对应于因子，其余的对应于响应变量，在本例中是不同的物种)和几行对应于不同处理级别和物种丰度的行。一个非常小的版本如下所示：
> TEST<-data.table(Time=c("0","0","0","7","7","7","12"),
Zone=c("1","1",
浏览 1提问于2014-02-19得票数 0
回答已采纳
1回答
我想测试一个猜想，即某个算法的性能取决于c违反三角形不等式的程度，其中c是被算法求解的实例的距离矩阵。
我的问题是，修正距离矩阵c的好方法是什么，这样我就可以控制它违反三角形不等式的程度？
到目前为止，我得出的最好结果是这样的:对于参数p和q，随机选择c元素的p%，对于每个选定的元素，乘以1+u，其中u是从0，q中均匀选择的。
我想我们可以做得更好。有什么想法吗？
浏览 1提问于2016-11-28得票数 0
回答已采纳
1回答
实时求逆矩阵
、、、
我一直在开发一个实时工作的C语言控制软件。软件实现了被控系统的离散状态空间观测器等。对于观察者的实现，需要计算4x4维的矩阵的逆。逆矩阵计算必须每隔50微秒进行一次，值得一提的是，在这段时间内还将进行其他非常耗时的计算。因此，逆矩阵计算所消耗的时间必须远远少于50微秒。还有必要说明的是，所使用的DSP没有支持浮点运算的ALU。 我一直在寻找一些有效的方法来做到这一点。我的一个想法是准备计算矩阵4x4的行列式的通用公式和计算4x4矩阵的伴随矩阵的通用公式，然后根据下面给定的公式计算矩阵的逆矩阵。 ? 你对这种方法有什么看法？
浏览 46提问于2020-09-10得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
