【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))