通过预先计算需要哪些元素，然后使用逻辑索引一次获取所有元素，可以使代码运行得更快：

% Make symmetric matrix
matrixsize = 625;
z = rand(matrixsize);
z = z + z';

% Compute unique elements
uniqueelements = tril(true(size(z,1)));

% Pull unique elements from matrix
c = z(uniqueelements);
% Transpose to turn c from a column into a row.
c = c';

然后将每一个放入您的行中：

numMatrices = 571;
X = zeros(numMatrices,numel(c1));
for n = 1:numMatrices
    X(n,:) = % However you're dynamically referencing your c variables
end

如果按照Luis的评论将矩阵存储在一个三维数组中，仍然可以一次性执行所有操作，并且可以免费获得X数组：

matrixsize = 625;
numMatrices = 571;

zs = rand(matrixsize,matrixsize,numMatrices);
zs = zs + pagetranspose(zs);

uniqueelements = tril(true(size(zs,1)));
numElementsPerMatrix = sum(uniqueelements,"all");
uniqueelements = repmat(uniqueelements,1,1,size(zs,3));

cs = reshape(zs(uniqueelements),numElementsPerMatrix,size(zs,3));
cs = cs';

注意：使用tril而不是triu是因为MATLAB以列为主的顺序从逻辑索引返回元素，而原始循环以行为主的顺序遍历矩阵。