矩阵乘方求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

想做一个矩阵乘方的程序。有两个序列a<-c(1,2,3), b<-c(1,2),现在想做 c=a^b, 最后得到一个3*2矩阵c=(1,2,3, 1,4,9)，程序如下

a=c(1,2,3)
b=c(1,2)
d<- matrix(, nrow=3, ncol=2)
for (i in 1:3){
d=a^b
}
d后面是中括号i，不知为什么打不上


但是总是得到number of items to replace is not a multiple of replacement length的错误，不知道该怎么修改程序，请大家指点，谢谢啦。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：程序

向量的长度不匹配

fanyonghui 发表于 2015-12-15 20:04
向量的长度不匹配

有没有什么解决的办法呢

我在我的计算机上试了一下, a^b 就出错了. 得到一个1$\times 3$的向量. d是一个矩阵.

fanyonghui 发表于 2015-12-15 20:15
我在我的计算机上试了一下, a^b 就出错了. 得到一个1$\times 3$的向量. d是一个矩阵.

有一列的运算结果被覆盖了，没有输出。这个程序有问题但是不知怎么改好，能不能把你的贴出来看一下呢

a=c(1,2,3)
b=c(1,2)
d <- matrix(, nrow=3, ncol=2)
for (i in 1:2){
  d[,i]=a^b[i]
}

蓝色 发表于 2015-12-15 22:01
a=c(1,2,3)
b=c(1,2)
d

谢谢蓝色老师

1. a <- c(1,2,3)
   
2. b <- c(1,2)
   
3. d <- sapply(b, function(x) a ^ x)

复制代码

> a=c(1,2,3)
> b=c(1,2)
> c <- outer(a,b,"^")
> c
     [,1] [,2]
[1,]    1    1
[2,]    2    4
[3,]    3    9