平均上、下三角矩阵
是一种特殊的方阵，其中上（下）三角矩阵的元素在主对角线上方（下方）全为0，而主对角线及其上（下）方的元素可以是任意值。这种矩阵结构在数学和计算机科学中有广泛的应用。
优势：
- 空间效率高：由于上（下）三角矩阵中大部分元素为0，因此存储空间较小，节省了存储资源。
- 计算效率高：对于上（下）三角矩阵，许多运算可以通过简单的索引计算得出，从而减少了计算量。
- 适用于特定问题：上（下）三角矩阵常用于解线性方程组、矩阵分解、数值计算等领域，特别适用于具有稀疏性质的问题。
应用场景：
- 线性方程组求解：上（下）三角矩阵可以简化线性方程组的求解过程，提高计算效率。
- 矩阵分解：上（下）三角矩阵可以用于矩阵LU分解、Cholesky分解等数值计算中。
- 图论算法：在图论算法中，上（下）三角矩阵可以用于表示图的邻接矩阵，从而进行图的遍历、最短路径等操作。
推荐的腾讯云相关产品： 腾讯云提供了多种云计算相关产品，以下是一些推荐的产品：
- 云服务器（ECS）：提供弹性的虚拟云服务器，可根据需求灵活调整配置。
- 云数据库MySQL版（CDB）：提供高性能、可扩展的MySQL数据库服务，支持自动备份和容灾。
- 云存储（COS）：提供安全可靠的对象存储服务，适用于存储和处理各种类型的数据。
- 人工智能平台（AI Lab）：提供丰富的人工智能算法和模型，支持图像识别、语音识别等应用。
- 物联网套件（IoT Hub）：提供全面的物联网解决方案，支持设备接入、数据管理和应用开发。
以上是腾讯云相关产品的简介，更详细的信息可以参考腾讯云官方网站：https://cloud.tencent.com/
相关·内容
2回答
我有来自三个市场的每日回报(GLD，间谍和USO)。我的目标是在130天的滚动基础上，从相关矩阵中计算出平均成对相关性。
我的出发点是：
import numpy as np
import pandas as pd
import os as os
import pandas.io.data as web
import datetime as datetime
from pandas.io.data import DataReader
stocks = ['spy', 'gld', 'uso']
start = datetime.datetime
浏览 2提问于2016-05-25得票数 6
回答已采纳
2回答
无法填充我的loop循环
、、、
因此，我正在研究R的矩阵，如下所示：
diff_0
SubPop0-1, SubPop1-1, SubPop2-1, SubPop3-1, SubPop4-1,
SubPop0-1, NA NA NA NA NA
SubPop1-1, 0.003403100 NA NA NA NA
SubPop2-1, 0.005481177 -0.002070277
浏览 4提问于2017-11-24得票数 0
回答已采纳
1回答
假设我有一个由均值定义的二维高斯矩阵的混合物，协方差矩阵。混合中每一个高斯的重量是相同的。现在我要计算概率密度函数的值。所以很明显，我会计算混合物中每一个高斯人的数值，并取平均值。但是，我的混合物含有许多高斯(> 10^4)，所以我担心这可能在数值上是不稳定的，而且可能也很慢，当我必须计算每个PDF中的指数时。解决这个问题的好办法是什么？
浏览 0提问于2016-10-25得票数 2
回答已采纳
2回答
内存有效的平均成对距离
、、、
我知道scipy.spatial.distance.pdist函数，以及如何从结果矩阵/ndarray计算平均值。 >>> x = np.random.rand(10000, 2)
>>> y = pdist(x, metric='euclidean')
>>> y.mean()
0.5214255824176626 在上面的示例中，y变得非常大(几乎是输入数组的2,500倍)： >>> y.shape
(49995000,)
>>> from sys import getsizeof
浏览 30提问于2019-05-03得票数 1
回答已采纳
2回答
我需要计算n×n矩阵中每个非对角元素的平均值。上下三角形是多余的。下面是我目前使用的代码：
A <- replicate(500, rnorm(500))
sapply(1:(nrow(A)-1), function(x) mean(A[row(A) == (col(A) - x)]))
这似乎是可行的，但在更大的矩阵中却不能很好地扩展。我拥有的不是很大，大约是2-5000^2，但即使是1000^2，花费的时间也比我想的要长：
A <- replicate(1000, rnorm(1000))
system.time(sapply(1:(nrow(A)-1), function(
浏览 5提问于2012-12-17得票数 10
回答已采纳
访问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回答
我在下三角矩阵的高斯消元法上遇到了麻烦，我不能想象这里的循环是如何工作的。我试着向后运行循环，但没有帮助。现在我得到的只是上三角矩阵的高斯消元法。
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
1回答
我需要从下列文件中为图的边指定权重：
L.肖和S. Boyd“图Laplacian特征值的凸优化”的“分布平均的快速线性迭代”
我有图的邻接矩阵(50x50矩阵)，有512个非零值.
我也有一个256乘1向量的最优权重。
对于我所使用的软件，我需要一个50×50矩阵，其中边(i，j)的权重在邻接矩阵的相关位置(与边(j，i)相反的符号)。
我的尝试是在下面，但我无法使它发挥作用。
function weights = construct_weight_mtx(weight_list, Adj)
weights = zeros(size(Adj));
positions = find(Adj)
浏览 4提问于2016-04-04得票数 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
回答已采纳
我有一个上对角线三角形矩阵和一个下对角线三角形矩阵，我想合并这两个矩阵，但也要添加一个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
回答已采纳
2回答
我想生成一个随机的实对称方阵，其条目均匀分布在0和1之间。我的尝试是：a = rand(5); b = a + a.'
我担心的是，虽然根据文档，矩阵a是均匀分布的，但是矩阵b可能不是，因为两个随机数的平均值可能与原始数不同。
我尝试使用hist(a); hist(b)，但不确定如何解释生成的图形。编辑:根据Oli矩阵b不再均匀分布，有没有办法做到这一点？
浏览 2提问于2012-03-17得票数 7
回答已采纳
我想用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
回答已采纳
我有一个总共包含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
1回答
相关矩阵的滚动均值
、、、
我必须计算相关矩阵的滚动平均值。是我干的
df_corr = df.rolling(window=3).corr()
但现在我想要
2017-06-05 00:00:00+02:00 NaN
2017-06-06 00:00:00+02:00 NaN
2017-06-07 00:00:00+02:00 0.86
...
其中我计算右上矩阵的平均值(不包括对角线)。
浏览 2提问于2017-07-04得票数 2
回答已采纳
3回答
我正在寻找一种方法来执行对称稀疏矩阵-矩阵乘法:X=A，其中稀疏矩阵A以前存储在CSR3格式(上三角)，而矩阵B是密集非对称矩阵。在MKL库中是否有一个例程可以这样做？还是它们都需要CSR格式的完全稀疏矩阵(从CSR格式获取句柄)而不是三角矩阵(我构建了三角矩阵，因为我需要在MKL Pardiso中使用它)？我知道mkl_sparse_d_mv(.)例程，但是我找不到从对称稀疏矩阵获得稀疏矩阵句柄的方法，这个对称稀疏矩阵以前是像CSR格式的上三角矩阵一样存储的。谢谢你，丹妮尔
浏览 7提问于2020-06-24得票数 1
我有一个对称矩阵，我想把它转换成R中的上三角/下三角矩阵，有没有这样做的方法？
我不能使用upper.tri和lower.tri做到这一点。使用它们可以得到一个矩阵，其中的项可以是TRUE或FALSE。
浏览 2提问于2014-10-15得票数 28
回答已采纳
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回答
我有一个非常大的矩阵，需要一些计算。由于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
1回答
将不对称方阵变换为对称方阵
、、、、
我有一个矩阵格式的数据存储在一个CSV文件中，看起来是这样的，
我想让这个6*6矩阵成为一个对称矩阵，就像这样，
如何使用python (或matlab)将n乘n矩阵(方阵)转换为对称矩阵？或者有其他工具可以这样做吗？
请给我任何建议，谢谢！
浏览 0提问于2014-05-05得票数 0
回答已采纳
我试图用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
回答已采纳
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
回答已采纳
2回答
我有一个矩阵A (Mx2)和一个矩阵B (Nx2)，使用M>=N。我必须计算垂直欧几里德距离来测试A和B之间的相似性。关于图像，单个欧几里德距离是每个pi和p'i之间的线段，一旦计算了每个垂直距离，我就可以计算平均值或和来获得单个参数。你知道做这件事的matlab例程或者这样的脚本吗？
浏览 0提问于2012-04-29得票数 2
回答已采纳
1回答
我很难理解这段代码中的错误是什么：
我有一堆矩阵，我想取每个矩阵的上三角部分，把它放到一个向量中，用它做一些事情，然后把结果映射回来。下面是代码：
%%
n=10;
m=3;
% generate a random 'stack of matrices'
bar=randn(n,n,m);
% index the upper triangular part
inds=triu(true(n,n));
% linearize
bar_lin=permute(bar,[3 1 2]);
bar_lin=bar_lin(:,inds);
% de-linearize
foo=
浏览 0提问于2017-06-30得票数 0
回答已采纳
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
回答已采纳
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
回答已采纳
2回答
问题很简单，但是在Pytorch中有没有一种从现有矩阵创建上三角矩阵的原生方法？我正在考虑使用一个掩模，但即使这样也需要创建上三角矩阵。
浏览 133提问于2019-11-21得票数 0
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
回答已采纳
尝试实现前向代换法，即以L为下三角矩阵，x，b为向量的解Lx =b问题的求解过程。
这是一项简单的任务：
def tri_solve(L,b):
n = len(b)
x = np.zeros(n)
x[0] = b[0]/L[0,0];
for i in range(1,n):
comp = 0;
for k in range(0,i):
index = L[i,k]
preSolution = x[k]
comp = comp + index * pre
浏览 0提问于2022-03-02得票数 0
我有一个名为" variables“的矩阵，它包含9个变量(9列)。我已经用下面的代码获得了成对相关矩阵：
matrix.cor <- cor(variables, method="kendall", use="pairwise")
现在，我想要获得作为所考虑变量数量的函数的平均成对相关性。我的意思是，2个变量，3个变量，4个变量的所有可能相关性的平均值...最多9个变量，以查看添加变量的效果。我有这个R代码(摘自一篇包含更多因子和列的文章)，但它运行得不好，我只获得了9个变量的平均值。
pairwisecor.df = ddply(data,c(
浏览 9提问于2017-11-01得票数 2
1回答
4种观点中6点的三角剖分
、、、
我对一个物体有四个不同的视图，每一个都有相应的K和M摄像机矩阵(内参数和外部参数)。
我有6个3D点，我在每个2D视图上标记(所以242d点)。
我想要做的是-把那些2D点投影到3D中，然后在空间中得到我的6个3D点。
当你有两个视图和它们的矩阵时，大多数的论文都展示了如何做到这一点。不过，我有4种看法。我在Python中使用OpenCV。
作为第一种方法，也许我可以使用两个视图的所有组合来进行三角剖分，所以我得到了4选择2=6组3D点，然后对所有的三角点进行平均。
是否有更好的办法更可取？你认为我应该用OpenCV以外的东西吗？(但仍然必须是Python )
浏览 2提问于2018-03-05得票数 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
回答已采纳
2回答
非零值的Numpy均值
、、、
我有一个大小为N*M的矩阵，我想找出每一行的平均值。这些值是从1到5，没有任何值的条目被设置为0。然而，当我想用下面的方法计算平均值时，它给出了错误的平均值，因为它还计算了值为0的条目。
matrix_row_mean= matrix.mean(axis=1)
我如何才能得到仅非零值的平均值？
浏览 0提问于2016-07-23得票数 23
回答已采纳
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回答
通过对一些示例数据使用pairwise.t.test，我可以得到P-值的矩阵：
attach(airquality)
pairwise.t.test(Ozone, Month)
这意味着：
Pairwise comparisons using t tests with pooled SD
data: Ozone and Month
May Jun Jul Aug
Jun 1.00000 - - -
Jul 0.00026 0.05113 - -
Aug 0.00019 0.04987
浏览 2提问于2015-03-02得票数 1
回答已采纳
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回答
我在TensorFlow中有一个N维的一维向量，
如何构造成对平方差的和？
示例
输入向量
[1,2,3]
输出 6
计算为
(1-2)^2+(1-3)^2+(2-3)^2.
如果输入为N-dim向量l，则输出应为sigma_{i，j}((L)^2).
添加了问题：如果我有一个2d矩阵，并且希望对矩阵的每一行执行相同的处理，然后对所有行的结果进行平均，我如何才能做到呢？非常感谢！
浏览 1提问于2018-06-25得票数 3
回答已采纳
1回答
径向功率谱分析
、、、
首先，我需要绘制这个径向平均光谱1-D与波数|k|的关系图。但是早期我使用两个矩阵来计算|k|，所以我有一个表示|k|的矩阵。所以我也要对这个矩阵进行径向平均？
其次-如果我在使用FFT2函数之前从我的数据中删除了一个平均值，我希望看到的只是值的振幅的差异，但实际上我注意到径向平均谱形状的差异。
我正试着用它来分析磁场数据。谢谢你的帮助！！
浏览 4提问于2013-08-04得票数 0
2回答
你知道在numpy中有像hermitian矩阵类这样的东西吗？我想优化矩阵计算，比如
B=U*A* U.H
，其中A(因此，B)是厄米特的。在不指定的情况下，将计算B的所有矩阵元素。事实上，它应该能够在这里节省大约2倍的系数。我错过了什么吗？
我需要的方法应该取A的上/下三角形，U的完整矩阵，并返回B的上/下三角形。
浏览 1提问于2015-11-06得票数 2
2回答
我需要确定一个矩阵，该矩阵表示在3D中对齐两个已知全等三角形所需的变换。我之前尝试过询问有关使用的问题，但我遇到了很多问题，所以我现在正在考虑尝试不同的想法。
谁能告诉我如何创建一个矩阵来表示将一个三角形(DEF)移动到与另一个三角形(ABC)相同的位置和方向所需的平移和旋转，其中A位于原点？
浏览 0提问于2011-10-11得票数 2
回答已采纳
1回答
是否有一种方法仅用一个变量(第一个加载矩阵)在C++中确定给定矩阵的行列式，而在下一个递归函数中仅使用该矩阵的一个引用？
如何利用矩阵中元素的坐标来确定给定矩阵的子矩阵的行列式，而不把它们作为矩阵，只使用第一矩阵中的元素及其坐标？可以使用递归完成吗?还是不应该使用递归？
浏览 1提问于2018-05-30得票数 1
回答已采纳
我有一个现有的三角剖分(顶点x和y以及连通矩阵tri)，我想在现有的三角剖分上应用delaunayTriangulation类的PointLocation方法(就像旧版本的MATLAB中过时的tsearch函数)。
但是，PointLocation方法显然需要一个delaunayTriangulation实例作为输入。delaunayTriangulation类似乎总是执行自己的三角剖分过程，在给定顶点x和y的情况下，这将导致与现有的tri矩阵不同的连通矩阵。
是否有方法将PointLocation (或类似于tsearch)应用于我现有的三角测量？我有Matlab 2013 a。
浏览 19提问于2015-01-06得票数 1
回答已采纳
设一个顶点，它是三角形的一部分，也是一个四边形的一部分。
据我所知，这个顶点的法线是四边形和三角形的法线的平均值。
三角形是在四边形之前画的。我什么时候应该调用glNormal，用什么向量？
我应该调用glNormal 2次，每次使用相同的向量(平均法向量)吗？
上一次用平均法向量绘制顶点时，我应该调用glNormal吗？
浏览 3提问于2014-04-26得票数 1
回答已采纳
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回答
我有一个向量，其中元素的顺序很重要，例如：
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
回答已采纳
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
回答已采纳
1回答
我正在尝试计算矩阵中每对行的总和。假设我有一个m x n矩阵，比如 [[1,2,3],
[4,5,6],
[7,8,9]] 我想创建一个所有行对的求和矩阵。因此，对于上面的矩阵，我们希望 [[5,7,9],
[8,10,12],
[11,13,15]] 一般来说，我认为新的矩阵将是(m general 2) x n。 import torch
x = torch.tensor([[1,2,3], [4,5,6], [7,8,9]])
y = x[None] + x[:, None]
torch.cat((y[0, 1:3, :], y[1, 2:3, :])) 它手动创建了我正
浏览 14提问于2021-07-26得票数 0
拉格朗日多项式需要这个。我很好奇，如果不使用for循环，如何才能做到这一点。目前的代码如下：
tj = 1:n;
ti = zeros(n,n-1);
for i = 1:n
ti(i,:) = tj([1:i-1, i+1:end]);
end
我的tj不仅仅是一个1:n向量，但这并不重要。当这个for循环完成工作时，我宁愿使用一些矩阵操作。我试图寻找一些适当的矩阵与它相乘，但到目前为止没有运气。
浏览 8提问于2016-02-05得票数 2
回答已采纳
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
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
