R中的矩阵代数

Question

我有:x = 1 x 999布朗运动矢量y = 1 x 999布朗运动矢量(我可以很容易地模拟这些)

我想设置一个名为z的1000x1000矩阵,如下所示: - 第一行和第二列将是零 - 每个其他元素将是x和y的元素元素,即z的第二列的第二行将是x [ 1]*y [1]等直到第1000列的第1000行将是x [999]*y [999](并且再举一个例子,第4列的第3行将是x [3]*y [4] ]

我该怎么做呢？

Answer 1

您正在寻找outer:

x <- 1:3
y <- 2:4
cbind(0, rbind(0, outer(x, y)))
#      [,1] [,2] [,3] [,4]
# [1,]    0    0    0    0
# [2,]    0    2    3    4
# [3,]    0    4    6    8
# [4,]    0    6    9   12

如果你想绘制z每一对x和y值,你可能会发现使用它更方便

to.plot <- expand.grid(x=c(0, x), y=c(0, y))
to.plot$z = to.plot$x * to.plot$y
to.plot
#    x y  z
# 1  0 0  0
# 2  1 0  0
# 3  2 0  0
# 4  3 0  0
# 5  0 2  0
# 6  1 2  2
# 7  2 2  4
# 8  3 2  6
# 9  0 3  0
# 10 1 3  3
# 11 2 3  6
# 12 3 3  9
# 13 0 4  0
# 14 1 4  4
# 15 2 4  8
# 16 3 4 12

然后你可以用以下东西绘图:

library(scatterplot3d)
scatterplot3d(to.plot$x, to.plot$y, to.plot$z)