计算R堆积矩阵的逆的有效方法

Question

我有以下输入：

• R矩阵S_1,...,S_R所有这些都是n-by-n且可逆的。它们存储在单元格数组S{r}，r=1,...,R。
中
• R矩阵x_1,...,x_R所有这些都是n-by-k。它们存储在单元格数组x{r}，r=1,...,R
中

我需要以下输出：

• (n*R)-by-k矩阵B，条目为B = (S_1\x_1,...,S_R\x_R)。

计算此问题的一种方法是通过forloop：

% Preallocate resources
B = zeros(n*R,k);
% Loop:
for rr= 1:R
     B( (rr-1)*n + 1: (rr-1)*n + n,:) =  S{rr}\x{rr};
end

这是一种更有效的方法来执行此计算。特别是，有没有办法避免forloop？