确实比较简单 能把完整的程序贴出来不 因为不知道你贴的F11是什么东东

S11=[];S21=[];S31=[];S41=[];
S12=[];S22=[];S32=[];S42=[];
l=ceil(rand(1)*(2491-90))
for j=1:41
a=[F11{j}(l:l+90)]';
S11=[S11,a];
b=[F21{j}(l:l+90)]';
S21=[S21,b];
c=[F31{j}(l:l+90)]';
S31=[S31,c];
d=[F41{j}(l:l+90)]';
S41=[S41,d];
e=[F12{j}(l:l+90)]';
S12=[S12,e];
f=[F22{j}(l:l+90)]';
S22=[S22,f];
g=[F32{j}(l:l+90)]';
S32=[S32,g];
h=[F42{j}(l:l+90)]';
S42=[S42,h];
end
其中F11-F42是八个<1*41 cell>的元胞数组，每个胞元都是一个[53x47 double]的矩阵。

这是我写的，数据不多还好，多的话就难办了：S11=[];S21=[];S31=[];S41=[];
S12=[];S22=[];S32=[];S42=[];
l=ceil(rand(1)*(2491-90))
for j=1:41
a=[F11{j}(l:l+90)]';
S11=[S11,a];
b=[F21{j}(l:l+90)]';
S21=[S21,b];
c=[F31{j}(l:l+90)]';
S31=[S31,c];
d=[F41{j}(l:l+90)]';
S41=[S41,d];
e=[F12{j}(l:l+90)]';
S12=[S12,e];
f=[F22{j}(l:l+90)]';
S22=[S22,f];
g=[F32{j}(l:l+90)]';
S32=[S32,g];
h=[F42{j}(l:l+90)]';
S42=[S42,h];
end
其中F11-F42是八个<1*41 cell>的元胞数组，每个胞元都是一个[53x47 double]的矩阵。

能不能用多元数组解决这个问题呢 比如说建立一个41*53*47的数值 根据你的逻辑写计算，然后对对应维度计算平均值

两种办法。
1， 所有结果一股脑放到一个cell数组里面
for h=1:33
   result{h}=rand(3);
end

2，用eval对命名的变量赋值
for h=1:4
for k=1:5
  eval(sprintf('s%d%d=rand(3);',h,k));
end
end