[heinz_2008]

搞定了！！！用三维数组！！！

Y(:,:,1)=eye(4);
for i=1:T
X(i)=i;
Y(:,:,i+1)=Y(:,:,i)*M;
Z(i)=V*Y(:,:,i+1)*S;
end
plot(X,Z,'b');