link管理
链接快照平台
- 输入网页链接，自动生成快照
- 标签化管理网页链接
相关文章推荐
|
睡不着的跑步鞋 · MTK核心板|MTK开发板|MTK模块|4G ...· 2 月前 · |
|
|
近视的仙人掌 · 毕马威：银行业客户体验白皮书-新意义、新内涵 ...· 3 月前 · |
|
|
另类的砖头 · 7款强大的Javascript网格插件推荐 ...· 4 月前 · |
|
|
逼格高的酱牛肉 · Plugin to set scale ...· 4 月前 · |
|
|
豪气的感冒药 · 2024年最新款【5X6.5尺四季涼被】-銀 ...· 4 月前 · |
link管理
›
在矩阵上循环到R中的矩阵_在R中的for循环中调用矩阵元素_R中对角矩阵的逆矩阵 -
| for循环 矩阵乘法 对角矩阵 矩阵 |
| https://cloud.tencent.com.cn/developer/information/%E5%9C%A8%E7%9F%A9%E9%98%B5%E4%B8%8A%E5%BE%AA%E7%8E%AF%E5%88%B0R%E4%B8%AD%E7%9A%84%E7%9F%A9%E9%98%B5-ask |
|
飘逸的小摩托
7 月前 |
精选内容/技术社群/优惠产品,
尽在小程序
立即前往
在矩阵上循环到R中的矩阵
是指在R语言中对矩阵进行循环操作。R是一种流行的编程语言，特别适用于数据分析和统计建模。下面是关于在矩阵上循环到R中的矩阵的完善且全面的答案：
概念： 在R中，矩阵是一个二维的数据结构，由行和列组成，可以存储相同类型的数据。在矩阵上循环是指对矩阵中的每个元素进行遍历和操作。
分类： 在矩阵上循环可以分为两种方式：按行循环和按列循环。按行循环是逐行遍历矩阵，按列循环是逐列遍历矩阵。
优势： 在矩阵上循环可以方便地对矩阵中的每个元素进行操作和计算，提高数据处理的效率和准确性。同时，R语言提供了丰富的矩阵操作函数和向量化运算，可以简化矩阵上的循环操作。
应用场景： 在矩阵上循环在数据分析和统计建模中广泛应用。例如，可以使用循环来计算矩阵的和、平均值、方差等统计指标，或者对矩阵中的每个元素进行条件判断和转换操作。
推荐的腾讯云相关产品和产品介绍链接地址： 腾讯云提供了云服务器、云数据库、云存储等一系列云计算产品，可以满足各种应用场景的需求。以下是一些相关产品和介绍链接地址：
- 云服务器（ECS）：提供弹性计算能力，支持多种操作系统和应用场景。详情请参考：https://cloud.tencent.com/product/cvm
- 云数据库MySQL版（CDB）：提供高性能、可扩展的关系型数据库服务。详情请参考：https://cloud.tencent.com/product/cdb_mysql
- 云存储对象存储（COS）：提供安全可靠的对象存储服务，适用于图片、视频、文档等大规模数据存储。详情请参考：https://cloud.tencent.com/product/cos
总结： 在R中，通过循环可以对矩阵进行遍历和操作，提高数据处理的效率和准确性。腾讯云提供了一系列云计算产品，可以满足各种应用场景的需求。
页面内容是否对你有帮助？
有帮助
没帮助
相关· 内容
1
回答
问题如下：
我有一套20个方程式：
r1 = s1*h1_1 + s2*h1_2 + ... s20*h1_20
r2 = s1*h2_1 + ...
r20 = s1*h20_1 + ...
其中r，s和h是矩阵，*表示点态乘积。它可以用矩阵形式R= H*S重写，我想为S解这个方程，所以我需要计算inv(H)*R，但是如果H的每个元素都是矩阵，我如何计算inv(H)？我不能简单地将这些较小的矩阵串联成一个更大的矩阵H，然后再将其反转--因为这将给出一个与之不同的结果，例如用符号值反演矩阵H，然后用较小的矩阵代替这些符号值(因为方程组中存在点态积)。
到目前为止，我想出了一个解决
浏览
0
提问于2015-07-23
得票数
4
回答已采纳
2
回答
如何从C中的列主矩阵中提取向量？
、
、
来自MATLAB背景，我经常使用花哨的矩阵操作命令(如vec = matrix(:,1) )来提取matrix的第一列作为向量。
将一些代码移植到C上，需要与FORTRAN和MATLAB接口，现在我可以用列的主要顺序将矩阵存储在一维数组中。
所以基本上，我使用的是宏
#define SUB2IND_2D(s1, s2, i1, i2) (s1)*(i2) + (i1)
还有这些循环
for(size_t r=0; r<ROWS; ++r)
for(size_t c=0; c<COLS; ++c)
size_t index = SUB2IND_
浏览
4
提问于2014-04-01
得票数
3
回答已采纳
1
回答
目前，我正忙于一种新的因素分析。我创建了以下二进制运算符：
%M%从R数据中计算正定对称矩阵
%E%计算公共度估计，并将这些估计放在方阵M的主对角线上。
%F%一次从新的平方矩阵中提取一个因子。
%R%从%F%的因子负载中复制一个新的平方矩阵，并从原始矩阵中减去它。
这是一个周期。步骤2到4对残差重复，%F%提取第二个因子，%R%再现新的矩阵等。在这种方式下，完整的过程(步骤2-4，2-4)依次应用于越来越小的矩阵。
我的问题是:我如何在R中编程这些迭代？显然，在每个迭代中，步骤4是下一个迭代中从步骤2到步骤4的新周期的输入。
我用的一些R码：
res= a%
浏览
3
提问于2014-10-31
得票数
0
2
回答
我想为两个矩阵/数据帧之间的所有行组合计算距离度量。
结果将是具有单元i，j的矩阵，其对应于由应用于第一矩阵的行i和第二矩阵的行j的函数所给出的结果。下面是一个示例，用一个示例函数来说明我想要对for循环做什么。
x<-matrix(rnorm(30),10,3) ## Example data
y<-matrix(rnorm(12),4,3)
results<-matrix(NA,nrow(x),nrow(y))
for (i in 1:nrow(x)){
for (j in 1:nrow(y)){
r1<-x[i,]
r2<-y[j,
浏览
0
提问于2012-05-26
得票数
5
1
回答
Numpy乘法矩阵(两个以上矩阵)
、
、
、
如果我需要在5个矩阵之间进行乘法(考虑到矩阵的良好形状才能进行乘法)，np.dot不起作用。对于ex，在matlab中你可以很容易地写出：
Rf = T^-1 * S^-1 * R * S * T
浏览
2
提问于2020-02-11
得票数
0
1
回答
正在尝试编写一个按行对矩阵排序的函数。我可以写一些东西来循环遍历一个值向量上的值，但不能增加复杂性来使其循环到某个矩阵上。
sww = function(x){
n <- length(x)
for(i in 1:(n-1)){
for (j in (i+1):n) {
if(x[i] > x[j]){
tmp = x[i]; x[i] = x[j]; x[j] = tmp
return(x)
有人知道如何让它在整个矩阵上循环吗？
编
浏览
6
提问于2020-02-10
得票数
0
回答已采纳
1
回答
用Matlab编码实现矩阵向量
、
、
、
、
嘿，我正在对一个信号使用FIR滤波器。我想让hf变量成为一个矩阵，包含变量HF中L的值。如下所示：
HF = hf1、hf2、hf3、hf4、hf5；
这样我就可以在子图中使用HF矩阵，我可以看到每个hf向量的图。
正如你所看到的，我得到了编写HF = [HF, hf];的技巧，它应该实现HF到矩阵hf的向量。但是我不能让它工作。
下面是我使用的代码：
w0 = 9*pi/80;
for L = [5,10,50,100,1000]
n = 1:L;
dum = (sin(w0.*n))./(pi.*n);
h = [fliplr(dum),w0/pi,dum];
%
浏览
0
提问于2014-04-11
得票数
0
1
回答
操作系统: Linux Centos 6.5 / Centos 7
机器：“特殊”2.5兆字节物理内存(是的，真的)
编程语言:r
Problem：在R中，实践(似乎是)在循环之前分配一个与最终结果一样大的矩阵。在这种情况下，最终的矩阵大小在启动时不确定。
因此，在R中创建一个扩展矩阵的最干净的方法是什么？
关键的一点是，我希望在扩展时避免内存副本的开销。如果这是唯一的方法，那是可以的，但不是首选的。
重要:我非常熟悉内核级编程、内核级内存管理、页面处理系统等。是R的新手，学习R. (很有趣，也很有用)在这个例子中，我特别尝试创建一个场景来演示一个矩阵的内存扩展，而不需要预先分配内存，或者返回
浏览
1
提问于2015-05-31
得票数
0
我有很多2-x-2矩阵S1，S2，...，SN，在每个矩阵上，我想执行一个左右矩阵乘法，就像R*S*R^T中一样，其中R也是一个2-x-2矩阵。显然，我可以只使用for循环来编写，但我预计对于MATLAB中的大N，它会非常慢。在不使用for循环的情况下，有没有一种简单有效的方法来实现这一点呢？提前感谢！
浏览
2
提问于2012-11-23
得票数
1
回答已采纳
1
回答
使用oct2py对多个(10MM+)大小的数据调用corrcoef.m，返回R，P矩阵，为ML算法生成训练集。昨天，我的工作没什么问题。今早从顶部运行脚本，返回一个相同的测试集，通过oct2py传递给Octave。
我将被遣返：
Oct2PyError:八度计算错误:错误: isnan:未定义单元格错误:调用: corrcoef在第152行，第5列CorrCoefScript在第1行，第7列
首先，集合中没有null/nan值。事实上，根本就没有零。在任何列中都没有一致性，以至于在相关系数计算中没有标准偏差被返回。这是一种精确的声音。
其次，当我通过图形用户界面将测试集加载到Octav
浏览
1
提问于2019-03-22
得票数
1
回答已采纳
a=18;
b=22;
for i=1:10
r1 = randi([18 22],1,1)
name= (b-a).*rand(r1,2) + a
现在保存此所有名称值矩阵在结果矩阵中，每个循环行大小中的所有属性值的不是固定的，而是coiumn为2。
浏览
0
提问于2017-04-09
得票数
0
1
回答
我在R中有一个方阵M，它的所有条目都是实数，介于0.5到1.9之间。我想要建立邻接矩阵的条件是，当每个元素小于0.6时，该元素应该替换为零，其他wise元素应该替换1。这是我要对seq(0.5, 1.9, 0.01)中所有141个阈值做的，这样我就可以得到141个邻接矩阵。我怎么能拿到这个？我如何在R中保存或打印所有这些矩阵？任何帮助都将不胜感激。请容忍我在R中知识贫乏:-)
浏览
2
提问于2015-02-20
得票数
1
回答已采纳
3
回答
为什么打印"B“比打印"#”慢得多？
、
、
、
、
我生成了两个1000 x 1000矩阵
第一矩阵：O和#。
第二矩阵：O和B。
使用以下代码，第一个矩阵需要8.52秒才能完成：
Random r = new Random();
for (int i = 0; i < 1000; i++) {
for (int j = 0; j < 1000; j++) {
if(r.nextInt(4) == 0) {
System.out.print("O");
} else {
System.out.print("#");
浏览
6
提问于2014-02-21
得票数
2951
回答已采纳
2
回答
检查计算孤岛时的异常
、
、
、
当涉及到矩阵遍历和任何类型的搜索时，我发现我经常需要验证一个潜在的项目，以添加到我的队列或递归。
假设我们看的是这个矩阵：
var matrixTest = [
[1,1,0,0,0],
[1,1,0,0,0],
[0,0,1,0,0],
[0,0,0,1,1]
在本例中，我迭代所有项，并在遇到1时执行BFS (标记为已访问)。
我在BFS函数中有一个子例程，在给定坐标r和c的情况下，它会验证从坐标r和c指向的所有可能路径(r+1，r-1，c+1，c-1)。
问题出在这段代码上
if(!visitMap[r+1][c] &&am
浏览
16
提问于2016-08-07
得票数
1
回答已采纳
2
回答
我能用Rcpp加速我的R代码吗？
、
、
、
、
我定义了一个R函数，它包含一个矩阵、一个向量和一个参数a。我需要为a的不同值计算函数的结果。在R中编写代码很简单，但当矩阵“大”且参数值的数目很大时，则非常慢。
我可以在R中定义函数并在Rcpp中执行for循环吗？
它能加快计算速度吗？
foo函数在R中的一个最小示例是
f <- function(X,y,a){
p = ncol(X)
res = (crossprod(X) + a*diag(1,p))%*%crossprod(X,y)
set.seed(0)
X <- matrix(rnorm(50*5),50,5)
y <- rnorm(50)
a &
浏览
0
提问于2018-07-22
得票数
2
回答已采纳
1
回答
矩阵或嵌套列表的部分列和
、
、
、
、
下面列出的代码是试图获得另一个矩阵的部分列和的矩阵，因此得到的矩阵的行row[r]是原始矩阵从row=0到row=r的部分列和。
例如，给定
A = [[0,0,0],
[4,5,6],
[7,8,9],
[10,11,12]]
B=[[0,0,0],
[4,5,6],
[11,13,15],
[21,24,27]]
是否允许我消除以下代码中的for-循环，并允许我使用纯列表理解？，您认为列表理解比循环或map/lambda的计算效率更高，前提是我正在处理的实际矩阵相对较大？
我的当前代码如下：
import numpy as np
#
浏览
1
提问于2020-02-12
得票数
2
回答已采纳
1
回答
在经过一些计算之后，我正在思考如何从输入矩阵的每一行中写入多个文件。我现在使用的代码如下所示：
akl <- function(dii) {
ddi <- as.matrix(dii)
m <- rowMeans(ddi)
M <- mean(m) # mean(ddi) == mean(m)
r <- sweep(ddi, 1, m)
b <- sweep(r, 2, m)
return(b + M)
require(plyr)
akl.list <- llply(1:nrow(aa
浏览
0
提问于2013-01-29
得票数
0
回答已采纳
1
回答
给定一个numpy数组M，我想要计算这个数组的每两个矩阵组合的矩阵乘积M[i] @ M[j].T。在对该矩阵(乘积)应用了一些操作之后，我希望将结果存储在[i,j]位置的另一个矩阵中。有没有一种方法可以在不重复两个嵌套循环的情况下快速计算这个值？
也就是说，我想避免的(因为这实际上需要几个小时)是：
import numpy as np
M = np.random.rand(7000,3,3)
r = np.zeros((len(M), len(M)))
for i in range(len(r)):
for j in range(len(r[0])):
n = M
浏览
1
提问于2019-06-01
得票数
0
回答已采纳
1
回答
在调用下面的递归函数时，我遇到了内部错误。
下面的函数检查矩阵表中的引用id (父id的一种类型)，如果引用id在矩阵表中可用，则执行应该返回该id，但如果在矩阵表中找不到该id，则再次从旧引用id(存储在$ref_id中)的成员表中获取引用id(enroller_id)，并再次调用该函数，然后继续调用，直到获得特定id可用或不在矩阵中为止，最后如果找不到，则返回默认id($default_ref_id)。
但是我不知道为什么我在那里收到内部服务器错误，我也尝试过将max_execution_time增加到30到50或100，但无法修复。
请，这将是伟大的，如果你能帮助我，我已经被困在这个问题
浏览
1
提问于2015-03-17
得票数
1
1
回答
我正在使用Matlab和我想利用我的2个GPU，我有一个大的矩阵，不适合1 GPU，但一半的矩阵适合1 GPU。所以我想把矩阵切成两半，让我的两个GPU中的每一个工作在矩阵的一半上。我所做的是
parpool('local',gpuDeviceCount)
gpuDevice
dp = 0.00001;
R = zeros(1,2);
parfor k=1:1:2
if k==1
M = gpuArray([dp:2*dp:10])
M = gpuArray([2*dp:2*dp:
浏览
3
提问于2020-04-14
得票数
0
1
回答
我尝试在Theano中实现非负矩阵分解。更详细地说，我试图找到两个矩阵L和R，以便它们的乘积L x R尽可能精确地表示给定矩阵M。
为了找到L和R矩阵，我使用反向传播。在某种程度上，我已经注意到L和R中的值可能是负值(当然，没有什么可以阻止back支持这样做)。我试图通过在反向传播步骤之后添加以下行来纠正这种行为：
self.L.set_value(T.abs_(self.L).eval())
self.R.set_value(T.abs_(self.R).eval())
在那之后，我的程序变得更慢了。
我做错了什么吗？我是否以错误的方式更新了张量的值？有更快的方法吗？
按照评论中的要
浏览
1
提问于2016-12-05
得票数
0
回答已采纳
1
回答
如何从for循环matlab中存储一系列向量
、
、
、
、
我有一个for循环，它生成一个向量。我想把这些向量存储在一个矩阵中。通常我会这样做：
for r=1:100
vec=[x:y]+r;
mat(:,r)=vec
但这不起作用，因为我有这样的东西：
dr=10/20
for r=1:dr:20
vec=[x;y]+r;
我现在如何将向量存储在矩阵中？因为我不能使用r作为列索引，因为r的值通常不是整数。
浏览
5
提问于2014-07-22
得票数
0
回答已采纳
2
回答
将Excel数据导入R
、
、
、
我正在处理一个由261 x 10矩阵组成的excel文件。该矩阵由10只股票从2010年到2015年的周回报组成。因此，我有10个变量(股票)，每个变量有261个观察值(每周回报)。对于我的硕士论文，我必须在我的矩阵上应用Rüschendorf和Puccetti (2012)开发的“重排算法”。我不会深入讨论这个概念的理论方面的更多细节。问题是，我下载了一个能够在R中执行重排算法的包。我对它进行了测试，它工作得很好。
实际上，我唯一需要知道的就是如何将我的excel矩阵导入到R中，以便能够对其执行重排算法。我可以通过使用R中的矩阵编程公式对矩阵的每个元素进行编码，将矩阵重写为R(手动)：
A
浏览
5
提问于2015-07-28
得票数
1
1
回答
我试着在R中创建一个循环，在每2行中用1替换我的矩阵0，但我只能创建一个数字为1的矩阵，我不知道如何快速而正确地实现它！
repeat{
while(i<3){
if(l<7) {
r<-rbind(r,1)
l<-l+1
}else{
r<-cbind(r,1)
i<-i+1
开始示例矩阵
浏览
2
提问于2013-08-29
得票数
1
1
回答
下面是从R运行的C代码，用于将矩阵的每一行与向量进行比较。相同值的数目存储在两列矩阵的第一列中。
我知道在R中可以很容易地做到这一点(就像检查结果一样)，但是对于更复杂的用例来说，这是第一步。
当未使用openmp时，它工作正常。当使用openmp时，它给出了相关(0.99)，但结果不一致。
Question1:我做错什么了？
Question2:我使用双for循环来填充输出矩阵(ret)中的零。什么是更好的解决方案？
此外，在包中使用代码时也会观察到不一致。我试图使用inline使代码可复制，但它不识别openmp语句(我试图在cfunction参数中包含'omp.h'，.)
浏览
1
提问于2019-02-13
得票数
0
1
回答
我现在必须学习PageRank，我写了这样的算法：
我已经确定了while循环内部的复杂性，我认为它是O(n^2)。但是，接下来我被时间循环本身的复杂性困住了，它本身就是由增量决定的。Delta是R在迭代L1中的i范数和R在迭代i+1上的差。
是否有任何方法来决定这样的复杂性？
编辑:更多解释：
R是秩向量。
epsilon是用户给出的值(我们希望epsilon很小才能很好地估计PageRank，但不要太小，需要很长时间才能计算出来)
creerMatricePageRank()为我们绘制邻接矩阵，并在混合中加入秩源向量E。
A是邻接矩阵。
浏览
1
提问于2017-11-24
得票数
2
1
回答
可以帮助我将可除矩阵的大小转化为较小的矩阵，用元素求出从6x6矩阵到2x3矩阵的和平均值i，e，其中r1c1、r1c2、r2c1、r2c2、r3c1、r3c2等元素的平均值在2x3矩阵的r1c1上移动。
浏览
4
提问于2015-11-16
得票数
0
回答已采纳
1
回答
高斯消元并行性
、
、
、
我已经成功地在CUDA中实现了一个用于高斯消除的单线程程序，并希望实现并行性。到目前为止，并行代码看起来像这样：
__global__ void ParallelGaussian(double* A)
int index = threadIdx.x;
int stride = blockDim.x;
if (index < ROWS) //Skip additional threads
for (unsigned int r = index; r < ROWS; r += stride)
浏览
27
提问于2017-12-16
得票数
0
1
回答
将特定公式应用于R中的大型矩阵/数据帧，以使我的数据规范化
、
、
、
、
我有一个很大的“距离矩阵”(实际上是R中的170x170数据帧)，例如：
A B C
A 0.198395022 0.314012433 0.32704998
B 0.314012433 0.262514533 0.318539233
C 0.32704998 0.318539233 0.211224133
我正在尝试应用一个特定的公式(我已经有了)，将这种变化纳入0-1的范围，这是我的统计建模所需要的。我期望在整个数据框架中获得类似的结果(在应用公式时，期望输出)：
A B
浏览
4
提问于2020-01-09
得票数
0
1
回答
我想写一个函数，使用memcpy将两个矩阵复制到一个矩阵中。这就像Matlab中的C= A；B。因为，我想对所有类型的数据使用这个函数，这意味着整型和浮点型矩阵，所以我使用了void指针。
void Repmatrix(void *M3, void *M1, int R1, int C1, void *M2, int R2, int C2, size_t t)
if (RC == 'R')
int i;
for (i = 0; i < R1*C1; i++)
me
浏览
3
提问于2014-03-09
得票数
0
3
回答
处理列表对象
、
我在R中有一个list对象，它包含更多的列表，每个列表包含三个向量。生成三个矩阵的最快方法是什么?首先是第一个矩阵，第二个矩阵是行，第二个矩阵是行，第三个矩阵是第三个矩阵。例如，考虑到：
metalist <- list(list(c(1,1),c(11,11),c("a","a")),
list(c(2,2),c(22,22),c("b","b")),
list(c(3,3),c(33,33),c("c","c")))
我想
浏览
13
提问于2016-03-15
得票数
4
回答已采纳
2
回答
查感知Jacobi 2D微基准测试的混淆结果
、
、
、
、
我的目标是实现一个低延迟版本的jacobi 2d模板。默认版本在下面的代码片段中描述：
for(auto i = 1; i < N - 1; ++i) {
for(auto j = 1; j < N - 1; ++j) {
*B[i][j] = 0.2 * (*A[i][j] + *A[i][j-1] + *A[i][1+j] + *A[1+i][j] + *A[i-1][j]);
for(auto i = 1; i < N - 1; ++i) {
for(auto j = 1; j < N - 1; ++j) {
浏览
14
提问于2022-03-23
得票数
1
2
回答
我在R中有一个很大的(35000 X 3)矩阵的列表，我想把它们组合成一个矩阵，但它大约有10亿行长，超过了。
包允许更大的矩阵，但支持rbind将多个矩阵放在一起。
是否有其他包或技术支持从较小的矩阵创建非常大的矩阵？
在你问之前，这不是一个RAM问题，只是一个R限制，即使是在64位R上。
浏览
5
提问于2018-05-25
得票数
2
1
回答
MPI allReduce迭代
、
、
我在使用MPI_allReduce时遇到了一个问题，我必须在每次迭代后更新矩阵R_T，R_T_original，因为矩阵R_T在每个过程中都是相同的，但是矩阵mid_L，mid_B，mid_A对于每个过程都是不同的。在我的例子中，矩阵A，B，L有3行，我有3个进程，每个进程都有一行矩阵(A，L，B)，变量no zero有3个不同的值(P0-->no zero=2，P1-->nozero=2，P2--> nozero=3)，在我的调试输出中，我一直收到结果，直到i=1，程序进入索引阶段。我认为进程P2还需要计算迭代i=2，而进程P0和P1已经终止了迭代i。我该如何解决这个问题？
浏览
3
提问于2020-05-21
得票数
0
2
回答
我希望从我的矩阵中删除一个星号(*)，并将该矩阵写入文本文件，其余的元素将在没有空格或任何类型的分隔符的情况下相互连接。我写了这段代码
for b = 1 : length(out7num_r7_nt_back)
if ~(out7num_r7_nt_back(b) == '*')
out7num_r7_back(b) = '';
disp(out7num_r7_nt_back);
dlmwrite('my_data.txt',out7num_r7_nt_back, '');
我收到
浏览
2
提问于2012-11-12
得票数
3
1
回答
代码片段遍历一个一维矩阵。(N是矩阵的大小)。
for (i=0; i< N; i++) // outer loop for Rows
当我在处理器模拟器上运行这段代码来测量TAGE精度时，我意识到随着阵列(N)的尺寸增加，TAGE的精度也会增加。
原因是什么？
浏览
10
提问于2021-01-21
得票数
0
回答已采纳
4
回答
递增2-d java数组的问题
、
、
、
、
三天来，我一直在尝试解决这个问题，但我还是搞不懂。其思想是创建一个矩阵，用一个矩阵进行增量。程序要求行和列的大小，然后创建矩阵。
我举了一个预期结果的例子：
1 2 3 4
5 6 7 8
我得到的是：
1 1 1 1
1 1 1 1
这是我的密码：
public static void main (String[] args) throws IOException {
BufferedReader in = new BufferedReader(new InputStreamReader
浏览
6
提问于2013-07-24
得票数
1
回答已采纳
2
回答
matlab中的余弦相似内建函数
、
、
我想在matlab中计算矩阵不同行之间的余弦相似度。我用matlab编写了以下代码：
for i = 1:n_row
for j = i:n_row
S2(i,j) = dot(S1(i,:), S1(j,:)) / (norm_r(i) * norm_r(j));
S2(j,i) = S2(i,j);
矩阵S1为11000*11000，代码执行非常耗时。那么，我想知道在matlab中是否有函数比上面的代码更快地计算矩阵行之间的余弦相似度？
浏览
1
提问于2018-01-04
得票数
9
回答已采纳
1
回答
替换同一维矩阵列表中的for循环以获得更高效率
、
、
、
、
我有一个相同维度的矩阵列表。在所有矩阵中，大多数值都是零。 这些非零值根据矩阵在列表中的位置而变化，用于替换的值存储在与矩阵列表相同长度的向量中。 因此，我所做的就是使用for循环替换非零值。然而，作为R的新手，我已经认识到这可能不是最有效的方法，特别是因为我有很多矩阵和相当多的非零概率(所有这些实际上都在某个函数中)。 背景是一个自行开发的马尔可夫模型，目前我希望完全控制发生的事情(即不使用任何软件包)，并且矩阵实际上是年龄和状态相关的转换矩阵。 那么问题来了:有没有办法替换for循环，比如apply？我尝试了几个小时，在论坛上搜索了几个小时，但似乎找不到针对这种特定情况的解决方案。谢谢你
浏览
17
提问于2021-08-30
得票数
3
2
回答
我必须生成一些数学运算的几次迭代，每次得到长度为100的向量。
在运行for循环之后，我希望将n次迭代的每个向量存储到一个list或data.frame中。
基本上，将我的结果(长度为100的向量)保存到要转换为矩阵或直接转换为矩阵的列表中。
for (i in 1:10){
r <- rnorm(100,0,1)
在每一次迭代中，r生成长度为100的向量。最终想要的输出是一个包含100行和10列的data.frame，每个列代表每个生成的迭代。
如何建立这样的矩阵或data.frame？
浏览
0
提问于2019-07-31
得票数
1
回答已采纳
1
回答
仿真学习调用函数
、
、
、
我花了好几个小时思考下面的问题。我正在运行一个模拟研究，我想在模拟研究之外定义函数，以便能够在我的代码末尾调用这些函数。
此示例说明了问题，但不可复制(下面您将找到该问题的可复制示例)。我使用"metafor“包进行元分析。
我想使用我在最终模拟代码之外定义的以下函数：
mat <- matrix(NA, nrow = 8, ncol = 3)
funtr.stu <- function(i) {
for (y in 1:8) {
mat[y,i] <- tr[[y]]$k0
return(mat)
"tr“是一个列表，由8次分析
浏览
3
提问于2013-04-02
得票数
0
回答已采纳
1
回答
我正在使用R中的doParallel/Foreach后端并行运行一个进程，我注册了一组由20个核心组成的集群，并运行了大约100次。我把一个矩阵传递给并行进程的每一次迭代，在子过程中，我用它自己的行的随机样本替换这个矩阵。我想知道的是:对于由同一个子进程处理的后续迭代，这种修改是否会持续？例如，当子进程1完成第一次迭代时，它是用原始矩阵开始第二次迭代，还是用随机样本开始第二次迭代？
一个最小的例子是：
library(doParallel)
X <- matrix(1:400, ncol=4)
cl<-makeCluster(2)
clusterExpo
浏览
6
提问于2016-10-11
得票数
1
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
回答已采纳
2
回答
我正在尝试编写一个R脚本，作为第一步，它为输入矩阵的每一行计算dist()和其他内容，然后，作为脚本的第二步，使用在第一步中获得的每对输出矩阵进行另一次计算。我的问题是我不能“保存”第一步得到的所有矩阵。谁能告诉我一个好的策略？
我的代码如下所示：
n<- nrow (aa)
output <- matrix (0, n, n)
for (i in 1:n)
for (j in i:n)
akl<- function (dii){
ddi<- as.matrix (dii)
m<
浏览
0
提问于2013-01-10
得票数
0
回答已采纳
在尝试使用时，我不太适应Matlab的语法，这种语法需要迭代更新包含嵌入在R3中的曲面上的点的矩阵，因此，它由3行和n列组成，第一个元素给出为A:,1，最后一个元素A:，n。
最初，两点将被选择在XY平面沿一条线。用函数f(x，y)将XY平面上穿过A：、1和A:，n的线段上的多个点投影到曲面上，这样我们就得到了初始矩阵，其中第一行包含X轴中的若干坐标；第二行是Y轴中的坐标；第三行是f(x，y)的每个点的曲线高度。每一列都是一点。
因此，如果前两列是A:,1和A:,2，我想得到一个更新的矩阵，其中第一列不变，A:,1，第二列是A:,1和A:,2的内部平均值。换句话说，列(x)的第一项是第一列和第
浏览
0
提问于2020-02-13
得票数
0
回答已采纳
1
回答
应用时移
、
考虑数组r，维数N，m，其中N是时间索引，m是网络中的元素数。考虑到这些 m 元素相互连接，但元素i在元素j中的影响需要一定时间才能发生，即延迟矩阵delayMat给出的每个连接之间都存在延迟，其维数为m，m。
这么说，我有以下问题:我想用数组r的值创建一个矩阵m，m，但是在delayMat之后使用适当的时间筛选(延迟)。这样做的一种方法是：
delay = n - delayMat
drate = np.zeros([m,m])
for i in range(m):
for j in range(m):
drate[j,i]=r[delay[j,i],i]
其中，n
浏览
4
提问于2020-03-23
得票数
2
2
回答
我有一个192 x 3的矩阵，顺序是(192 X 3)：
order(:，1)和order(:，2)都包含1- 16的重复值，order(:，3)包含1和2的重复值。我需要打乱矩阵，同时防止最后一列中相同值的三个以上的任何重复，因此order(:，3)不应该显示超过3个1或2的重复。
这就是我所拥有的，它在较小版本的矩阵上工作得很好，但似乎在较大的矩阵上遇到了问题：
not_good = true;
while not_good
not_good = false;
order = Shuffle(order);
% returns an array of 1s
浏览
0
提问于2012-09-21
得票数
1
回答已采纳
4
回答
从1到52的随机数矩阵不重复，但总是重复1个数。
、
、
、
我用C编写了一个项目，我试图创建一个从1到52的随机数矩阵，不重复，但是总是重复一个数字！
void baralhar(int b[]){
int x,y,r;
for(x=0;x<53;x++){
r=rand() % 52+1;
for(y=0;y<=x;y++){
if(r==b[y])
r=rand() % 52+1;
b[x]=r;
产出: 49 2 3 23 15 50 29 12 33 37 6 21 9 14
浏览
7
提问于2014-01-10
得票数
2
回答已采纳
1
回答
如何构造这个块三对角(稀疏)矩阵？
、
、
、
在R或C++中是否有快速填充(稀疏)矩阵的方法：
A, B, 0, 0, 0
C, A, B, 0, 0
0, C, A, B, 0
0, 0, C, A, B
0, 0, 0, C, A
其中A，B，C是5x5矩阵，0是零的5x5矩阵。
实际上，我使用的矩阵是按列排列的数百到数千行。在R中，我知道可以使用rbind和cbind，但这是一个有点乏味和昂贵的解决方案。
更新:如何使用此矩阵
推荐文章