1.2260 思路就是把三维矩阵转为一个元胞序列，元胞里的每个元素就是原始矩阵的每一层，然后利用cellfun函数进行批量计算。
c=reshape(vertcat(a{:})*vertcat(b{:})',[size(a)size(b)]);