如何在Matlab循环中利用变量名片段优化重复代码？

核心思路 #

放弃使用不推荐的eval，通过结构体统一管理变量+多层循环遍历分组/变量名的方式，大幅减少重复代码，同时保证代码的可读性与性能。

推荐实现方案（更规范易扩展） #

先将分散的输入变量整合到一个大结构体中，再通过循环遍历分组（sd/base）和角度变量（psi_rad/theta_rad）来完成逻辑：

pp = 4;

% 1. 把所有输入变量整合到一个大结构体（避免动态变量名，规范管理）
input_data.sd.psi_rad = psi_rad_sd;
input_data.sd.theta_rad = theta_rad_sd;
input_data.base.psi_rad = psi_rad_base;
input_data.base.theta_rad = theta_rad_base;

% 2. 预分配输出结构体（提升运行效率，MATLAB推荐操作）
corresponding = struct( ...
    'sd', struct('t', cell(1,pp), 'psi_rad', cell(1,pp), 'theta_rad', cell(1,pp)), ...
    'base', struct('t', cell(1,pp), 'psi_rad', cell(1,pp), 'theta_rad', cell(1,pp)) ...
);

% 3. 多层循环处理逻辑
for p = 1:pp
    % 遍历sd和base分组
    for group = {'sd', 'base'}
        group_name = group{1};
        % 获取当前参与者、当前分组下的角度结构体
        psi_struct = input_data.(group_name).psi_rad(p);
        theta_struct = input_data.(group_name).theta_rad(p);
        
        % 计算对应时间戳
        common_t = intersect(double(psi_struct.time), double(theta_struct.time));
        corresponding.(group_name)(p).t = common_t;
        
        % 提取对应时间戳的psi_rad数据
        corresponding.(group_name)(p).psi_rad = psi_struct.value(ismember(psi_struct.time, common_t));
        % 提取对应时间戳的theta_rad数据
        corresponding.(group_name)(p).theta_rad = theta_struct.value(ismember(theta_struct.time, common_t));
    end
end

% 可选：拆分出独立的两个结构体（匹配你的需求）
corresponding_sd = corresponding.sd;
corresponding_base = corresponding.base;

方案优势 #

• 完全规避eval，代码更安全、易维护
• 扩展成本极低：新增角度变量只需在input_data中添加对应字段；新增分组只需在循环列表中加字符串
• 预分配结构体避免动态扩容的性能损耗
• 变量集中管理，逻辑更清晰

另一种实现（无需整合大结构体） #

如果不想调整输入变量的存储方式，可通过getfield和assignin动态获取/赋值变量，但不如上面的方案直观：

pp = 4;
groups = {'sd', 'base'};
angle_vars = {'psi_rad', 'theta_rad'};

% 预分配输出结构体
corresponding_sd = struct('t', cell(1,pp), 'psi_rad', cell(1,pp), 'theta_rad', cell(1,pp));
corresponding_base = struct('t', cell(1,pp), 'psi_rad', cell(1,pp), 'theta_rad', cell(1,pp));

for p = 1:pp
    for group_idx = 1:length(groups)
        group_name = groups{group_idx};
        % 获取当前分组下的psi和theta结构体
        psi_struct = getfield(evalin('base', sprintf('psi_rad_%s', group_name)), p);
        theta_struct = getfield(evalin('base', sprintf('theta_rad_%s', group_name)), p);
        
        % 计算公共时间戳
        common_t = intersect(double(psi_struct.time), double(theta_struct.time));
        % 更新输出结构体的t字段
        output_struct_name = sprintf('corresponding_%s', group_name);
        output_struct = evalin('base', output_struct_name);
        output_struct(p).t = common_t;
        
        % 遍历角度变量处理数据
        for var_idx = 1:length(angle_vars)
            var_name = angle_vars{var_idx};
            angle_struct = getfield(evalin('base', sprintf('%s_%s', var_name, group_name)), p);
            output_struct(p).(var_name) = angle_struct.value(ismember(angle_struct.time, common_t));
        end
        
        % 将更新后的结构体存回工作区
        assignin('base', output_struct_name, output_struct);
    end
end

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