这个很简单吧
假设固定矩阵是a，100列的矩阵是b
result<-list()
for (i in 1:100){
result[[i]]<-a*b[,i]#a*b的矩阵乘法和b*a是不一样的，看你自己的需要
}