R中的矩阵计算

R中的矩阵计算,r,matrix,plot,R,Matrix,Plot,我有3个向量:a(101乘45),b(101乘45)和k(1乘45)被操纵到a+b*k^T,其中a是rep(a,45),b是rep(b,45)和k是T(矩阵(rep(k,45),ncol=45))我有matlab背景,把一切都看成是矩阵,这个观点对R来说似乎不是最好的。我很难适应,并且已经创建了以下脚本,我相信可以大大改进 nof=45 A = matrix(rep(a,nof),ncol=nof) B = matrix(rep(b,nof),ncol=nof) K = t(matrix(rep

我有3个向量:a(101乘45),b(101乘45)和k(1乘45)被操纵到a+b*k^T,其中a是
rep(a,45)
,b是
rep(b,45)
和k是
T(矩阵(rep(k,45),ncol=45))
我有matlab背景,把一切都看成是矩阵,这个观点对R来说似乎不是最好的。我很难适应,并且已经创建了以下脚本,我相信可以大大改进

nof=45
A = matrix(rep(a,nof),ncol=nof)
B = matrix(rep(b,nof),ncol=nof)
K = t(matrix(rep(k,nof),ncol=nof))
mf <- A+B%*%K
plot(mf)
nof=45
A=矩阵(代表(A,nof),ncol=nof)
B=矩阵(代表(B,nof),ncol=nof)
K=t(矩阵(代表(K,nof),ncol=nof))

mf请确保您的代码是。请包括
a
b
k
,以便您的代码可以运行。我发现代码是可读的;快捷方式可以是
a请注意,Stack Exchange有一个站点codereview.stackexchange.com,您可以在该站点发布运行的代码并获得有关编码风格的反馈。这样的问题更适合那里。如果通过包含
a
b
k
来更新代码以运行,我可以迁移它。