如何提取对称矩阵非重复元素为行向量并优化批量处理？

我拥有571个625×625的对称矩阵，命名为z1、z2……z571。需要将每个矩阵的非重复元素提取为行向量，例如从z1提取得到行向量c1，从z2提取得到c2，以此类推。我当前用于从z1生成c1的代码如下：

c1=zeros(1,195625);
z1=zeros(625);
for i=1:625
    for j=1:625
        if (i<=j)
            z1(i,j)=c1(1,n);
            n=n+1;
        end
    end
end

第一个问题：这段代码运行耗时较长，是否有提速方法？
第二个问题：若要从z1到z571分别生成c1到c571，并将这些行向量组合成矩阵X，该如何操作？

问题一：提取/填充操作的提速方案 #

首先注意到你提供的代码逻辑和需求反向：代码是把c1元素填充到z1上三角区域，而你实际需求是从z1提取非重复元素生成c1。针对提取操作，完全可以抛弃循环，用Matlab内置向量化操作实现极速提速：

提取非重复元素（生成c1）的最优代码 #

% 生成上三角（含对角线）的逻辑索引，直接提取元素
idx = triu(true(625));
c1 = z1(idx);

提速原理 #

Matlab的嵌套循环因解释执行特性开销极大，而triu等内置函数是底层优化的预编译代码，配合逻辑索引的向量化提取，速度能提升几个数量级。

如果你的原始需求是将c1填充到z1的上三角区域，同样可以用向量化操作替代循环：

idx = triu(true(625));
z1(idx) = c1;

问题二：批量处理生成矩阵X #

推荐方案：先整合为三维数组 #

避免使用z1、z2...z571这类分散变量，先把所有矩阵整合到一个三维数组（维度为625×625×571），再批量提取：

% 1. 合并所有矩阵为三维数组Z
Z = cat(3, z1, z2, z3, ..., z571);

% 2. 生成上三角索引
idx = triu(true(625));

% 3. 批量提取并转置为矩阵X（X为571行，195625列）
X = squeeze(reshape(Z(idx,:), 195625, 571))';

兼容分散变量的方案（效率略低） #

如果无法提前合并，可通过循环配合eval处理分散变量：

% 初始化结果矩阵X
X = zeros(571, 195625);
idx = triu(true(625));

% 循环处理每个矩阵
for k = 1:571
    % 动态获取z1到z571变量
    current_z = eval(['z', num2str(k)]);
    % 提取非重复元素存入X的第k行
    X(k,:) = current_z(idx);
end

内容的提问来源于stack exchange，提问作者mahdokht