【R语言 矩阵相乘】R语言矩阵相乘100次

【D1 D2】2*1
【T1 T2】1*2

要求D1和D2随机的变动, 矩阵相乘100次

rm(list=ls())
gc()

options(scipen = 2000)


##################写成函数###########3


#################定义TT矩阵(1*2)
TT <- matrix(c(1,3),1,2)

DD<- matrix(c(1,2),2,1)
result1 <- DD %*% TT


m1=result1
######################设定随机取的整数范围
x <- 1:100


m=data.frame()
#################变换DD矩阵(2*1)##############
for (i in 2:100){
    D2<- matrix(c(sample(x,1,replace=TRUE),sample(x,1,replace=TRUE)),2,1)
    # print(D2)
    result <- D2%*% TT
    print(result)
    m <- rbind(m,result)
    result1 <- result  %*% result1

}

(finally_result <- result1)
(m_all <- rbind(m,m1))