优化代码CPT3=zeros([1,21^7],"single"); for p3=1:21 for p4=((p3-1)*21+1):p3*21 for p5=((p4-1)*21+1):p4*21 for p6=((p5-1)*21+1):p5*21 for p7=((p6-1)*21+1):p6*21 for p8=((p7-1)*21+1):p7*21 A=[p3,p4-(p3-1)*21,p5-(p4-1)*21,p6-(p5-1)*21,p7-(p6-1)*21,p8-(p7-1)*21]; h=min(A); CPT3(1,(h-1)*21^6+p8)=1; end end end end end end

你可以尝试使用向量化操作来优化这段代码，这样可以减少循环次数，提高代码效率。

可以先创建一个21 x 21 x 21 x 21 x 21 x 21的数组，然后使用reshape函数将其变形为一个 1 x (21^6) 的向量。接着，使用 repmat 函数将 p3~p8 扩展为 1 x (21^6) 的向量，然后进行一些简单的运算，最终将结果赋值给 CPT3。

下面是优化后的代码：

CPT3 = zeros([1,21^7],'single');
A = reshape(1:21^6,21,21,21,21,21);
A = reshape(A,[],1);

p3 = repmat(1:21,1,21^5);
p4 = repmat(1:21^2,1,21^3);
p5 = repmat(1:21^3,1,21^2);
p6 = repmat(1:21^4,1,21);
p7 = repmat(1:21^5,1,21);
p8 = 1:21^6;

A1 = p4 - (p3-1)*21;
A2 = p5 - (p4-1)*21;
A3 = p6 - (p5-1)*21;
A4 = p7 - (p6-1)*21;
A5 = p8 - (p7-1)*21;
h = min([p3;A1;A2;A3;A4;A5]);

CPT3(1,(h-1)*21^6+p8) = 1;

这样，你就可以用更少的代码行数和更快的速度得到相同的结果。