我有一个长（~7000）的cellarray，每个都包含一个3x3矩阵。我想用相同的三维向量a对这些矩阵中的每个矩阵进行预乘和后乘。此操作可表示为：
[a1 a2 a3]'*M*[a1 a2 a3] ,
其中a = [a1 a2 a3]
我的问题是：如何快速实施？我尝试过一个简单的for循环，结果比在我们拥有的3x3矩阵中重塑矢量a并将其转换为单元数组后使用cellfun要快得多。有什么想法吗？
编辑：也许包含3x3矩阵的单元格数组可以用另一种形式表示，因此可以用完全矢量化的方式计算操作。