我需要做一个以矩阵为参数的迭代,其中一部分大概是这样:
x0=matrix(......)
y0=matrix(......)
for(i in 1:10)
{
m[i]=t(x[i-1])%*%y[i-1]%*%t(y[i-1])%*%x[i-1]
x[i]=3*x[i-1]
y[i]=2*y[i-1]
i=i+1
}
我想表达的是:x[i]是迭代中的第i个矩阵,以此来做循环,但是R中却判定x[i]是矩阵x的第i个分量。
所以想请教下高手们,这程序应该怎么写?
不知道我说清楚了没。。。先谢过了!