如何用循环为不同.lvm文件生成独立Octave绘图？

问题根源分析 #

你遇到的所有图都复用第一个文件数据的问题，核心原因是config_horizontale函数里硬编码了i=1——不管主循环遍历到第几个文件，函数都只会取三维数组里第1层（也就是第一个文件）的数据来绘图。另外，主循环里的N1m赋值写法也存在索引混乱的隐患。

修复步骤&完整代码 #

我帮你调整了代码，主要做了这些优化：

• 把主循环里的文件索引i作为参数传给绘图函数，让函数明确知道当前要处理哪个文件的数据
• 修正函数内部的索引逻辑，用传入的索引获取对应文件的均值
• 优化主循环里的数组赋值，确保每个文件的数据正确存储到三维数组的对应位置
• 给每个绘图窗口加上文件名标题，方便你区分不同文件的结果

主程序代码

myfiles = dir('*.lvm'); % 获取所有.lvm文件列表
% 预分配三维数组，避免循环中动态扩容（可选但能提升运行效率）
file_count = length(myfiles);
n = zeros(size(dlmread(myfiles(1).name, '', 2, 0)), file_count); % 根据第一个文件的尺寸预分配
N1m = zeros(size([n(:,2,1) n(:,3,1) n(:,4,1) n(:,5,1)]), file_count);

%% 循环读取文件并处理
for i = 1:file_count
    current_file = myfiles(i).name; 
    % 读取文件：跳过前2行，从第0列（即第一列）开始读
    n(:,:,i) = dlmread(current_file, '', 2, 0); 
    % 提取当前文件的第2到第5列数据
    CM1N1 = n(:,2,i); 
    CM2N1 = n(:,3,i); 
    CM3N1 = n(:,4,i); 
    CM4N1 = n(:,5,i); 
    % 将当前文件的4列数据存入N1m的第i层
    N1m(:,:,i) = [CM1N1 CM2N1 CM3N1 CM4N1]; 
    
    % 调用绘图函数，传入当前文件索引、文件名
    [Nmean1m] = config_horizontale(N1m, i, current_file); 
end

修正后的config_horizontale函数

function [Nmean1m] = config_horizontale(N1m, file_idx, file_name)
    % 计算当前文件（第file_idx层）的均值
    Nmean1m = mean(N1m(:,:,file_idx)); 
    xM = 20:20:80; 
    
    % 创建独立的绘图窗口
    figure
    % 绘制当前文件的均值数据
    plot(xM, Nmean1m, '-+')
    % 添加标题和坐标轴标签，提升可读性
    title(['文件：', file_name])
    xlabel('X位置')
    ylabel('均值数据')
    grid on % 可选：添加网格线，更易观察数据趋势
end

测试说明 #

把你提供的示例数据保存为多个.lvm文件（比如data1.lvm、data2.lvm等），运行上面的代码后，每个文件都会生成独立的绘图窗口，标题显示对应的文件名，绘图数据也会精准对应各自文件的均值结果。

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